Skip to main content

Evaluation of China’s Hubei control strategy for COVID-19 epidemic: an observational study



To fight against COVID-19, many policymakers are wavering on stricter public health interventions. Examining the different strategies both in and out of China’s Hubei province, which contained the epidemic in late February 2020, could yield valuable guidance for the management of future pandemics. This study assessed the response process and estimated the time-varying effects of the Hubei control strategy. Analysis of these strategies provides insights for the design and implementation of future policy interventions.


We retrospectively compared the spread and control of COVID-19 between China’s Hubei (excluding Wuhan) and non-Hubei areas using data that includes case reports, human mobility, and public health interventions from 1 January to 29 February 2020. Static and dynamic risk assessment models were developed to statistically investigate the effects of the Hubei control strategy on the virus case growth after adjusting importation risk and policy response timing with the non-Hubei strategy as a control.


The analysis detected much higher but differential importation risk in Hubei. The response timing largely coincided with the importation risk in non-Hubei areas, but Hubei areas showed an opposite pattern. Rather than a specific intervention assessment, a comprehensive comparison showed that the Hubei control strategy implemented severe interventions characterized by unprecedentedly strict and ‘monitored’ self-quarantine at home, while the non-Hubei strategy included physical distancing measures to reduce contact among individuals within or between populations. In contrast with the non-Hubei control strategy, the Hubei strategy showed a much higher, non-linear and gradually diminishing protective effect with at least 3 times fewer cases.


A risk-based control strategy was crucial to the design of an effective response to the COVID-19 outbreak. Our study demonstrates that the stricter Hubei strategy achieves a stronger controlling effect compared to other strategies. These findings highlight the health benefits and policy impacts of precise and differentiated strategies informed by constant monitoring of outbreak risk.

Peer Review reports


The world is combating the ongoing COVID-19 pandemic, for which potential therapeutics and vaccines are still being investigated. Accordingly, great expectations are placed on non-pharmaceutical public health interventions to contain the epidemic [1, 2]. This has fueled interest in exploring their effectiveness in epidemic control. So far, the potential effects of anti-contagion policies were usually estimated and reported through process-based epidemiological simulations [3,4,5,6,7,8]. A limited number of control studies based on current and past epidemics focused on the static evaluation of individual or partial interventions [2, 9]. These approaches suffer from a few limitations. First, the information about how long policies should be maintained is obscure or unavailable. Second, the transmission risk in policy scenarios and the interactions between interventions in policy packages are neglected. Therefore, current studies cannot generate an efficient assessment of the time-varying effects of the interventions.

In response to the outbreak of COVID-19 in December 2019, a series of public health interventions has been employed by China’s national, provincial, or municipal governments effectively curbing the epidemic by the end of February. These interventions were triggered by launching the public health emergency response [10, 11]. Besides Wuhan, other cities in Hubei province were recognized as key areas of epidemic growth with city lockdown as milestones for COVID-19 control. Different control strategies were implemented inside and outside of Hubei. This concrete difference in policy interventions ultimately led to different levels of effectiveness. A careful and comprehensive comparison of these two control strategies may help to identify whether or when these interventions should be deployed, consolidated, or relaxed. Such guidance would fill an urgent need in the governmental decision-making process, especially as COVID-19 is still running rampant and policymakers in many countries are wavering on stricter interventions.

Population movement from Wuhan during the 2020 Chinese New Year mass migration constituted an importation risk for the travelers’ destinations [12,13,14]. In this study, we retrospectively compared the spread and control of COVID-19 in and out of Hubei using data on case reports, human mobility, and public health interventions during the period from 1 January to 29 February 2020. We first characterized the association between the importation risk of COVID-19 and the policy-making process in each study prefecture. Then we documented the difference in control strategies and measures in and out of Hubei province across mainland China. Finally, we developed static and dynamic models to quantify the time-varying effects of Hubei and non-Hubei control strategies.


Data collection

To explore the role of importation risk on control strategy decisions across mainland China and ascertain their impacts, we collected a wide array of data. These data sets included virus case data, control measures implementation details, potentially related population and economic data, and travel data from 16 cities in Hubei province (excluding Wuhan) and 30 provinces outside Hubei. They originated from 46 prefectures in total and ranged in time from the beginning of the COVID-19 epidemic to the time of its fundamental containment (1 January to 29 February 2020). The first three types of data were extracted from local official websites (cases, control measures, and economic/population data). Travel data were retrieved through Baidu Qianxi platform, which derives its data from Location-Based Services. Specifically, Daily Baidu Mobility Indexes (dBMIs) of population outflow from Wuhan to each study prefecture were obtained.

We calculated the total aggregate population outflow from Wuhan between 1 to 26 January 2020 to track movements from Wuhan to each study prefecture before its closure on 23 January 2020 and measure their risks of COVID-19 importation. This data was quantified using the sum of dBMI values, \({x}_{1i}={\sum }_{t=1\mathrm{ Jan}}^{26\mathrm{ Jan}}{\mathrm{dBMI}}_{ti}\) for prefecture i and day t. We collected data to model the response timing in COVID-19 control strategies. Specifically, for non-Hubei provinces we collected the dates of the Level One emergency response policy (Additional file 2: Table S1). For the cities in the Hubei province, we collected the dates on which the cities initiated shutdown procedures (Additional file 3: Table S2). Additional file 1: S1 provides detailed information about data preparation.

Data analysis

To analyze the effect generated by the control strategy, we first introduced static models and then dynamic models, which were extended from Jia et al. [12]. The static models generate a cross-sectional analysis of the effect of the control strategy on daily infections and the dynamic models investigate the time-varying effect. Only data after implementing COVID-19 control strategies (26 January) were included for statistical modelling. Because our dynamic models were developed based on a sigmoidal growth pattern of cases, case data with abnormal fluctuations on the epidemic curves were first pre-processed. A variety of events in the early phases of the epidemic lead to abnormal spikes in the number of cases. As more information about the virus surfaced, health professionals amended the diagnosis criteria for the disease (especially for cities in Hubei, Additional file 4: Table S3). Additionally, bulk reporting of jail infections (as happened in Shangdong and Zhejiang) caused other spikes. These anomalies were assumed to be infected several days before and should have been diagnosed and reported earlier than their actual report date. A presumed ought-to-be-reported date was generated for each of them through random assignment. Additional file 1: S2.1 provides a mathematical description for this process.

The static models based on the gravity model [15] have been developed to characterize the effect of population outflow from Wuhan on infections in other prefectures [12]. We extended it to statistically and cross-sectionally investigate the role of the outbreak control strategies. The saturated static model has the following multiplicative exponential form:

$${y}_{i}=c{\cdot e}^{{\sum }_{k}^{3}{\beta }_{k}{x}_{ki}}{e}^{{\lambda }_{1}\cdot {\mathrm{I}}_{Hubei}+{\lambda }_{2}\cdot {D}_{response}+{{\lambda }_{3}\cdot \mathrm{I}}_{Hubei}*{D}_{response}}$$

where \({y}_{i}\) is the cumulative number of confirmed cases in prefecture i; \({x}_{1i}\) is the aggregate outflow from Wuhan between 1 to 26 January to prefecture i, as described above; \({x}_{2i}\) is per capita GDP; \({x}_{3i}\) is the population density; c and \({\beta }_{k}\) are parameters to estimate. \({\mathrm{I}}_{Hubei}\) is an indicator function with \({\mathrm{I}}_{Hubei}=1\) for implementing Hubei strategy, otherwise \({\mathrm{I}}_{Hubei}=0\); \({D}_{response}\) is the response timing (23, 24, or 25 January); \({\mathrm{I}}_{Hubei}*{D}_{response}\) denotes their interaction; \({\lambda }_{j}\) s are the parameters.

A dynamic model under the Cox proportional hazards framework replaces the constant parameters in the static model with a time-varying hazard function \({h}_{0}\left(t\right)\) to model a sigmoidal growth of COVID-19 cases:

$$h\left(t|{x}_{i}\right)={{h}_{0}\left(t\right)e}^{{\sum }_{k}^{3}{\beta }_{k}{x}_{ki}} {e}^{{\lambda }_{1}\cdot {\mathrm{I}}_{Hubei}+{\lambda }_{2}\cdot {D}_{response}+{\lambda }_{3}\cdot {\mathrm{I}}_{Hubei}*{D}_{response}}$$

where \(h\left(t|{x}_{i}\right)\) is the hazard function describing cumulative number of confirmed cases at time t given other variables \({x}_{i}=\{{x}_{1i},{x}_{2i},{x}_{3i}\}\). Here, the logistic function with parameters \(\alpha\), \(\gamma\) and \(\omega\) was utilized as the sigmoidal growth approximation [16, 17]:

$${h}_{0}\left(t\right)=\frac{\alpha }{1+{e}^{-\gamma t+\omega }}$$

To further improve the model fitting, \({\lambda }_{j}\)s are allowed to be time dependent, that is, \({\lambda }_{j}={\lambda }_{j}(t)\). They are empirically determined by the results from the static models.

The R package minpack.lm [18] with a nonlinear least-squares Levenberg–Marquardt (LM) algorithm was used for model fitting and parameter estimation. Models were evaluated using the Bayesian Information Criterion (BIC) and R2. Additional file 1: S2.2 and S2.3 provide detailed methodologies for this section.

After statistical modelling, the indices of control effectiveness were calculated by exploiting the integral of the differences between predicted and actual case number:

$$IC{E_i} = \sum\limits_{t = 26 {\rm{Jan}}}^{29 {\rm{Feb}}} {\left[ {h\left( {t|{x_i}} \right) - \hat h\left( {t|{x_i}} \right)} \right].}$$

The normalized \({ICE}_{i}\)s were used for final measurements of control effectiveness.


Importation risk and response timing in each prefecture

In Wuhan, the 2020 Chinese New Year travel rush was interrupted by its lockdown on 23 January and the fluxes sharply declined to almost no movement since 27 January (Fig. 1a). The total aggregate amount of people who entered any study prefecture from Wuhan during this observation period (1–26 January) was used to measure their importation risk (Fig. 1b). More prefectures in Hubei are colored in dark red, indicating high risk of case importation. Comparison using the Wilcoxon Rank Sum test showed that the importation risk to other parts of Hubei is much higher than that to outside provinces with P < 0.001 (Additional file 7: Figure S1). Regional differences in risk, both in and out of Hubei, were also observed.

Fig. 1
figure 1

Geographical distribution of the risk of COVID-19 importation to each prefecture and their response timing. a Fluxes of population outflow from Wuhan, represented by dBMIs, during period of 1 January to 29 February in 2019 and 2020. The vertical lines in grey and red represent the dates of Wuhan lockdown and Chinese Lunar New Year, respectively. dBMI daily Baidu Mobility Index. b Spatial distribution of COVID-19 importation risk, measured by their aggregate population inflow from Wuhan until 26 January 2020. Provincial inflow from Wuhan was plotted at the upper left corner and municipal inflow in Hubei plotted at the lower right corner. The black area is Wuhan city. c Governments’ response timing in each prefecture. The dates to raise provincial public health response to Level One for COVID-19 control were plotted using different colors (upper left). Similarly, time to execute city shutdown in Hubei was plotted (lower right). d The importation risk distribution grouped by governments’ response. Samples with insufficient size at the response timing were excluded, such as Qinghai, Tibet and Xiangyang

Since 23 January 2020, all the provinces across China successively launched or raised their major public health emergency response to Level One, the highest level (Fig. 1c). Following Wuhan’s outbound traffic closure on 23 January, all other cities in Hubei subsequently announced their shutdown (Fig. 1c). The decline timelines on human mobility had a high consistency with the dates to execute the interventions (Additional file 8: Figure S2). To represent their respective response timing for COVID-19 control, we used the dates of the Level One response launch in non-Hubei provinces and the date of shutdown procedures in Hubei cities (excluding Wuhan). Figure 1d showed that the time to trigger COVID-19 control outside of Hubei was generally consistent with the distribution of the local importation risk with P < 0.001 using Jonckheere–Terpstra test. The time to trigger control measures was contrary to the importation risk with P = 0.042 for cities in Hubei other than Wuhan.

Control strategy in Hubei and non-Hubei regions

The major public health emergency response triggered an array of actions by provincial and/or local governments. The main policy instruments deployed included: (1) travel restrictions, (2) case finding and contact tracing, (3) isolation and management of infected individuals and exposed contacts, (4) social distancing, and (5) closed-off community management (Table 1). Despite many similar interventions implemented both in Hubei and non-Hubei areas, important differences remained in control strategies. The control strategy in Hubei required almost all people stay under ‘monitored’ self-quarantine at home. The concrete measures involved unprecedentedly strict closure and traffic restrictions, the ‘monitored’ stay-at-home order, extreme social distancing without any public or business activities, and complete closed-off community management. In contrast, other provinces allowed work resumption from 10 February after application and approval. All resumed enterprises were encouraged to work from home, and businesses’ employees were encouraged to telecommute. The limited enterprises approved to work on-site were required to take rigorous measures to prevent gatherings and cross-infection.

Table 1 Comparison between Hubei and non-Hubei control strategy

COVID-19 epidemic trend and association with control strategy

As of 29 February 2020, a total of 30,702 cases of COVID-19 were confirmed in mainland China excluding Wuhan (Fig. 2a). Among of them, 17,785 (58.9%) cases were reported in Hubei (excluding Wuhan) and 12,917 (42.1%) outside Hubei. These small numbers of cases, relative to tens of millions of people in Hubei or more than 1 billion people outside the province, suggest that the control policy had been working for the epidemic. Specifically, most prefectures experienced a typical sigmoidal growth in cases. This growth was characterized by an accelerated increase in the beginning and a flat period after mid-February (Additional file 9: Figure S3). As expected, Hubei experienced a more serious epidemic characterized by more rapid growth in infections. Abnormal fluctuations were detected in Shandong, Zhejiang, and cities in Hubei, as described in “Methods” section.

Fig. 2
figure 2

COVID-19 progress and its association with importation risk and control strategy. a The epidemic curves of COVID-19 in Hubei and non-Hubei area by 29 February 2020. b Relationship over time between the number of confirmed cases (cumulative until 29 February 2020) and total population inflow (up to 26 January 2020) from Wuhan, both on a logarithm scale. c The relationship between the log-transformed importation risk (the total population outflow from Wuhan up to 26 January 2020) and the log-transformed number of confirmed cases by prefectures on 29 January 2020. Circles are prefectures in Hubei; rectangles are prefectures outside Hubei; and the point sizes are proportional to the population density of the prefecture. The linear fitting is done for overall (black), Hubei (red) and non-Hubei (cyan) data. d The distribution of confirmed cases on 29 January 2020, grouped by governments’ response including response timing and response strategy on a logarithm scale. Samples with insufficient size at the response timing were excluded, such as Qinghai, Tibet and Xiangyang

Despite the strong association between case progress and importation risk (Fig. 2b), the daily total number of cases in the Hubei area tended to split apart from those in provinces outside Hubei. The overall fit line for all study prefectures highlights this split between the Hubei and non-Hubei areas taken separately (Fig. 2c). Further analysis confirmed this split pattern was maintained over time (Additional file 10: Figure S4a, b). The infectious cases grouped by governments’ responses showed a similar V-shape distribution to their importation risk (Fig. 2d and Additional file 10: Figure S4c, d). These data reveal the potential effects of the intervention control strategies.

The static and dynamic models for quantifying the effect of control strategy

The preprocessing of case data with abnormally reported dates generated smooth epidemic curves (Additional file 11: Figure S5). In subsequent modelling, the cumulative Wuhan population inflow was always included because of its strong correlation with new daily infections. Each prefecture’s GDP and population density were excluded because statistical tests provided no evidence of correlation (P > 0.05).

The static models including aggregate inflow population from Wuhan (before 26 January) and the control strategy as independent variables generated the consistently negative estimates of the coefficients on the Hubei strategy (Additional file 5: Table S4). This finding indicates that the implementation of the Hubei strategy was a continuous protective factor (all P < 0.05) in contrast with the non-Hubei strategy. Its protective effects declined over time and remained at a stable level in late February (Additional file 12: Figure S6a). It is noteworthy and significant that the R2 increased more at the accelerated growth period of COVID-19 cases after introducing the control strategy (Additional file 12: Figure S6b). This observation implies it contributed more during the earlier dates of the epidemic. Subsequently, the response timing was added but no overall effect was found (all P > 0.05). When the interaction between the control strategy and response timing was added, the effect of response speed on the epidemic became statistically significant before early February (P < 0.05) and all had negative slope estimates (Additional file 5: Table S4). This suggests that the response timing in non-Hubei area may be a ‘risk’ factor during the earlier stage in the implementation of the control strategy. Early implementation of control measures triggers better detection of actual cases, and this should not be considered a ‘risk’ from a policy perspective. This point is further addressed in “Discussion”.

The dynamic model using two variables of total population outflow from Wuhan (during 1 to 26 January) to each prefecture and the control strategy showed R2 = 0.910 (Fig. 3a); and the inclusion of the response timing increased R2 to 0.922 (Additional file 6: Table S5). According to the features of the coefficients from the static models, we also introduced a quadratic function of time t for the variable of the control strategy and a truncated function of time t with the cut-off T for the response timing. As expected, the BIC improved and R2 further increased to 0.938 (Additional file 6: Table S5). The truncation date T was fixed on 10 February based on the BIC statistics.

Fig. 3
figure 3

Estimation of time-varying effect of Hubei control strategy. a The fitted performance of our dynamic model (see dynamic model I in Additional file 1: Section S2.3). b Change of estimated relative risk over time, comparing between Hubei control strategy and non-Hubei control strategy (as reference). c Comparison of index of control effectiveness between all study prefectures

Finally, the time-varying effect of the control strategy estimated by the dynamic model with time-dependent effects is shown in Fig. 3b. We used the non-Hubei strategy for COVID-19 control enacted on 23 January as a reference over time. The response timing was unexpectedly shown to be a weak risk factor at about 2 weeks after its implementation. In contrast, the Hubei control strategy showed a very strong protective effect and the protective effect rapidly declined if the implementation was delayed. Precisely, the Hubei control strategy showed a 4 times greater protective effect than the non-Hubei strategy on 26 January if both were taken on 23 January. Despite a narrowing effect, the Hubei area achieved about 3 times fewer cases than non-Hubei areas in late February. The 2-day delay of the Hubei strategy finally narrowed down the effect to 2 times fewer cases. The marginal analysis showed that the Hubei control strategy had the marginal effect curve equal to its effect trajectory enacted on 23 January (Fig. 3b). Overall, the Hubei strategy for COVID-19 control would achieve 3 times fewer cases than the non-Hubei strategy.

Evaluation of control effectiveness

The differences in the growth trends between predicted and observed cases can be used to benchmark the control of COVID-19 for the prefectures (Additional file 13: Figure S7 and Additional file 14: Figure S8), providing an under-performing or over-performing order. We used the integral of the differences over time to create a total index for control effectiveness. After sorting these indexes, we can identify a list of under-performers and over-performers (Fig. 3c). Indeed, most prefectures in Hubei were prone to be over-performers because of their strong control measures. Prefectures like Zhejiang, Xiaogan and Suizhou were observed as under-performers.


The COVID-19 epidemic broke out in December 2019. China employed almost all the available infectious disease control tools at an unprecedented scale [1] and successfully controlled the epidemic by the end of February 2020. Although importation risks were considered in policy making, the risk-based differentiated need for COVID-19 control in and out of Hubei was inadequately assessed (Fig. 1). The whole of Hubei province was defined as a key area of the epidemic, and a Hubei-specific strategy for COVID-19 control was executed. Other areas followed a non-Hubei control strategy.

The Hubei-specific control strategy consisted of a series of strict stay-at-home measures. We compared non-pharmaceutical interventions between the Hubei and non-Hubei control strategies. The Hubei measures served to enforce adherence to self-quarantine at home directives (Table 1). In contrast, the control strategy executed outside Hubei was typically a series of physical distancing measures (e.g. social distancing) which have been experimentally proven to be effective in delaying and reducing the height of the peak and median epidemic size [7, 19].

Previous studies have documented the strong linear relationship between population movement from Wuhan and the number of infections [20,21,22]. This study has revealed the effect of applying different control strategies in and out of Hubei. Despite the higher importation risk in Hubei (excluding Wuhan), the epidemic was under control 1 week after COVID-19 containment outside Hubei (Fig. 2a). Both the static and dynamic models showed that the Hubei strategy was a very strong protective factor in contrast with the non-Hubei strategy. This protective effect proved to be time varying. As expected, the higher controlling effect appeared in the earlier dates of the outbreak when the spread was more severe. The protective effect later narrowed gradually and then stabilized (Fig. 3b). By mid-February, the Hubei control strategy obtained 3 times fewer COVID-19 cases compared to the non-Hubei strategy. It’s worth noting that the estimated effects are likely to be conservative because of the introduction of ‘clinically diagnosed cases’ into the analysis, an expanded case definition specifically for Hubei province due to the insufficient testing capacity in the 5th version of the diagnosis and treatment guidelines [23, 24].

The government response is also of interest. Through the statistical modelling, a weak ‘protective’ effect was found for the later implementation of non-Hubei control strategy, and this protective effect disappeared about 2 weeks after implementing the measure (Fig. 3b). This result differs from a previous study [8]. This seemingly unreasonable result has a remarkable interpretation. The Level One response required suspected and confirmed cases to be isolated and reported immediately [6]. During the earliest phase of the epidemic, the virus was diffused, and the response triggered aggressive case and contact identification. This earlier response means that more cases were found earlier and isolated. Therefore, the accelerated response being a ‘risk’ factor is an artifact of improved tracing and disease detection.

Some limitations in our study have been recognized. First, we focused on province-level data outside Hubei unlike the city-level data available for the Hubei area. It is known that the importation risk and daily amount of infections in most cities outside Hubei was incompatible with those in the Hubei area. These differences could confound the effects discovered for the Hubei control strategy. Second, random assignment for sharp discontinuities on cumulative case curves was used to smooth the data. These corrections did not include other considerations of disease progression, such as the incubation period. Early phase fluctuations in case data for the epidemic were at least partly driven by changes in case definitions and adapting health system processes. These issues made it necessary to simplify the assignment process. Third, the cumulative population movement from the epicenter using location-based service data is a simplified estimation of the importation risk at each prefecture. Risk considering case underestimation and reporting delay through probability or mechanistic modelling [25,26,27] may be another choice. However, extra data collection and complicated calculation may defer the decision making about the epidemic control.


A risk-based control strategy would improve the effective response for COVID-19 control. Our study shows that the stricter Hubei strategy can achieve better control effectiveness than distancing-focused non-Hubei strategies. These findings highlight the policy impacts and health benefits of precise and differentiated strategies informed by constant monitoring of outbreak risk.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its additional information files. They are all publicly available. Specifically, control measures (Table 1, Additional file 2: Table S1 and Additional file 3: Table S2), case definitions (Additional file 4: Table S3), case data (Additional file 15: Case data), and population and economic data (Additional file 16: Population and economic data) were extracted from local official websites, and travel data (Additional file 17: Travel data) were retrieved through Baidu Qianxi platform (



Coronavirus disease 2019


Daily Baidu Mobility Index


  1. Wilder-Smith A, Freedman DO. Isolation, quarantine, social distancing and community containment: pivotal role for old-style public health measures in the novel coronavirus (2019-nCoV) outbreak. J Travel Med. 2020;27(2):taaa020.

    Article  Google Scholar 

  2. Hsiang S, Allen D, Annan-Phan S, Bell K, Bolliger I, Chong T, et al. The effect of large-scale anti-contagion policies on the COVID-19 pandemic. Nature. 2020;584(7820):262–7.

    Article  CAS  Google Scholar 

  3. Kraemer MUG, Yang CH, Gutierrez B, Wu CH, Klein B, Pigott DM, et al. The effect of human mobility and control measures on the COVID-19 epidemic in China. Science. 2020;368(6490):493–7.

    Article  CAS  Google Scholar 

  4. Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science. 2020;368(6489):395–400.

    Article  CAS  Google Scholar 

  5. Flaxman S, Mishra S, Gandy A, Unwin HJT, Mellan TA, Coupland H, et al. Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature. 2020;584(7820):257–61.

    Article  CAS  Google Scholar 

  6. Hellewell J, Abbott S, Gimma A, Bosse NI, Jarvis CI, Russell TW, et al. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. Lancet Glob Health. 2020;8(4):e488–96.

    Article  Google Scholar 

  7. Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. Lancet Public Health. 2020;5(5):e261–70.

    Article  Google Scholar 

  8. Tian H, Liu Y, Li Y, Wu CH, Chen B, Kraemer MUG, et al. An investigation of transmission control measures during the first 50 days of the COVID-19 epidemic in China. Science. 2020;368(6491):638–42.

    Article  CAS  Google Scholar 

  9. Hatchett RJ, Mecher CE, Lipsitch M. Public health interventions and epidemic intensity during the 1918 influenza pandemic. Proc Natl Acad Sci USA. 2007;104(18):7582–7.

    Article  CAS  Google Scholar 

  10. China’s State Council Information Office. Fighting COVID-19: China in action. The State Council Information Office of the People’s Republic of China. 2020. Accessed 15 Sept 2000.

  11. World Health Organization. Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19). World Health Organization. 2020. Accessed 15 Sept 2020.

  12. Jia JS, Lu X, Yuan Y, Xu G, Jia J, Christakis NA. Population flow drives spatio-temporal distribution of COVID-19 in China. Nature. 2020;382(7812):389–94.

    Article  Google Scholar 

  13. Liu K, Ai S, Song S, Zhu G, Tian F, Li H, et al. Population movement, city closure in Wuhan and geographical expansion of the 2019-nCoV pneumonia infection in China in January 2020. Clin Infect Dis. 2020;71(16):2045–51.

    Article  CAS  Google Scholar 

  14. Zhong A, Guo S, Chen T. Correlation between travellers departing from Wuhan before the Spring Festival and subsequent spread of COVID-19 to all provinces in China. J Travel Med. 2020;27(3):taaa036.

    Article  Google Scholar 

  15. Zipf GK. The p1 p2/d hypothesis: on the intercity movement of persons. Am Sociol Rev. 1946;11(6):677–86.

    Article  Google Scholar 

  16. Burger R, Chowell G, Lara-Diiaz LY. Comparative analysis of phenomenological growth models applied to epidemic outbreaks. Math Biosci Eng. 2019;16(5):4250–73.

    Article  Google Scholar 

  17. Roosa K, Lee Y, Luo R, Kirpich A, Rothenberg R, Hyman JM, et al. Short-term forecasts of the COVID-19 epidemic in Guangdong and Zhejiang, China: February 13–23, 2020. J Clin Med. 2020;9(2):596.

    Article  Google Scholar 

  18. R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. R Core Team. 2020. Accessed 15 Feb 2020.

  19. Thu TPB, Ngoc PNH, Hai NM, Tuan LA. Effect of the social distancing measures on the spread of COVID-19 in 10 highly infected countries. Sci Total Environ. 2020;2020(742):140430.

    Article  Google Scholar 

  20. Colizza V, Barrat A, Barthelemy M, Vespignani A. The role of the airline transportation network in the prediction and predictability of global epidemics. Proc Natl Acad Sci USA. 2006;103(7):2015–20.

    Article  CAS  Google Scholar 

  21. Balcan D, Colizza V, Goncalves B, Hu H, Ramasco JJ, Vespignani A. Multiscale mobility networks and the spatial spreading of infectious diseases. Proc Natl Acad Sci USA. 2009;106(51):21484–9.

    Article  CAS  Google Scholar 

  22. Halloran ME, Vespignani A, Bharti N, Feldstein LR, Alexander KA, Ferrari M, et al. Ebola: mobility data. Science. 2014;346(6208):433.

    Article  CAS  Google Scholar 

  23. Lin L, Li TS. Interpretation of “Guidelines for the Diagnosis and Treatment of Novel Coronavirus (2019-nCoV) Infection by the National Health Commission (Trial Version 5).” Zhonghua Yi Xue Za Zhi. 2020;100(11):805–7.

    CAS  PubMed  Google Scholar 

  24. Tsang TK, Wu P, Lin Y, Lau EHY, Leung GM, Cowling BJ. Effect of changing case definitions for COVID-19 on the epidemic curve and transmission parameters in mainland China: a modelling study. Lancet Public Health. 2020;5(5):E289–96.

    Article  Google Scholar 

  25. Du Z, Wang L, Cauchemez S, Xu X, Wang X, Cowling BJ, et al. Risk for transportation of coronavirus disease from Wuhan to other cities in China. Emerg Infect Dis. 2020;26(5):1049–52.

    Article  CAS  Google Scholar 

  26. Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. Lancet. 2020;395(10225):689–97.

    Article  CAS  Google Scholar 

  27. Xu XK, Wang L, Pei S. Multiscale mobility explains differential associations between the gross domestic product and COVID-19 transmission in Chinese cities. J Travel Med. 2021;28:taaa36.

    Article  Google Scholar 

Download references


The authors would like to thank Dr. Paul W. Bible of Marian University in Indianapolis Indiana for his helpful comments.


This work was supported by Guangdong Basic and Applied Basic Research Foundation (2021A1515011591), Guangdong Medical Science and Technology Research Foundation (A2021104) and Sun Yat-sen University 2020 Emergency Funding for COVID-19 Prevention and Control (No grand number provided). The funders played no role in study design, data collection, and analysis, decision on publishing, or preparation of the manuscript.

Author information

Authors and Affiliations



YL and YH conceived and designed this study. FZ and MJ collected the data. YL and ZD developed the model and ran the analysis. JL, JG and DY advised on model development. YL and FZ drafted the manuscript. JL, JG and DY revised the manuscript and gave scientific comments. SG and YH provided critical revision for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yuantao Hao.

Ethics declarations

Ethics approval and consent to participate

This is an observational study using publicly available aggregate data without any individual information. No ethical approvals were required.

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.

Supplementary Information

Additional file 1.

Text. Supplementary methods for data preparation and analysis.

Additional file 2: Table S1.

Time to execute Level One public health emergency response.

Additional file 3: Table S2.

Time to execute city shutdown in Hubei.

Additional file 4: Table S3.

Definitions of suspected and confirmed cases of COVID-19 in six versions implemented before February 29, 2020.

Additional file 5: Table S4.

Results from the static models.

Additional file 6: Table S5.

Parameter estimation for the dynamic models.

Additional file 7: Figure S1.

Comparison of COVID-19 importation risk between Hubei and non-Hubei areas. The risk of COVID-19 importation to each prefecture is defined as their total population inflow from Wuhan until 26 January 2020.

Additional file 8: Figure S2.

The response timeline of COVID-19 control on human mobility. (a) Provincial population outflow outside Hubei. The flow of each province was represented by the responding BMIs of their capital city. (b) Municipal population outflow in Hubei excluding Wuhan. (c and d) The decline timelines of the outflow BMIs were compared after they were scaled to the values on 22 January. Legends for provincial population flow outside Hubei (a and c) and municipal population flow in Hubei (b and d) are shared, and respectively displayed in (a) and (b). All the colors here followed the response timelines shown in Fig. 1c. BMI: Baidu Mobility Index.

Additional file 9: Figure S3.

The epidemic curves of COVID-19 in each prefecture by 29 February 2020. Each circular bar on the polar coordinate system represents daily number of reported cases in the prefecture. The grey and red circular lines in bold indicate the dates of Wuhan lockdown and Chinese Lunar New Year. The polar coordinates also show the number of prefectures with COVID-19 case report in parentheses. (a) Provincial epidemic curves outside Hubei. (b) Municipal epidemic curves in Hubei other than Wuhan.

Additional file 10: Figure S4.

Associations between number of COVID-19 cases, importation risk and control strategy. (a-b) The relationship between the log-transformed importation risk (the total population outflow from Wuhan up to 26 January 2020) and the log-transformed number of confirmed cases by prefectures on 9 February 2020 (a) and 19 February 2020 (b). Circles are prefectures in Hubei; rectangles are prefectures outside Hubei; and the point sizes are proportional to the population density of the prefecture. The linear fitting is done for overall (black), Hubei (red) and non-Hubei (cyan) data. (c-d) The distribution of confirmed cases on 9 February 2020 (c) and 19 February 2020 (d), grouped by governments’ response including response timing and response strategy in a logarithm scale. Samples with insufficient size at the response timing were excluded, such as Qinghai, Tibet and Xiangyang.

Additional file 11: Figure S5.

Change of the epidemic curves after re-assignment of report date for cases with abnormal fluctuations. They included Shandong and Zhejiang jail cases intensively reported on 20 February 2020 and clinically diagnosed cases in Hubei area due to the amendment of the diagnosis and treatment program of the COVID-19.

Additional file 12: Figure S6.

Results generated by static models. (a) Change of relative risk over time, generated by the static model II. (b) Change of R2 over time, compared between static models.

Additional file 13: Figure S7.

Predicted versus actual case growth in the prefecture outside Hubei.

Additional file 14: Figure S8.

Predicted versus actual case growth in the prefecture in Hubei other than Wuhan.

Additional file 15.

The original case data.

Additional file 16.

The original population and economic data.

Additional file 17.

The original travel data.

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Zheng, F., Du, Z. et al. Evaluation of China’s Hubei control strategy for COVID-19 epidemic: an observational study. BMC Infect Dis 21, 820 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: