Abstract

Coronavirus disease (COVID-19) cases and COVID-19-related deaths have been increasing worldwide since the outbreak in 2019. Before the mass vaccination campaign for COVID-19, the main methods for COVID-19 control in China were mass testing and quarantine. Based on the transmission mechanism of COVID-19, we constructed a dynamic model for COVID-19 transmission in two typical regions: Beijing and Xinjiang. We calculated the basic reproduction number , proved the global stability of COVID-19 transmission via the Lyapunov function technique, and introduced the final size. We assessed the effectiveness of mass testing and quarantine. Sensitivity analysis indicated that the more the people were tested per day, the larger is the quarantine proportionality coefficient, the earlier the source location was determined, and the better is the controlling effect. In addition, it was more effective to increase the coefficient of quarantine if the population density in the region was low. To eliminate the pandemic, the government has to expand testing and quarantine, requiring a large amount of continuous manpower, material, and financial resources. Therefore, new control measures should be developed.

1. Introduction

The sudden and unexpected outbreak of the coronavirus disease (COVID-19) pandemic is the most serious global public health emergency and has spread to more than 200 countries. COVID-19 is a contagious respiratory illness caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is a zoonotic virus. The symptoms of COVID-19 include fever or chills, cough, dyspnea, difficulty in breathing, fatigue, headache, nasal congestion or runny nose, muscle or body aches, sore throat, loss of smell or taste, nausea or vomiting, and diarrhea. Notably, some people become infected but do not develop any symptoms or feel unwell. Infectious droplets are released when someone with COVID-19 sneezes, coughs, or talks, which can settle in the mouth or nose of people who are nearby or possibly be inhaled into the lungs. COVID-19 is thought to spread mainly through close person-to-person contact (within approximately 6 feet). People who are infected but do not show symptoms may also spread the virus to others [1, 2].

In December 2019, COVID-19 was first reported in Wuhan, the capital of Hubei Province in central China. Early in the outbreak, many patients were reported to have a link with a large seafood and live animal market [3, 4]. On January 23, 2020, the Wuhan shutdown was implemented to limit the movement of people in and out of the city. Meanwhile, suspected and confirmed cases had been isolated, public transport had been suspended, schools and entertainment venues had been closed, public gatherings had been banned, and health checks and information had been widely disseminated [5, 6]. Through the lockdown, the spread of the virus was controlled. Residents returned to their normal lives 7 months after the city’s strict lockdown was lifted.

In June and July 2020, COVID-19 affected Beijing and Xinjiang Uygur Autonomous Region (for Xinjiang), respectively, giving rise to the second wave of infections; however, this wave was localized within the regions and did not cause widespread transmission. The two regions curbed the rapid spread of COVID-19 cases by adopting strict prevention and control measures. The details are as follows: on June 11, 2020, a new locally transmitted case was reported in Beijing. The Xinfadi wholesale market, which is the largest food wholesale market in Beijing, was identified as a possible source of infection within 22 hours. Most cases were linked to the market. The market was shut down, high-risk individuals were identified, and nearby communities were subjected to lockdown. The government ensured an increase in COVID-19 testing, including testing of suspected persons, and followed up on their status. With the implementation of targeted measures, Beijing controlled the surge in cases in less than a month. The number of daily infections returned to zero by July 7. The market was reopened to the public on August 16. On July 15, 2020, a 24-year-old woman from Urumqi (the capital of Xinjiang) tested positive for SARS-CoV-2, making her the first locally transmitted COVID-19 case in Xinjiang since February. Based on epidemiological investigations and laboratory testing, this round of outbreak was associated with a gathering, and most confirmed cases were reported in Urumqi. No newly confirmed or asymptomatic COVID-19 cases have been reported in Xinjiang since August 16. Before the mass vaccination campaign of COVID-19, there was successful management of the outbreak in Wuhan, which offered a valuable experience to those in Beijing and Xinjiang. Expanding the testing of the population and asking people from medium- and high-risk areas or who had close contact with confirmed or asymptomatic cases to quarantine themselves gated widespread community spread or spread across regions. These measures have effectively reduced the spread of the virus and prevented mass infections, thereby helping consolidate the efforts of COVID-19 control and prevention. These two outbreaks were quickly curbed. However, the time of pandemic control in Xinjiang was longer than that in Beijing.

With the massive global pandemic, it is natural that the outbreak emerges accidentally in other areas worldwide. There are geographical, economic, cultural, health care, and population density differences between people living in Beijing and Xinjiang. For instance, the population density in Beijing is about 1313 people per square kilometer, whereas it is 14.98 people per square kilometer in Xinjiang. Based on the initial spread of the infection, in Beijing, the source of the infection, that is, the market, was identified within 22 hours after the first case was confirmed. In Xinjiang, the epidemiological analysis showed that the virus source of this outbreak was from the same source of infection identified 13 days later. The other regions can learn from these two regions’ experiences in containing the spread of the virus. Quantifying the effectiveness of these control measures is of crucial importance for cities preparing for rapid responses in the event of another outbreak; if asymptomatic cases are identified quickly, the source of the infection can be controlled and the infection rate can be reduced. Therefore, massive screening of asymptomatic cases and quarantine of susceptible people in communities are necessary.

The dynamic model is a useful method for identifying the key factors affecting disease outbreaks and evaluating the intervention strategies. Based on the transmission mechanism, the model can dynamically predict future trends according to current information. The quantities used often include the basic or effective reproduction number and the final epidemic size [711]. A mathematical model is used to estimate the degree of the COVID-19 epidemic in cities worldwide [1216]. The SEIR model is the most widely adopted model to study the COVID-19 outbreak in China and other countries [1720]. Based on the SEIR model, some researchers have increased the number of states that have to follow the quarantine process [6, 21]. We analyzed the transmission of COVID-19 in two typical regions: Beijing and Xinjiang. Between the two regions, the outbreak occurred simultaneously and similar control measures were adopted, though there were significant differences, for instance, in the duration of confirming the source, the ability of mass testing, and proportional coefficient of quarantine.

In this study, we considered a dynamic model with intervention measures. The model is introduced in Section 2. In Section 3, the basic reproduction number, , a threshold quantity for the stability of equilibria, was calculated. Using the Lyapunov function, we proved that the disease-free equilibrium was globally asymptotically stable if , whereas an endemic equilibrium was globally asymptotically stable if . The final size was derived. In Section 4, we estimate the parameter values using Markov Chain Monte Carlo (MCMC) simulations. In Section 5, we explore the influence of mass testing and quarantine in Beijing and Xinjiang. The final size was evaluated for COVID-19 in Beijing and Xinjiang.

2. Model Formulation

The model was used to estimate the COVID-19 epidemic in Beijing and Xinjiang. We divided the total population () into five groups: susceptible (), unfound infected (), hospitalized (), removed (, removed group includes recovered and death populations), and quarantine (). Susceptible people were infected by asymptomatic infected people. Through a latent period, some infected people developed symptoms, went to the hospital, and were diagnosed, while others remained asymptomatic. Some asymptomatic people who were identified by nucleic acid testing were taken to the hospital. Therefore, hospitalized patients included symptomatic and asymptomatic individuals. Patients hospitalized were isolated and treated; thus, we considered that the susceptible people could not be infected by the hospitalized patients. Other asymptomatic people not diagnosed in the crowd recovered without treatment. Recovered patients included the treated and untreated patients. To minimize risks, susceptible people from medium- and high-risk areas or who had close contact with confirmed or asymptomatic cases were tested. If their results were negative on nucleic acid testing, they were prohibited from leaving the outbreak regions and quarantined. Therefore, the susceptible individuals were quarantined in proportion to the number of new cases hospitalized. Those quarantined remained till the outbreak was over in their regions. Hence, people quarantined did move toward susceptible people.

Based on the previous assumptions, the dynamic model was described by the following ordinary differential equations:where , , and , where is the infection coefficient. denotes the discovery rate of infectious people. and are the recovery rates. is the rate of natural births and deaths, where is the daily proportion of nucleic acid testing. With contact tracing, the uninfected individuals, who got in contact with new hospitalized cases, move to the quarantined compartment at a rate of . All parameters were nonnegative. A flow diagram of system (1) is shown in Figure 1. Additionally, the total population was included . Adding all five equations in system (1), we obtain , for all , which implies for all . Let the feasible region for system (1) be

And the initial condition is

Theorem 1. The region is positively invariant with respect to system (1).

Proof. To show the positivity of each individual, we first show that . We multiply both sides of the second equation of system (1) by and obtainWe obtainSince , we obtain for all time . Now to show that , we multiply both sides of the third equation of system (1) with and obtainWe may writeon integrating from to and getSince and for all times , we have . We multiply both sides of the fourth equation of system (1) by . By integrating from to , we obtainBecause , and for all times we have . From the fifth equation of system (1) and , we obtainWe multiply both sides of the above equation with and getOn integrating from to , we obtain thatSince , and for all times , we have . Similarly, to show that , we multiply both sides of the first equation of system (1) by and obtainwhich can be rewritten asOn integrating from to , we obtain thatand thus for all times .
. Therefore, the region is positively invariant. The solution enters in an infinite time or asymptotically approaches . The region attracts all solutions in [22].

3. Dynamical Analysis

For the convenience of calculation, adding the equations in system (1), it follows that is a first integral, and introducing the densities , , , , , and , system (1) becomes

3.1. Global Stability of the Disease-Free Equilibrium (DFE)

By solving the equations of system (16) at a steady state, we obtain a disease-free equilibrium (DFE) , where . We can evaluate the basic reproduction number for system (16) following the next-generation matrix formulated by van den Driessche and Watmough [23]. For the epidemic model, shows that an average of each infected individual infects more than one individual. If , the disease dies out. On the other hand, we may expect the disease to spread in the community if . The basic reproduction number for system (16) is

Theorem 2. The DFE of system (32) is locally asymptotically stable if and unstable if .

Proof. By linearizing system (16) around , the local stability of the DFE solution can be examined. Let and denote the vector field of the system in (16). Jacobian matrix associated with this :The characteristic equation for the above matrix isWe can see that all the eigenvalues have negative real parts if and only if which shows that the DFE is locally asymptotically stable. The DFE is unstable if .

Theorem 3. The DFE of system (16) is globally asymptotically stable in whenever .

Proof. Consider the following Lyapunov function:Then, if ,where if and only if . The largest compact invariant set in is a singleton . The global stability of follows from LaSalle’s invariance principle [2426]. This establishes the theorem.

3.2. Global Stability of the Endemic Equilibrium

By solving the equations of system (16) in the steady state, there is an endemic equilibrium if , where

Now, we prove that is globally asymptotically stable in the interior of . We establish the following result.

Theorem 4. The unique endemic equilibrium is locally asymptotically stable if .

Proof. By linearizing system (16) about , we obtain the Jacobian matrix:The characteristic equation for the above matrix iswhereHence, all the eigenvalues are negative real parts if , which shows that is locally asymptotically stable.

Theorem 5. The endemic equilibrium is globally asymptotically stable whenever .

Proof. We define the following Lyapunov function:where . By calculating the time derivative of along the solution of system (16), we obtainfor all . Therefore, is globally asymptotically stable.
Using the previous conclusion,We deriveSimilarly, we obtainAccording to LaSalle’s invariance principle [27, 28], the endemic equilibrium is unique and globally asymptotically stable in the interior of if . The theorem is proven.

3.3. Final Size

The final size of the epidemic is the number of individuals ultimately becoming infected, which is the total number of cases during the entire outbreak [29, 30]. In analyzing system (16), we adopt the conventions that, for an arbitrary continuous function with nonnegative components:

Due to the small-time span, we considered natural birth and death to be zero.

Letting in system (16), it becomes

The disease basic reproductive number is

We observe thatfor all , while we know that and ; hence, and exist, and we have , and . Adding , and equations to system (32), we obtain

Integration of the equation with respect to from to gives

Adding and equations to system (32),

From , we obtain that

An epidemic ends when no infection is found; that is, . The integration of the above equation yields

We now derive the equation for the final epidemic size, defined by . This yields an implicit relation:

4. Estimation of Epidemiological Parameters

We obtained data on COVID-19 cases in Beijing and Xinjiang from the National Health Commission of The People’s Republic of China. The data included the cumulative and new infections of confirmed, suspected, death, cured, and asymptomatic cases; the cumulative number of close contacts; and hospitalization time [31].

In China, the per-day rate of natural birth and death is . Suppose that the recovery rate of asymptomatic cases is . Beijing is a metropolitan area with a population of 21.536 million (). From June 11 to July, 9,335 confirmed local cases and 49 asymptomatic cases had been reported, and no deaths were reported of the cases in Beijing. According to the data in Beijing, the mean time from infection to being tested positive for the virus was 10.8 days , and the treatment time was 10 days . The testing capacity was approximately 0.3 million people per day . The initial values on June 11 in Beijing were and .

The population in Xinjiang is 24.87 million (). Since the outbreak of the virus on July 15, 828 confirmed COVID-19 cases and 391 asymptomatic cases have been reported in Xinjiang, mainly in Urumqi. No newly confirmed COVID-19 case had been reported in Xinjiang since August 16. The treatment time was 11 days . From July 15 to 19, about 0.24 million people who had close contact with infected people and asymptomatic carriers had got the nucleic acid test, where the testing capacity was 0.06 million people per day (). The mean time from infection to being tested positive for the virus was 14 days (). On July 18, a citywide, free nucleic acid test was launched in Xinjiang. The testing capacity was up to 0.5 million people per day () since July 20, the meantime from getting infected to being tested positive for the virus was shortened (). The initial values on June 11 in Xinjiang were , , and .

Based on the mathematical model, the cumulative numbers of confirmed cases and identified asymptomatic cases were determined in Beijing and Xinjiang. We used the MCMC method to fit the model and adopted an adaptive Metropolis–Hastings algorithm to carry out the MCMC procedure. The algorithm was run for 30,000 iterations with a burn-in of the first 10,000 iterations, and the Geweke convergence diagnostic method was employed to assess the convergence of chains [32, 33]. We could acquire the parameter values for and , and the initial values for and in the two regions. The mean values and 95% confidence intervals (95% CI) are presented in Table 1.

The infection coefficient in Beijing () is lower than that in Xinjiang (). in Beijing is greater than that in Xinjiang because the population density in Beijing (approximately 1313 people per square kilometer) is higher than that in Xinjiang (approximately 14.98 people square kilometer). in Beijing is lower than that in Xinjiang because the source of infection in Beijing was identified quicker than that in Xinjiang. The time evolution of the infected cases and comparison with the cumulative number of confirmed cases and identified asymptomatic cases are shown in Figure 2. Evidently, the theoretical prediction is nearly in agreement with the actual data, which also validates the accuracy of the model. The final sizes of the pandemic in Beijing and Xinjiang were , and , , respectively. The final scale in Beijing is smaller.

5. Spread of COVID-19 and the Effectiveness of Interventions in Beijing and Xinjiang

The basic reproduction numbers in Beijing and Xinjiang were 0.3783 and 1.0882, respectively. The basic reproductive number in Beijing was less than one, and it was verified that the disease dies out. The basic reproductive number in Xinjiang was greater than one; however, with the increasing expansion of the nucleic acid testing, the outbreak was curbed. To draft policies for control measures, we analyzed the cumulative number of confirmed and identified asymptomatic COVID-19 cases for different values of (Figure 3). The analysis showed that, to control the spread of COVID-19 (the smaller the final size, the shorter the duration), a larger coefficient of quarantine is necessary. In addition, as the population density in Xinjiang is less than that in Beijing, it is effective to increase the coefficient of quarantine in Xinjiang. We provided a sensitivity analysis of the basic production number with the parameters and in Beijing and Xinjiang (Figure 4). Apparently, improving the testing capacity can decrease the basic reproduction number . When the testing capacity reaches a certain level, it will make the basic reproduction number less than one (Figure 4). In Xinjiang, the basic reproduction number was equal to if the other parameters are fixed. This implies that the infection could not persist over time in Xinjiang if the number of people tested was greater than 415,600 per day. Without mass vaccination campaigns, mass testing may be an effective prevention measure. Early detection of new cases is a useful control strategy for COVID-19.

6. Conclusion and Discussion

In this study, we constructed a mathematical model for COVID-19 transmission in Beijing and Xinjiang in 2020. There are differences in geographical, economic, cultural, health care, and lifestyles of people living in the two regions. The population in Xinjiang is similar to that in Beijing; however, the population density in Xinjiang is 100 times less than that in Beijing. The outbreak occurred in Beijing and then in Xinjiang. Similar measures (mass testing and quarantine) were adopted in the two regions, though the time to confirm the source, the ability of mass testing, and the proportional coefficient of quarantine were different. To explore effective control and prevention measures for COVID-19, we proposed a deterministic model with mass testing and quarantine to describe the spread of COVID-19 in China. First, the basic reproduction number of the model was provided. Then, it was found that the model had two nonnegative equilibria: the DFE and endemic equilibrium. Through the analysis of the model, it was found that the global behavior of the system was completely determined by the size of the basic reproduction number , that is, the DFE was globally asymptotically stable if , while an endemic equilibrium existed uniquely and was globally asymptotically stable if . Based on mathematical analysis, data fitting, and sensitivity analysis, we revealed the effects of mass testing and quarantine. Without the mass vaccination campaign of COVID-19, it was found that the higher the mass testing is, the earlier the source of infection was determined for effectively controlling the spread. In addition, a larger coefficient of quarantine was necessary to control the spread of COVID-19. It was more effective in increasing the coefficient of quarantine if the population density of the region was less. If the government implements mass testing and quarantine, it must devote a large amount of continuous manpower, material, and financial resources to eliminate the epidemic. With the development of the vaccine for COVID-19, more people are vaccinated in a community, which makes it possible for eliminating the spread of COVID-19.

In the future, we plan to study the mass vaccination campaign for COVID-19 in the whole of China.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request (email: [email protected]).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (61873154), Shanxi Key Laboratory (201705D111006), and Shanxi Scientific and Technology Innovation Team (201705D15111172) Grant.