## ABSTRACT

Recent outbreak of coronavirus disease 2019 (COVID-19) in China has lead a global pandemic around the world. For controlling COVID-19 outbreaks, most countries take two typical intervention strategies: suppression approach like immediately lock-downing cities at epicentre of and mitigation; or mitigation approach that slows down but not stopping epidemic of COVID-19 for reducing peak healthcare demand. Both strategies have their apparent merits and limitations; it becomes extremely hard to conduct one intervention strategy as the most feasible way to certain country. Targeting at this problem, this paper conducts a feasibility study by defining a mathematical model named SEMCR that access the effectiveness of mitigation, suppression and hybrid interventions for controlling COVID-19 outbreaks in London and Wuhan. The model first extends traditional SEIR (Susceptible-Exposed-Infectious-Recovered) model by considering one key fact: there is a direct link between Exposed and Recovered population. Then, it defines threshold and parameters to classify two stages of COVID control: active contain by isolation of cases and contacts, passive contain by suppression or mitigation. The model was fitted and evaluated with public dataset containing daily number of confirmed active cases including Wuhan, London, Hubei province and the UK during January, 2020 and March 2020. The simulated results show that 1) Suppression taken in Wuhan significantly reduces total exposed and infectious populations, but it has to be consistently maintained at least a period of 90 days (by the middle of April 2020) where its success relies on sufficient support from other parts of China. 2) In London, it is possible to take a hybrid intervention of suppression and mitigation for every 2 or 3 weeks over a longer period. While the total infectious populations will be doubled, the economic performance of London will be less affected due to limited mobility constraints. 3) Both in Wuhan and London cases, one important issue of fitting practical data is that there are a large portion (like 42.5% in Wuhan) of self-recovered population who are asymptomatic or mild symptomatic. These people might think they have been healthy at home because they did not go to hospital for COVID-19 tests. Early release of intervention intensity might increase a risk of the second breakout. One limitation of our model is that its prediction of infections and deaths depends on a parameter estimation of intervention intensity that presented by average-number contacts with susceptible individuals as infectious individuals in a certain region. We assume that each intervention has the same effect on the reproduction number in different regions over time. The practical effectiveness of implementing intervention intensity might be varied with respect to cultures or other issues of certain county.

## 1 INTRODUCTION

Throughout human history, Infectious diseases (ID), also known as transmissible diseases or communicable diseases, are considered as serious threats to global public health and economics [1]. From the 1918 influenza pandemic in Spain resulting in nearly 50 million deaths in 1920s, to recent ongoing global outbreak of coronavirus disease 2019 (COVID-19) from Wuhan in China killing over 11 thousands people in all over the world [2], infectious disease is a leading contributor to significant mortality and causes huge losses to society as well as personal family burden. Among a variety of factors leading to emergence and breakout of ID, one primary issue is human mobility where in the cities with developed transportation systems, pathogens can be spread to large geographic space within a short period of time. For instance, the ongoing global epidemic outbreak of COVID-19 has spread COVID-19 to at least 146 countries and territories on six continents in 2 months, infecting Wuhan Cases. As shown in Fig.1 with Wuhan as a stimulated case using data from [6], we demonstrate performance of taking different intervention strategies: a) No interventions: the peak of more than 274 thousands people. Thus in the early stage, many researchers have been working in traditional ID propagation models [3-7] like SIR, SEIR,.et.al, for understanding COVID-19 transmission with human mobility and predicting outbreak process of epidemics. Meanwhile, as realizing a long period of this battle against COVID-19, many of them recently focus on intervention strategies [8-10] that balance a trade-off between limited human mobility and potential economic loss in COVID-19 control. It poses an important research area that explores how and when to take what level of interventions in light of multiple natures and capabilities of countries.

In traditional compartmental models paradigm in epidemiology, SIR (Susceptible-Infectious-Recovered) [3] and SEIR (Susceptible –Exposure-Infectious-Recovered) [4] are two popular approaches to stimulate and predict how infectious disease is transmitted from human to human. These two models have defined several variables that represent the number of people in each compartment at a particular time. As implied by the variable function of time, these models are dynamic to reflect the changes and fluctuations of these numbers in each compartment over time. As for the investigation of public health interventions for COVID-19 control in Wuhan, Zhong, et.al [6] introduced a modified SEIR model in prediction of the epidemics trend of COVID-19 in China, where the results show that under strong public intervention of “lockdown Hubei”, the epidemic of COVID-19 in China would achieve peak by late February and gradually decline by the end of April 2020. Some other extended models [8] [12] are also proposed for predicting the epidemics of COVID-19 in Wuhan and give some similar forecasts. While above methods demonstrate good performance in prediction of COVID-19 breakout by taking strong public intervention, also named as suppression strategy [13] that aims to reverse epidemic growth, one important challenge is that taking suppression strategy only is to treat disease controls as single-objective optimisation of reducing the overall infectious populations as soon as possible, and require strategic consistency in a long term. In real-world, taking public health intervention strategies is actually a multiple-objective optimisation problem including economic loss and society impacts. Thus, most countries have taken different intervention strategies, like enhanced surveillance and isolation to affected individuals in Singapore [14], four-stage response plan of the UK [15], mitigation approaches [13] and even multiple interventions taken in many EU countries. Due to the fact that standalone intervention strategy has apparent merits and limitations, it becomes highly necessary to study the feasibility of intervention strategies to certain country in light of its multiple natures and capabilities.

Targeting at this problem, this paper conducts a feasibility study that analyses and compares mitigation and suppression intervention strategies for controlling COVID-19 outbreaks in London and daily infectious populations is up to 2.5 million, but will be completed in 100 days. The epidemics will last a relatively shorter period of 90 days, but lead to more death.

b) Contain phase: taking 90% effectiveness of surveillance and isolation from the second day of confirmed case potentially enables controlling a new outbreak of COVID-19, but it needs to be followed a long time. c). Suppression intervention from the 32^{nd} day: the peak of daily infectious populations greatly reduces to 30 thousands, but this intervention needs to be followed at least 120 days. Nearly 3 months suppression will lead to high economic loss even crisis. d). Mitigation intervention from the 32^{nd} day: the peak of daily infectious populations increases to 350 thousands, but the period of taking intervention extends to 150 days. It implies a more death but less economic loss compared to suppression intervention. e). Hybrid intervention of taking both suppression and intervention every 2 weeks: the epidemics of COVID-19 appears a long-term multimodal trend where the peaks of daily infectious populations are about 40-60 thousands. This might lead to less daily critical cases and offer more time to hospital for releasing their resources.

Above initial analysis demonstrates the complexity and difficulty of how and when to take what level of interventions. In this paper, we propose a mathematical model: SEMCR to analysis this challenge. The model extends traditional SEIR (Susceptible-Exposed-Infectious-Recovered) model by considering one key fact: there is a direct link between Exposed and Recovered population. Then, it defines threshold and parameters to classify two stages of COVID control: active contain by isolation of cases and contacts, passive contain by suppression or mitigation. The model was fitted and evaluated with public dataset containing daily number of confirmed active cases including Wuhan, London, Hubei province and the UK during January, 2020 and March 2020. For each point, we design and set up experimental protocols for comparison and exploration, highlighting following multi-fold contributions:

Suppression taken in Wuhan significantly reduces total exposed and infectious populations, but it has to be consistently maintained at least a period of 90 days (by the middle of April 2020) where its success relies on sufficient support from other parts of China.

In London, it is possible to take a hybrid intervention of suppression and mitigation for every 2 or 3 weeks over a longer period. While the total infectious populations will be doubled, the economic performance of London will be less affected due to limited mobility constraints.

Both in Wuhan and London cases, one important issue of fitting practical data is that there are a large portion (like 42.5% in Wuhan) of self-recovered population who are asymptomatic or mild symptomatic. These people might think they have been healthy at home because they did not go to hospital for COVID-19 tests. Early release of intervention intensity might increase a risk of the second breakout.

One limitation of our model is that its prediction of infections and deaths depends on a parameter estimation of intervention intensity that presented by average-number contacts with susceptible individuals as infectious individuals in a certain region. We assume that each intervention has the same effect on the reproduction number in different regions over time. The practical effectiveness of implementing intervention intensity might be varied with respect to cultures or other issues of certain county.

The remainder of this paper is arranged as follows. Section 2 introduces an extended SEIR model for controlling COVID-19 breakout. In the Section 3, the materials and implementation of experiment are reported. Section 4 provides detailed experimental evaluation and discussion. The conclusion and future directions are given in Section 5.

## 2 METHODOLOGY

### 2.1 Problem formulation of COVID-19 breakout

We modified the original SEIR model to account for a dynamic Susceptible [S], Exposed [E], Infectious [I] and Recovered [R] or Dead [D] population’s state by extending two components: Mild [M] cases and Critical [C] cases. Conceptually, the modified modal is shown in Fig.2. Here, we assume that susceptible population *E* is initial susceptible population of certain region; and *β* represents the effectiveness of intervention (strict surveillance and isolation) in Contain phase. If effectiveness of intervention in contain phase is not sufficiently strong, susceptible individuals may contract the disease with a given rate when in contact with a portion of exposed population (asymptomatic but infectious) *E*. After a incubation period, the exposed individuals become the infectious population *I* (sympotomatic) at a ratio *κ*. Notably, infectious population starts from mild cases *M* to ciritial cases *C* at a ratio *a*. Finally, a portion *d* of critical cases lead to deaths; the rest of infectious population will be recovered.

There are two enhanced features in our model in comparison to popular SEIR models [6] [8] [12]. The first one is a straightforward relationship between Exposed and Recovered population. It is based on our observation that in the early breakout of COVID-19, some portion of exposed people may have no obvious symptoms or only develop as mild cases, but they might be not able to get an official test due to lack of testing kits. This group of populations might be self-recovered in some days, but will not realize they have actually been infected.

The second feature in our model is that we separate infectious population into mild and critical cases in light of their symptoms. It mainly concerns a phenomena of relatively higher mortality rate in the early breakout of COVID-19 in Wuhan after taking immediate suppression intervention on 23^{rd} Jan 2020. Most likely, suppression intervention has led to fear in Wuhan citizens and chaos in hospital systems, as a result of delaying rescue of many critical cases.

Therefore, by introducing above two features in our model, it will be helpful for evaluating positive and negative effects of different interventions. If we assumed the overall population of a certain region is *N*, the number of days is *t*, the dynamic transmissions of each components of our model are defined as follow:

Regarding Mild cases, Critical cases and Death, the dynamic transmission is as below:

Lastly, we define a benchmark in SEMCR model to reflect the strength of intervention over time, as M_{t}. It is presented by average number of contacts per person per day in a region.

### 2.2 Implementation of dynamic transmission of SEMCR

In order to apply the SEMCR model into practical cases, it needs to estimate the defined parameters in last section, including β, σ and γ^{1}, γ^{2}, a, b, where β is the product of the people exposed to each day by confirmed infected people (k) and the probability of transmission (b) when exposed (i.e., β= kb) and σ is the incubation rate which is the rate of latent individuals becoming symptomatic (average duration of incubation is 1/σ).

According to recent report [8], the incubation period of COVID-19 has been reported to be between 2 to 14 days, we chose the midpoint of 7 days. γ is the average rate of recovery or death in infected populations. Using epidemic data from [6], we modeled the SEMCR model to determine the probability of transmission (b) which was used to derive β and the probability of recovery or death (γ). The number of people who stay susceptible in each region is similar to that of its total resident population. Other transmission parameters are estimated with early prediction of Hubei cases in [6] on January 23 2020 using Monte Carlo simulation, as shown in the Table.1

Notably, as for the strength of intervention M, it is closely related to the population density in a region. We use a benchmark reported in [6] that assumes Hubei province with no intervention as M = 15, and after suppression intervention, M reduced to 3. When applying SEMCR model into other simulated cases, M should be initialized according to the population density and human mobility in these places. Also, after taking any kind of interventions, the change of M should follow a reasonable decline or increase over few days, not immediately occur at the second day.

Following previous assumptions, the implementation of dynamic transmission of SEMCR model follows steps as below:

### 2.3 Model evaluation protocol

In order to utilise our proposed SEMCR model into practical cases, we design an evaluation protocol to access multiple effects of taking different intervention strategies to control breakout of COVID-19 in 4 typical cases, including Hubei province, Wuhan city, the UK and London, as shown in Fig.3.

The first stage is initial parameters estimation using COVID-19 data from four cases: Hubei, Wuhan, UK and London. In this stage, a preliminary qualitative assessment of each case is performed, by comparing their similarity and dissimilarity on area, transportation, population density, migration flows, date of the first confirmed case, etc. We would determine value of initial parameters in the SEMCR including N, M_{t}, and date of the first confirmed case. Notably, in Wuhan, the date of the first confirmed case is not officially released. The work [8] estimates the first confirmed case is by the end of Nov 2019; and Zhong [11] points out the first confirmed case in Wuhan is on 23nd December 2019. Here, we take the same settings of first confirmed case on 23nd December 2019.

The next step is to estimate and normalize other parameters in the model. Assuming that COVID-19 has similar transmission ratio and incubation rate in all four cases, we use parameter values fitted from [11], where incubation rate is 1/7; the rate of transmission for the I to S is 0.157; the rate of transmission for the E to S is 0.787. As for estimation of other parameters, we follow the COVID-19 official report from WHO [15], including the proportion of Mild, Severe and Critical cases, the probability of death, etc.

Thirdly, how to take intervention strategies needs to be evaluated by tuning parameters in SEMCR model. The key tuning operation is to adjust the level of M_{t} over a period. For instance, we assume that no intervention strategies result in unaltered internal mobility of a region, taking suppression strategy in Wuhan means a reduction of M to 3. But in other cases with larger area, it is extremely difficult to take a complete suppression strategy. So the reduction of M will be relatively adjusted to 4 or 5.

Final stage, we perform quantitative analysis of effectiveness of different intervention strategies, including: strict surveillance and isolation, suppression strategy, mitigation strategy, and multiple intervention. The evaluation metric of cross-validation is employed to evaluate the performance of COVID-19 progression model. The final two evaluation indicators are the length of intervention and the peak time. The length of intervention is calculated due the date that confirmed cases are nearly clear to zero.

### 2.4 Data collection

The most recent epidemiological data based on daily COVID-19 outbreak numbers reported by the National Health Commission of China were retrieved. Also regarding the daily update from world meter, we record the number of confirmed cases and death each day in four cases. In order to stimulate four cases, we require the exact confirmed cases in the first 2-5 days to initialise parameters of our model.

## 3 EXPERIMENTS

### 3.1 Evaluation of cases with no interventions

We first simulate four cases that predict the breakout of COVID-19 without taking any interventions. The initial populations are given as London (9.3 million), Wuhan (14.18 million), UK (66.49 million) and Hubei (58.9 million). The parameter M representing average number of contacts per person per day is given as 15 to London, Wuhan and Hubei; 12 to the UK. The simulation results are given in Figure.4.

The results show that in the peak time, there will be up to 100K, 200K, 450K and 900K exposed population that are infection but no symptoms at London, Wuhan, UK and Hubei. This implies that: 1) the total infectious population of these four cases will be 7.68 million in London, 12.15 million in Wuhan, 50.45 million in Hubei and 47.98 million in the UK. 2) The total death of these four cases will be 761K in London, 120K in Wuhan, 0.493 million in Hubei and 0.475 million in the UK. It equals to about over 80% of total population of each region will be infectious, with the morality rate nearly 1%. The above figures show that without intervention, the breakout of COVID-19 will lead to huge infectious population and deaths. The main reason is that COVID-19 is estimated as relatively high production number R_{0} equals to 3 [14], where the transmission ratio β_{2} from Susceptible to Exposed is up to 78%. Thus, in some regions with high migration and dense population, it will easily lead to a breakout.

However, one notable issue is that the initialisation of parameter M seems to impact on the occurrence of peak time and the length of overall period. In the cases of London, Wuhan and Hubei with M = 15, their peak time are all roughly on the 60^{th} day from the date of first confirmed case; but UK with M = 12, their peak time is delayed to the 100^{th} day from the date of first confirmed case. That means to regions with similar total population, the low population density will reduce overall infectious population, but delay the peak time of breakout as a result of longer period.

### 3.2 Effectiveness of surveillance and isolation

As for the strategy of taking surveillance and isolation intervention in the contain phase, recent study [15] has developed a stochastic transmission model parameterised to the COVID-19 outbreak. It proves that highly effective contact tracing and case isolation can control a new outbreak of COVID-19 within 3 months, where for a production number R_{0} of 3.5, more than 90% contacts had to be traced. Thus, we plan to transfer this finding as a tuning parameter β_{2} to evaluate if the outbreak of COVID-19 could be controlled.

Considering that 90% of contacts had to be traced, it implies that the surveillance and isolation is effective to scale down the group of contacts. In other words, we stimulate a situation that from the second day of receiving confirmed case, only 5% of the overall population will be possible to contact the infectious ones. Then, the β_{2} = 0.05. The results are shown in the Fig.5. The results show that in all four cases, the outbreak of COVID-19 has been successfully controlled, and has no peak time. There will be only 44, 56, 42 and 120 people that are infection and finally recovered at London, Wuhan, UK and Hubei. The overall period of COVID outbreak is less than 100 days. This finding is same as the proof in [15].

Notably, the simulated situation of Wuhan, London and Hubei are exactly same, it is mainly because we initialise the same value of confirmed cases as 3, and M = 15 in all cases. When taking highly effective surveillance and isolation, the transmission of COVID-19 has been limited and controlled within a small group of population. The population difference will not affect the total infectious ones. But in the UK case with M = 12, low population density will limit the effectiveness of surveillance and isolation, as a result of more infectious population. It implies another fact that to these countries with low population density, it is challenging to take high quality of intervention like surveillance and isolation.

Then, we evaluate when the outbreak of COVID-19 could occur by adjusting the value of parameter β_{2} in the UK case. As shown in Table.2, we record the total recovered population, and if there is a peak of breakout. As the increased β_{2}, the total recovered population is dramatically increased and generates a peak as β_{2} = 0.075. That means that if we cannot guarantee at least 92.5% of UK population are not contacted by infectious people, it is indispensable to have an outbreak of COVID-19. Meanwhile, if the isolation intervention is less ineffective; the more population will be infectious; and the peak day will be brought forwarded. Thus, in the early stage of control COVID-19, taking high effective surveillance and isolation in any regions is necessary to avoid the outbreak of COVID-19.

### 3.3 Effectiveness of Suppression Intervention

The suppression strategy is widely recognised as the most effective solution to reduce the infectious population, where it was taken in Wuhan on 23nd January 2020. In the work [11], Zhong has reported that taking suppression strategy has successfully limits the overall infectious population in Hubei on 22^{nd} February 2020 to 50K. And if the suppression strategy was taken one week earlier, this figure would reduce to 18K. Thus, we would like to simulate the situation of taking suppression strategy in London. The parameter M is given as 15 to 4 in London, in comparison to Wuhan from 15 to 3. The simulation results are given in Figure.6.

The results show that by taking suppression strategy, the cases in Wuhan and London appear a quite similar trend that daily exposed and infectious population will be greatly reduced. The outbreak of COVID-19 will be controlled by the 100^{th} days, and can be nearly ended by the 150^{th} days. The difference is that the daily infectious population of London will be nearly double to the ones in Wuhan. It is probably because we stimulate the date of taking suppression strategy in Wuhan (the 32^{nd} day) is 3 days earlier than London (the 35^{th} day). Another possible reason is that considering the impact of culture difference, the value of M is only limited to 4 in London, but not as lower as 3 in Wuhan. In fact, the suppression in Wuhan is actually applied to limit mobility in community level, which is hard to be followed by London.

We also consider another two situations of taking suppression intervention 1 week earlier or 1 week later in London, as shown in Fig.6. The results show that the overall infectious population will be greatly reduced to 35 thousands if taking actions one week earlier; oppositely it will increase to 1.5 million infection if taking actions one week later.

However, two important issues of taking suppression strategy are the intervention length and the peak time. In order to effectively control the breakout, the intervention length of taking actions at above three timing points all require at least 100 days. It potentially leads to significant side effects to economic and society, including job loss, mental health, etc. Also, the peak time will arrive earlier than its nature transmission if taking actions one week earlier. This might increase a shorter time to government for preparing sufficient resources, and causing more difficulties in control the breakout of COVID-19.

### 3.4 Effectiveness of Mitigation Intervention

The mitigation strategy is recently highlighted by researchers from Imperial College [16] [17], where this strategy was initially taken by UK government. The aim of mitigation is to use other strategies like drugs to help individuals so that not to interrupt transmission completely, but to reduce the health impacts of an epidemic. In this cases, population immunity builds up through the epidemic, leading to an eventual rapid decline in case numbers and transmission dropping to low levels. Thus, we simulate the situation of the UK by taking mitigation strategy with different level of strengths. The parameter M representing average number of contacts per person per day is given as 12 to the UK, and gradually reduced to 10 and 8 from the 32^{th} day of first confirmed case. The simulation results are given in Figure.7.

The results show that firstly taking mitigation strategy is effective to reduce daily infectious population, further lead to huge reduction of total infectious population. By reducing the parameter of M from 12 to 10 or 8, the peak of daily infectious population is reduce from 4.5 million to 2.5 million or 1 million. On the other hand, the peak time of infectious population will be delayed as taking mitigation intervention, where the peak dates of infectious population at M = 12, 10 and 8 are roughly on the 80^{th}, 110^{th} and 170^{th} days. The overall period of breakout will be extended from 180 days to 200 or 250 days.

Thus, above simulated figures appear a similar trend as findings from [17]. Taking mitigation intervention in the UK is capable of reducing the impact of an epidemic by flattening the curve, reducing peak incidence and overall death. While total infectious population may increase over a longer period, the final mortality ratio may be minimised at the end. But as similar as taking suppression strategy, the mitigation interventions need to remain in place for as much of the epidemic period as possible. However, the timing of introducing this mitigation intervention is important, where too early execution may allow transmission to return once they are lifted and sufficient “herd immunity” has not been developed. Thus, it is necessary to balance the timing of introduction with the scale of disruption imposed and the likely period over which the interventions can be maintained.

### 3.5 Effectiveness of hybrid intervention

In terms of above discussion, the effectiveness of taking any one intervention (either suppression or mitigation) is likely to be limited. It is highly necessary to consider the possibility of taking multiple interventions to be combined to have a substantial impact on social and economic cost reduction. Thus, we have simulated one simple situation of taking multiple strategy in London from the 35 days, by giving a hybrid of suppression and mitigation strategies every two weeks. The M is given as a pattern of 6-4.5-3-4.5-6 in London, where it means mitigation and suppression strategies are taken in an every two weeks roll. For minimizing side effects of taking interventions on human mobility, the application of first mitigation to reduce M from 12 to 6 spends 2 days. The simulation results are given in Figure.8.

The results show that the epidemic appears a multi-modal decline trend over 500 days. The first peak of infectious population occurs on the 50^{th} days with 1000 infections after taking suppression intervention to reduce M from 12 to 3. After two weeks, mitigation strategies are taken so that the second peak of infectious population rises up to 750 infections. The total infectious population will be 530K over 500 days; the deaths will be limited to 520.

Apparently, taking multiple intervention in the UK is capable of reducing the impact of an epidemic by fluctuating the curve, reducing overall infections and death. While the total period will be extended, final mortality ratio may be minimised at the end. But the longer period of limiting human mobility might increase economic risks and reduce employment ratio. There will be plenty of choices to taking multiple interventions through adjusting the strength and length of intervention. The consequence possibility show a similar multimodal curve but with different peak incidence.

## 4 DISCUSSION AND LIMITATION

As we mentioned in section 1, the question on how and when to take what level of interventions to control an epidemic is highly challenging, particularly in light of multiple natures and capabilities of countries. In many cases, it is even hard to evaluate effectiveness until the end of epidemic, and there is always controversy on taking any one intervention. However, our findings contribute to several useful suggestions on controlling COVID-19 breakout.

The first point is that highly effective surveillance and isolation strategy is necessary to control an epidemic in early stage. In theory, if this strategy can be executed in excellent level, there will be no huge breakout later on. The case in Wuhan is a sample that have not taken any effective surveillance and isolation in the early stage of COVID-19, resulting of huge difficulties on late interventions. Also, considering the area, transportation, migration flows and population density of a region, most countries cannot achieve excellent level of isolation in contain phase, probably only in fine or good level. The breakout of an epidemic is avoidable. But considering its low social and economic cost, this strategy is still a very cost-effective option.

Secondly, the cases in Wuhan and London has demonstrated high effectiveness of suppression strategy to reduce the overall infection. Probably in the UK or similar countries, suppression will minimally require a combination of social distancing of the entire population, home isolation of cases and household quarantine of their family members. But its practical effectiveness is not possible to achieve as same as Wuhan, as the success of suppression strategy in Wuhan is based on locking down human mobility to community level and sufficient resource support from other cities or provinces in China. If there are no sufficient external support, it will be risky to take intensive suppression to entire country due to huge impacts on its economics. Also, taking such intensive intervention to control an epidemic will need to be maintained until vaccine released (up to 12 months or more). If intensive interventions are relaxed at any time points, the transmission will quickly rebound. This is more like a multi-modal curve when taking multi-intervention strategies.

Thirdly, we also find out that while COVID-19 is estimated as a high production rate (R = 3), experimental evaluation results show that in either Wuhan or London cases fitted with real data in the last 6 weeks, high percentage of exposed or infectious population (at least 42.% of infectious population in Wuhan) are actually self-recovered. These people may have no or mild symptoms but been not checked as confirmed cases. This is one important issue that Zhong’s SERI model [11] has ignored. It will answer the model [8] predicts practical infectious population in Wuhan that ten times over figures in [11]. Similarly, it could explain the estimation of practical morality ratio can be varied in [11] and [8].

Lastly, one limitation of our model is that its prediction of infections and deaths depends on a parameter estimation of intervention intensity that presented by average-number contacts with susceptible individuals as infectious individuals in a certain region. We assume that each intervention has the same effect on the reproduction number in different regions over time. The practical effectiveness of implementing intervention intensity might be varied with respect to cultures or many other issues of certain county. As implementing hybrid intervention, the policy needs to be very specific and well-estimated at each day according to the number of confirmed cases, deaths, morality ratio, health resources, etc.

## 5 CONCLUSION

This paper conducts a feasibility study by defining a mathematical model named SEMCR that analyses and compares mitigation and suppression intervention strategies for controlling COVID-19 outbreaks in London and Wuhan Cases. The model was fitted and evaluated with public dataset containing daily number of confirmed active cases including Wuhan, London, Hubei province and the UK. The experimental findings show that the optimal timing of interventions differs between suppression and mitigation strategies, as well as depending on the definition of optimal. In future, our model could be extended to investigate how to optimise the timing and strength of intervention to reduce COVID-19 morality and healthcare demand in mobile application.

## Data Availability

The date in this manuscript are from public dataset, which has been referenced in the manuscript.