Risk assessment of the step-by-step return-to-work policy in Beijing following the COVID-19 epidemic peak.

Novel coronavirus (COVID-19) is a new strain of coronavirus first identified in Wuhan, China. As the virus spread worldwide causing a global pandemic, China reduced transmission at considerable social and economic cost. Post-lockdown, resuming work safely, that is, while avoiding a second epidemic outbreak, is a major challenge. Exacerbating this challenge, Beijing hosts many residents and workers with origins elsewhere, making it a relatively high-risk region in which to resume work. Nevertheless, the step-by-step approach taken by Beijing appears to have been effective so far. To learn from the epidemic progression and return-to-work measures undertaken in Beijing, and to inform efforts to avoid a second outbreak of COVID-19, we simulated the epidemiological progression of COVID-19 in Beijing under the real scenario of multiple stages of resuming work. A new epidemic transmission model was developed from a modified SEIR model for SARS, tailored to the situation of Beijing and fitted using multi-source data. Because of strong spatial heterogeneity amongst the population, socio-economic factors and medical capacity of Beijing, the risk assessment was undertaken spatiotemporally with respect to each district of Beijing. The epidemic simulation confirmed that the policy of resuming work step-by step, as implemented in Beijing, was sufficient to avoid a recurrence of the epidemic. Moreover, because of the structure of the model, the simulation provided insights into the specific factors at play at different stages of resuming work, allowing district-specific recommendations to be made with respect to monitoring at different stages of resuming work . As such, this research provides important lessons for other cities and regions dealing with outbreaks of COVID-19 and implementing return-to-work policies. © Springer-Verlag GmbH Germany, part of Springer Nature 2020.

Introduction

Novel coronavirus (COVID-19) is a new strain of coronavirus with a high transmission ability that has spread widely since January 2020 (North China Morning Post 2020; Zhu et al. 2020; Li et al. 2020; Sujath et al. 2020). It is recognized as the most serious respiratory virus since the 1918 H1N1 influenza pandemic in terms of public health threat (Dorigatti et al. 2020). Since 31 December 2019, and as of 11 August 2020, 20,075,600 cases of COVID-19 (in accordance with the applied case definitions in the affected countries) were reported, including 736,372 deaths (European Centre for Disease Prevention and Control 2020). With the virus has a foothold on every continent, the pandemic threat has already been realized (World Health Organization 2020a; Chen et al. 2020), and many countries have published serval pandemic-prevention measures thereafter in response, such as lockdown. However, in the post-pandemic world, the reopening of society is in strong need, and the corresponding risk assessment is necessary for megacities to carefully lift the issued measures avoiding the ongoing outbreaks of COVID-19.

As COVID-19 has spread worldwide, China has eased transmission of the virus at considerable social and economic cost. We focus on Beijing, China because it provides an informative case showing progress from pandemic outbreak to government curbs on the pandemic, and further to the resumption of work. In fact, Beijing started to resume work on 10 February 2020 (Xinhua net 2020a), and production levels were back to 88.9% as early as 15 March 2020 (Xinhua net 2020b). As one of the largest cities in China, Beijing has 21.5 million residents, of which 7.6 million are in-migrants (domestic migrants) (Beijing Municipal Bureau Statistics 2019). This flux of migrant workers makes Beijing a high-risk region in terms of resuming work. In addition, because of the large distances between workplaces and residences in the city and the possibility of asymptomatic transmission of the virus, the ubiquity of commuting raises transmission possibilities in Beijing (Chinazzi et al. 2020; Chen et al. 2020). To begin economic recovery while maintaining a low risk of a second outbreak of COVID-19, step-by-step measures to return to work were taken by the Beijing government. The resumption of work in Beijing could bring significant risks if not well controlled. Several existing studies showed that influenza spread has an intrinsic spatio-temporal behavior, and the diffusion of infected cases in a given area is likely to be uneven and determined by the distribution of susceptible hosts (Qian et al. 2009; Moustakas 2017; Moss et al. 2019; Amorós et al. 2020). For example, individuals in schools have many more contacts with other people than those who stay at home most of the time, so regions with more schools are at higher risk. In addition, mode of employment and the scale of one’s workplace also have a great impact on the probability of infection (Aschwanden 2004). Controlling for the number of susceptible hosts, k is larger for companies operating within more confined spaces than, for example, hotels and some public places, such as parks.

To analyze the goal of resuming social and economic operations while avoiding a re-outbreak of COVID-19, we conducted a simulation analysis of the risk of pandemic outbreak for Beijing at different stages of resuming work. The basic transmission model we employed is the modified SEIR model proposed by Lipsitch et al. (2003), revised according to the circumstances of the pandemic in Beijing represented using multi-source data; that is, population (residents and in-migrants), socio-economic factors and medical capacity. Although the economy is not the focus of pandemic models, it is appropriate to consider some economic impacts (Enserink and Kupferschmidt 2020), since transmitted infectious diseases can cause extensive morbidity, mortality and economic burden (Choi et al. 2008; Christakos et al. 2017). Heterogeneity exists across the districts of Beijing in the population distribution and local socio-economic factors (Song et al. 2019; Brunton et al. 2017). At different stages of resuming work, the business functions in operation vary (Beijing Youth Daily 2020), leading to different situations in different districts (i.e. different outbreak risk). Districts with large populations will not necessarily become high-risk regions, because the majority of residents will stay home for much of time if the district has few open workplaces. Thus, it is unreasonable to employ a single simulation model in developing pandemic-prevention measures for all districts (Yin et al. 2019; Chien and Chen 2020). From this perspective, we study the risk of resuming work for each district based on the spatially heterogeneous characteristics of the population, socio-economic factors and medical capacity. In particular, two parameters of the transmission model are evaluated locally: the daily contact number per capita and the mean daily rate at which infectious cases (who are not in quarantine) are detected and isolated. Each is evaluated using district-specific socio-economic factors and the medical capacity associated with the local population. Then, we produce a spatial–temporal risk assessment based on the district-specific parameters and propose relevant pandemic-prevention measures. The overall aim is to provide valuable insights for other cities or regions anticipating the journey from pandemic outbreak to transmission prevention to resumption of socio-economic activities.

Methodology and materials

Revised SEIR model

COVID-19 is a disease caused by the virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Langousis and Carsteanu 2020; Sivakumar 2020) which is genetically related to the coronavirus responsible for the SARS outbreak of 2003 (World Health Organization 2020b). To model the epidemiological progression of COVID-19, the transmission model for SARS (Lipsitch et al. 2003) was adapted. The transmission model for SARS is a modification of the standard SEIR model (Anderson and May 1991). The standard SEIR model tracks susceptible, exposed (infected but not yet infectious), infectious and recovered individuals in the compartments S, E, I and R, respectively, whereas the Lipsitch transmission model further incorporates quarantine measures into the SEIR model. Specifically, susceptible individuals are divided into ‘susceptible’ and ‘susceptible but isolated’ (compartments and ), exposed individuals are divided into ‘exposed’ and ‘exposed but isolated’ (compartments and ) and infectious individuals are divided into ‘those who have not been isolated’ (), and ‘those who have been hospitalized’ (). Based on the evidence for indirect transmission resulting from transmission by asymptomatic infected persons and pre-symptomatic transmission (Cai et al. 2020; Tong et al. 2020; Bai et al. 2020; Nishiura et al. 2020; Rothe et al. 2020), we further modified the model by assuming that a fraction of all infectious persons are asymptomatic ( of all infectious persons are symptomatic), that a proportion of people who make contact with an asymptomatic person are infected, and a proportion of people who make contact with a symptomatic person are infected. In this situation, the probability of a susceptible person becoming infected by contacting an infectious person is , and the probability of a susceptible who makes contact not becoming infected is . For simplicity, we use the notation to represent and to represent . In addition, we model the situation in which Beijing faces the return of considerable numbers of people from elsewhere in the country; thus, two parameters and are incorporated into the model, representing imported susceptible people and imported exposed people, respectively. That is,

It should be noted that the superscript in and indicates the stage of resuming work. For different stages of resuming work, there will be different populations of imported migrants. In this study, the resumption rate is reflected by the corresponding proportion of imported migrants. Moreover, during the specific stage of resuming work, the daily contact number per capita is determined by the corresponding reopened socio-economic factors, denoted by . That is the parameter in the revised SEIR model varies with the resumption rate, because reopened businesses at different stages of resumption have different abilities to attract people. At the initial stage of resuming work, there were only essential services reopened. To reflect the ability of these reopened essential services attracting people, we used the density of POI regarding the essential services to characterize the value of the parameter for this initial stage. The mean daily rate at which infectious cases are detected and isolated is determined by the medical capacity which is denoted by . This parameter is constant because medical resources were at full operation in the initial stages of resumption. The other parameters in Eq. (1) are independent of resumption rate, and are also considered constant in this study. Figure 1 illustrate the mechanism for our transmission model.

Fig. 1

Schematic representation of the revised transmission model

The classic SEIR model is accepted widely as a basic model to describe the properties of COVID-19 transmission. Hou et al. (2020) and He et al. (2020a) employed a mixed SEIR compartmental model assuming that susceptible individuals contacting an exposed individual would probably become infected. Yang et al. (2020) modified the original SEIR model, introducing move-in and move-out parameters, because the initial outbreak of COVID-19 in China occurred during the Spring Festival. At that time, hundreds of millions of people moved between the provinces of China. Despite these modifications some properties of COVID-19 transmission, particularly pandemic-prevention measures, were not considered. Therefore, here, a quarantined component is added into the classic SEIR model as a fifth element in addition to the four elements: S (susceptible), E (exposed), I (infectious) and R (recovered) (Lopez and Rodó 2020). Furthermore, He et al. (2002b) added another element, hospitalized, representing the part of the infectious population that cannot contact other people. Infectious individuals still active in society and potentially spreading the virus should be distinguished from hospitalized infectious individuals. Wei et al. (2020) considered a more detailed scenario, where asymptomatic infected persons are introduced. In contrast, our model incorporates the above considerations into the classic SEIR model with some simplifications providing a comprehensive characterization of the properties of COVID-19 transmission in Beijing.

Estimating the reproductive number

In the absence of quarantine, the reproductive number for this model is given by , where is the mean duration of infectiousness, . The effect of quarantining a fraction q of contacts of infectious persons is to multiply this expression by (Lipsitch et al. 2003). On the other hand, the reproductive number can be calculated deterministically by the formula , where is the ratio of the infectious period to the serial interval. is the sum of the mean infectious and mean latent periods, and is the exponential growth rate of the cumulative number of cases in the pandemic (Lipsitch et al. 2003).

Data collection and preprocessing

Pandemic data for Beijing were collected from the website of the Beijing Municipal Health Commission (http://wjw.beijing.gov.cn/xwzx_20031/wnxw/) including daily confirmed cases, daily recovered cases and daily death cases for Beijing Municipality, as well as the total confirmed cases for each district of Beijing Municipality. The Beijing pandemic was first reported on 20 January 2020. To initialize the model, we assume that there were only exposed individuals at that time as a result of passenger inflow from Wuhan direct to Beijing. After the travel ban from Wuhan, the growth of the pandemic in Beijing can be viewed as based on internal progression within Beijing without further imported cases (Tian et al. 2020).

With Beijing reopened, socio-economic factors and medical capacity within the different districts of Beijing play important roles in the evolution of the epidemic. To better describe the evolution of the epidemic in Beijing, we collected multi-source data with which to fit the revised SEIR model. Socio-economic factors and medical capacity are described by the densities (pcs per square kilometer) of socio-economic activity places (SEAPs) and medical services, respectively, within each district of Beijing. The density of SEAPs for each district is measured by the density of POIs (points of interest) collected from the Gaode Map API (https://lbs.amap.com/api/webservice/guide/api/search) including 10 major categories: Accommodation Services, Governmental Organization, Financial and Insurance Services, Commercial Business, Enterprises, Shopping, Food and Beverages, Transportation Services, Science/Culture and Education Services, and Sports and Recreation. Medical Services per se is a major category of the POI data. Thus, the corresponding densities were used directly to describe the medical capacity for each district.

With respect to SEAPs, each major category consists of related granular categories, each of which has different population flows. Therefore, we assigned weights to the granular categories to calculate the density of each major category, specifically as the weighted average of the densities of the corresponding granular categories (see Table 5 in “Appendix”). For example, the major category Accommodation Services includes two granular categories, Hotel and Hostel. Because Hotel had a smaller population flow than Hostel, the weights of Hotel and Hostel were set as 0.4 and 0.6, respectively. It should be noted that the weights here represent the relative, rather than the absolute, magnitude. The densities of the other nine major categories were obtained in the same way. Next, the density of SEAPs was obtained by the weighted average of the densities of the 10 major categories, wherein the weights of the major categories were pre-defined according to their relative population flows also (see Table 5 in “Appendix”). Medical capacity was quantified by the density of medical services consisting of four granular categories in a similar way, but with the weights for the four granular categories of medical services pre-defined based on their theme relevance to the epidemic (see Table 5 in “Appendix”). For example, the largest weight was assigned to the granular category, Special Hospital, because one of the themes of it is infectious disease. The POI dataset was generated in 2017, while the acquisition time was March 2020. Population numbers of both residents and in-migrants for each district of Beijing were collected from the Beijing Statistical Yearbook in 2019 (Beijing Municipal Bureau Statistics 2019).

We employed non-linear least squares for parameter estimation by minimizing the distance between the predicted , , dead and corresponding observed numbers from 20 January 2020 to 20 March 2020. With respect to the parameters, the duration of quarantine, the average time of progression from latent infection to infectious, the proportion of asymptomatic cases and the per capita death rate were determined according to relevant reports (World Health Organization 2020c; Lauer et al. 2020; Mizumoto et al. 2020; Xia et al. 2020). The imported exposed population, the daily number of contacts per capita, the fraction of all persons contacted by an infectious person who are successfully quarantined, the per capita recovery rate, the transmission probability and the mean daily rate at which infectious cases are detected and isolated are estimated and given in Table 1. For each district of Beijing, we assumed that all parameters were the same as for the greater Beijing Municipality, except the daily contact number , detection rate and quarantine fraction . The first two parameters are considered to reflect the spatial heterogeneity of people’s behavior and the medical capacity of each district. Districts with more companies and schools, for example, tend to have larger daily contact numbers. To reduce the uncertainty of the model, we employed socio-economic factors and medical capacity in the evaluation of the parameters and , respectively. The socio-economic factors and medical capacity were determined by the POI data. That is the parameters and were determined by the densities of SEAPs and medical services regarding the whole Beijing Municipality. Then, the quarantine fraction of each district was estimated based on local pandemic data with a fixed daily contact number and detection rate . Under the circumstances, we set the local resident population as the susceptible population for each district.

Table 1

Initial conditions and parameters for the Beijing pandemic

Symbol	Parameter	Value (α=0.05)	Source
S_Q	Susceptible but quarantined	0	Assumption
E	Exposed-infected but not yet infectious	70	Assumption
E_Q	Exposed but quarantined	0	Assumption
I_U	Infectious who are not yet detected and isolated	0	Assumption
1/r_q	Duration of quarantine	14	Reference
1/r	Average time of progression from latent infection to infectious	5	Reference
m	The per capita death rate	0.0009	Reference
k	The daily number of contacts per capita	5.4668	POI
w	The mean daily rate at which infectious cases are detected and isolated	0.8579	POI
q	The fraction of all person contacted by an infectious person are successfully quarantined	0.99±0.0847	Estimated
v	The per capita recovery rate	0.032±0.0029	Estimated
b	The probability of transmission per contact between a susceptible and an infectious person	0.6±0.1626	Estimated

For simplicity, the districts of Beijing Municipality are named as follows: DC (Dongcheng District); XC (Xicheng District); CY (Chaoyang District); FT (Fengtai District); SJS (Shijingshan District); HD (Haidian District); MTG (Mentougou District); FS (Fangshan District); TZ (Tongzhou District); SY (Shunyi District); CP (Changping District); DX (Daxing District); HR (Huairou District); PG (Pinggu District); MY (Miyun District); YQ (Yanqing District).

Results

Chronology of resumption of socio-economic activities

To slow the pandemic, migrant residents returned to Beijing batch-by-batch in terms of their profession. Based on the timeline of resuming work (see Fig. 2), we simplify the situation into six scenarios with resumption rates from 50 to 100%:

Fig. 2

The timeline for resuming work in Beijing. EDZ (Economic Development Zone)

Resumption of work at 50% by 20 February 2020, including Accommodation Services, Governmental Organization, Financial and Insurance Services, Commercial Business, and Enterprises;

Resumption of work at 60% by 1 March 2020. In addition to the five types of workplace mentioned above, Shopping was added, reflecting reports that around 90% of shops and markets reopened on 23 February;

Resumption of work at 70% by 10 March 2020. In this scenario, we assume no extra categories of business resume;

Resumption of work at 80% by 20 March 2020. Food and Beverages were added to the resuming businesses, since it is reported that around half of the restaurants in the Chaoyang District reopened by 23 March 2020;

Resumption of work at 90% by 30 March 2020. Transportation Services were added to the resuming businesses;

Resumption of work at 100% by 10 April 2020. In this scenario, we assume the resumption of all business categories, where Science/Culture and Education Services and Sports and Recreation are the final two categories reopened. These two public services are more confined to places with high population mobility.

At different stages of resuming work, the factors contributing to the contact number are selected based on the process in Fig. 2, leading to a heterogeneous situation (see Table 2). At a certain stage of resuming work, the parameter k is the weighted average of the densities of the selected categories of POI data within each district. The baseline represents the circumstance under the Spring Festival, when local residents stay at home with only three major categories of POI opened (i.e. Governmental Organization, Commercial Business, and Accommodation Services). At this stage, the contact number in each district is contributed by the corresponding density of open business functions and housing units. For example, within Dongcheng District, the value of parameter k is the weighted average of the densities of Governmental organization, Commercial Business, and Accommodation Services, where the weights are pre-defined (see Table 5 in “Appendix”). Then, for each stage of resuming work, the parameter k for each district is calculated in the same way with corresponding reopened major categories. The parameter w represents the medical capacity to characterize the mean daily rate at which infectious cases are detected and isolated. In this research, medical capacity was captured by the density of Medical Services which can also be estimated from the POI data. It should be noted that Medical Services have been operating at full capacity through the period of this study. Therefore, the parameter w is a constant amid resuming work. As in Sect. 2.3, the parameter w for each district equals to the density of Medical Services within each district.

Table 2

Two district-specific SEIR model parameters at different stages of resuming work

District	k							w
	Baseline	50%	60%	70%	80%	90%	100%
DC	6.814	17.500	18.485	18.485	22.611	22.625	25.294	0.694
XC	6.166	14.355	15.431	15.431	18.929	18.939	21.503	0.702
CY	1.204	4.612	4.897	4.897	6.357	6.361	7.264	0.130
FT	0.982	3.094	3.369	3.369	4.295	4.300	4.853	0.115
SJS	1.170	3.000	3.180	3.180	3.908	3.912	4.625	0.123
HD	1.135	3.887	4.061	4.061	4.982	4.985	6.128	0.131
MTG	0.052	0.080	0.083	0.083	0.100	0.100	0.117	0.003
FS	0.073	0.166	0.180	0.180	0.229	0.230	0.276	0.012
TZ	0.213	0.744	0.791	0.791	0.999	1.000	1.163	0.020
SY	0.174	0.461	0.493	0.493	0.612	0.613	0.696	0.021
CP	0.172	0.485	0.525	0.525	0.711	0.712	0.881	0.025
DX	0.200	0.799	0.841	0.841	1.043	1.043	1.171	0.024
HR	0.057	0.105	0.111	0.111	0.134	0.134	0.156	0.006
PG	0.093	0.185	0.197	0.197	0.229	0.230	0.277	0.009
MY	0.035	0.061	0.065	0.065	0.077	0.088	0.100	0.004
YQ	0.023	0.033	0.035	0.035	0.044	0.050	0.058	0.002

The impact of resuming work on virus transmission

Under each scenario, the estimated total confirmed cases from 20 February to 10 April 2020 are distributed unevenly due to the heterogeneity of the multi-source data used and the step-by-step return-to-work policies. After resuming work, the pandemic increases for a short period in most districts, with the central region showing the most significant increase because of the high density of enterprises and transportation links.

To characterize the real situation more closely, we introduced more parameters to the model. The model fit for the Beijing pandemic is satisfactory (see Fig. 3a), although it should be noted that the uncertainty for some parameters is large because imported cases are not considered by the modified SEIR model. The simulated total number of confirmed cases revealed Chaoyang District to be the highest risk region, while Fengtai District produced the most cases up to 1 Aug. 2020 (see Fig. 3b). Besides, the large error associated with Fengtai District in 24 Jun. is also caused by imported cases. In fact, the pandemic in Beijing appears to have a short duration. In early March 2020, the total confirmed cases within each district gradually becomes stable, and the new confirmed cases thereafter are all related to imported cases.

Fig. 3

a–c Comparison between the fitted model and observed data with respect to Beijing. Id represents infectious individuals who are detected and isolated, and R represents recovered individuals. d Validation of the estimated total number of confirmed cases with respect to each district in Beijing. The error is calculated into absolute value to represent the size of it

Different from the total number of confirmed cases, the number of active cases is required to assess the current risk of COVID-19. Active cases refer to undetected infectious individuals, , who can still come into contact with susceptible individuals, such that these susceptible people may become infected. The total number of confirmed cases consists of existing confirmed, recovered and dead cases, where clearly the latter two bring no risk to the susceptible population. The simulated evolution of active cases of COVID-19 for each district of Beijing under different stages of return-to-work is demonstrated in Fig. 4a–d. Different from a simple cumulative process of total confirmed cases, confirmed active cases behave more dynamically within the process of resuming work. For example, the number of active cases in Huairou District increases at the start of resuming work, and gradually decreases in the post-resuming work stage. At the different stages of resuming work, the risk comes from different resuming business categories, influencing the spatial distribution of k and the infectious population. Enterprises and factories reopened first, as they were viewed as imperative for the domestic economy, followed by markets, the catering industry and public transportation. About 60% of economic operations resumed by 1 March 2020, such that attention should be paid to towns with the most densely located shopping services. Work resumed by 80% by 20 March, such that towns with the most densely located food and beverage operations should also be monitored. When transportation returned from lockdown to full capacity, towns with more transportation services should be included for monitoring. Finally, all operations resumed by 10 April 2020, such that science/culture and education services, especially schools in Haidian District and Chaoyang District, should be included for close monitoring. In the post-COVID-19 future in Beijing, we find that the epidemic can be contained (see Fig. 4e–g) based on our simulation with existing epidemic-prevention measures published in Beijing. Figure 4g illustrates that there are five districts with moderate risk, seven districts with mild risk and four districts with low risk on 30 May. Southern Beijing poses a higher risk than northern Beijing, because of a high density of SEAPs.

Fig. 4

Spatio-temporal prediction of the number of active cases in Beijing before and after resuming work. a–d active cases during resuming work at the 50%, 70%, 90% and 100% stages, respectively. e–g active cases in the post-pandemic future with society in full operation

Reproductive number is an intrinsic property which characterizes the evolution of an epidemic. Table 3 shows that the reproductive number in each district was generally lower than 1 during the resumption of work, estimated by the district-specific transmission model for quarantine. Amongst the top 10 provinces and cities with the largest numbers of total confirmed cases, Beijing moved downward in the rank list quickly (Kraemer et al. 2020). However, the risk assessment provided by the revised SEIR model provides only an empirical estimate. The small reproductive numbers for all the districts of Beijing represent that the epidemic was dying at different stages of resuming work with the step-by-step return-to-work policy in Beijing. In fact, the deterministic reproductive number estimated by the pandemic data represents a similar risk pattern. The deterministic reproductive number decreases to null, because there is no evidence that the epidemic remains active in Beijing except through imported active cases. Besides, the deterministic and model-based reproductive numbers cannot be directly compared by the size of values. The minimum deterministic reproductive number is 1 if there are any new cases. In contrast, the model-based reproductive number can close to zero when there still are new cases. It should be noted that the model-based reproductive number tends to increase with progress towards resuming work because of the increasing contact number, but the reproductive numbers are all lower than 1 because of the quarantine, which means that the epidemic is still under control. Nevertheless, the model-based reproductive number varies in the districts of Beijing because of the heterogeneous medical capacity, although the cumulative number of cases in each district is almost static.

Table 3

District-specific reproductive numbers at different stages of resuming work

District	The reproductive number R
	Deterministic						Model-based
	50%	60%	70%	80%	90%	100%	50%	60%	70%	80%	90%	100%
DC	1.47	1	1	Null	Null	Null	0.28	0.29	0.29	0.36	0.36	0.40
XC	1.77	Null	Null	Null	Null	Null	0.22	0.24	0.24	0.30	0.30	0.34
CY	1.67	2.33	Null	Null	Null	Null	0.32	0.34	0.34	0.44	0.44	0.51
FT	1.61	1.36	Null	Null	Null	Null	0.24	0.26	0.26	0.33	0.33	0.37
SJS	1.43	Null	Null	Null	Null	Null	0.22	0.23	0.23	0.29	0.29	0.34
HD	1.68	1	Null	1	Null	Null	0.27	0.28	0.28	0.35	0.35	0.43
MTG	1.17	Null	Null	Null	Null	Null	0.02	0.03	0.03	0.03	0.03	0.04
FS	1.44	Null	Null	Null	Null	Null	0.04	0.04	0.04	0.06	0.06	0.07
TZ	1.47	Null	Null	Null	Null	Null	0.16	0.17	0.17	0.21	0.21	0.24
SY	1.37	Null	Null	Null	Null	Null	0.10	0.10	0.10	0.13	0.13	0.15
CP	1.54	Null	Null	Null	Null	Null	0.09	0.10	0.10	0.14	0.14	0.17
DX	1.58	Null	Null	Null	Null	Null	0.16	0.17	0.17	0.21	0.21	0.23
HR	1.31	Null	Null	Null	Null	Null	0.03	0.03	0.03	0.04	0.04	0.04
PG	Null	Null	Null	Null	Null	Null	0.20	0.22	0.22	0.25	0.25	0.30
MY	1.31	Null	Null	Null	Null	Null	0.02	0.02	0.02	0.02	0.03	0.03
YQ	1	Null	Null	Null	Null	Null	0.01	0.01	0.01	0.01	0.02	0.02

The deterministic R is null means there is no case identified during the related period

Discussion and conclusion

Faced with the COVID-19 outbreak, a variety of different stringent measures were taken in a short time by the Chinese government (Tian et al. 2020a, b; Ferretti et al. 2020). From 23 January, a total of 31 provinces implemented a public health level 1 response, shutting down businesses and industry, cancelling mass gatherings and restricting movement. Medical resources were also sent from other regions of China to support Wuhan. As the capital city of China, Beijing implemented even stricter measures than other provinces, which were successful in preventing and controlling the COVID-19 infection.

Overall, the model-based reproductive number suggests that the quarantine and social distancing measures were effective in curbing the epidemic, while the deterministic reproductive number indicates that the epidemic is not yet dismissed in Beijing and areas with lesser medical capacity should be monitored particularly closely. To support economic recovery while minimizing disease-related risk from the outbreak, stringent measures should be taken by local governments in different regions. Regions with high density public transportation services are suggested to strengthen public guidance, control the number of people entering stations, increase the use of temperature detection devices in stations and clean train and bus carriages before each trip. Regions with densely distributed indoor public places, such as hotels and restaurants, should pay special attention to air circulation, and surface cleaning and disinfection. In enterprises and Commercial Business, monitoring of individuals’ temperatures, strict control of the numbers of workers in confined spaces and installation of windows that can be opened for regular ventilation is encouraged (Wu et al. 2020; Hellewell et al. 2020).

While there has been encouraging news that no new infections inside China are being reported, other parts of the world are becoming the new epicenters of the outbreak. Clearly, China is at a different stage in the COVID-19 pandemic cycle than many other major economies. At this time, the world’s mega-cities such as New York, Paris and London are coping with the impacts of the first wave of the pandemic. The experience of Beijing should be of great significance to other countries where the pandemic curve is still rising, especially developing regions with poor health services (Simon et al. 2005).

Efforts towards economic recovery will be required for all the cities suffering from the pandemic, especially mega-cities with dense industrial and populated areas. Evidence from Beijing shows that step-by-step measures can support economic recovery while avoiding a second pandemic outbreak. Beijing started resuming work on 10 February, when the number of infected cases was still rising, so only some major projects were resumed at that time. About ten days later, by February 19, some 61% of companies and industries resumed work, as the number of confirmed cases in Beijing reached 400 and stabilized. Meanwhile, some markets and shops reopened, to meet the basic requirements of daily life for people in Beijing. By 15 March, about 89% of companies and industries resumed work. Thus, the rate of resuming work in Beijing was about 10% every ten days. The model simulation undertaken here demonstrates that this pace of resuming work can be relatively effective, allowing control of the number of infections while restarting the economy as soon as possible. We hope that cities in other countries find this research useful in implementing their own policies to resume economic activity in the future.

Table 4

Top 10 (from top to bottom) towns with the most densely distributed POIs

Food and beverages	Enterprise	Shopping	Transportation	Financial and insurance	Science/culture and education	Commercial house	Sports and recreation	Government	Accommodation
Chaoyangmen, DC	Jianwai,CY	Chongwenmenwai,DC	Jinrongjie,XC	Chaoyangmen, DC	Xiangheyuan,CY	Chaoyangmen, DC	Sanlitun,CY	Tiyulu,DC	Chongwenmenwai,DC
Jianwai,CY	Chaoyangmen, DC	Hujialou, CY	Dongzhimen, DC	Jianwai,CY	Zhongguancun, HD	Andingmen, DC	Chongwenmenwai,DC	Jinrongjie,XC	Dazhanlan,XC
Dazhanlan,XC	Dongzhimen, DC	Zhanlanlu, XC	Chongwenmenwai,DC	Jinrongjie,XC	Jianwai,CY	Dongzhimen, DC	Jianwai,CY	Longtan,DC	Jingsong,CY
Andingmen,DC	Hujialou, CY	Xichanganjie, XC	Jianguomen,DC	Hujialou, CY	Beixiaguan, HD	Jinrongjie,XC	Hujialou, CY	Chongwenmenwai,DC	Qianmen,DC
Sanlitun,CY	Zhongguancun, HD	Dazhanlan,XC	Chaowai,CY	Chongwenmenwai,DC	Chaoyangmen, DC	Pinggu town,PG	Chaowai,CY	Niujie,XC	Jianguomen,DC
Chaowai,CY	Chaowai,CY	Dahongmen,FT	Guanganmenwai,XC	Chaowai,CY	Beitaipingzhuang, HD	Beixinqiao, DC	Dongzhimen, DC	Yuetan,XC	Jianwai,CY
Chongwenmenwai,DC	Chongwenmenwai,DC	Chaowai,CY	Fangzhuang,FT	Dongzhimen, DC	Chongwenmenwai,DC	Jiaodaokou, DC	Andingmen, DC	Donghuashi, DC	Anzhen,CY
Beixinqiao,DC	Haidian, HD	Chunshu,XC	Chunshu,XC	Beixinqiao,DC	Dongzhimen, DC	Chaowai,CY	Chaoyangmen, DC	Beixinqiao, DC	Dongzhimen, DC
Jiaodaokou,CY	Shangdi, CY	Andingmen,DC	Zhongguancun, HD	Dongsi,DC	Yanyuan, HD	Hujialou, CY	Donghuashi, DC	Chaoyangmen, DC	Jiaodaokou, DC
Jiuxianqiao,CY	Beixinqiao, DC	Wangsiying,CY	Jianwai,CY	Jianguomen,DC	Haidian, HD	Xinjiekou,XC	Jiaodaokou, DC	Beixiaguan, HD	Shuangjing,CY

Table 5

The weight for each mid categories and big categories respectively

Big category	Weight	Mid category	Weight
Food & Beverages	0–1	Chinese Food Restaurant	0–2
		Foreign Food Restaurant	0–1
		Fast Food Restaurant	0–1
		Leisure Food Restaurant	0–1
		Coffee House	0–1
		Tea House	0–1
		Icecream Shop	0–1
		Bakery	0–1
		Dessert House	0–1
Shopping	0–05	Shopping Plaza	0–15
		Convenience Store	0–05
		Home Electronics Hypermarket	0–05
		Supermarket	0–05
		Plants & Pet Market	0–05
		Home Building Materials Market	0–05
		Comprehensive Market	0–15
		Stationary Store	0–05
		Sports Store	0–05
		Commercial Street	0–15
		Clothing Store	0–05
		Franchise Store	0–05
		Special Trade House	0–05
		Personal Care Items Shop	0–05
Sports & Recreation	0–05	Sports Stadium	0–3
		Golf Related	0–1
		Recreation Center	0–3
		Holiday & Nursing Resort	0–05
		Recreation Place	0–15
		Theatre & Cinema	0–1
Accommodation Service	0.1	Hotel	0–4
		Hostel	0–6
Commercial Business	0–1	Industrial Park	0–4
		Building	0–5
		Residential Area	0–1
Governmental Organization	0–13	Social Groups	0–2
		Governmental Organization	0–5
		Foreign Organization	0–1
		Democratic Party	0–2
Science/Culture & Education Service	0–14	School	0–5
		Research Institution	0–1
		Training Institution	0–3
		Driving School	0–1
Transportation Service	0–17	Airport Related	0–2
		Railway Station	0–2
		Coach Station	0–2
		Subway Station	0–1
		Light Rail Station	0–1
		Bus Station	0–1
		Commuter Bus Station	0–1
Finance & Insurance Service	0–03	Finance & Insurance Service Institution	0–15
		Bank	0–15
		ATM	0–1
		Insurance Company	0–2
		Securities Company	0–2
		Finance Company	0–2
Enterprises	0–13	Enterprises	0–1
		Famous Enterprise	0–3
		Farming,Forestry,Animal Husbandry and Fishery Base	0–1
		Company	0–4
		Factory	0–1
Medical Service	1*	Hospital	0–3
		Special Hospital	0–4
		Clinic	0–2
		Emergency Center	0–1

*Note that only one category contributes to the parameter , thus the weight is 1

Table 6

The top-level categories contributing to the contact number of each district at different stages of resuming work

Big category	The level of resuming work (%)
	baseline	50	60	70	80	90	100
Food & Beverages					√	√	√
Shopping			√	√	√	√	√
Sports & Recreation							√
Accommodation Service	√	√	√	√	√	√	√
Commercial House	√	√	√	√	√	√	√
Governmental Organization & Social Group	√	√	√	√	√	√	√
Science/Culture & Education Service							√
Transportation Service						√	√
Finance & Insurance Service		√	√	√	√	√	√
Enterprises		√	√	√	√	√	√