








key: cord-355267-ndzgxk0k
authors: Kassa, Semu M.; Njagarah, John B.H.; Terefe, Yibeltal A.
title: Analysis of the mitigation strategies for COVID-19: from mathematical modelling perspective
date: 2020-06-05
journal: Chaos Solitons Fractals
DOI: 10.1016/j.chaos.2020.109968
sha: 
doc_id: 355267
cord_uid: ndzgxk0k

In this article, a mathematical model for the transmission of COVID-19 disease is formulated and analysed. It is shown that the model exhibits a backward bifurcation at [Formula: see text] when recovered individuals do not develop a permanent immunity for the disease. In the absence of reinfection, it is proved that the model is without backward bifurcation and the disease free equilibrium is globally asymptotically stable for [Formula: see text]. By using available data, the model is validated and parameter values are estimated. The sensitivity of the value of [Formula: see text] to changes in any of the parameter values involved in its formula is analysed. Moreover, various mitigation strategies are investigated using the proposed model and it is observed that the asymptomatic infectious group of individuals may play the major role in the re-emergence of the disease in the future. Therefore, it is recommended that in the absence of vaccination, countries need to develop capacities to detect and isolate at least 30% of the asymptomatic infectious group of individuals while treating in isolation at least 50% of symptomatic patients to control the disease.

A novel coronavirus, named a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2; previously known as 2019-nCoV), was identified as the infectious agent causing an outbreak of viral pneumonia in Wuhan, China, in December 2019 [47] . The World Health Organization (WHO) medical team codenamed the new outbreak caused by SARS-CoV-2 as "coronavirus disease 2019 (COVID- 19) ". The infection is in the same category as the Severe Acute Respiratory Syndrome (SARS) which emerged in Southern China in 2002, spreading to up to 30 Countries, with a total of 8,098 cases and claiming 774 lives [35] . COVID-19 is also in the same category as the Middle East Respiratory Syndrome (MERS) which was first identified in Saudi Arabia in 2012, and ended up spreading to 27 countries around the world, reaching a total of 2,519 cases confirmed and claiming up to 866 lives [45] .

Since January 2020, an increasing number of cases confirmed to be infected with COVID-19 have been detected outside Wuhan, and currently the infection has spread all over the world. As of May 28, 2020, (09:09 GMT), COVID-19 had affected all continents including island nations (213 countries and territories as well as 2 international conveyances), with the total number of cumulative infections globally standing at 5,807,012 cases and 357,800 deaths [52] , and the numbers still increasing.

The major portal of entry of the virus into the body is the tissue lining the T-zones of the face (including the nose, eyes and mouth). The infection is characterised by loss of the sense of smell (a condition referred to as hyposmia/anosmia), taste and poor appetite. Although, such conditions have been observed in COVID-19 patients, many carriers of the infection may not show any severe symptoms like fever and cough but have hyposmia, loss of taste and loss of appetite.

Whereas knowledge of the virus dynamics and host response are essential for formulating strategies for antiviral treatment, vaccination, and epidemiological control of COVID-19, estimation of changes in transmission over time can provide insights into the epidemiological situation and help to identify whether public health control measures are having a measurable effect [5, 39] . The analysis from mathematical models may assist decision makers to estimate the risk and forecast the potential spread of the disease in the population. Understanding the transmission dynamics of the infection is crucial in designing alternative intervention strategies [29] . In general, by approaching infectious Kassa) mechanisms and adherence to the advice given by public health agencies. In particular, the people in some parts of Asia have fully embraced the measures while many in other parts of the world were very much hesitant to use them. For example, wearing a face mask everyday in public appearances is like a ritual in most of the countries in Southeast Asia, while the same is considered as a bad gesture in many of the other parts of the world [50] (even if it is now becoming a "new normal" also everywhere in the globe). One key difference between these societies and the people in the West is that, the communities in the Southeast Asia have experienced similar disease outbreaks before and the memories are still fresh and painful [50] . That means, recent history of a similar event plays a role in behavioural change of the population especially at the beginning of the outbreak in addition to the perceived threat from the disease. Therefore, in this paper, we consider a mathematical model that takes into account 1. the transmission dynamics of COVID-19 similar to the SIRS model, 2. the contribution of the asymptomatic infectious individuals in the transmission dynamics of the disease in the population, 3. the effect of indirect transmission of the disease through the environment, 4. behavioural change of individuals in the society to apply self-protective measures, and 5. the intensity of historical events from recent similar outbreaks.

By analysing the proposed basic mathematical model, the effect of each of these factors is investigated in terms of their contribution to the control strategies of the disease. Moreover, the use of isolation or quarantining and strict social distancing measures are also considered as mitigation strategies, and a comparative study is made for different scenarios.

The layout of this paper is as follows: The model is described and formulated in the next section. Its qualitative analysis is presented in Section 3. Estimation of the parameters and the sensitivity analysis of the reproduction number of the model with respect to involved parameters are discussed in Section 4. Numerical simulations of the model with some assumed intervention scenarios are also presented and analysed in this same section. Concluding remarks of the study are given in Sections 5.

In this section, we present a mathematical model for the transmission dynamics of COVID-19 which spreads in a population. The susceptible individuals can be infected through either direct contact with infectious individuals or indirect contact with the novel coronavirus infected environment. The population under consideration is grouped into disjoint compartments. Individuals who are susceptible to the disease and without formal awareness about the prevention mechanisms or who did not decide to use any protective mechanisms are grouped in the class. Individuals who are susceptible but are aware of and decide to apply any of the existing protective mechanisms after receiving public health information on how to protect themselves from the novel coronavirus infection are placed in the class. COVID-19 infected individuals who are asymptomatic and symptomatic are grouped in classes and , respectively. Some studies consider the asymptomatic class as the "exposed" class (see for instance [9] ). But since the individuals in this group are known to be infectious and some of them also recover from the disease without going through the group [8] , we used a name "carrier" to avoid confusion. The class contains the recovered individuals from COVID-19. Finally, denotes the amount of the novel coronavirus pathogen that contaminates the environment due to shedding by COVID-19 infectious individuals. In the analysis of the model, we intentionally excluded the actual exposed class for mathematical simplicity. However, a 5 days incubation period is taken into consideration in the numerical simulation part of this paper.

By combining the direct and indirect ways of transmission, the force of infection will have the form

where = + + + + and is the concentration of the novel coronavirus in the environment which increases 50% chance of triggering the disease transmission. The proposed flow diagram for the transmission dynamics of COVID-19 is depicted in Figure 1 while the description of each of the state variables is given in Table 1 . Pathogen concentration on contaminated surfaces or objects in the environment. Table 1 Description of the model variables The dynamics of the pandemic is described by using the following system of differential equations (see Table 2 for the description of the involved parameters):

where

with 0 is the value of the force of infection corresponding to the threshold infectivity in which individuals start reacting swiftly (that means, the point at which the behaviour change function changes its concavity). We append the following nonnegative initial conditions to the system (2):

(0) = 0 , (0) = 0 , (0) = 0 , (0) = 0 , (0) = 0 , and (0) = 0 . Rate of disease transmission from the environment The pathogen concentration in the environment that yields 50% of chance for a susceptible individual to catch the viral infection from the environment Modification parameter (transmission of relative to ) Rate of loosing immunity after recovery Natural death rate Fraction of recovered individuals moving into the class after loosing immunity Table 2 Description of the model parameters

In this section, we study the quantitative and qualitative analysis of the model system Eq. (2).

We begin by determining the biologically feasible set for the model (2) . The following theorem implies that the solutions of (2) are nonnegative and bounded from above, provided that the initial conditions are nonnegative. 

Proof: The proof of this Theorem is outlined in Appendix A.

To determine the equilibrium solutions, we set the right-hand-side of Eq. (2) equal to zero and obtain 

Then, the disease-free equilibrium (DFE) is found to be

The basic reproduction number, which is very important for the qualitative analysis of the model, is determined here below by using the method of the next generation matrix used in [11, 42] . For the model under consideration, using the notation = ( , , ), we have the vector functions

with 1 = + + and 2 = + + represent the rates at which the disease compartments increase and decrease in size due to the infection, respectively. Then the next generation matrix is

where

The basic reproduction number denoted by  0 is defined as the average number of secondary cases produced in a completely susceptible population by a typical infected individual during its entire period of being infectious [11, 42] . Mathematically,  0 is the spectral radius of in Eq. (7) and after further simplification, we obtain

The next result is a direct application of Theorem 2 in [42] . (2) is locally asymptotically stable whenever  0 < 1 and unstable if  0 > 1.

The epidemiological implication of Theorem 3.2 is that the transmission of COVID-19 can be controlled by forcing the dynamics through its parameter values so that  0 < 1 if the initial total numbers in each of the subpopulation involved in Eq.(2) are in the basin of attraction of  0 . To ensure that elimination of the disease is independent of the initial size of the subpopulation, the disease-free equilibrium must be globally asymptotically stable when  0 < 1. This is what we present here below.

where = 1 1

The proof of Theorem 3.3 is carried out using the center manifold theory in [6] and is given in Appendix B.

In the above theorem (Theorem 3.3), if the parameter describing the waning of immunity, , is zero, we can observe that both and are negative. Hence, Eq. (9) fails to be true. In this case, we give below a Theorem which asserts the global stability of the DFE of the model.

The disease-free equilibrium of system (2) in the case when = 0 is globally asymptotically stable for  0 < 1.

To prove the theorem for the case = 0, we use Kamgang-Sallet Stability Theorem stated in [24] . Let = ( 1 , 2 ) with 1 = ( , , ) ∈ R 3 and 2 = ( , , ) ∈ R 3 . Then the system (2) can be written aṡ

and Analysis of the mitigation strategies for COVID-19

We show that the five sufficient conditions of Kamgang-Sallet Theorem (in [24] ) are satisfied as follows.

1. The system (2) is a dynamical system on Ω. This is proved in Theorem 3.1.

2. The equilibrium * 1 is GAS for the subsysteṁ 1 = 1 ( 1 , 0)( 1 − * 1 ). This is obvious from the structure of the involved matrix. 3. The matrix 2 ( ) is Metzler (i.e., all the off-diagonal elements are nonnegative) and irreducible for any given ∈ Ω. Again this is straight forward from the formulation of matrix 2 ( ). 4. There exists an upper-bound matrix̄ 2 for the set

This can also be verified by checking the eigenvalues of̄ 2 under the condition of  0 < 1.

Hence, by the Kamgang-Sallet Stability Theorem [24] , the disease-free equilibrium is globally asymptotically stable for  0 < 1.

□ The reality behind Theorem 3.4 is that, if immunity is permanent ( = 0), coronavirus will be effectively controlled in the community by reducing  0 effectively to a value less than unity.

The novel coronavirus being a new strain of coronaviruses, information about the dynamics of the infection is still evolving. Biological studies of parameter values describing the vital dynamics of the infection are still ongoing as more laboratory tests become available. Although some studies have been done on the early dynamics of the disease most especially on data from Wuhan, extensive reading reveals that some of the disease dynamics parameters are highly variable and some processes are not fully explored. In this work, we use new cases data from Hubei Province of China extracted from WHO situation reports 1-57 [47] , i.e. for the period January 21, 2020 to March 17, 2020. We ought to fit the proposed model to the extracted data and estimate the unknown parameters.

The total population of Hubei province was estimated as 59.2 million. The life expectancy of Hubei province varies depending on the area of dwelling (i.e urban or rural) as well as gender [48] . For urban dwellers, the average life expectancy is estimated to be 75.68 years (with an average being 73.72 years for men and 77.79 years for women). The life expectancy of China of the year 2019 was estimated to be 76.79 years whereas that for the year 2020 is estimated as 76.96 years [32] . Owing to the negligible difference in the provincial and Countrywide value, we use the Country life expectancy for the year 2019 which gives an average mortality rate of = 3.57 × 10 −5 per day. The recruitment rate is thus given as = × 0 , where 0 is taken to be the total population size, 59.2 million.

The average time period taken for symptoms to appear after exposure is observed to vary considerably, with ranges between 2-14 days [7] , 2-24 days [51] , and some outliers going to up to 27 days. The observed median incubation period was nearly 5 days [16] .

The time to recovery from the onset of symptoms varies depending on the seriousness of the infection. Individuals presenting mild illness were observed to recover in an average period of 2 weeks while those presenting serious/critical illness recovering in about 3 to 6 weeks. For our parameter estimation and simulations, we consider a nominal value of 0.0476 day −1 (corresponding to 3 weeks) estimated from an interval (0.0238, 0.0714). We note that a patient is considered recovered: (1) if two swab tests taken in a time interval of at least 24 hours both test negative, (2) if the time taken after the end of respiratory symptoms and fever is at least 72 hours.

The waning of the immunity after recovery is estimated to range between 4 months to 1 year, which gives an interval for as (0.00274, 0.00824) day −1 . For our simulation, we consider a nominal value of = 0.00274 day −1 (approximately one year). We propose that at the end of the epidemic, at least 20 − 65% of the recovering population will learn from the experiences during the infection and even when the acquired immunity wanes, such individuals will become susceptible individuals with past history/knowledge of the disease.

The rate of recovery for the symptomatic individuals ( ) in Wuhan varied considerably but majority of individuals who recovered from the virus were discharged from hospital after 2 1 2 weeks [4] . However, the patients in Wenzhou-China stayed in hospital for 27 days (0.037 per day) on average. In the model fitted on the early trends data from Wuhan-China [28] , the recovery rate obtained for symptomatic cases was 0.0897 per day (accounting for 10 days to recovery). The rate of recovery ( ) for carrier individuals is expected to be higher [9] .

According to the WHO situation report [46] , it is estimated that up to 80% of COVID-19 cases are asymptomatic or show mild symptoms, 15% show severe symptoms and up-to 5% end up with critical infection and require oxygen or a ventilator. The proportion of individuals who do not show symptoms or have mild symptoms can be as high as 94% [18] . For our model fitting, we use a range of (0.65, 0.86) for both and with selected initial values within the prescribed interval. Although Hubei province was put on a lockdown on January 23, 2020, the first major decline in the number of new confirmed cases was only observed on February 20, 2020 (Situation report 31 [47] ), approximately 1 month after the lockdown was imposed. From February 14, 2020, the method of identification of new cases was revised to include both cases confirmed through laboratory tests and clinical observations. As such, there was an observed spike in the number of new cases on February 14, 2020 to 4823 compared to 1508 cases (reported on February 13, 2020) and 2420 cases (reported February 15, 2020).

Applying the above described set of assumptions in the bounds for some of the parameters, we optimized the model output to fit the daily new cases data reported from the Hubei province, China. The parameter values for which the model best fits to the incidence data are given in Table 3 . Figure 2 shows the plot of the reported new-cases data together with the incidence of the disease obtained from the model. As we can observe from the graph, the model Analysis of the mitigation strategies for COVID-19 slightly overestimates the reported data except for the two highest points. In addition, since our model does not assume any control measure at this stage while the reported data after the 31st day may represent the effect of the strict lockdown measure taken by the authorities, the parameters estimated give a good result. We note however that, in some recent work [37, 36] , the use of fractional-order calculus is recommended to get better data fit. When we calculate the value of  0 from Eq. (8) using the estimated parameters given in Table 3 , we obtain  0 ≈ 2.91, which is within the range of values reported in [10, 53] .

We examine the sensitivity of  0 to variations in parameter values and establish the significance of the sensitivity indices. We used the Latin hypercube Sampling (LHS) scheme, which is an efficient stratified Monte Carlo sampling that allows for simultaneous sampling of the multi-dimensional parameter space [19] . For each run, 1000 simulations were done and Partial Rank Correlation Coefficients (PRCCs) [1] calculated between each of the selected input parameters and the disease threshold. The PRCCs indicate the degree of effect each parameter has on the outcome, which in this case is the disease threshold. The sign of the PRCC identifies the specific qualitative relationship between the input parameter and the output variable. The positive value of the PRCC of the variables implies that when the value of the input parameter increases, the future number of cases will also increase. On the other hand, processes underlying the parameters with negative PRCCs have a potential to contain of the number of cases when enhanced. The results of sensitivity analysis are indicated in Figure 3 (a) and the box plot (Figure 3(b) ) gives the five-number summary for the computed disease threshold value from the sampled parameter space.

The processes described by parameters 1 , 2 , and with the greatest positive PRCCs have the greatest potential of worsening pandemic when increased. On the other hand, parameters (ℎ and ) with negative PRCCs have the greatest potential in helping contain the infection when maximised. In this respect, we note that increasing social/physical distancing directly reduces 1 as this lowers the likelihood of a susceptible individual getting in contact with a potentially infected individual. In addition, practising good hygiene (such as regularly washing hands, using sanitisers to disinfect the infected environment and avoiding touching the T-zones of the face) is associated with lowering the likelihood of contracting the virus from infected surfaces. Anything contrary to the above increases the likelihood of getting the infection through the two aforementioned routes. We further note that practising good hygiene also involves the infected individuals reducing the shedding of the virus into the environment. It is evident from the results in figure Table 3 . The processes underlying the parameters 1 , 2 , and have the greatest potential of making the epidemic worse if increased, whereas processes described by and ℎ have the greatest potential of containing the epidemic when enhanced. Table 4 that reducing the rate at which the virus is shed into the environment is significant in reducing the severity of the problem. From the five number summary of the results in Figure 3 (b), the lower quartile of the computed values of  0 is about 2, the median around 2.9 and the upper quartile of about 4. The obtained value of  0 is within the range of 3.11(95% , 2.39 − 4.13) obtained in the early studies in [38] . We note that for a selected combination of underlying processes much higher values of  0 can be obtained, which is an indication of possible worsening of the situation. In a similar way, we observe that for particular underlying processes (a selected combination of parameters) the value of  0 can be reduced to values below one.

As indicated in [1] , we note that although some parameters in the model may have very small magnitudes of PRCCs (non-monotonically related to the disease threshold output), they may still produce sizeable changes in the disease burden. To identify the most important parameters in containing or aggravating the epidemic, we computed p-values for the simulated parameters using Fisher's Transformation [1] . We note that the computed PRCCs are bounded between the interval [−1, 1]. For this matter, some sampling distribution of variables that are highly correlated is skewed. The Fisher's Transformation ( ) = 0.5 log 1+ 1− is used to transform the skew distribution to a normal distribution and then compute p-values for each of the parameters based on the PRCCs [1] . The PRCCs for the parameters together with their corresponding p-values are indicated in Table 4 . We carry out pairwise comparison of the significant parameters (whose p-values are less than 0.05, see Table 4 ) to ascertain whether the processes described by the compared parameters are different. We computed the p-values for the different pairs of significant parameters while accounting for the false discovery rate (FDR) adjustment and the results are given in Table 5 . The major question posed at this point is: Are the different pairs of significant parameters different after FDR adjustment? Based on the FDR adjusted p-values in Table 5 , the compared pair of parameters are rendered to be different if their p-value is less than 0.05 and not different otherwise. We summarise our results in Table 6 , where "TRUE" indicates that the compared parameters are significantly different and "FALSE" indicating otherwise. Table 6 Are the parameters different after FDR adjustment? We observe that the more sensitive parameters are also significantly different (see Table 6 ) except for the 1 − pair which may not necessarily be related.

We examine effect of variation of the sensitive parameters on the reproduction number ( 0 ). The results of the variation of parameters with more negative PRCCs are indicated in the bar graphs in Figure (4) .

From Figure 4 , it is evident that the decay of the virus from the environment (Figure 4(b) ) which can be accelerated by disinfecting surfaces reduces the value of  0 and consequently the disease burden. In addition, we observe that an increased proportion of individuals with knowledge of similar infections from the past that are practising selfprotection and preventive measures (see Figure 4 (a)) is important in slowing down the infection at the initial stage. Such proportions of individuals would normally have knowledge about prevention and control mechanisms of the infection just at the onset of the disease.

We observe in Figure 5 that the increase in person-to-person contact, 1 (Figure 5(a) ), poor personal hygiene, 2 ( Figure 5(b) ), and the rate of shedding of the virus into the environment by both carriers (Figure 5(c) ) and symptomatic individuals ( Figure 5(d) ) increase the value of  0 and therefore the disease burden. It is evident that the most effective way of containing the infection is by minimizing contact, which is why in most cases imposing a lockdown becomes an effective way of slowing the spread of the infection. In addition, good hygiene practices by all individuals are twofold: (1) avoiding touching surfaces, always washing hands with soap and water, or using alcohol based hand sanitizer, which reduced the likelihood of contracting the pathogen from the environment; (2) those who are sick with symptoms like cough and flu, ought to use masks, when they cough or sneeze, must do so in a sanitary tissue which must then be properly disposed off. We also note that hygienic practices without social/physical distancing may not significantly slow down the infection.

In summary, we observe that it is possible to reduce the value of  0 to a value less than unity by reducing only the value of 1 below 0.1 (see Figure 5(a) ). This observation is in direct agreement with mitigation approaches that are aimed at minimising human-to-human contact (such as social distancing and imposing a lockdown). Therefore, the Table 3 paremeter 1 is more influential in the model and can also play a significant role in eradication of the disease. The other parameters (see Figures 4, 5(b) , 5(c) and 5(d)) may reduce the value of  0 significantly when applied in combination but not as independent mitigation processes.

There are various intervention mechanisms for COVID-19 that are being implemented in different parts of the world. The strategies differ from country to country depending on the scientific information available to decision makers. To investigate the outcomes of the mitigation strategies, we include the isolated and/or quarantined classes to the model system (2) . We assume that the individuals in the asymptomatic class ( ) are detected at a rate of and placed in isolation 1 class, while the individuals in the symptomatic class ( ) are identified at a rate of and quarantined in 2 class. Moreover, the system is formulated as a mixed-delay system of differential equations, where the time delay is assumed to account for the incubation period of the disease. Then the system becomes 

where the force of infection is now modified to

with representing the average percentage of contacts reduced due to the social distancing measures and representing the total average rate (or percentage) of disinfecting the environment. In the modified model system (12) , the parameters Analysis of the mitigation strategies for COVID-19 Table 3 .

For the simulation purpose of this study, we considered five different cases or scenarios of how to apply the interventions. The strategies described in each of the cases below are in addition to the awareness creation for voluntary self-protective mechanisms which are widely communicated through various media outlets. Here, we assume that the average effectiveness of the self-protective measures is 15% (as estimated from the data and reported in Table 3) , and the individuals who decided to use any one of them are strict in following the appropriate rules. Here below, we consider each of the scenarios for mitigation strategies case by case.

In this case, we assume that about 40% of the symptomatic infectious individuals and only 1% of the asymptomatic infectious individuals are detected and quarantined. This scenario is based on the assumption that among the people in the class only about 40% show "above mild" symptoms and hence visit health care facilities, while the remaining individuals in this class (nearly 60% of them) remain at home or at large in the society. Then, through contact tracing mechanisms corresponding the hospitalized individuals, some people will be traced and tested, thereby about 1% of the total asymptomatic individuals can be detected.

A similar scenario is being applied currently in some sub-Saharan African countries.

The time profile in Figure 6 shows the situation described in Case 1. From this graph we can observe that the infection stabilizes around its endemic equilibrium, which is nearly at 5000 cases. (This number depends on the initial conditions and the demographic variables of the population under study.) This shows that the disease persists in the population.

Analysis of the mitigation strategies for COVID-19 Figure 6 : Dynamics of the disease with no additional intervention is applied.

In this case, we assume that strict and longer time (6 weeks) of social distancing rules are enforced by the government nearly 4 weeks after the first positive case of COVID-19 is reported in the community.

We assume for simulation purpose that the implementation of the intervention strategy is divided into the following 4 time phases.

Phase 1: The first phase in this case, is 24 days long (measured starting from the time the first positive case of COVID-19 is reported). During this phase, because of lack of information and the nature of the infection, we assume (as in Case 1) that only 40% of symptomatic infectious individuals and 1% of the asymptomatic infectious individuals are detected and quarantined. phase 2: The second phase is assumed to last for 6 weeks (42 days). During this period, it is also presumed that; * 80% of the symptomatic class and 30% of the asymptomatic class are detected and quarantined, * a mandatory strict social distancing rule is imposed, which is assumed to have a 70% reduction of effective contacts of individuals in the society, * environmental disinfection is widely carried out, which is assumed to result in a 50% reduction in the rate of infection from the environment, and to contribute about the same percent impact in increasing the rate of decay of the pathogen from the environment.

Phase 3: The third phase is assumed to be 4 weeks (28 days) long, and is characterised by the partial lifting of the 'lockdown' imposed in Phase 2. During this period, it is assumed further that; * 70% of the symptomatic class and 25% of the asymptomatic class are detected and quarantined, * a relaxed social distancing rule is exercised, which is assumed to have a 25% reduction of effective contacts of individuals in the society, * environmental disinfection is partially carried out, which is assumed to have an impact of reducing the rate of infection from the environment by 30% and increasing the rate of decay of the pathogen from the environment by the same 30%.

Phase 4: The last and fourth phase is the time when the social distancing rule is fully lifted. Due to the lesson learnt from the previous phases, we assume that the following interventions will continue during this period as well; * 70% of the class and 10% of the class are detected and quarantined, Analysis of the mitigation strategies for COVID-19 possible to significantly reduce the infection to a level where it cannot be a public threat. Otherwise, any lower proportion of this effort will imply the resurgence of a second wave of the infection in the community. Therefore, to contain COVID-19 in every given community, public health authorities need to work more on the detection and quarantining of the asymptomatic infectious individuals.

In this case, we assume that there is no lockdown imposed but only a large number of testing is applied to detect and quarantine a larger proportion of infected cases.

If it is possible to intensify the effort of tracing the asymptomatic infectious individuals and be able to quarantine at least 35% of them continuously and effectively, our simulation shows that there is a possibility for the disease to be contained without imposing the strict lockdown rule on the total population. The plot in Figure 10 shows the time profile of the count of the infected groups while about 50% of the individuals from class are effectively quarantined (for example inside appropriate health facilities). We can observe that this intervention strategy can also produce the required result in containing the outbreak as some countries (like South Korea) are currently following this pattern.

In this case, we assume that the length of the lockdown period is nearly twice to the scenario in Case 2. However, the effort in detecting the asymptomatic infectious individuals is kept at the minimum level. This scenario is more applicable in highly resource constrained countries as the current cost of testing is high. In this case, it is assumed that the length of the duration of each phase (except for Phase 2) is the same as that given in Case 2. More still, it is supposed that; 1. the conditions in Phase 1 remain the same, 2. Phase 2 lasts 11 weeks with 50% of the symptomatic individuals and 5% of the asymptomatic individuals detected and quarantined. Moreover, strict social distancing rules are still in place with an effect of reducing 70% of human contacts and 50% of environmental variables, 3. Phase 3 lasts 4 weeks (the same as in Case 2), but with 50% of individuals in the class and 10% of individuals in the class detected and quarantined, while partial social distancing rules are in place with an effect of reducing 25% of human contacts and 25% of environmental variables, 4. Phase 4 continues with detecting and quarantining 50% of members in the class and 10% of members in the class, while the other mandatory intervention are lifted.

The time profile of the infection following the scenario in Case 5 is plotted in Figure 11 .

The simulation for this scenario shows that even if we increase the length of lockdown period to 11 weeks (like it was practised in the Hubei province, China) the disease may re-emerge after some period of time. However, the heights of the peaks in the subsequent waves of the disease are much lower than the first peaks. Therefore, once again, unless the authorities apply some kind of strict contact tracing mechanism and conduct enough testing to detect and isolate up to 30% of the asymptomatic infectious individuals, the disease persists in the community with multiple subsequent waves.

In general, from the simulations, we can observe that in all of the above scenarios a transition from one phase to the other intervention phase is characterised by a surge in new cases. However, the number will eventually go down if the intervention in the immediate next phase is effective, or else the disease re-emerges in the population.

Analysis of the mitigation strategies for COVID-19 Figure 11 : Dynamics of the disease with at most 10% of the population in the class and at least 50% of the class are detected and quarantined just after Phase 1 period, with strict social distancing rule imposed for 11 weeks.

We presented a mathematical model for the dynamics of COVID-19 whose first cases were reported in December 2019 in Wuhan-China. The model incorporates a behaviour change function to account for the proportion of individuals who decided to use any of the self-protective measures and adhered to them. In addition, it also considers a proportion of individuals with a history/knowledge of similar infections from the past and practice necessary protective measures right from the onset of the epidemic. The model also accounts for asymptomatic carriers of the infection as well as the concentration of the pathogen within the environment. The basic properties of the model including well-posedness, the disease free equilibrium and its stability, model basic reproduction number as well as the existence of backward bifurcation were examined. To estimate the parameter values, the model was fitted to the data on daily new cases reported in WHO situation reports 1-57 [47] , which accounts for the period from January 21, 2020 to March 17, 2020. From the nominal values from the data fitting, we obtained a reproduction number,  0 ≈ 2.9 (2.1-4) which compares well with the values of  0 obtained in other researches, for instance, (2.24 − 3.58) [55] and 3.11(95%CI, 2.39 − 4.13) [38] . From our sensitivity analysis simulations, we observed that for some given parameter combinations, the value of  0 can be reduced to below 1, and similarly for values much higher than 4.

We observed that if the recovering individuals do so with permanent immunity ( = 0), then reducing the reproduction number to a value below unity is enough to contain the infection. On the other hand, if recovering individuals do so with temporary immunity ( ≠ 0), the proposed model exhibits backward bifurcation, which implies that reducing the value of  0 below 1 is not enough to contain the infection.

By applying the Latin Hypercube sampling scheme, we observed that if the disease is to be easily contained, measures such as; physical/social distancing (which reduces the rate of disease transmission directly from human to human), improved personal hygiene (which reduces the rate of disease transmission from the environment to humans), and minimal shedding of the pathogen into the environment by both asymptomatic and symptomatic individuals, have the greatest potential of slowing the epidemic when enhanced. We further observed that increased decay of the pathogen from the environment (achieved by disinfecting surfaces) alone is less significant in reducing/curbing the number of cases.

We further observed that having high numbers of people with knowledge from previous similar infections, that are practising the prescribed self-protective measures can delay/slow down the otherwise potentially explosive outbreak. Consequently, the daily number of cases is kept at low manageable levels. In addition, increasing the average effectiveness of the self-protective measures and adherence to such measures is vital in realising low peaks of the number of cases. Furthermore, due to the absence of vaccination or any approved medication, developing capacity to detect carrier groups is very important. From our results, it is recommended that countries should develop capacities to identify and quarantine at least 30% of carriers as well as at least 50% of symptomatic cases if the infection is to be controlled. Our model predicts a possible resurgence of the number of cases, if the asymptomatic cases are still many by the time disease spread curbs/lockdown measures are lifted. In addition, we observe from simulations that although disease spread curbs (such as a lockdown measure) may be imposed, their real impact on the number of new cases may only be realised after approximately 21 days, and the reduction (when it appears) could be sharp in the case of a strict lockdown measure with high impact in reducing effective contact between individuals in the population.

When providing the mitigation strategies, we did not account for the delay between the actual incidence and the point when cases have been confirmed since the actual parameters describing such a delay are not known. In health systems where testing of suspected cases is done after individuals show symptoms or on demand, it is likely to have a big gap between the actual incidence and confirmation of new cases. The impact of the delay between actual incidence and confirmation of case can be explored in future work. In addition, we assumed that all individuals who recover, do so with the same level of immunity. However, this may not necessarily be the case since immunity of individuals is affected by a number of factors, including age, cortisol levels and nutrition among others. The impact of differentiated levels of immunity on the disease dynamics and potential resurgence of the epidemic can be explored in the future when relevant data becomes available. Our model did not include the possibility of vaccination or treatment. We, however acknowledge their importance in controlling the infection. Therefore, optimal control of the infection in the presence of these mitigation strategies can be explored in future works when relevant data becomes available. Since the disease has been observed to affect age groups differently, it is plausible to consider age-structured models to better understand the effect of the disease in the respective age groups. 

Likewise, from the fourth equation of (2), we obtain ′ = (1 − ) + (1 − )(1 − ) + − ( + + ) ≥ −( + + ) .

Solving (15) 

Similarly, using the last two equations of Eq. (2), we have

and

because ( ), ( ), ( ), and ( ) are nonnegative for ≥ 0. Solving Eqs. (17) and (18) gives

and

respectively. Thus, any solution of Eq. (2) is nonnegative for ≥ 0 and any initial condition in Ω. Finally, the total number of the population ( ) at time is governed by

Thus, for the initial data 0 ≤ (0) ≤ , by Gronwall inequality, we obtain

Moreover, for the environmental variable , we have

because ( ) and ( ) are less than for all ≥ 0. Applying again the Gronwall inequality, for 0 ≤ (0) ≤ ( + ) , leads into

Combining the above two steps and Theorem 2.1.5 in [11] for the existence of unique bounded solution, we infer that any solution of Eq. (2) is nonnegative and bounded. Hence, Eq. (2) defines a dynamical system on Ω.

□

Sensitivity and uncertainty analysis of complex models of disease transmission: An HIV model as an example

Effects of information-dependent vaccination behavior on coronavirus outbreak: insights from a SIRI model

Vertically Transmitted Disease: Models and Dynamics

A day-by-day breakdown of coronavirus symptoms shows how the disease COVID-19 goes from bad to worse

Temporal changes in ebola transmission in sierra leone and implications for control requirements: a real-time modelling study

Dynamical models of tuberculosis and their application

Symptoms of coronavirus

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster

A mathematical model for simulating the phase-based transmissibility of a novel coronavirus

Vital surveillances: The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19)-China 2020

Mathematical Epidemiology of Infectious Diseases

Optimal COVID-19 epidemic control until vaccine deployment. medR iv: preprint

Vaccinating behaviour, information, and the dynamics of SIR vaccine preventable diseases

Endemic disease, awareness, and local behavioural response

Modelling the influence of human behaviour on the spread of infectious diseases: a review

Clinical characteristics of 2019 novel coronavirus infection in china

Multiple sources and routes of information transmission: implications for epidemic dynamics

COVID-19: What proportion are asymptomatic?

Sampling and statistical analyses tools (SaSAT) for computational modelling

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

Epidemiology of coronavirus respiratory infections

Mathematical modeling of the spread of the coronavirus disease 2019 (COVID-19) taking into account the undetected infections. the case of China

Analysis of an epidemic model with awareness decay on regular random networks

Computation of threshold conditions for epidemiological models and global stability of the disease-free equilibrium (DFE)

Epidemiological models with prevalence dependent endogenous self-protection measure

The impact of self-protective measures in the optimal interventions for controlling infectious diseases of human population

Effect of negligence and length of time delay in spontaneous behavioural changes for the response to epidemics

Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative

Early dynamics of transmission and control of COVID-19: a mathematical modelling study

Modeling the role of information and limited optimal treatment on disease prevalence

New coronavirus outbreak: Framing questions for pandemic prevention

China Life Expectancy

Modeling the Interplay between Human Behavior and the Spread of Infectious Diseases

Mathematical modeling of COVID-19 transmission dynamics with a case study of wuhan

SARS (Severe Acute Respiratory Syndrome)

Fractional derivatives applied to MSEIR problems: Comparative study with real world data

Modeling chickenpox disease with fractional derivatives: From caputo to atangana-baleanu

Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions

Transmission dynamics of the etiological agent of SARS in Hong Kong: impact of public health interventions

Mathematical modeling of epidemic diseases: A case study of the COVID-19 coronavirus. aR ive

Dynamical Systems and Numerical Analysis

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1

Statistical physics of vaccination

Middle East Respiratory Syndrome: MERS situation update

Coronavirus disease 2019 (COVID-19): Situation Report-46

Coronavirus disease (COVID-2019) situation reports

Total population, life expectancy and ethicity of hubei

Virological assessment of hospitalized cases of coronavirus disease

Coronavirus: Why some countries wear face masks and others don't

Coronavirus Incubation Period

COVID-19 Coronavirus Pandemic

Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in wuhan, china: a modelling study

Duration of antibody responses after severe acute respiratory syndrome

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

The 

Analysis of the mitigation strategies for COVID-19 Figure 7 : Dynamics of the disease with a mandatory 6-weeks lockdown and a 4-weeks of partial social distancing is imposed as described in Case 2 * environmental disinfection is partially carried out, which is assumed to have an impact of reducing the rate of infection from the environment by 20% and increasing the rate of decay of the pathogen from the environment by the same 20%.The time profile of the disease dynamics after implementing the above described interventions strategy is plotted in Figure 7 . The figure shows that the count of the infected individuals decreases down to nearly zero in Phases 2 and 3, but the disease resurges back into the society soon after. However, the peak of the second wave looks to be much smaller than the first one. That means, the intervention mechanisms described in the above 4 phases of this case are not enough to contain the disease, and unless some additional intervention mechanisms are developed the disease persists in the society.

In this case, we assume that early action with shorter time social distancing rule is applied. In this scenario, it is assumed that the interventions described in Case 2 started half way through the time that Phase 2 was implemented in Case 2. That means, the implementation of the interventions described in the four phases of Case 2 is assumed to be followed, but the length of the time in Phases 1 and 2 is reduced as described below. Otherwise, all the details of the interventions in Case 2 are kept the same. The time profile for this set of interventions is given in Figure 8 . The general behaviour of the graph in Figure 8 is the same as that of Figure  7 . However, this strategy has an advantage in significantly reducing the height of the first peak. The height of the subsequent peaks are found to be the same unless some additional measures are taken after Phase 3.Unfortunately, the strategies in both of the above two scenarios (Case 2 and Case 3) do not help to fully contain the disease once it spreads in the population. As it can be seen from Figures 7 and 8 , another wave of outbreaks of the disease emerge at a later stage. Here, we can see that the asymptomatic infectious individuals play the greater role in becoming the major source for the second wave. Therefore, if there is a possibility to track and detect people with asymptomatic infection, and if they can be effectively quarantined for the required period of time, then there is a possibility for the disease to be contained. As it can be observed from Figure 9 , if we can increase the rate of detecting and quarantining the asymptomatic individuals to a proportion of about 30%, it is Appendix A. Proof of Theorem 3.1The proof of Theorem 3.1 is outlined here below based on the following two steps. First, we show that all solutions of Eq.(2) are nonnegative as required in [3, 41] . To show that the state variables and of the model are positive for all ≥ 0, we use proof by contradiction. We suppose that a trajectory crosses one of the positive cones at times 1 or 2 such that:Using the first equation of Eq. (2), the first assumption leads towhich contradicts the first assumption that ′ ( 1 ) < 0. Thus, ( ) remains positive for all ≥ 0. Here, 1 is chosen so that our point to be on the positive axis of ( ) so that ( 1 ) is positive.Using the second equation of Eq. (2), ′ ( 2 ) = ℎ + + > 0, which also contradicts the assumption ′ ( 2 ) < 0. Hence, ( ) remains positive for all ≥ 0. Based on the third equation of Eq. (2),Analysis of the mitigation strategies for COVID-19

Proof: The theorem is the direct application of Theorem 4.1 in [6] . To check the existence of backward bifurcation of the model Eq. (2) at  0 = 1, we use the center manifold theorem [6] . For this purpose, we introduce the following change of variables.so that = 1 + 2 + 3 + 4 + 5 , = 1 ( 4 + 3 ) + 2 6 6 + , and ( ) = 0 + .Moreover, by using the vector notation = ( 1 , 2 , 3 , 4 , 5 , 6 ) , the model Eq.(2) can be written in the form ′ ( ) = = ( 1 , 2 , 3 , 4 , 5 , 6 ) as follows: where, 1 = + + , 2 = + + , and 3 = + . When  0 = 1 and 1 is considered as a bifurcation parameter, from (8) we get Further more, 1 < * 1 if and only if  0 < 1 and 1 > * 1 whenever  0 > 1. The Jacobian of the system (26) at the associated ( 0 ) isThe transformed system Eq. (26), with 1 = * 1 , has a non-hyperbolic equilibrium point such that the linear system has a simple eigenvalue with zero real part and all other eigenvalues have negative real parts. Hence, the centre manifold theory [6] can be used to analyse the dynamics of the model Eq. (26) near 1 = * 1 . By using the notation in [6] , the following computations are carried out.The right-eigenvectorAnalysis of the mitigation strategies for COVID-19 associated with the zero eigenvalue of ( 0 ) such that ( 0 ). = 0 at 1 = * 1 is given byof ( * ) such that The right-eigenvector and the left-eigenvector need to satisfy the condition . = 1.The bifurcation coefficient at the DFE ( 0 ) is given by Analysis of the mitigation strategies for COVID-19Thus, the bifurcation coefficient , can be positive for the right choice of the parametric values that satisfy the condition in Eq. (9) .The second bifurcation coefficient is given by Clearly, > 0 because and are positive. When = 0, and in (30) are negative and in (32) is negative as well. Hence, by Theorem 4.1 in [6] , the model will not exhibit a backward bifurcation at  0 = 1.

☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Semu Mitiku Kassa, Hatson J. B. Njagarah, Yibeltal A. Terefe.

