Main

After implementing the ‘zero-COVID’ policy for more than two years, China has recently begun to adjust its coronavirus disease 2019 (COVID-19) response strategies, notably by announcing the 20 measures on 1 November and further the 10 measures on 7 December 2022 (refs. 1,2). Although increasing vaccination coverage among the elderly and protection of high-risk groups were emphasized, the 20 measures included restricting testing coverage, cutting quarantine period for close contacts or inbound travelers and suspending the tracing of secondary contacts1. The 10 measures further prohibited region-wide mass testing, and home isolation or quarantine was allowed2. Since then, Omicron infections have spread rapidly in major cities in China, including Beijing, where the predominant SARS-CoV-2 variant, Omicron BF.7, has put great pressure on the healthcare system since early December3,4.

Regular mass testing and intensive contact tracing were suspended in late November, and nucleic acid testing has been conducted on a voluntary basis thereafter5. As such, the daily number of confirmed cases reported by Beijing Municipal Health Commission (http://wjw.beijing.gov.cn/) was no longer an accurate reflection of the epidemic, making it difficult to assess the transmission dynamics. Here, we tracked the effective reproduction number of Omicron BF.7 in Beijing in November–December 2022 (refs. 6,7) using our previous epidemic nowcast framework, which combined real-time mobility data and case data in the disease transmission models8.

Effective reproduction number (R t) of Omicron in Beijing

We categorized the COVID-19 response adjustments in Beijing into three stages:

  1. i.

    Stage 1 (1–11 November): Zero-COVID strategy was strictly implemented with regular mass testing, intensive contact tracing and lockdown of residential buildings or communities once polymerase chain reaction (PCR)-positive infections and their close contacts were traced.

  2. ii.

    Stage 2 (12–25 November): Although mass testing and contact tracing were maintained after the announcement of 20 measures on 11 November, lockdown was limited to the residential buildings or only the floors or units where people with PCR-positive infections lived.

  3. iii.

    Stage 3 (after 25 November): The requirements of regular mass testing, intensive contact tracing and lockdown were gradually relaxed and finally suspended on 30 November. Nucleic acid testing is conducted on a voluntary basis after the announcement of 10 measures on 7 December.

We parameterized the disease transmission model with the daily number of passengers from Beijing Mass Transit Railway (MTR) and Beijing Subway, which collectively manage 22 of 27 rapid transit lines in Beijing (Fig. 1). We fitted the model to two data streams to estimate the effective reproduction number Rt: (i) the daily number of symptomatic cases reported to Beijing Municipal Health Commission between 1 and 11 November (Stage 1), when the testing and reporting behavior were relatively stable and (ii) the proportion of participants who self-reported to have ever tested positive by PCR or rapid antigen test (RAT) since 1 November, based on Weibo and WeChat online polls conducted via convenience sampling between 10 and 22 December (Fig. 1 and Extended Data Table 1). See Methods for details.

Fig. 1: Data used in the inference.
figure 1

a, The daily number of reported symptomatic cases between 1 and 18 November (orange), the daily number of reported asymptomatic cases between 1 and 18 November (blue) and the proportion of participants of Weibo or WeChat polls who reported to have ever tested positive (PCR or RAT) between 10 and 22 December (green). The details of Weibo or WeChat polls are provided in Extended Data Table 1. The dots and bars in green indicate the mean and 95% confidence interval (CI) of the binomial proportions of the positive participants in the polls. b, The daily number of passengers of Beijing (BJ) subway lines managed by Beijing MTR (five lines, https://www.mtr.bj.cn/) and Beijing Subway (17 lines, https://www.bjsubway.com/). c, The daily number of subway passengers relative to that on 1 November. The daily number of passengers are provided in Extended Data Table 2.

Within one week of the announcement of 20 measures, Rt increased from 1.04 (95% CrI: 0.84–1.29) on 11 November to 3.44 (95% CrI: 2.82–4.14) on 18 November (Fig. 2). In response to the increasing number of cases, public health and social measures (PHSMs) were implemented: residents were urged to stay home over the weekend of 19–20 November; 95% of staff were advised to work from home during the week of 21–25 November and kindergartens and primary and secondary schools were closed on 21 November. Consequently, mobility decreased, and Rt dropped to 0.99 (95% CrI: 0.83–1.16) on 27 November.

Fig. 2: Estimated effective Rt, daily incidence and cumulative IAR in Beijing.
figure 2

a, Estimated Rt between 1 November and 22 December. Lines indicate the maximum a posteriori probability (MAP) estimates and shades indicate 95% CrI. b, Estimated daily number of infections between 1 November and 22 December. Lines indicate MAP estimates and shades indicate 95% CrI. c, Estimated cumulative infection attack rates between 1 November and 22 December. Lines indicate MAP estimates and shades indicate 95% CrI.

The surge of infections saturated the capacity of PCR testing and quarantine facility in late November. The requirement of regular mass testing, intensive contact tracing and lockdown were gradually relaxed and finally suspended with the announcement of 10 measures. PHSMs were relaxed and consequently Rt increased to 2.43 (95% CrI: 2.11–2.52) on 7 December (Fig. 2), which was substantially higher than Rt of 1.9 under similar PHSMs in the early stage of Hong Kong’s Omicron wave in February–March 20229.

Omicron infections spread rapidly again starting from early December, and many symptomatic individuals and their close contacts self-isolated10. Within one week after the announcement of 10 measures, Beijing’s mobility dropped to low levels, and we estimated that Rt dropped below 1 on 11 December and to 0.72 (95% CrI: 0.30–0.93) on 14 December. The mobility level started to increase in the week of 17–23 December, and Rt increased to 1.09 (95% CrI: 0.39–1.96) on 19 December. We estimated that Rt dropped below 1 again on 21 December due to depletion of the susceptible population (Fig. 2).

Daily incidence and cumulative infection attack rate

We estimated the daily incidence and cumulative infection attack rate from the fitted model accordingly (Fig. 2). On 30 November, when regular mass testing was suspended, we estimated the daily number of infections was 94,272 (95% CrI: 52,270–160,042) in Beijing. The daily incidence increased rapidly since then and peaked on 11 December, with an estimated daily number of infections of 1.03 million (95% CrI: 0.61–1.49; that is, 4.7% of the population).

We estimated that the cumulative IAR was 43.1% (95% CrI: 25.6–60.9) on 14 December. The daily incidence slightly decreased between 11 and 16 December but started to increase again on 17 December coincident with increased mobility. We estimated that the cumulative IAR reached 75.7% (95% CrI: 60.7–84.4) on 22 December. Assuming that the mobility and population mixing level would remain at the same levels after 22 December, we estimated the cumulative IAR would reach 92.3% (95% CrI: 91.4–93.1) on 31 January 2023 (Extended Data Fig. 1).

In the base-case scenario above, we assumed each participant of the online polls underwent multiple RATs and the ascertainment probability of previous infection was 1 after 11 November (Fig. 2). As such, the actual IAR was likely to be slightly higher than our base-case estimates. As a sensitivity analysis, we included the ascertainment probability in the inference framework (Extended Data Fig. 2). The resulting estimate of ascertainment probability was 81.5% (95% CrI: 78.9–84.2). The corresponding estimated Rt was slightly higher (for example, 3.92 (95% CrI: 3.19–4.73) versus 3.44 (95% CrI: 2.82–4.14) on 18 November and 3.05 (95% CrI: 2.63–3.24) versus 2.43 (95% CrI: 2.11–2.52) on 7 December). The estimated cumulative IAR was higher accordingly, reaching 55.1% (95% CrI: 34.3–73.5) on 14 December and 88.5% (95% CrI: 80.1–93.0) on 22 December. Assuming that the mobility and population mixing level would remain at the same levels after 22 December, we estimated the cumulative IAR would reach 96.2% (95% CrI: 96.1–96.3) on 31 January 2023.

Discussion

Our study tracked Rt of the Omicron BF.7 variant in Beijing by fitting an epidemic transmission model parameterized with mobility data to early-stage case count and recent survey data on cumulative incidence8. When mass testing and confirmation of all cases become impossible during the surge of infections, it is important to continuously monitor infection prevalence in the community through various surveillance programs, such as REACT-type studies in the United Kingdom and Hong Kong11,12, wastewater surveillance13 and serological surveillance14. Given the rollout of the second booster vaccines, data from such surveillance studies could also inform the assessment of vaccine effectiveness in real time12,14.

By 22 December, we estimated that the Omicron epidemic had peaked in Beijing and 76% of the Beijing residents had been infected (Fig. 2). Assuming no changes in PHSMs or population behavior, we anticipated that the number of infections would start dropping toward the end of December and the pressure on the healthcare system would be alleviated in January 2023 (Extended Data Fig. 1). However, vaccination should be ramped up in the coming months in anticipation of new waves that might arise due to importation of different Omicron variants from increased population movement between Beijing and other provinces during the Spring Festival in January 2023.

Our study has several limitations. First, passenger statistics of Beijing MTR and Beijing Subway were the only real-time mobility data we found. These passengers only account for about 95% of Beijing’s subway volume and might not be representative of the population using other transportation means. Second, the subway mobility data might not be the best proxy for changes in contacts and hence Rt. During the week of 8–14 December, most of the population stayed home or self-quarantined and only made contact with their household members15. Thus, it is difficult to accurately estimate the rapid changes in population mixing and PHSM effectiveness. Third, our estimation relied on IARs inferred from the results of Weibo and WeChat online polls. As most active Weibo or WeChat users are between 18 and 60 years old, the poll results might not be representative of all age groups. Self-reported infection status or test results can also introduce bias, and the ascertainment probability of infection was uncertain. Future refinements of our framework should include more reliable data to estimate IARs, such as age-specific seroprevalence data14

Given that COVID-19 response strategies have been adjusted nationwide, the number of new infections is expected to surge across China. The epidemics in other major cities and counties might be similar to that of Beijing with multiple incidence peaks, especially when population mixing starts to increase after the first peak which is further exacerbated by the upcoming mass population movement during the annual Spring Festival. Surveillance programs should be rapidly set up to monitor the spread and evolution of SARS-CoV-2 infections; further work should be done to track the transmissibility, incidence and IAR of the epidemics.

Methods

Data

Ethics statement

The study was based on anonymized data that are publicly available. No ethical approval was required for use of the data.

Mobility data from Beijing Subway and Beijing MTR

Beijing Mass Transit Railway Operation Corporation Limited (https://www.bjsubway.com/, Beijing Subway) and Beijing MTR Corporation Limited (https://www.mtr.bj.cn/, Beijing MTR) operate 22 major lines out of 27 rapid transit lines in Beijing, which serve about 95% of the subway passengers daily. Both Beijing Subway and Beijing MTR publish the recorded daily number of trips regularly (https://weibo.com/bjsubway and https://weibo.com/bjmtr) and we obtained the data between 1 November and 22 December 2022, from their Weibo accounts (Extended Data Table 2).

Data from online polls conducted via Sina Weibo or Tencent WeChat

Since the announcement of 10 measures on 7 December 2022, verified Weibo and WeChat users have conducted online polls to track the proportion of participants who self-reported to be ever positive by PCR or RAT since 1 November, in view of the rapid increase of Omicron infections around them.

Sina Weibo has a verification policy to confirm the identity of a personal user or organization, similar to Twitter’s account verification (https://kefu.weibo.com/faqdetail?id=37). On Weibo, verified users can create Weibo online polls (https://kefu.weibo.com/faqdetail?id=20533). Similar to Twitter polls, the default duration of Weibo polls is 1 day and the poll results are instantly displayed. Also, the poll creators are able to stratify the poll results by provinces inferred from the IP addresses where the participants voted.

Similarly, verified Tencent WeChat users can register as a public account, which enables them to interact with subscribers and provide them with services. WeChat online polls can be conducted in public accounts with results stratified by provinces inferred from IP addresses (https://kf.qq.com/faq/161221FN7fi6161221NRBJzq.html).

In this study, we included data from 22 cross-sectional online polls conducted on Weibo or WeChat by verified users whose polls were read more than 100,000 times nationwide within the day posted (Extended Data Table 1). Results of these polls were publicly available, and there were two to four options in most of the polls; for example:

  1. i.

    Two-option polls (such as https://weibo.com/chinaetfs, 11 December)

    1. a.

      Option 1: I have tested positive in this wave (recovered or not yet recovered)

    2. b.

      Option 2: I have not tested positive in this wave

  2. ii.

    Three-option polls (such as https://weibo.com/u/2987102112, December 10–14)

    1. a.

      Option 1: I have tested positive and recovered in this wave

    2. b.

      Option 2: I have tested positive but not yet recovered in this wave

    3. c.

      Option 3: I have not tested positive in this wave

  3. iii.

    Four-option polls (such as https://weibo.com/u/7126731879, 10 December)

    1. a.

      Option 1: I have not tested positive

    2. b.

      Option 2: I have tested positive and recovered

    3. c.

      Option 3: I have tested positive but not yet recovered

    4. d.

      Option 4: I have tested positive twice

Transmission model

We used our previous age-structured Susceptible-Exposed-Infectious-Removed (SEIR) model to simulate the transmission of Omicron in Beijing16:

$$\begin{array}{*{20}{l}} {\frac{{dS_a\left( t \right)}}{{dt}}} \hfill & = \hfill & { - S_a\left( t \right)\mathop {\sum}\limits_{b = 1}^m {\mathop {\sum}\limits_{i = 1}^k {\beta _{ab}\left( t \right)I_{b,i}\left( t \right)} } } \hfill \\ {\frac{{dE_a\left( t \right)}}{{dt}}} \hfill & = \hfill & { - \gamma _EE_a\left( t \right) + S_a\left( t \right)\mathop {\sum}\limits_{b = 1}^m {\mathop {\sum}\limits_{i = 1}^k {\beta _{ab}\left( t \right)I_{b,i}\left( t \right)} } } \hfill \\ {\frac{{dI_{a,1}\left( t \right)}}{{dt}}} \hfill & = \hfill & { - \gamma _II_{a,1}\left( t \right) + \gamma _EE_a\left( t \right)} \hfill \\ {\frac{{dI_{a,2}\left( t \right)}}{{dt}}} \hfill & = \hfill & { - \gamma _II_{a,2}\left( t \right) + \gamma _II_{a,1}\left( t \right)} \hfill \\ {} \hfill & {} \hfill & \ldots \hfill \\ {\frac{{dI_{a,k}\left( t \right)}}{{dt}}} \hfill & = \hfill & { - \gamma _II_{a,k}\left( t \right) + \gamma _II_{a,k - 1}\left( t \right)} \hfill \\ {\frac{{dR_a\left( t \right)}}{{dt}}} \hfill & = \hfill & {\gamma _II_{a,k}\left( t \right)} \hfill \\ {N_a} \hfill & = \hfill & {S_a\left( t \right) + E_a\left( t \right) + \mathop {\sum}\limits_{i = 1}^k {I_{a,i}\left( t \right) + R_a\left( t \right)} } ,\hfill \end{array}$$

where m was the number of age groups in the population; k=4 was the number of infectious states of the model; Sa(t), Ea(t) and Ra(t) were the number of susceptible, exposed and removed individuals in age group a at time t; Ia,i(t) was the number of individuals in the i-th infectious state of age group a at time t; Na was the total number of people in age group a; 1/γE was the duration from exposure to become infectious; and k/γI was the duration of being infectious. The mean generation time was calculated as \(T_{GT} = 1/\gamma _E + \frac{{k + 1}}{{2\gamma _I}}\), following the formulation of Svensson17.

Following our previous framework to parameterize SARS-CoV-2 transmission models with mobility data8, we formulated the average rate at which an individual in age group a made infectious contacts with age group b at time t (that is, βab(t) as

$$\beta _{ab}\left( t \right) = c_{ab}\gamma \left( t \right)\left( {1 - e^{ - g\left( t \right)}} \right)$$

where cab was the contact rates between age group a and age group b from Mistry et al18 and g(t) was the normalized number of passengers of subway lines managed by Beijing Subway and Beijing MTR on day t (such that g(t)=1 on 1 November 2022).γ(t) was the scaling factor for the mobility data proxy g(t). Given the changes in PHSMs during the weeks of 12–25 November, we assumed γ(t) was γ1 between 1 and 11 November, increased linearly between 12 and 25 November from γ1 to γ2, and remained at γ2 after 25 November. γ1 and γ2 were estimated in the parameter inference (Extended Data Tables 3, 4).

The time-varying next generation matrix for this SEIR model was:

$$\begin{array}{*{20}{c}} {NGM\left( t \right)} & = & {T_g\left[ {\begin{array}{*{20}{c}} {\frac{{\beta _{11}\left( t \right)S_1\left( t \right)}}{{N_1\left( t \right)}}} & \cdots & {\frac{{\beta _{1m}\left( t \right)S_1\left( t \right)}}{{N_1\left( t \right)}}} \\ \vdots & \ddots & \vdots \\ {\frac{{\beta _{m1}\left( t \right)S_m\left( t \right)}}{{N_m\left( t \right)}}} & \cdots & {\frac{{\beta _{mm}\left( t \right)S_m\left( t \right)}}{{N_m\left( t \right)}}} \end{array}} \right]} \end{array}$$

where Tg was the mean generation time.

The effective reproduction number Rt corresponded to the dominant eigenvalue of NGM(t)19,20. The incidence rate of infection and reported onsets at time t were calculated as follows:

$$A_{infection}\left( t \right) = \mathop {\sum}\limits_a {S_a\left( t \right)\pi _a\left( t \right)}$$
$$A_{onset}\left( t \right) = p_{report}\mathop {\int}\limits_0^t {f_{incubation}\left( {t - u} \right)\left( {\mathop {\sum }\limits_a A_{a,infection}(u)} \right)du},$$

where preport was the proportion of infections ascertained and reported as symptomatic cases by Beijing Municipal Health Commission (http://wjw.beijing.gov.cn/). Similarly, the cumulative IAR at time t was calculated as follows:

$$IAR\left( t \right) = \mathop {\int}\limits_0^v {A_{infection}(u)du}.$$

We assumed that the epidemic was seeded by M local infections on 1 November 2022.

Inference

The set of parameters that were subject to statistical inference, which we denoted by θ, included: (i) the seed size M, (ii) the scaling factors γ1 and γ2 and (iii) the ascertainment proportion preport between 1 and 11 November (Extended Data Table 3). We estimated θ from (i) the daily number of symptomatic cases reported to Beijing Municipal Health Commission between 1 November and 11 (Fig. 1) and (ii) the proportion of participants who reported to be positive by PCR or RAT since 1 November in Weibo or WeChat online polls conducted between 10 and 22 December (Extended Data Table 1).

The likelihood function L(θ) is a product of the two components L1(θ) and L2(θ). L1(θ) was formulated as below:

$$L_1\left( \theta \right) = \mathop {\prod }\limits_t \left( {\frac{{A_{onset}\left( t \right)!}}{{n_{case}\left( t \right)!\left( {A_{onset}\left( t \right) - n_{case}\left( t \right)} \right)!}}p_{report}^{n_{case}\left( t \right)}\left( {1 - p_{report}} \right)^{A_{onset}\left( t \right) - n_{case}\left( t \right)}} \right),$$

where ncase(t) was the daily number of symptomatic cases confirmed and reported by Beijing Municipal Health Commission between 1 and 11 November and Aonset(t) was daily number of infections from the model convoluted by the incubation period distribution, assuming no onset-to-confirmation delay under mass testing and intensive contact tracing between 1 and 11 November. Similarly, assuming each participant of the online polls underwent multiple RATs and the ascertainment probability of infection was 100% after 1 November, L2(θ) was formulated as:

$$L_2\left( \theta \right) = \mathop {\prod }\limits_i \left( {\frac{{n_{pol,i}!}}{{n_{pos,i}!\left( {n_{pol,i} - n_{pos,i}} \right)!}}IAR\left( {t_{pol,i}} \right)^{n_{pos},{{{\mathrm{i}}}}}\left( {1 - IAR\left( {t_{pol,i}} \right)} \right)^{n_{pol,i} - n_{pos,i}}} \right),$$

where tpol,i was the time or date when Weibo or WeChat online poll i was carried out, npol,i was the number of participants of Weibo or WeChat online poll i and npos,i was the number of participants of Weibo or WeChat online poll i who reported to have ever tested positive by PCR or RAT since 1 November.

However, the sensitivity of RAT ranges between 70% and 80% for Omicron and thus the ascertainment probability of infection could be less than 100% after 1 November21. In a sensitivity analysis, we assumed that the ascertainment probability of infection (pascertain) was not 100%, and estimated pascertain and other parameters with the updated L2(θ) as:

$$\begin{array}{l}L_2\left( \theta \right) = \mathop {\prod }\limits_i \Big( \frac{{n_{pol,i}!}}{{n_{pos,i}!\left( {n_{pol,i} - n_{pos,i}} \right)!}}\left( {p_{ascertain}IAR\left( {t_{pol,i}} \right)} \right)^{n_{pos,i}}\\\left( {1 - p_{ascertain}IAR\left( {t_{pol,i}} \right)} \right)^{n_{pol,i} - n_{pos,i}} \Big),\end{array}$$

L(θ) was formulated as follows:

$$L\left( \theta \right) = L_1\left( \theta \right) \times L_2\left( \theta \right).$$

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.