Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The impact of multi-level interventions on the second-wave SARS-CoV-2 transmission in China

  • Yuanchen He,

    Roles Data curation, Formal analysis, Visualization, Writing – original draft

    Affiliation School of Public Health, Fudan University, Key Laboratory of Public Health Safety, Ministry of Education, Shanghai, China

  • Yinzi Chen,

    Roles Conceptualization, Validation, Writing – review & editing

    Affiliation School of Public Health, Fudan University, Key Laboratory of Public Health Safety, Ministry of Education, Shanghai, China

  • Lin Yang,

    Roles Methodology, Writing – review & editing

    Affiliation School of Nursing, Hong Kong Polytechnic University, Hong Kong, China

  • Ying Zhou,

    Roles Writing – review & editing

    Affiliation School of Public Health, Shenzhen University, Health Science Center, Shenzhen, China

  • Run Ye ,

    Roles Supervision, Writing – review & editing

    erinwang@fudan.edu.cn (XW); yerun85@163.com (RY)

    Affiliation Department of Tropical Diseases, Navy Medical University, Shanghai, China

  • Xiling Wang

    Roles Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing

    erinwang@fudan.edu.cn (XW); yerun85@163.com (RY)

    Affiliations School of Public Health, Fudan University, Key Laboratory of Public Health Safety, Ministry of Education, Shanghai, China, Shanghai Key Laboratory of Meteorology and Health, Shanghai, China

Abstract

Background

A re-emergence of COVID-19 occurred in the northeast of China in early 2021. Different levels of non-pharmaceutical interventions, from mass testing to city-level lockdown, were implemented to contain the transmission of SARS-CoV-2. Our study is aimed to evaluate the impact of multi-level control measures on the second-wave SARS-CoV-2 transmission in the most affected cities in China.

Methods

Five cities with over 100 reported COVID-19 cases within one month from Dec 2020 to Feb 2021 were included in our analysis. We fitted the exponential growth model to estimate basic reproduction number (R0), and used a Bayesian approach to assess the dynamics of the time-varying reproduction number (Rt). We fitted linear regression lines on Rt estimates for comparing the decline rates of Rt across cities, and the slopes were tested by analysis of covariance. The effect of non-pharmaceutical interventions (NPIs) was quantified by relative Rt reduction and statistically compared by analysis of variance.

Results

A total of 2,609 COVID-19 cases were analyzed in this study. We estimated that R0 all exceeded 1, with the highest value of 3.63 (1.36, 8.53) in Haerbin and the lowest value of 2.45 (1.44, 3.98) in Shijiazhuang. Downward trends of Rt were found in all cities, and the starting time of Rt < 1 was around the 12th day of the first local COVID-19 cases. Statistical tests on regression slopes of Rt and effect of NPIs both showed no significant difference across five cities (P = 0.126 and 0.157).

Conclusion

Timely implemented NPIs could control the transmission of SARS-CoV-2 with low-intensity measures for places where population immunity has not been established.

Introduction

The emergence of coronavirus disease 2019 (COVID-19), which was first detected in December 2019 in Wuhan, China, has resulted in over 188 million cases and 4 million deaths around the world as of 16 July 2021 [1]. The first-wave of epidemic has been contained effectively by a series of stringent public health interventions in mainland China, especially the decisive city-level lockdown [2, 3]. The absence of herd immunity [46] and emerging SARS-CoV-2 variants still posed an upcoming threat of re-emergence to the whole society. After relaxing the interventions and resuming a normal life, the northeast cities of China experienced the second-wave of COVID-19 epidemic from Dec 2020 to Feb 2021. The first confirmed local case was reported on 2 Jan 2021 in Shijiazhuang city, Hebei province, and was likely originated from Russian strains suggested by genomic sequencing [7]. The epidemic spread rapidly to over 10 neighboring cities in Heilongjiang, Jilin and Liaoning provinces, within only 10 days.

As the COVID-19 vaccine campaign was not launched in mainland China until early 2021, non-pharmaceutical interventions (NPIs) remained the main way to control the second-wave of COVID-19. The effects of NPIs on controlling COVID-19 transmission have been widely confirmed in the first- and second-wave [812]. A study from Italy by Manica et al. [13] has estimated that the three-tiered restriction system, which was taken by Italian government to counter the national second-wave, effectively lowered the reproduction number below 1 in 85 out of 107 provinces. While some researchers believed that without strengthening interventions, it would be exceedingly difficult to contain a COVID-19 resurgence [14, 15] even as vaccination is coming into effect. China has adopted different levels of control measures including mass testing, social distancing and lockdown in cities since 6 Jan 2021 (Fig 1), but there is limited literature so far to investigate whether early implementing NPIs could avoid strong-intensity control measures. In this study, we aimed to assess the impact of multi-level interventions on SARS-CoV-2 transmission in the worst-hit cities during the second wave of COVID-19 in China.

thumbnail
Fig 1. Epidemic curves and non-pharmaceutical interventions implemented in five worst-hit cities in China, Jan-Feb 2021.

(A)-(E) The daily number of COVID-19 infections and the timings of specific non-pharmaceutical interventions in Shijiazhuang, Suihua, Tonghua, Haerbin and Changchun cities.

https://doi.org/10.1371/journal.pone.0274590.g001

Materials and methods

City selection and case definition

We restricted analyses to prefecture-level cities that reported the total number of confirmed COVID-19 cases and asymptomatic SARS-CoV-2 infected individuals over 100 within a consecutive month during Dec 2020 to Feb 2021. The definition of symptomatic and asymptomatic COVID-19 cases (details in S1 File) followed the Protocol on Prevention and Control of Novel Coronavirus Pneumonia (7th Edition) [16] issued on Sept 11, 2020 by Chinese Center for Disease Control and Prevention (China CDC). All SARS-CoV-2 infected individuals would be quarantined and treated by municipal designated hospitals.

Data collection

The daily cumulative number of confirmed COVID-19 cases and asymptomatic cases were extracted from the official websites of national, provincial, and municipal health commissions. Individual information on COVID-19 cases, including diagnosis date, type of SARS-CoV-2 infected individuals (symptomatic or asymptomatic) and source of infection (imported from other cities or local transmission) was collected from open access municipal official press release platforms.

The NPIs taken by local governments were collected from two publicly official sources: the press conferences on the novel coronavirus and authorized governmental news media. Following the protocol [16] from China CDC and guidelines from WHO [17], we summarized 8 types of main interventions (as defined in S1 File), including testing (mass/for the potential exposed), school closure, home quarantine, within-city traffic restriction, environmental sanitization, social distancing and lockdown. Further, we reclassified the interventions into three levels according to their influence on transmission rate. Based on previous studies [18, 19], we classified mass testing and potential exposed group testing as “Mild” NPIs, city lockdown as “Strong” NPIs and the rest of interventions as “Moderate” NPIs.

Statistical analysis

Considering that less than 0.02‰ of the people in the selected provinces had been confirmed as COVID-19 cases before second-wave (the Note 3 of S1 File), basic reproduction numbers of selected cities, denoted by R0, were estimated through the exponential growth rate r as in Formulation (1). (1) where M is the moment generating function of the discretized generation time distribution [20]. The appropriate expression of M is determined by the shape of generation time distribution. In the main analysis, the generation time (i.e., the lag between infection of the primary case and of their secondary cases) was assumed to follow a gamma distribution with a mean (SD) of 7.5(3.4) days [21]. Alternative scenarios about the values of the serial interval were explored in sensitivity analyses, one with mean (SD) of 5.7(1.8) days [22] and the other with mean (SD) of 5.0(1.7) days [23]. Simple moving average (n = 7) was used to acquire smoother epidemic curves and fit exponential growth. The time period used to fit the exponential growth was selected based on the highest R-squared values.

Time-varying reproduction number (Rt) was estimated by the Cori method [24] using the daily number of confirmed and asymptomatic cases by their diagnosis date (the Note 4 of S1 File). Using the infection-to-diagnosis delay distribution, the time series of daily number of diagnosis was deconvoluted to the time series of daily number of new infections [25, 26]. This estimation method specifically distinguishes between local cases and imported cases [27] to better estimate community transmission. To fulfil the criteria for when to start Rt estimation required by the method (listed in the Note 5 and 6 of S1 File), we started the calculation from the 7th day of the first local cases (including confirmed and asymptomatic cases). We assumed the relative infectiousness of the asymptomatic and pre-symptomatic were 50% and 100% of the symptomatic transmission, and set the time window τ to 7 days. Simple linear regression lines were fitted to daily Rt values from the start of estimation to the first day of Rt < 1 in each city, and analysis of covariance (ANCOVA) was used to test the statistical difference of regression slopes of Rt. To evaluate the impact of NPIs on controlling COVID-19 transmission, we also used the following Formulation (2) to quantify the relative changes of Rt [28]. (2) where ε is the controlling effect indicator, t0 and t1 is the time when Rt estimation onset and the minimum time when the mean Rt < 1. represents the corresponding Rt value at time t0. The indicators were statistically compared by analysis of variance (ANOVA). For comparing the timeliness of taking measures across the cities, we also calculated the delay from the first case confirmation to the first intervention implementation and the cumulative cases when imposing the first intervention respectively. Statistical analyses and visualization were all performed using the R software (version 4.1.0), the ANCOVA and ANOVA tests were implemented via stats and ggpmisc packages. Statistical significance was considered to be P < 0.05.

Results

City description

During Dec 2020 to Feb 2021, there were 5 cities reporting over 100 COVID-19 cases within a month (Fig 1). For Shijiazhuang, Suihua and Tonghua cities, the total numbers of confirmed and asymptomatic cases were 1043, 844 and 320 respectively, followed by 291 and 111 in Haerbin and Changchun cities. The second-wave of COVID-19 was first detected in Shijiazhuang city, Hebei Province on 2 Jan 2021 and then spread to the cities of neighboring provinces (e.g., Jilin and Heilongjiang Provinces). Other four cities (Suihua, Tonghua, Haerbin and Changchun city) began reporting local cases on Jan 10, 12, 15 and 11. We included 1351 asymptomatic cases in the analyses, among whom 658 (48.7%) were re-diagnosed as confirmed COVID-19 cases later. Shijiazhuang city had the most SARS-CoV-2 infected individuals (n = 1043) and the longest epidemic period (44 days) in the epidemic.

Epidemic curves and NPIs

Starting from the end of 2020, the second-wave of COVID-19 grew rapidly in the northeast of China, especially in Shijiazhuang, Suihua and Tonghua cities. Haerbin and Changchun cities implemented NPIs on the same day of reporting the first case, while the other three cities delayed actions by 1–4 days as shown in Table 1. The cumulative cases were 233, 8, 13, 3 and 7 respectively in the cities when the first intervention was implemented. In the epidemic, the municipal CDC implemented mass testing for case investigation and contact tracing. Each city implemented its most stringent control measures within a week from the start of local epidemic.

thumbnail
Table 1. The description of COVID-19 epidemic before imposing NPIs in five worst-hit cities in China, Jan-Feb 2021.

https://doi.org/10.1371/journal.pone.0274590.t001

Notably, the three worst-hit cities, Shijiazhuang, Suihua and Tonghua, took “Strong” measures (city lockdown) over a month (implemented in Jan 7- Feb 8, Jan 12- Feb 22 and Jan 18- Feb 20 respectively) to contain the epidemic. The cumulative cases in these three cities were 299, 12 and 122 respectively on the first day of city lockdown. These three cities kept multiple rounds of mass testing during the lockdown period, Shijiazhuang also implemented environmental sanitization, within-city traffic restriction and potential exposed group testing. In addition, Suihua city adopted environmental sanitization on Jan 20 as Shijiazhuang did on Jan 17.

Haerbin and Changchun cities only adopted “Mild” and “Moderate” control measures instead. The highest city-level measures were initiated when local cases accumulated to 14 and 8 respectively, which specifically referred to within-city traffic restriction and environmental sanitization. Besides, the testing strategies in the two cities were both transferred from potential exposed group testing to mass testing during the “Moderate” period.

Basic reproduction number (R0) and time-varying reproduction number (Rt)

The start of the best fitted time period lagged 4–9 days behind the start of epidemic time (curves after moving average were shown in S1 Fig of S1 File). All cities’ R0 values were well above 1 (Table 2), with the highest in Haerbin (3.63, 95% CI: 1.36–8.53) and the lowest in Shijiazhuang (2.45, 95% CI: 1.44–3.98). Different generation time distributions showed similar results of R0 with main analysis, which validated robustness of R0 (in the S2 Fig of S1 File).

thumbnail
Table 2. Basic reproduction number (R0) estimated by the exponential growth method.

https://doi.org/10.1371/journal.pone.0274590.t002

The Rt values of Shijiazhuang, Suihua, Tonghua and Haerbin city showed a continuous downward trend (Fig 2), while those of Changchun city showed a small rise on January 20 followed by a steady decline. The starting time of Rt < 1 was not the same calendar time between these 5 cities, but around the day 11 to 13 of reporting local COVID-19 cases. The daily estimated values of Rt in all cities are shown in S1 File. When different levels of interventions were taken at the same time, we only retain the highest level of NPIs in each city for visualization (Fig 2). The results revealed that the highest level of NPIs had been implemented before or exactly the day we began to estimate effective reproduction number, and Rt kept falling since the implementation. Fitting the regression lines with Rt of each city, the slopes were not statistically different (P = 0.126, Fig 3). Comparing the effectiveness of NPIs as shown in Table 3, Haerbin city had the highest effect among five cities though without city lockdown, while Shijiazhuang city had the lowest effect with the most kinds of NPIs implemented. The effectiveness of NPIs were not significantly different as the P value was above 0.05.

thumbnail
Fig 2. Temporal dynamics of SARS-CoV-2 transmission in five worst-hit cities in China, Jan-Feb 2021.

The light, medium and dark grey regions show periods of mild, moderate and strong NPIs respectively. The red line represents the first day of Rt < 1 in each city. D represents the days from the initial identification of local cases to Rt below the threshold. (A) Estimated Rt values from Jan 8 to Feb 6 in Shijiazhuang city. (B) From Jan 16 to Feb 11 in Suihua city. (C) From Jan 18 to Feb 8 in Tonghua city. (D) From Jan 21 to Feb 8 in Haerbin city. (E) From Jan 17 to Feb 6 in Changchun city. Rt = time-varying reproduction number.

https://doi.org/10.1371/journal.pone.0274590.g002

thumbnail
Fig 3. The scatter plots and fitted regression lines of Rt estimations in five worst-hit cities in China.

The “red” lines represent the cities with “Strong” level NPIs, the “blue” lines represent the cities without “Strong” level NPIs. The equations of linear regression (top right) were fitted from the Rt values from the start of calculation to the first day of Rt < 1 in each city. The P value below the equations were the p-value for the hypothesis test based on the analysis of covariance (ANCOVA) evaluating whether the decline in transmission rate was different across cities. Rt = time-varying reproduction number.

https://doi.org/10.1371/journal.pone.0274590.g003

thumbnail
Table 3. The effectiveness of NPIs in five worst-hit cities in China.

https://doi.org/10.1371/journal.pone.0274590.t003

Discussion

Here we retrospectively assessed the transmission potentials and epidemic evolution of COVID-19 using Rt in northeast of mainland China, the worst affected area in the second-wave epidemic. Three of the cities included in the study adopted “Strong” level NPIs, while the other two only implemented “Mild” and “Moderate” level NPIs. The Rt in all five cities fell below 1 within 2 weeks from the first case confirmation. The decline rate of Rt and effectiveness of NPIs were not significantly different across cities.

The estimates of R0, ranging from 2 to 4, showed no significant difference among selected cities. Reviewing the first wave, a narrow range of R0, between two and three, were frequently estimated from data on confirmed cases [2931]. Specifically in some studies that located in Wuhan, China [21, 3235], range of R0 estimates reported from at least 2.20 to 2.71. The contact patterns after the first wave couldn’t be comparable with the pre-pandemic level, even with a few restriction policies. The mean number of contacts only reached about 50% or less of the pre-pandemic level [3638]. Considering the transmissibility of SARS-CoV-2, which usually measured by R0, was significantly influenced by population contact level, R0 in the first- and second- waves should be compared with caution and required more information.

The rapid decline of Rt suggests the effectiveness of non-pharmaceutical interventions implemented in five cities. Owing to stringent NPIs implemented in the early phase, Rt has been below the epidemiological threshold as of Jan 26, 2021 in all cities. The comparable Rt slopes suggested that with only “Mild” and “Moderate” measures, Haerbin and Changchun cities still successfully drove Rt below 1 as in other three cities. The prominent effectiveness enlightened that if NPIs are taken as soon as possible, as Haerbin and Changchun city began to implement city-wide measures when there were few sporadic cases, the highest level of interventions (i.e., city lockdown) may not be needed, thus reducing the impact on socioeconomics. The heterogeneity of sociodemographic and geographic characteristics in five cities may confound the association between Rt and effectiveness of NPIs, but the confounding effect might be minor. Geographically, the five cities are in the Hebei, Heilongjiang and Jilin provinces, concentrated in northeast China. For the age structure, the proportion of population aged 0–14, 15–59, and 60+ ranged from 10–19%, 62–67% and 18–25%, respectively [3942]. Economically, these cities mainly develop heavy industry with the per capita GRP in 2020 being 52.961, 30.257, 34.338, 51.597 and 77.634 thousand respectively in Shijiazhuang, Suihua, Tonghua, Haerbin and Changchun cities, reported by the official Statistical Yearbook [4345]. In the analysis, we classified interventions that besides testing and lockdown policies, rest of the interventions were all at the “Moderate” level. The individual effects of these interventions were hard to quantify in reality partly because of close intervals in implementation [19] and partly because of the limited impact on Rt with a single control measure [46]. The concurrent interventions might obscure the association between Rt and the impact of NPI at a specific level, which required further exploration.

A world research found even with concurrent lockdown policies, school closure and social distancing still worked on reducing COVID-19 infections. But the within-city traffic restriction had no such effect [47]. Another systematic review also found that no evidence on the effectiveness of public transport closure [48]. This counterintuitive result might be due to the interventions imposed before the traffic restriction had reduced human mobility already, as shown in Shijiazhuang and Haerbin cities. These research enlightened that effects of interventions in the “Moderate” level might vary and need more detailed exploration.

The study presented here is subject to several limitations. First, the study focused primarily on population-based rather than individual-based interventions, but targeted test, trace and isolation strategy might be also effective to prevent further COVID-19 transmission [49]. On the other hand, the success of individual-based interventions relies highly on the response speed of disease surveillance system in the early stages [23, 50]. As this was a pronounced challenge for places where community transmission had occurred, Shijiazhuang, Suihua and Tonghua implemented mass testing as the first intervention to quickly identify exposed population. Second, we extracted the parameters of generation time distribution from published literatures. As the virus mutating, the epidemiological characteristics of first-wave might not be suitable for this epidemic in mainland China. But sensitivity analyses of R0 by changing the generation time reveled the robustness of our estimates. Third, to our knowledge, an asymptomatic individual infected 102 new cases across three cities through crowded indoor lectures in this epidemic. Considering the source of superspreading events was from Suihua city and the efficacy of contact tracing inevitably wouldn’t reach 100% [51], even only lagging for 1 day to impose NPIs, Suihua still take city lockdown to contain the epidemic. The key point is that reproduction number only represents average number of secondary cases in a region [52]. This average ignores the 80/20 rule [53], which also means that superspreading events are not detected by R0 or Rt and hard to capture on a population level.

Conclusion

In conclusion, our results provided an overview and quantify the effects of multi-level interventions of second-wave COVID-19 epidemic in mainland China. The findings suggested a fast downward trend of the COVID-19 outbreak in the northeast provinces, indicating that non-pharmaceutical interventions still played an integral part in stopping the epidemic as in the first-wave, and decisive enforcement might even have dramatic effects. Our study indicated that timely implemented non-pharmaceutical interventions could control the transmission of SARS-CoV-2 with low-intensity measures for places where population immunity have not been established.

Supporting information

S2 File. The daily number of reported cases in five cities.

https://doi.org/10.1371/journal.pone.0274590.s002

(XLSX)

References

  1. 1. Johns Hopkins University. COVID-19 Dashboard 2021 [cited 2021 July 16]. Available from: https://coronavirus.jhu.edu/map.html.
  2. 2. Pan A, Liu L, Wang C, Guo H, Hao X, Wang Q, et al. Association of Public Health Interventions With the Epidemiology of the COVID-19 Outbreak in Wuhan, China. JAMA. 2020;323(19):1915–23. pmid:32275295
  3. 3. Leung K, Wu JT, Liu D, Leung GM. First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment. Lancet. 2020;395(10233):1382–93. pmid:32277878
  4. 4. Aleta A, Martín-Corral D, Pastore Y Piontti A, Ajelli M, Litvinova M, Chinazzi M, et al. Modelling the impact of testing, contact tracing and household quarantine on second waves of COVID-19. Nat Hum Behav. 2020;4(9):964–71. pmid:32759985
  5. 5. Xu S, Li Y. Beware of the second wave of COVID-19. The Lancet. 2020;395(10233):1321–2. pmid:32277876
  6. 6. López L, Rodó X. The end of social confinement and COVID-19 re-emergence risk. Nat Hum Behav. 2020;4(7):746–55. pmid:32572175
  7. 7. Shunxiang Q, Xiang Z, Peter H, Nankun L, George FG, Yang S, et al. Two Reemergent Cases of COVID-19—Hebei Province, China, January 2, 2021. China CDC Weekly. 2021;3(2):25–7. pmid:34594899
  8. 8. Bisanzio D, Reithinger R, Alqunaibet A, Almudarra S, Alsukait RF, Dong D, et al. Estimating the effect of non-pharmaceutical interventions to mitigate COVID-19 spread in Saudi Arabia. BMC Med. 2022;20(1):51. pmid:35125108
  9. 9. Chung PC, Chan TC. Impact of physical distancing policy on reducing transmission of SARS-CoV-2 globally: Perspective from government’s response and residents’ compliance. PLoS One. 2021;16(8):e0255873. pmid:34375342
  10. 10. Gountas I, Quattrocchi A, Mamais I, Tsioutis C, Christaki E, Fokianos K, et al. Effect of public health interventions during the first epidemic wave of COVID-19 in Cyprus: a modelling study. BMC Public Health. 2021;21(1):1898. pmid:34666740
  11. 11. Li Y, Campbell H, Kulkarni D, Harpur A, Nundy M, Wang X, et al. The temporal association of introducing and lifting non-pharmaceutical interventions with the time-varying reproduction number (R) of SARS-CoV-2: a modelling study across 131 countries. Lancet Infect Dis. 2021;21(2):193–202. pmid:33729915
  12. 12. Sharma M, Mindermann S, Rogers-Smith C, Leech G, Snodin B, Ahuja J, et al. Understanding the effectiveness of government interventions against the resurgence of COVID-19 in Europe. Nat Commun. 2021;12(1):5820. pmid:34611158
  13. 13. Manica M, Guzzetta G, Riccardo F, Valenti A, Poletti P, Marziano V, et al. Impact of tiered restrictions on human activities and the epidemiology of the second wave of COVID-19 in Italy. Nat Commun. 2021;12(1):4570. pmid:34315899
  14. 14. Ge Y, Chen Z, Handel A, Martinez L, Xiao Q, Li C, et al. The impact of social distancing, contact tracing, and case isolation interventions to suppress the COVID-19 epidemic: A modeling study. Epidemics. 2021;36:100483. pmid:34284227
  15. 15. Lemey P, Ruktanonchai N, Hong SL, Colizza V, Poletto C, Van den Broeck F, et al. Untangling introductions and persistence in COVID-19 resurgence in Europe. Nature. 2021;595(7869):713–7. pmid:34192736
  16. 16. China CDC. COVID-19 Prevention and Control Protocol. 2020 [cited 2021 September 23]. Available from: http://www.nhc.gov.cn/jkj/s3577/202009/318683cbfaee4191aee29cd774b19d8d/files/f9ea38ce2c2d4352bf61ab0feada439f.pdf.
  17. 17. World Health Organization. Readiness for influenza during the COVID-19 pandemic 2020 [cited 2021 July 7]. Available from: https://www.who.int/publications/i/item/WHO-2019-nCoV-Influenza-readiness-COVID-19-2020.1.
  18. 18. Haug N, Geyrhofer L, Londei A, Dervic E, Desvars-Larrive A, Loreto V, et al. Ranking the effectiveness of worldwide COVID-19 government interventions. Nat Hum Behav. 2020;4(12):1303–12. pmid:33199859
  19. 19. 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. pmid:32512579
  20. 20. Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc Biol Sci. 2007;274(1609):599–604. pmid:17476782
  21. 21. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N Engl J Med. 2020;382(13):1199–207. pmid:31995857
  22. 22. Hu S, Wang W, Wang Y, Litvinova M, Luo K, Ren L, et al. Infectivity, susceptibility, and risk factors associated with SARS-CoV-2 transmission under intensive contact tracing in Hunan, China. Nat Commun. 2021;12(1):1533. pmid:33750783
  23. 23. Ferretti L, Wymant C, Kendall M, Zhao L, Nurtay A, Abeler-Dörner L, et al. Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing. Science. 2020;368(6491). pmid:32234805
  24. 24. Cori A, Ferguson NM, Fraser C, Cauchemez S. A new framework and software to estimate time-varying reproduction numbers during epidemics. American journal of epidemiology. 2013;178(9):1505–12. pmid:24043437
  25. 25. Leung K, Wu JT, Leung GM. Real-time tracking and prediction of COVID-19 infection using digital proxies of population mobility and mixing. Nat Commun. 2021;12(1):1501. pmid:33686075
  26. 26. Li B, Deng A, Li K, Hu Y, Li Z, Shi Y, et al. Viral infection and transmission in a large, well-traced outbreak caused by the SARS-CoV-2 Delta variant. Nat Commun. 2022;13(1):460. pmid:35075154
  27. 27. Thompson RN, Stockwin JE, van Gaalen RD, Polonsky JA, Kamvar ZN, Demarsh PA, et al. Improved inference of time-varying reproduction numbers during infectious disease outbreaks. Epidemics. 2019;29:100356. pmid:31624039
  28. 28. Yi H, Kaifa W, Wendi W. Analysis of transmissibility of COVID-19 and regional differences in disease control. Acta Mathematicae Applicatae Sinica. 2020;43(02):227–37. (in Chinese).
  29. 29. Kucharski AJ, Russell TW, Diamond C, Liu Y, Edmunds J, Funk S, et al. Early dynamics of transmission and control of COVID-19: a mathematical modelling study. Lancet Infect Dis. 2020;20(5):553–8. pmid:32171059
  30. 30. Wu W, Bai R, Li D, Feng A, Xu A, Lv J. Preliminary prediction of the epidemic trend of coronavirus disease 2019 in Guangdong Province (in Chinese). Jinan Univ (Nat Sci Med Ed). 2020;41(02):181–5.
  31. 31. Zhao S, Musa SS, Lin Q, Ran J, Yang G, Wang W, et al. Estimating the Unreported Number of Novel Coronavirus (2019-nCoV) Cases in China in the First Half of January 2020: A Data-Driven Modelling Analysis of the Early Outbreak. J Clin Med. 2020;9(2). pmid:32024089
  32. 32. Wang LP, Wang J, Zhao HY, Shi YY, Wang K, Wu P, et al. Modelling and assessing the effects of medical resources on transmission of novel coronavirus (COVID-19) in Wuhan, China. Math Biosci Eng. 2020;17(4):2936–49. pmid:32987508
  33. 33. Tuite AR, Fisman DN. Reporting, Epidemic Growth, and Reproduction Numbers for the 2019 Novel Coronavirus (2019-nCoV) Epidemic. Ann Intern Med. 2020;172(8):567–8. pmid:32023340
  34. 34. Imai N, Cori A, Dorigatti I, Baguelin M, Donnelly CA, Riley S, et al. Report 3: Transmissibility of 2019-nCoV. Imp Coll London; 2020.
  35. 35. 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. pmid:32014114
  36. 36. Mobility, commuting, and contact patterns across the United States during the COVID-19 outbreak. 2020 [cited 2022 June 26]. Available from: https://covid19.gleamproject.org/mobility.
  37. 37. Gimma A, Munday JD, Wong KLM, Coletti P, van Zandvoort K, Prem K, et al. Changes in social contacts in England during the COVID-19 pandemic between March 2020 and March 2021 as measured by the CoMix survey: A repeated cross-sectional study. PLoS Med. 2022;19(3):e1003907. pmid:35231023
  38. 38. Zhang J, Litvinova M, Liang Y, Zheng W, Shi H, Vespignani A, et al. The impact of relaxing interventions on human contact patterns and SARS-CoV-2 transmission in China. Sci Adv. 2021;7(19). pmid:33962957
  39. 39. Jilin Statistical Bureau. The 7th national population census bulletin of Jilin province [cited 2022 August 19]. Available from: http://tjj.jl.gov.cn/tjsj/qwfb/202105/t20210524_8079098.html. http://tjj.jl.gov.cn/tjsj/tjgb/pcjqtgb/202105/t20210524_8079042.html.
  40. 40. Shijiazhuang Statistical Bureau. The 7th national population census bulletin of Shijiazhuang city [cited 2022 August 19]. Available from: http://tjj.sjz.gov.cn/col/1584345215439/2021/05/31/1622426676587.html.
  41. 41. Suihua Statistical Bureau. The 7th national population census bulletin of Suihua city [cited 2022 August 19]. Available from: https://suihua.dbw.cn/system/2021/06/01/058654474.shtml.
  42. 42. Haerbin Statistical Bureau. The 7th national population census bulletin of Haerbin city [cited 2022 August 19]. Available from: http://www.harbin.gov.cn/art/2021/5/31/art_25924_1117402.html.
  43. 43. Hebei Statistical Bureau. Hebei Statistical Yearbook (2021) [cited 2022 August 19]. Available from: http://tjj.hebei.gov.cn/hetj/tjnj/2021/zk/indexch.htm.
  44. 44. Heilongjiang Statistical Bureau. Heilongjiang Statistical Yearbook (2021) [cited 2022 August 19]. Available from: http://tjj.hlj.gov.cn/tjjnianjian/2021/zk/indexch.htm.
  45. 45. Jilin Statistical Bureau. Jilin Statistical Yearbook (2021) [cited 2022 August 19]. Available from: http://tjj.jl.gov.cn/tjsj/tjnj/2021/ml/indexc.htm.
  46. 46. Davies NG, Kucharski AJ, Eggo RM, Gimma A, Edmunds WJ. Effects of non-pharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study. Lancet Public Health. 2020;5(7):e375–e85. pmid:32502389
  47. 47. Askitas N, Tatsiramos K, Verheyden B. Estimating worldwide effects of non-pharmaceutical interventions on COVID-19 incidence and population mobility patterns using a multiple-event study. Sci Rep. 2021;11(1):1972. pmid:33479325
  48. 48. Mendez-Brito A, El Bcheraoui C, Pozo-Martin F. Systematic review of empirical studies comparing the effectiveness of non-pharmaceutical interventions against COVID-19. J Infect. 2021;83(3):281–93. pmid:34161818
  49. 49. Keeling MJ, Hollingsworth TD, Read JM. Efficacy of contact tracing for the containment of the 2019 novel coronavirus (COVID-19). Journal of Epidemiology and Community Health. 2020;74(10):861–6. pmid:32576605
  50. 50. Griffin S. Covid-19: Test and trace programmes are important but no silver bullet, say scientists. BMJ. 2020;369:m2151. pmid:32467147
  51. 51. Contreras S, Dehning J, Loidolt M, Zierenberg J, Spitzner FP, Urrea-Quintero JH, et al. The challenges of containing SARS-CoV-2 via test-trace-and-isolate. Nat Commun. 2021;12(1):378. pmid:33452267
  52. 52. Adam D. A guide to R—the pandemic’s misunderstood metric. Nature. 2020;583(7816):346–8. pmid:32620883
  53. 53. Galvani AP, May RM. Epidemiology: dimensions of superspreading. Nature. 2005;438(7066):293–5. pmid:16292292