 Research article
 Open Access
 Published:
Modeling of suppression and mitigation interventions in the COVID19 epidemics
BMC Public Health volume 21, Article number: 723 (2021)
Abstract
Background
The global spread of the COVID19 pandemic has become the most fundamental threat to human health. In the absence of vaccines and effective therapeutical solutions, nonpharmaceutic intervention has become a major way for controlling the epidemic. Gentle mitigation interventions are able to slow down the epidemic but not to halt it well. While strict suppression interventions are efficient for controlling the epidemic, longterm measures are likely to have negative impacts on economics and people’s daily live. Hence, dynamically balancing suppression and mitigation interventions plays a fundamental role in manipulating the epidemic curve.
Methods
We collected data of the number of infections for several countries during the COVID19 pandemics and found a clear phenomenon of periodic waves of infection. Based on the observation, by connecting the infection level with the medical resources and a tolerance parameter, we propose a mathematical model to understand impacts of combining intervention measures on the epidemic dynamics.
Results
Depending on the parameters of the medical resources, tolerance level, and the starting time of interventions, the combined intervention measure dynamically changes with the infection level, resulting in a periodic wave of infections controlled below an accepted level. The study reveals that, (a) with an immediate, strict suppression, the numbers of infections and deaths are well controlled with a significant reduction in a very short time period; (b) an appropriate, dynamical combination of suppression and mitigation may find a feasible way in reducing the impacts of epidemic on people’s live and economics.
Conclusions
While the assumption of interventions deployed with a cycle of period in the model is limited and unrealistic, the phenomenon of periodic waves of infections in reality is captured by our model. These results provide helpful insights for policymakers to dynamically deploy an appropriate intervention strategy to effectively battle against the COVID19.
Background
The COVID19 pandemic has become a major global threat for human lives. In the absence of vaccines, effective treatment, and with limited knowledge of the virus [1–4], nonpharmaceutic interventions have been adopted to slow down the disease propagation. Nations around the world implemented a number of containment policies aimed at mitigating the epidemic. With the progression of the epidemic, individuals have improved their awareness of infection and changed their behavior to reduce their risk of infection by wearing facemasks and washing hands frequently [5]. Policies such as lockdown of the city [6, 7], travel restrictions [8–10], school closure [11], quarantine [12, 13] or stayathome [14], social distancing [15–19], bans on gatherings of more than a number of people, tracking individuals who are potentially infected [20, 21] have been implemented to reduce the contact rate and halt the epidemic. Some countries like Singapore used contact tracing to efficiently slow down the epidemic, while other countries such as the UK, opted to herd immunity and then changed to a strict lockdown.
Policymakers are confronted with difficult choices for controlling the epidemic. On one hand, strict measures on suppressing the epidemic can save people’s lives from deaths, while likely increasing the risk of economical losses; on the other hand, gentle mitigation interventions can reduce negative economical impacts but risk people’s life. Hence, it is necessary to estimate the epidemic dynamics in order to implement efficient, economical interventions accordingly [22, 23]. Many studies have developed mathematical models to evaluate the role of restriction measures on the dynamics of the COVID19 pandemic [13, 23–28]. Most of previous studies focused on modeling and estimating the transmission rate β to reflect the impact of interventions [20, 29]. For instance, a twostep control strategy relating suppression and mitigation was proposed [13, 23] based on modeling. Optimal control of the COVID19 from the point view of economics was also studied in [25]. The impacts of several interventions including case isolation, voluntary home quarantine, social distancing, and closure of schools have been explored in Ref. [26]. It shows that the most effective combination of interventions in reducing deaths and intensive care unit (ICU) capacity is a combination of case isolation, home quarantine, and social distancing. In a recent study [30], an adaptive intervention strategy was proposed for the control of the COVID19 epidemic by considering individuals’ mobility between cities. The authors assumed that intervention measures change the transmission rate with a controlled parameter and the floating population. In [24], strategies of intervention combination of suppression and mitigation were compared among 16 countries, while the cycle of the intervention is set at a fixed period. Although intensive interventions can effectively reduce the transmission rate, the epidemic may rebound if interventions are relaxed [26].
Since the epidemic is highly dynamic, rapidly changing with the increase of infection cases, appropriate interventions should be responsive to the outbreak and change with the epidemic dynamics accordingly. To capture this aspect, in this work, we explore a standard epidemiological model modified by considering two types of interventions, i.e., suppression and mitigation, which are dynamically implemented by considering factors, such as the city’s medical resources and its tolerance for infection. We consider that mitigation interventions are a combination of measures, such as general social distancing measures, hygiene rules, case isolation, shielding of vulnerable groups, school closures, etc, targeting the basic reproduction number approaching but larger than 1, while suppression interventions are taken as additional measures of strict physical distancing, including lockdowns, targeting the basic reproduction number less than 1. We propose a combined strategy of suppression and mitigation to control the disease propagation. Instead of setting the cycle of the strategies as a fixed period, the system is allowed to dynamically adjust the control strategy depending on the infection level. By doing so, the epidemic is under control below an acceptable level and a strict suppression intervention is not necessarily deployed during the period.
Methods
In order to understand the impact of different interventions on controlling the epidemic, we firstly collect the data from data source DXY [31] and observe the curves of existing infections in different countries, as shown in Fig. 1. We clearly see that in some countries, the number of existing infections shows a typical, nonuniform periodic wave with different peaks and periods. For instance, in Thailand (Fig. 1a), the number of existing infections reaches the peak around 1500 cases on April, 1st, 2020, and decreases to 100 cases soon before May, 25th, 2020. Then, it follows a period wave of infections less than 200 cases. Recently, it increases to 500 cases again. Also in Qatar (Fig. 1b), during the period from July 12th, 2020 to Sep. 11th, 2020, the number of infections was controlled less than 4000 cases and then it decreases again from Nov. 1st, 2020. Similar phenomena can also be found in Tajikstan and Vietnam (Fig. 1c and d). Although the periods of the waves and the peak values of infections may be different, such a phenomenon of period waves indeed reflects various interventions deployed by nations to battle against the COVID19 pandemic.
In the following, we perform a modeling study to understand interventions for the appearance of periodic waves of infections. Since the COVID19 can cause infections with no symptoms [32] and severely infected individuals have to be hospitalized for treatment, we model them in a population based on the classical susceptibleexposedinfectedrecovered (SEIR) epidemiological model, where individuals belong to eight states: susceptible (S), exposed (E), asymptomatic (A), symptomatic (I), severe infection (P), hospitalized (H), dead (D), and recovered (R). See Fig. 2. The total population size is denoted as N=S+E+A+I+P+H+D+R. At time t, susceptible individuals get infected by having contacts (asymptomatic, symptomatic or severe) infectious individuals and move to exposed state (E) at a rate β(t), where β(t) varies with time due to the involvment of interventions. The infectiousness of asymptomatic and severe infected individuals compared with symptomatic infected individuals are adjusted by factors θ_{a} and θ_{p}, respectively. A fraction of exposed individuals (f_{I}) moves to symptomatic infectious state (I) at rate δ, while the remainder of them (1−f_{I}) moves to asymptomatic infectious class (A) at the same rate δ. Thus, the \(\frac {1}{\delta }\) refers to as incubation period. The infectious periods of symptomatic and asymptomatic infections are denoted as \(\frac {1}{\gamma _{I}}\) and \(\frac {1}{\gamma _{A}}\), respectively. A fraction of symptomatic infections (f_{P}) develops to be severely infected by moving to severe infection (P) at rate γ_{I}, while the remainder of symptomatic infections (1−f_{P}) move to recovered state (R) at rate γ_{I}. Severe infections will be hospitalized for treatment at rate ω. Thus, ω can be taken as detection rate for infections. A fraction of patients (f_{D}) die at rate μ and the remainder of them recover. The dynamics of the model is given by,
where N is the total population size.
The basic reproduction number R_{0} is usually estimated in the early stage of the outbreak as a constant value. With the introduction of interventions, the transmission rate β is timevarying. Consequently, the reproduction number changes with time, often named as effective reproduciton number and expressed as R_{t}. By linearing the system (Eq. (1)) at the disease free state (S,E,A,I,P,H,D,R)=(S_{0},0,0,0,0,0,0,0) and setting the vector v=(E,A,I,P)^{T}, we obtain \(\dot {v} = (\mathcal {F}\mathcal {V})v\), where \(\mathcal {F}\) is the new infection rate in each class and \(\mathcal {V}\) is the transition rate for each class by transferring in or out of each class, given by
from which we can obtain the effective reproduction number of the model (Eq. (1)) as the spectral radius \(\rho (\mathcal {F}\mathcal {V}^{1})\) of the next generation matrix \(\mathcal {F}\mathcal {V}^{1}\). The effective reproduction number R_{t} is given by
Specifically, the basic reproduction number is recovered as \(\mathbf {R}_{0}=\beta _{0} \left (\frac {\theta _{a} (1f_{I})}{\gamma _{A}}+\frac {f_{I}}{\gamma _{I}}+\frac {\theta _{p} f_{I} f_{P}}{\omega }\right)\), where β_{0} is the transmission rate in the early stage. The difference between interventions of suppression and mitigation is the aim for affecting the epidemic. Suppression aims at halting the epidemic with extremely strict strategies to satisfy the condition R_{t}<1, while mitigation aims at slowing down the epidemic with relaxed control strategies, resulting in a reduced R_{t} larger than 1. A strict suppression is more likely to have negative impacts on economics and social lives, while a gentle mitigation may not be able to efficiently control the epidemic. Thus, how to implement an effective intervention strategy to balance the aims is an important and difficult issue to answer.
In order to devise efficient intervention strategies to control epidemics, one has to consider medical resources a city possesses, since infected cases require medical resources for treatment. Also, the timing of interventions being deployed plays a key role in containing epidemic. A delayed, relaxed intervention strategy might be uneffective for controlling the epidemic, since the accumulated existing hospitalized infections may overload the medical resources and further promote the spread of epidemic. Thus, at this stage, a strict suppression strategy is preferentially deployed. This is exactly the strategy some countries like UK adopted. On the other hand, if the existing hospitalized infections are not at risk for medical resources, a gentle, relaxed mitigation strategy would be sufficient enough to control the epidemic.
In the following, we assume that intervention strategies deployed by a city are closely related with its capacity of medical resources I_{m} and a tolerance level for epidemic I_{t}. At time t_{0}, if the total hospitalized cases during a period τ are higher than the capacity of medical resources I_{m} with factor c (0<c≤1), a strict suppression strategy will be implemented; otherwise, if the number of hospitalized cases is higher than some tolerance parameter I_{t}, a mitigation intervention will be deployed. Then, the effective reproduction number R_{t} is given by
where α_{s} and α_{m} are the mean intensity values associating with suppression and mitigation, respectively, followed by a given distribution [29]. Here, for simplicity, we assume that both α_{s} and α_{m} are constants satisfying α_{s}<α_{m}<1. The parameter I_{m} represents the capacity of medical resources, which may depend on the economical level. For instance, in the US, I_{m}=1.2 beds per thousand people. The parameter I_{t} represents the tolerance level for epidemic, e.g., \(I_{t}=\frac {100}{N}\) means that the tolerance of 100 infections in the population. The relationship between parameters I_{m} and I_{t} is tuned by the factor c. A value of c satisfying cI_{m}>I_{t} can capture the condition of suppression and mitigation, i.e., α_{s}<α_{m}. Given the parameters for controlling strength, i.e., α_{s} and α_{m}, and the starting time t_{0}, the role of interventions is captured by the evolving curve of R_{t} or β(t).
It is to note that the deployment of an intervention depends on the number of recently hospitalized cases \(\sum _{t(\tau 1)}^{t} H(t)\), which is related with the accumulated confirmed cases \(\sum _{t(\tau 1)}^{t} I(t)\). In reality, evaluating the exact value of I(t) is fundamental to understand the impacts of interventions on the epidemic. For the COVID19, since the real number of infections is generally unknown and only the number of confirmed cases is available, it is necessary to use some approaches to evaluate [29, 33]. Here, in the following, for the real data analysis, we apply a deconvolution method proposed in [29] to estimate the real number of infection cases from the reported infection cases as the initial values of the model. Then, we study the impact of interventions of suppression and mitigation on the epidemic.
Results
Modeling results
We first perform simulations on the proposed model in a host population with no interventions. Then, we compare the effects of interventions such as suppression, mitigation, and the combination of them with the proposed model. The total population is set as N=10^{6} and the fractions of exposed individuals and symptomatic infected individuals are E(0)=6×10^{−4} and I(0)=4×10^{−4}, respectively. The basic reproduction number is set as R_{0}=2.2 for the COVID19 as Ref. [34]. The incubation period is \(\frac {1}{\delta }=5.1\) days [35]. The fraction of symptomatic infections is f_{I}=0.82 as calculated in Ref. [36]. The infectious periods for symptomatic and asymptomatic infection are \(\frac {1}{\gamma _{I}}=\frac {1}{\gamma _{A}}=2.3\) days [24]. Since it has been found that there is no difference in the transmission rates between symptomatic and asymptomatic patients [32], we set θ_{a}=1. Usually, the infectiousness of severe infected individuals is higher than that of symptomatic infectious individuals. Here, we assume θ_{p}=1 for simplicity. The fraction of infections requiring hospitalization is set as f_{P}=0.3 and the delay between severe infection and hospitalization is \(\frac {1}{\omega }\)=2.7 days [24] and the hospital stay is assumed as \(\frac {1}{\mu }=8\) days. According to the number of daily deaths and the existing infections in WuHan, the fraction of hospitalized infections to die is calculated as f_{D}=0.065. Without specification, the factor constant c is set as c=0.8. We did simulations with several values of c and found that the choice of c does not change the main results of the present study. Most of the parameters are chosen according to the recently published results and some of them are assumed, as summarized in Table 1. Figure 3 illustrates the results for the baseline model, where no interventions are involved. We observe that the infected and hospitalized number reaches the peak at around 50 and 60 days, respectively, approaching 4.5% and 3% of the population. The fraction of deaths reaches 0.8% of the population, approximately 8000 people.
Effects of suppression intensity α _{s} and the starting time t _{0}
In this scenario, we assume that if the fraction of accumulated hospitalization cases during the past τ days, e.g., τ=7, \(\sum _{t7}^{t}H(t)\), is larger than the medical resources by a factor, cI_{m}, a suppression intervention will be implemented. We take it as a constant control intensity with the average value α_{s}=0.3,0.1 starting at time t_{0}=10,20,30. These parameters are chosen such that the effective reproduction number R_{t} is approximatively less than 1. In Fig. 4, we see that with the introduction of a suppression measure with α_{s}=0.3, the peak of infections dramatically reduces to a lower level. For instance, the epidemic peak is reduced by 4.5% to 0.7% at t_{0}=30. The earlier the intervention is deployed, the lower the peak of infections will be. If the intervention is deployed 20 days earlier, the infection could be further reduced less than 0.15%, see Fig. 4c. Earlier deployment of suppression strategies flatted the curve of dynamics with a significantly reduced peak value and the fraction of hospitalized cases remain at a lower level. The fraction of deaths is also reduced to a lower level due to the extremely efficient intervention of suppression. For instance, at t_{0}=20, the fraction of deaths is controlled less than 1‰ of the population. In addition, since the introduction of suppression intervention depends on the number of the accumulated hospitalized cases, both the existing infections and the hospitalized infections show a dynamical period of peaks. Earlier deployment of the intervention further brings a more frequent change.
With a more intensive suppression (α_{s}=0.1), the sharp reductions in infections and hospitalizations narrow the curve of the dynamics. In addition, a timely intensive suppression can reduce the deaths to 0.5‰ of the population at t_{0}=10. Even with a delayed intervention at t_{0}=30, the deaths can be reduced to approximately 1‰ of the population and the curve of hospitalizations decreases faster. Therefore, to save lives from epidemic, an efficient, strict suppression should be deployed as early as possible. The decreasing slopes of the epidemic curves also provide insightful information for evaluating the effect of different control measures.
Effect of mitigation intensity α _{m} and the starting time t _{0}
Since a suppression intervention is likely to bring negative effects on economics and social activities, a gentle mitigation intervention is feasible to gradually reduce the basic reproduction number R_{0}. If the ratio of the accumulated existing hospitalized cases during the past week is larger than some tolerance parameter I_{t}, the mitigation intervention is implemented. To compare different mitigation measures, we consider α_{m}=0.7,0.5, evaluating R_{t}=1.54,1.1∈[1,R_{0}], respectively. In Fig. 5, we see that with a very gentle mitigation α_{m}=0.7, the peak of infections reduces to approximately half of that with no interventions, and the ratio of deaths reduces to 6‰ of the population. Due to the limited role of mitigation intervention in controlling the epidemic, the timing of the intervention does not take obvious effects on the peaks of both infections and deaths. In addition, comparing with the sharp reduction in the curve of the epidemic under suppression intervention, the curve of the epidemic under mitigation further flattens and lasts for longer time. The dynamical period of the peaks of infections is not observed due to the gentle mitigation in slowing down the epidemic.
A more intensive mitigation intervention with α_{m}=0.5 further flatterns the curve of epidemic dynamics with a reduced peak value. The ratio of deaths is also dramatically reduced to less than half of that with no interventions, even starting at late time t_{0}=30. With a more intensive mitigation, the starting time of the mitigation measure plays a key role in changing the curve of the epidemic dynamics. Earlier interventions can efficiently suppress the epidemic spread and reduce the deaths.
A combined intervention of suppression and mitigation
From the above analysis, we find that with an appropriate suppression intervention, the epidemic can be efficiently controlled in short time, while with a gentle mitigation intervention, the epidemic can be slowed down to some degree. To further understand the effects of the two types of interventions on the epidemic, we investigate a strategy of combined intervention of suppression and mitigation. Instead of a manual adjustment of the cycle of the intervention period, we propose a dynamical intervention measure composed of suppression and mitigation depending on the infection level. If the accumulated number of the existing hospitalization cases during the last week is larger than the capacity of medical resources by a factor, then a strict suppression is deployed to suppress the epidemic spread in order to avoid overloading the medical resources. If it is larger than some accepted value, determined by policymakers, then a gentle, relaxed mitigation intervention is implemented. By doing so, the system can dynamically adjust the interventions depending on the existing hospitalization cases. Based on this assumption, we propose a strategy of combined intervention by setting parameters as α_{s}=0.3 (R_{t}=0.66<1) and α_{m}=0.7 (R_{t}=1.54>1) starting at different time t_{0}.
From Fig. 6, we see that when the intervention is deployed at delayed time t_{0}=30, due to the increased accumulated number of hospitalized infections, a suppression measure lasting for more than one and a half months is necessary in order to keep the epidemic under control, after which a relaxed mitigation measure can follow. Since the deployment of a relaxed mitigation measure may cause the increase of hospitalized infections again, then a strict suppression can be further imposed. Such a periodic iteration of suppression and mitigation proceeds with the hospitalized infections being controlled under a lower level. The period wave for the curve of epidemic dynamics is determined by the intervention time t_{0}, the medical resources I_{m} and the tolerance parameter I_{t}. A timely intervention switching between suppression and mitigation will lead to a significant reduction in both the hospitalized infections and deaths.
Our analysis is different from previous study [13], where a period of lockdown or quarantine is fixed manually. By dynamically alternating the intervention with the epidemic dynamics, it is possible to keep the infection under control while allowing a sustainable economy as well as normally social activities.
Real world data analysis
In the following, we explore the impact of intervention measures on epidemic by analyzing the data of Wuhan for the COVID19 epidemic. We are aware that the actual infection number is unknown and can only be inferred from other epidemiological observations (e.g., the daily confirmed cases). Such observations are lagging behind the infection events due to inevitable time delays between an individual being infected and reported (e.g., days for symptom onset). In Ref. [29], a deconvolution method is proposed to estimate the actual infection cases from the daily reported confirmed cases with the renewal process. Detailed descriptions of the method refer to Ref. [29]. In this work, we applied this method to infer the actual infection cases on Jan. 11th, 2020. The daily number of onset patients in Wuhan (i.e. daily confirmed cases) is obtained from the study by Pan et al. [38]. Then, we evaluate the infection cases on Jan. 11th as 284. Combined with the information of the ratio of symptomatic infection over asymptomatic infection \(\frac {f_{I}}{1f_{I}}\), the number of asymptomatic infected cases is evaluated as A(0)=62. By assuming that the incubation period is 5.1 days, we can get the number of exposed cases on Jan. 11th from the number of onset patients on Jan. 16th as E(0)=359. The numbers of recovered cases and death cases on Jan. 11th, 2020 are obtained from [31]. The total population of Wuhan is N=11081643. The basic reproduction number is estimated as R_{0}=2.2 [34] and accordingly, the transmission rate is calculated as β_{0}=0.7422. With the available data of beds number in Wuhan, the capacity of medical resource of Wuhan is calculated as I_{m}=0.0023 [37]. Other parameters are the same as in Table 1.
Figure 7a shows the epidemic dynamics with no intervention. We see that with no interventions being introduced, it will cause more than ninety thousand deaths in Wuhan. While in reality, as a part of the national emergency response, public transport were suspended, public gathering was banned in Wuhan [10, 39]. The measures taken in Wuhan dramatically reduced the death number less than four thousand by April 17th, 2020 and the peak of infections around thirtytwo thousand on Feb. 19th, 2020, as shown in Fig. 7b. It is to note that measures such as lockdown of Wuhan and the construction of mobile cabin hospitals play a key role in suppressing the epidemic, which reduce the reproduction number R_{t}<1 in very short time. To compare the model prediction with the real case, we also perform simulations on the proposed model with the real data. By setting t_{0}=12 on the lockdown day of Wuhan and α_{s}=0.3 in Fig. 7c, we see that the extremely strict intervention will immediately reduce the reproduction number R_{t} less than 1 and the deaths will reduce to less than two hundred and the arrival of the peak value will be one month earlier than expected. In reality, the intervention intensity α_{s} depends on several factors, such as medical resources and the construction of mobile cabin hospitals. The difficulty of obtaining these resources may hinder the control of the epidemic.
Next, we investigate what would happen if dynamical intervention were deployed in Wuhan, as shown in Fig. 8. Depending on the intervention time t_{0}, alternative interventions are dynamically deployed. For instance, with a delayed time t_{0}=40 (Fig. 8a), the number of infections will increase more than fifteen thousands, consequently, a strict suppression has to be deployed for more than 40 days in order to efficiently control the epidemic, after which a gentle mitigation and a strict suppression alternatively follow. If the intervention start by 10 days earlier (Fig. 8b), the peak value of infections can be reduced to less than fifteen thousands. Then, a shorter suppression followed by a shorter mitigation alternatively controls the epidemic. Finally, if the intervention were deployed at early time with t_{0}=12 (Fig. 8c), the arrival of the peak value of infections will be delayed and the epidemic is completely controlled periodically.
In order to understand the impacts of the detection rate on the epidemic, we carried out simulations with Wuhan data for different choices of ω as \(\omega =\frac {1}{2.7}, \frac {1}{8.1}\), respectively, in Fig. 9. With no interventions, we see that it will cause less than two hundred thousand severe infections and more than three hundred thousand hospitalized infections (Fig. 9a). With the delay of detection, e.g., \(\omega =\frac {1}{8.1}\), the number of severe infections is doubled while the number of the hospitalized infections is reduced. In addition, the arrival of the peak infection will be delayed for one month, as shown in Fig. 9d.
Next, we observe the impacts of the detection rate ω on the epidemic dynamics in Wuhan when the combined intervention strategy is involved at t_{0}=12, as shown in Fig. 9b and e. With a higher detection rate (\(\omega =\frac {1}{2.7}\)), the number of severe infections reaches less than five thousand and shows a periodical wave under control (Fig. 9b). In such a situation, the intervention measure of suppression frequently alternates with mitigation (Fig. 9c). With the decrease of the detection rate ω (e.g., \(\omega =\frac {1}{8.1}\)), severe infections are hard to be detected and their number is doubled to ten thousand (Fig. 9e). Since the hospitalized infections are decreased, consequently, a longer period of suppression and mitigation is necessary in order to keep the epidemic under control (Fig. 9f).
Discussion
So far, although several vaccines for the COVID19 are ready to deploy, nonpharmaceutical interventions are still effective actions to control the COVID19 pandemic. Therefore, how to implement appropriate interventions for the next stage is fundamentally important for the control of epidemic while sustaining normally social life. It is interesting to point out that extremely strict interventions, such as lockdown and quarantine measures, have substantial effects on changing the epidemic dynamics. With a strict suppression strategy, the number of deaths and that of infections can be reduced to a lower level while longterm measures may have negative impacts on economics and social life.
To solve such a dilemma, we propose a strategy of combined intervention of suppression and mitigation, which dynamically alternates according to the epidemic dynamics. The deployment of suppression interventions is assumed to be relevant with the capacity of medical resources, that is, if the accumulated number of hospitalized infections during a given period is close to the capacity of medical resources with a factor, then a strict suppression will be implemented to avoid overloading medical resources. If the number of hospitalized infections is more than some tolerance level, a relaxed mitigation, such as social distancing and hygiene measures, may be sufficient enough to control the epidemic. Depending on the tolerance parameter, the mitigation strategy can be dynamically switched on or off. Our study shows that early deployment of a suppression intervention will shift the peak value to an earlier date. While timely implement of a mitigation intervention will flatten the epidemic curve with a prolonged period. With a strategy of combined intervention of suppression and mitigation, the epidemic is contained with an acceptable level, where the two measures alternatively interchange with different periods. Such a dynamical deployment of interventions is able to keep the tradeoff between economics and epidemic, and take less negative impacts on social life. We believe that such an approach may be adopted as a strategy until the available of effective vaccine.
Our study has certain limitations as well. In reality, the periodic waves of infections in different nations show diverse features, e.g. peaks and periods, depending on how the prevention measures are deployed. Using a simplified Heaviside function alternating between two interventions is not sufficient enough to reflect such a complicated situation in reality. Instead, a nonlinear function is expected to offer a solution. Consequently, the present analysis and interpretation of the results are limited to a modeling study for such a phenomenon of periodic waves of infections.
Conclusions
By analysing the data for countries in battling against the COVID19 pandemic, we found a periodiclike wave in the number of infections in several countries, which indicates a clear relationship between the infection and the deployed interventions when facing the COVID19 epidemic. The present study explores the combined impacts of suppression and mitigation measures on the epidemic. The findings are as follows: (a) depending on a city’s capacity of medical resources and its tolerance on infectious population, a combined intervention of suppression and mitigation can efficiently reduce the peak of infections and negative effects on social lives and economics; (b) an immediate, strict suppression measure is highly efficient in reducing the number of infections and brings less losses on people’s lives, as occured in Wuhan, China; (c) a delayed intervention has to be accompanied by a longer suppression followed by a shorter mitigation alternatively, which is expected to balance the loss of people’s lives and that of economics for a longterm control of the COVID19.
Availability of data and materials
The number of existing confirmed cases in some countries is available in the DXYCOVID19Data Repository (https://github.com/BlankerL/DXYCOVID19Data). The data set about daily number of onset patients in Wuhan is available in the https://github.com/nianqiaoju/rt.
Declarations
Abbreviations
 COVID19:

Corona Virus Disease 2019
 R _{0} :

Basic reproduction number
 SEIR model:

SusceptibleExposedInfectedRecovered model
References
He D, Zhao S, Xu X, Lin Q, Zhuang Z, Cao P, Wang M, Lou Y, Xiao L, Wu Y. Low dispersion in the infectiousness of covid19 cases implies difficulty in control. BMC Public Health. 2020; 20(1):1558.
Imai N, Cori A, Dorigatti I, Baguelin M, Donnelly CA, Riley S, Ferguson NM. Report 3: Transmissibility of 2019nCoV. Imp Coll Lond. 2020. https://www.imperial.ac.uk/mrcglobalinfectiousdiseaseanalysis/covid19/report3transmissibilityofcovid19/. Accessed 25 Jan 2020.
Lin Q, Zhao S, Gao D, Lou Y, Yang S, S.Musa S, H.Wang M. A conceptual model for the coronavirus disease 2019(covid19) outbreak in wuhan, china with individual reaction and governmental action. Int J Infect Dis. 2020; 93:211–6.
Read JM, Bridgen JRE, Cummings DAT, Ho A, Jewell CP. Novel coronavirus 2019ncov: early estimation of epidemiological parameters and epidemic predictions. MedRxiv 2020:2020.01.23.20018549. 2020. https://doi.org/10.1101/2020.01.23.20018549.
Eikenberry SE, Mancuso M, Iboi E. To mask or not to mask: modeling the potential for face mask use by the general public to curtail the covid19 pandemic. Infect Dis Model. 2020; 5:293–308.
Carcione JM, Santos J. A simulation of a covid19 epidemic based on a deterministic seir model. Front Public Health. 2020; 8:230.
Shi WZ, Tong CZ, Zhang AS, Wang B, Shi ZC, Yao YP, Jia P. An extended weight kernel density estimation model forecasts covid19 onset risk and identifies spatiotemporal variations of lockdown effects in china. Commun Biol. 2021; 4(1):126.
Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, Pastore y Piontti A, Mu K, Rossi L, Sun K, Viboud C, Xiong X, Yu H, Halloran ME, Longini IM, Vespignani A. The effect of travel restrictions on the spread of the 2019 novel coronavirus(covid19) outbreak. Science. 2020; 368:395–400.
Jia JS, Lu X, Yuan Y, Xu G, Jia J, Christakis NA. Population flow drives spatiotemporal distribution of covid19 in china. Nature. 2020; 582:389–94.
Tian H, Liu Y, Li Y, Wu CH, Chen B, Kraemer MUG, Li B, Xu B, Yang Q, Wang BEA. An investigation of transmission control measures during the first 50 days of the covid19 epidemic in china. Science. 2020; 368:638–42.
Lee B, Hanley JP, Nowak S, Bates JH, HébertDufresne L. Modeling the impact of school reopening on sarscov2 transmission using contact structure data from shanghai. BMC Public Health. 2020; 20:1713.
Li D, Liu Z, Liu Q. Estimating the efficacy of quarantine and traffic blockage for the epidemic caused by 2019ncov: a simulation analysis. MedRxiv 2020:2020.02.14.20022913. https://doi.org/10.1101/2020.02.14.20022913.
PérezGarcía VM. Relaxing quarantine after an epidemic: A mathematical study of the spanish covid19 case. 2020. https://doi.org/10.13140/RG.2.2.36674.73929/1.
Medline A, Hayes L, Valdez K, Hayashi A, Vahedi F, Capell W, Sonnenberg J, Glick Z, Klausner JD. Evaluating the impact of stayathome orders on the time to reach the peak burden of covid19 cases and deaths: does timing matter?BMC Public Health. 2020; 20:1750.
Chen X, Zhang A, Wang H, Gallher A, Zhu X. Compliance and containment in social distancing: mathematical modeling of covid19 across townships. Int J Geogr Inf Sci. 2021; 35(3):446–65.
Du Z, Xu X, Wang L, Fox S, Cowling B, Galvani A, Meyers L. Effects of proactive social distancing on covid19 outbreaks in 58 cities, china. Emerg Infect Dis. 2020; 26:2267–9.
Falco D, Cioppa AD, Scafuri U, Tarantino E. Coronavirus covid19 spreading in italy: optimizing an epidemiological model with dynamic social distancing through differential evolution. arXiv preprint arXiv:2004.00553. 2020.
Lee SW, Yuh WT, Yang JM, Cho YS, Yoo IK, Koh HY, Marshall D, Oh D, ha EK, Han MY, Yon DK. Nationwide results of covid19 contact tracing in south korea: Individual participant data from an epidemiological survey. JMIR Med Inform. 2020; 8(8):e20992.
Price GN. The effect of social distancing on the spread of novel coronavirus: estimates from linked statelevel infection and american time use survey data. https://doi.org/10.13140/RG.2.2.30632.32007.
Ferretti L, Wymant C, Kendall M, Zhao L, Nurtay A, AbelerDorner L, Parker M, Bonsall D, Fraser C. Quantifying sarscov2 transmission suggests epidemic control with digital contact tracing. Science. 2020; 368(6491). https://doi.org/10.1126/science.abb6936.
Walker PGT, Whittaker C, Watson O. The global impact of covid19 and strategies for mitigation and suppression. Science. 2020; 369:413–22.
Gros C, Valenti R, Valenti K, Gros D. Strategies for controlling the medical and socioeconomic costs of the corona pandemic. arXiv preprint arXiv:2004.00493. 2020.
Hochberg ME. Importance of suppression and mitigation measures in managing covid19 outbreaks. arXiv preprint arXiv: 2005.03323. 2020.
Chowdhury R, Heng K, Shawon MSR, Goh G, Okonofua D, OchoaRosales C, GonzalezJaramillo V, Bhuiya A, Reidpath D, Prathapan S, Shahzad S, Althaus CL, GonzalezJaramillo N, Franco OH. Dynamic interventions to control covid19 pandemic: a multivariate prediction modelling study comparing 16 worldwide countries. Eur J Epidemiol. 2020; 35:389–99.
DjidjouDemasse R, Michalakis Y, Choisy M, Sofonea MT, Alizon S. Optimal covid19 epidemic control until vaccine deployment. MedRxiv 2020:04.02.20049189. https://doi.org/10.1101/2020.04.02.20049189.
Ferguson N, Laydon D, NedjatiGilani G. Impact of nonpharmaceutical interventions to reduce COVID19 mortality and healthcare demand. 2020. http://www.imperial.ac.uk/mrcglobalinfectiousdiseaseanalysis/covid19/report9impactofnpisoncovid19/. Accessed 16 Mar 2020.
Piccolomini EL, Zama F. Monitoring italian covid19 spread by an adaptive seird model. MedRxiv 2020:2020.04.03.20049734. https://doi.org/10.1101/2020.04.03.20049734.
Zu J, Li M, Li Z, Shen M, Xiao Y, Ji F. Epidemic Trend and Transmission Risk of SARSCov2 After Government Intervention in the Mainland of China: a Mathematical Model Study. 2020. https://doi.org/10.2139/ssrn.3539669.
Wang S, Yang X, Li L, Nadler P, Arcucci R, Huang Y, Teng ZZ, Guo YK. A bayesian updating scheme for pandemics: estimating the infection dynamics of covid19. IEEE Comput Intell Mag. 2020; 15:23.
Wu Z, Liao H, Vidmer A, Zhou M, Chen W. Covid19 plateau: a phenomenon of epidemic development under adaptive prevention strategies. arXiv preprint arXiv:2011.03376. 2020.
Blanker L. COVID19/2019nCoV Time Series Infection Data Warehouse. https://github.com/BlankerL/DXYCOVID19Data. Accessed 23 Jan 2021.
Yin G, Jin H. Comparison of transmissibility of coronavirus between symptomatic and asymptomatic patients: reanalysis of the ningbo covid19 data. JMIR Public Health Surveill. 2020; 6:19464.
Grassly NC, Fraser C. Mathematical models of infectious disease transmission. Nat Rev Microbiol. 2008; 6:477.
Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, Ren R, Leung K, Lau E, Wong J, Xing X, Xiang N, Wu Y, Li C, Chen Q, Li D, Liu T, Zhao J, Liu M, Tu W, Chen C, Jin L, Yang R, Wang Q, Zhou S, Wang R, Liu H, Luo Y, Liu Y, Shao G, Li H, Tao Z, Yang Y, Deng Z, Liu B, Ma Z, Zhang Y, Shi G, Lam T, Wu J, Gao G, Cowling B, Yang B, Leung G, Feng Z. Early transmission dynamics in wuhan, china, of novel coronavirusinfected pneumonia. N Engl J Med. 2020; 382:1199–207.
Lauer S, Grantz K, Bi Q, Jones F, Zheng Q, Meredith H, Azman A, Reich N, Lessler J. The incubation period of coronavirus disease 2019 (covid19) from publicly reported confirmed cases: Estimation and application. Ann Intern Med. 2020; 172:577–82.
Kenji M, Katsushi K, Alexander Z, Gerardo C. Estimating the asymptomatic proportion of coronavirus disease 2019 (covid19) cases on board the diamond princess cruise ship. Euro Surveillance. 2020; 25:1.
News P. How many beds are available in Wuhan to deal with COVID19?. 2020. https://baijiahao.baidu.com/s?id=1658385426681451867&wfr=spider&for=pc. Accessed 13 Feb 2020.
Pan A, Liu L, Wang C, Guo H, Hao X, Wang Q, Huang J, He N, Yu H, Lin X, Wei S, Wu T. Association of public health interventions with the epidemiology of the covid19 outbreak in wuhan, china. JAMA. 2020; 323:1915–23.
Zhang J, Litvinova M, Liang Y, Wang Y, Wang W, Zhao S, Wu Q, Merler S, Viboud C, Vespignani A, Ajelli M, Yu H. Changes in contact patterns shape the dynamics of the covid19 outbreak in china. Science. 2020; 368(6498):1481–6. https://doi.org/10.1126/science.abb8001.
Acknowledgments
Not applicable.
Funding
This work was supported by the National Key Research and Development Program of China under Grant No. 2017YFE0117500, the National Natural Science Foundation of China under Grant No. 91746203, and the Natural Science Foundation of Shanghai under Grant No. 20ZR1419000.
Author information
Authors and Affiliations
Contributions
YH: Conceptualization, Writing; ZX: Acquisition, Analysis; YG: Review and editing; BW: Methodology, Writing. All authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Han, Y., Xie, Z., Guo, Y. et al. Modeling of suppression and mitigation interventions in the COVID19 epidemics. BMC Public Health 21, 723 (2021). https://doi.org/10.1186/s12889021106636
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12889021106636
Keywords
 COVID19
 Basic reproduction number
 Interventions
 Suppression
 Mitigation