Summaries, Analysis and Simulations of Recent COVID-19 Epidemics in Mainland China ================================================================================== * Lequan Min ## Abstract **Background** Globally COVID-19 epidemics have caused tremendous disasters. China prevented effectively the spread of COVID-19 epidemics before 2022. Recently Omicron and Delta variants cause a surge in reported COVID-19 infections. **Methods** Using differential equations and real word data, this study modelings and simulates COVID-19 epidemic in mainland China, estimates transmission rates, recovery rates, and blocking rates to symptomatic and asymptomatic infections. The transmission rates and recovery rates of the foreign input COVID-19 infected individuals in mainland China have also been studied. **Results** The simulation results were in good agreement with the real word data. The recovery rates of the foreign input symptomatic and asymptomatic infected individuals are much higher than those of the mainland COVID-19 infected individuals. The blocking rates to symptomatic and asymptomatic mainland infections are lower than those of the previous epidemics in mainland China. The blocking rate implemented between March 24-31, 2022 may not prevent the rapid spreads of COVID-19 epidemics in mainland China. For the foreign input COVID-19 epidemics, the numbers of the current symptomatic individuals and the asymptomatic individuals charged in medical observations have decreased significantly after March 17 2022. **Conclusions** Need to implement more strict prevention and control strategies to prevent the spread of the COVID-19 epidemics in mainland China. Keeping the present prevention and therapy measures to foreign input COVID-19 infections can rapidly reduce the number of foreign input infected individuals to a very low level. Keywords * China COVID-19 * infection transmission rates * infection blocking rates * recovery rates * modeling * simulations ## 1 Introduction The 2019 coronavirus disease (COVID-19) has placed tremendous pressures onto the prevention, control, and healthcare systems worldwide. Many countries have experienced multiple outbreaks of the COVID-19 epidemic due to incomplete preventive measures and Omicron and Delta variants. As of 30 March 2022, there are more than 483.5 million confirmed cases of COVID-19 with more than 6.1 million deaths globally [1]. Modelling the dynamics of spread of disease can help people to understand the mechanism of epidemic diseases, formulate and evaluate prevention and control strategies, and predict tools for the spread or disappearance of an epidemic [2]. Since the outbreak of COVID-19 in Wuhan China, a large numbers of articles on modelings and predictions of COVID-19 epidemics have been published (for examples see [3–12]). China prevent effectively the spread of COVID-19 epidemics before Omicron and Delta variants appeared. Recently, reported numbers of symptomatic and asymptomatic COVID-19 infected individuals are increased rapidly. This paper summarizes, analyzes and simulates the recent COVID-19 epidemic in mainland China, estimates infection transmission rates, infection blocking rates, and preventive measures through modelings and numerical simulations. The rest of this paper is organized as follows. Section 2 introduces materials and methods. Subsection 3.1 modelings and simulates the dynamics of COVID-19 Epidemics in Mainland. subsection 3.2 modelings and simulates the dynamics of foreign input COVID-19 Epidemics in China. Subsection 3.3 discuses and compares the simulation results. Virtual simulations for China Mainland COVID-19 epidemic and foreign input COVID-19 epidemic are implemented in the Subsections 3.4 and 3.5, respectively. Concluding remarks are given in Section 4. ## 2 Materials and Methods The dataset of the China COVID-19 epidemics from December 31, 2021 to April 25, 2022 was collected and edited from the National Health Commission of the People’s Republic of China official website [13]. Using differential equation models stimulates the outcomes of the numbers of the current symptomatic individuals, the current asymptomatic individuals, the cumulative recovered symptomatic individuals and cumulative asymptomatic individuals discharged from observations. Equation parameters were determined by so-called minimization error square criterion described in references [12, 14, 15]. Using virtual simulations estimates outcomes of the spreads of COVID-19 epidemics in mainland China. Simulations and figure drawings were implemented via Matlab programs. ## 3 Analysis and Simulations of COVID-19 Epidemics in China ### 3.1 COVID-19 Epidemics in Mainland Figure 1 shows the outcomes of the numbers of the current symptomatic individuals (CSI) and the current asymptomatic individuals (CAI). Figure 2 shows the outcomes of the numbers of the cumulative recovered symptomatic individuals (CCSI) and the cumulative asymptomatic individuals (CCAI) discharged from observations. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/02/2022.04.03.22273225/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/05/02/2022.04.03.22273225/F1) Figure 1. Blue circles: outcome of the number of the current symptomatic individuals (CSI), blue line: outcome of the corresponding simulations of equation (1). Red circles: outcome of the numbers of the current asymptomatic individuals (CAI) charged in medical observations, red line: outcome of the corresponding simulations of equation (1). The lines colored by cyan, magenta, green and yellow correspond to the virtual simulation results of equation (1). See Section Mainland Epidemic Virtual Simulations for details. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/02/2022.04.03.22273225/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/05/02/2022.04.03.22273225/F2) Figure 2: Blue circles: outcome of the number of the current cumulative recovered symptomatic individuals (CCSI), blue line: outcome of the corresponding simulations of equation (1). Red circles: outcome of the numbers of the current cumulative asymptomatic individuals (CCAI) discharged in medical observations, red line: outcome of the corresponding simulations of equation (1). The lines colored by cyan, magenta, green and yellow correspond to the virtual simulation results of equation (1). See Section Mainland Epidemic Virtual Simulations for details. The recent COVID-19 epidemics in mainland China are still continuing. Although there existed a turning point of the current symptomatic infections on day 17 (January 17, 2021), after day 43 (February 12, 2022), numbers of the the current symptomatic and asymptomatic infections are increasing until reached to the second turning point on day 110 (April 20). For the current asymptomatic infections charged i medical observations, the number were increasing until reached recently to the first turning point on day 106 (April 16). In order to estimate numerically the transmission rates and blocking rates to symptomatic and asymptomatic infections, we need to set up mathematic models (similar to [12, 14, 15]) to simulate the dynamics of spread of infection disease. Because we modeling uniformly the epidemic situations appeared in different provinces and regions, the transmission rates *β**ij* of symptomatic infections and asymptomatic infections are changing over different transmission intervals. For the mainland epidemics over the *lth* transmission interval, the symptomatic infected individuals (*I*) and the asymptomatic infected individuals (*I**a*) infect the susceptible population (*S*) with the transmission rates of *β*11(*l*) and *β*21(*l*), respectively, making *S* become symptomatic infected individuals, and with the transmission rates of *β*12(*l*) and *β*22(*l*), respectively, making *S* become asymptomatic individuals. Then, a symptomatic individual is cured at a rate *κ*(*l*), an asymptomatic individual returns to normal at a rate *κ**a*(*l*). Here all parameters are positive numbers. Assume that the dynamics of an epidemic can be described by *m*-time intervals, which correspond different transmission rates, prevention and control measures, and medical effects. At *l*th time interval, the model has the form (similar to [12, 14, 15]): ![Formula][1] ![Formula][2] ![Formula][3] ![Formula][4] where Θ1(*i*) = (1 − *θ*1(*l*)) and Θ2(*l*) = (1 − *θ*2(*l*)) (*l* = 1, …, *m*) represent the blocking rates to symptomatic and asymptomatic infections, respectively. It can be assumed that the input transmissions can be divided into 16 time intervals (see solid points in Figs. 1 and 2). We need to determine the parameters of equation (1) for *l* = 1, 2, …, 16. Over the *lth* time interval [*t**l*−1, *t**l*], Denote ![Formula][5] ![Formula][6] Denote *I**c*(*t**l*) to be the number of the reported current symptomatic infected individuals, and *I**ca*(*t**l*) be the number of the reported current asymptomatic individuals charged in medical observations. Denote ![Graphic][7] to be the number of the reported current cumulative recovered symptomatic infected individuals, and *I**cra*(*t**l*) be the number of the reported current cumulative asymptomatic individuals discharged in medical observations. Using the minimization error square criterion: ![Formula][8] determines the *β**ij*(*l*)′*s, κ*(*l*)′*s, κ**a*(*l*)′*s*), *θ*1(*l*)′*s*, and *θ*2(*l*)′*s*. The calculated parameters are shown in Table 1. The corresponding simulation results of equation (1) are shown in Fig. 1 and Fig. 2. Observe that the simulation results of equation (1) were in good agreement with the data of the COVID-19 epidemics in mainland China (see the solid blue lines, the red lines in Fig. 1 and Fig. 2). View this table: [Table 1.](http://medrxiv.org/content/early/2022/05/02/2022.04.03.22273225/T1) Table 1. Equation parameters of the COVID epidemics in mainland China during 2021.12.31-2022.4.25. ### 3.2 Foreign input COVID-19 Epidemics in China Figure 3 shows outcomes of the numbers of the current symptomatic individuals (CSI) and the current asymptomatic individuals (CAI). Figure 4 shows outcomes of the numbers of the cumulative recovered symptomatic individuals (CCSI) and the cumulative asymptomatic individuals (CCAI) discharged from observations. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/02/2022.04.03.22273225/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/05/02/2022.04.03.22273225/F3) Figure 3: Blue circles: outcome of the number of the current symptomatic individuals (CSI), blue line: outcome of the corresponding simulations of equation (2). Red circles: outcome of the numbers of the current asymptomatic individuals (CAI) charged in medical observations, red line: outcome of the corresponding simulations of equation (2). The lines colored by cyan and magenta correspond to the virtual simulation results of equation (2). See Section Foreign Input Epidemic Virtual Simulations for details. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/05/02/2022.04.03.22273225/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/05/02/2022.04.03.22273225/F4) Figure 4: Blue circles: outcome of the number of the current cumulative recovered symptomatic individuals (CCSI), blue line: outcome of the corresponding simulations of equation (2). Red circles: outcome of the numbers of the current cumulative asymptomatic individuals (CCAI) discharged in medical observations, red line: outcome of the corresponding simulations of equation (2). The lines colored by cyan and magenta correspond to the virtual simulation results of equation (2). See Section Foreign Input Epidemic Virtual Simulations for details. The epidemics have reached the first turning point of the current symptomatic infected individuals on day 16. After day 40, the number of the current symptomatic individuals began to increase until on day 74 (March 15) to reach the second infection turning point. The number of the current asymptomatic individuals charged in medical observations reached the first infection turning point on day 35. After day 52, the number began to increase until on day 76 (March 17) to reach the second infection turning point. For the foreign input COVID-19 infected individuals, they were discovered immediately and no further transmissions generated each other after entering China. Therefore the model has simply the form ![Formula][9] ![Formula][10] ![Formula][11] ![Formula][12] where *β*11(*l*) and *β*22(*l*) represent input transmission rates of the symptomatic individuals and asymptomatic individuals over the *lth* time interval, respectively. It can be assumed that the input transmissions can be divided into 15 time intervals (see solid points in Fig. 3 and Fig. 4). We need to determine the parameters for equation (2) for *l* = 1, 2, …, 15 over the *lth* time interval [*t**l*−1, *t**l*]. Denote ![Formula][13] Denote *I**c*(*t**l*) to be the number of the reported current symptomatic infected individuals, and *I**ca*(*t**l*) be the number of the reported current asymptomatic individuals charged in medical observations. Denote ![Graphic][14] to be the number of the reported current cumulative recovered symptomatic infected individuals, and *I**cra*(*t**l*) be the number of the reported current cumulative asymptomatic individuals discharged from medical observations. Using the minimization error square criterion: ![Formula][15] determines the *β*11(*l*)′*s, β*22(*l*)′*s, κ*(*l*)′*s* and *κ**a*(*l*)′*s*. The calculated parameters are shown in Table 2. The corresponding simulation results of equation (2) are shown in Fig. 3 and Fig. 4. Observe that the simulation results of equation (2) were in good agreement with the data of the foreign input COVID-19 epidemics (see the solid blue lines and the red lines in Fig. 3 and Fig. 4). On the end points of the 15 investigated time-interval [*t**l*−1, *t**l*]′*s*, the simulated numbers and the actual reported numbers were approximate the same–errors were less than one person, respectively. View this table: [Table 2.](http://medrxiv.org/content/early/2022/05/02/2022.04.03.22273225/T2) Table 2. Equation parameters of the foreign input COVID epidemics in mainland China during 2021.12.31-2022.4.25. ### 3.3 Results and Discussions #### 3.3.1 Mainland Epidemics Recent China COVID-19 epidemics with both Omicron and Delta variants make more difficult to prevent spread of the diseases. On the end points of the twelve investigated time-interval [*t**l*−1, *t**l*]′*s*, that is on day 0, day 4, day 11, day 20, day 30, day 43, day 48, day 55, day 64, day 70 day 77 day 82, day 90, day 100, day 106, day 110 and day 115, we obtain the following results (also see Figs. 1 and 2). * On day 11, there was only 1 difference between the numbers of the reported and the simulated current symptomatic individuals. On the other end point days, the numbers of the reported and the simulated current symptomatic individuals were approximate the same. * On day 11 and day 20, there were 11 and 8 differences between the numbers of the reported and the simulated current asymptomatic individuals charged in observations. On the other end point days, the numbers of the reported and the simulated current asymptomatic individuals charged in observations were approximate the same. * On day 11, day 20, day 30, and day 43, there were only one or two differences between the numbers of the reported and the simulated current cumulative recovered symptomatic individuals. On the other end point days, the numbers of the reported and the simulated current cumulative recovered symptomatic individuals were approximate the same. * On day 11 and day 20, there were nine and three differences between the numbers of the reported and the simulated current cumulative asymptomatic individuals discharged from observations. On the other end point days, the numbers of the reported and the simulated current cumulative asymptomatic individuals discharged from observations were approximate the same. * The transmission rates of the symptomatic infections caused by the symptomatic individuals were increasing from day 21 (January 20) to day 70 (March 10), and then seems to stop (see *β*11(*l*)′*s* in Table 1). * The transmission rates of the asymptomatic infections caused by the symptomatic individuals have obviously increased after day 56 (February 24, see *β*12(*l*)′*s* in Table 1). * The transmission rates of the symptomatic infections caused by the asymptomatic individuals have obviously increased after day 55(see *β*21(*l*)′*s* in Table 1). * The transmission rates of the symptomatic infections caused by the asymptomatic individuals have obviously increased after day 55 (see *β*21(*l*)′*s* in Table 1). * During the first 70 days (before March 30), the transmission rates of the asymptomatic infections caused by the asymptomatic individuals are very low (see *β*22(*l*)′*s* in Table 1). After day 70, the transmission rates increased significantly. * The blocking rates Θ1(*l*)′*s* and Θ2(*l*)′*s* to symptomatic and asymptomatic infections were not hight. Even on day 90 (March 30), the blocking rates only reach about 87.6% and 39.0% (see Table 1), respectively. However, for the first and second epidemics in Beijing and the five wave epidemics in Shanghai, the blocking rates reached to over 95% in one month [14–16]. * The recovery rates *κ*(*l*) and *κ**a*(*l*) of the symptomatic individuals and asymptomatic individuals waved. The recovery rates reached first maximums during days 21-30 and days 31-43, respectively. In the last two weeks (After day 100, April 9), The recovery rates increased significantly (see Table 1). #### 3.3.2 Foreign input epidemics It seems that the foreign input COVID-19 infected individuals have been obtained good managements and therapies. On the end points of the 15 investigated time-interval [*t**l*−1, *t**l*]′*s*, that is on day 0, day 10, day 16, day 20, day 30, day 40, day 49, day 53, day 64, day 74 day 76 day 82 day 90, day 99, day 106, and day 115, we obtain the following results (also see Figs. 3 and 4). 1. On the end points of the 15 investigated time-interval [*t**l*−1, *t**l*]′*s*, the numbers of the reported and the simulated current symptomatic individuals were approximate the same. 2. On the end points of the 15 investigated time-interval [*t**l*−1, *t**l*]′*s*, the numbers of the reported and the simulated current asymptomatic individuals charged in medical observation were approximate the same. 3. On the end points of the 15 investigated time-interval [*t**l*−1, *t**l*]′*s*, the numbers of the reported and the simulated current cumulative recovered symptomatic individuals were approximate the same. 4. On the end points of the 15 investigated time-interval [*t**l*−1, *t**l*]′*s*, the numbers of the reported and the simulated current cumulative asymptomatic individuals discharged from medical observations were approximate the same. * The input transmission rates of the symptomatic infection individuals and the asymptomatic infection individuals waved. The maximal input transmission rates reached during days 50-64, and days 54-64, 100-115 respectively (see Table 2). * The recovery rates *κ*(*l*)′*s* of the symptomatic individuals were increasing during days 11-40 and days 54-82. The recovery rates *κ**a*(*l*)′*s* of the asymptomatic individuals were increasing during days 17-49, and days 65-90 (see Table 2). * After days 74 (March 14) and 76 (March 16, 2022), the turning point of the numbers of the current symptomatic individuals and the current asymptomatic individuals charged in observations appeared, respectively. #### 3.3.3 Comparing * The average recovery rate of the foreign input COVID-19 symptomatic infected individuals was about 0.071572. The average recovery rates of the mainland COVID-19 infected asymptomatic individuals was about 0.051501. The average recovery rate of the foreign input individuals was much higher than the one of the mainland people (see *κ*(*l*)′*s* in Table 1 and Table 2). * The average recovery rate of the foreign input COVID-19 asymptomatic infected individuals was about 0.06356. The average recovery rate for mainland COVID-19 infected asymptomatic individuals was about 0.03131. The average recovery rate of the foreign input individuals was much higher than the one of the mainland people (see *κ**a*(*l*)′*s* in Table 1 and Table 2). * The average input transmission rate of the foreign input symptomatic infection individuals was about 0.057037, which was much lower than the average transmission rate 0.17474 of the symptomatic infection causing by the mainland symptomatic individuals (see *β*11(*l*)′*s* in Table 1 and Table 2). * The average input transmission rate of the foreign input asymptomatic infection individuals was about 0.064133, which was much higher than the average transmission rate 0.034141. of the asymptomatic infection causing by the mainland symptomatic individuals (see *β*12(*l*)′*s* in Table 1 and *β*22(*l*)′*s* in Table 2). * The last 24 days’ prevention measures implemented for mainland COVID-19 infected individuals reduced significantly the spread speed of the symptomatic infections (see Fig. 1 and Table 1). However it needs more effective prevention measures such that the numbers of the infected symptomatic and asymptomatic individuals decrease more rapidly. * During days 83 to 115, the ratio of the average blocking rates of the asymptomatic infection and the symptomatic infection of the mainland epidemic was about 0.4958. Therefor it needs to increase largely the blocking rates to the asymptomatic infection. ## 4 Virtual Simulations ### 4.1 Mainland Epidemic Virtual Simulations Assume that after day 115 (April 24, 2022), it still keeps the transmission rates ![Graphic][16] *s*, the blocking rates (Θ1(16), Θ2(16)), the recovery rates *κ*(16), and *κ**a*(16) until day 145 (May 25, 2022). The simulation results of equation (1) are shown in Fig. 1 and Fig. 2 by cyan lines and magenta lines, respectively. Calculated results show that on day 145 (May 25), the numbers of current symptomatic and the asymptomatic infected individuals reach about 18329 and 1165535, respectively. The numbers of cumulative recovered symptomatic individuals and cumulative asymptomatic individuals discharged in medical observations reach about 140452 and 140452, respectively. Furthermore assume that after day 115, it still keeps the transmission rates ![Graphic][17], and recovery rates *κ*(16) and *κ**a*(16) but increases the blocking rates to (Θ1, Θ2) ≡ (99%, 99%) until day 145. The simulation results of equation (1) are shown in Fig. 1 and Fig. 2 by green lines and yellow lines, respectively. Calculated results show that on day 145, the numbers of current symptomatic and the asymptomatic infected individuals reach about 2564 and 22708, respectively. The numbers of cumulative recovered symptomatic individuals and cumulative asymptomatic individuals discharged in medical observations reach about 106872 and 524005 respectively. ### 4.2 Foreign Input Epidemic Virtual Simulations Assume that after day 115 (April 24, 2022), it still keeps the transmission rates *β*15(15), *β*22(11), the recovery rates *κ*(11), and *κ**a*(11) until day 145 (25 May, 2022). The simulation results of equation (2) are shown in Fig. 3 and Fig. 4 by cyan lines and magenta lines, respectively. Calculated results show that on day 145, the numbers of the current symptomatic and the asymptomatic infected individuals reduce to 93 and 233, respectively. The numbers of cumulative recovered symptomatic individuals and cumulative asymptomatic individuals discharged from medical observations reach about 8006 and 8470, respectively. ## 5 Concluding Remarks The main contributions of this paper are summarized as follows: 1. It is the first time to summary the COVID-19 epidemic from December 31 2021 to April 24, 2022 in mainland China. It shows a clear picture to prevent and control the spread of the COVID-19 China epidemics [13]. 2. It uses two models to simulate the recent China epidemics. The simulation results on the end points of transmission intervals were in good agreement with the real word data [13], in particular the case for foreign input infections. 3. The simulation results can provide possible interpretations and estimations of the prevention and control measures, and the effectiveness of the treatments. 4. The recovery rates of the foreign input symptomatic and asymptomatic infected individuals were much higher than those of the mainland COVID-19 infected individuals. 5. Virtual simulations suggest that * The evolution of the number of the current symptomatic individuals may be located in the region between the cyan line and the green line shown in Fig. 1 if the recovery rates and treatment efficacy keep the levels on April 24, 2022, depending on practical blocking rates to infections. Similarly, the evolutions of the number of the current asymptomatic individuals charged in medical observations may be located in the region between the magenta line and the yellow line shown in Fig. 1. * The evolution of the number of the current cumulative recovered symptomatic individuals may be located in the region between the cyan line and the green line shown in Fig. 2. The evolutions of the number of the current cumulative asymptomatic individuals discharged in medical observations may be located in the region between the magenta line and the yellow line shown in Fig. 2. * Although the infection turning points have appeared, it needs to rise the recovery rates *κ*(*l*)′*s, κ**a*(*i*)′*s* and the blocking rates Θ1(*l*)′*s* and Θ2(*i*)′*s* of the mainland COVID-19 infected individuals to reduce the spreads of the COVID-19 epidemics. * Keeping the measures implemented during days 107-115 (April 16-25, 2002) can rapidly reduce the number of foreign input infected individuals to a very low level in one month. * Different combinations of the eight parameters of Equation (1) may generate similar simulation results. Therefore need further study to obtain better parameter combinations to interpret COVID-19 epidemics. A recommendation is that the administration should at least maintain the prevention and control measures implemented 7 days after reaching the turning point [15,16]. Figure 1 and Table 1 suggest that the importance of the recommendation. It is not wise strategy to withdraw all prevention and control measures before the number of the all infected people have been cleared. 100% blocking the speed at which COVID-19 infection spreads is key Strategies for early clearance or reduction of epidemic spread possible [15, 16]. More strict prevention and control strategies are necessary to prevent the spread of COVID-19 with Omicron and Delta variations. It is expected that this research can provide better understanding, interpretation and leading the spread and control measures of epidemics. ## Data Availability All data produced are available online at [http://www.nhc.gov.cn/](http://www.nhc.gov.cn/). [http://www.nhc.gov.cn/](http://www.nhc.gov.cn/). ## Funding The author has not declared a specific grant for this research from any funding agency in the public, commercial or not for profit sectors. ## Conflict of Interest The author declares no potential conflict of interest. ## Data availability statement Data are available on reasonable request. Please email the author. ## Ethical Statement Not applicable/No human participants included. ## Footnotes * (1) Additional twenty days' epidemic data (from April 1 to 25 April， 2022) were simulated and analyzed in the new version. (2) In subsections 3.3.1 and 3.3.2, items (1)-(4) are additional discussions, respectively. * Received April 3, 2022. * Revision received April 29, 2022. * Accepted May 2, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1].WHO, “Coronavirus 2019,” March 30, (2022), [http://www.who.int/emergencies/diseases/novel-coronavirus-2019](http://www.who.int/emergencies/diseases/novel-coronavirus-2019). 2. [2].Daley D, Gani J, Epidemic Modeling: An Introduction. New York: Cambridge University Press, (2005). 3. [3].Zhao S, Chen H, “Modeling the epidemic dynamics and control of COVID-19 outbreak in China,” Quant Bio, no. 11, pp. 1–9, (2020). 4. [4].Yang Z, Zeng Z, Wang K, Wong SS, Liang W, et al., “Modified SEIR and AI prediction of the epidemics trend of covid-19 in China under public health interventions,” J Thorac Dis, vol. 12, pp. 165–174, (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.21037/jtd.2020.02.64&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F02%2F2022.04.03.22273225.atom) 5. [5].Zhao S, Linc Q, Ran J, Musa SS, Yang G, et al., “Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak,” Int J Infect Dis, vol. 92, pp. 214–217, (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.01.050&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32007643&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F02%2F2022.04.03.22273225.atom) 6. [6].Roosa K, Lee Y, Luo R, Kirpich A, Rothenberg R, et al., “Real-time forecasts of the COVID-19 epidemic in China from February 5th to February 24th, 2020,” Infect Dis Model, vol. 5, pp. 256–263, (2020). 7. [7].Tang B, Bragaazzi NL, Li Q, Tang S, Xiao Y, et al., “An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov),” Infect Dis Model, vol. 5, pp. 248–255, (2020). 8. [8].Li X, Li Z, “Study on transmission model of the novel coronavirus (in Chinese),” (2020), [http://www.paper.edu.cn/releasepaper/content/202002-42](http://www.paper.edu.cn/releasepaper/content/202002-42) (accessedJuly 3 2021). 9. [9].Chintalapudia N, Battinenia G, Sagaroa GG, Amenta F, “COVID-19 outbreak reproduction number estimations and forecasting in Marche, Italy,” Int J Infect Dis, vol. 96, pp. 327–333, (2020). 10. [10].Chang SL, Harding N, Zachreson C, Cliff OM, Prokopenko M, “Modelling transmission and control of the COVID-19 pandemic in australia,” Nat Commun, vol. 11, pp. 1–13, (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-13889-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31911652&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F02%2F2022.04.03.22273225.atom) 11. [11].Niu R Chan YC Wong EW Van Wyk MA, Chen G, “A stochastic SEIHR model for COVID-19 data fluctuations,” Nonlinear Dyn., vol. 106, pp. 1311–1323, (2021). 12. [12].Min L, “Summaries, analysis and simulations of covid-19 epidemics in shanghai,” J Infect Dis Ther, vol. 10, no. S1, p. 005, (2022). 13. [13].National Health Commsion of the People’s Republic of China, “Epidemic Bulletin (in Chinese),” 2021.12.31-2022.3.18, [http://www.nhc.gov.cn/](http://www.nhc.gov.cn/). 14. [14]. L. Min, “Modelling, Simulations and Analysis of the First and Second COVID-19 Epidemic in Beijing,” November 15, 2021, medRxiv preprint doi: [https://doi.org/10.1101/2021.07.04.21259205](https://doi.org/10.1101/2021.07.04.21259205). 15. [15]. L. Min, “Modelling, Simulations and Analysis of the First COVID-19 Epidemic in Shanghai,” April 7, (2021)., medRxiv preprint doi: [https://doi.org/10.1101/2021.06.20.21259203](https://doi.org/10.1101/2021.06.20.21259203). 16. [16]. L. Min, “Summaries, Analysis and Simulations of COVID-19 Epidemic in Shanghai,” January 12, 2022., medRxiv preprint doi: [https://doi.org/10.1101/2022.01.11.22269050](https://doi.org/10.1101/2022.01.11.22269050). [1]: /embed/graphic-3.gif [2]: /embed/graphic-4.gif [3]: /embed/graphic-5.gif [4]: /embed/graphic-6.gif [5]: /embed/graphic-7.gif [6]: /embed/graphic-8.gif [7]: /embed/inline-graphic-1.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-13.gif [10]: /embed/graphic-14.gif [11]: /embed/graphic-15.gif [12]: /embed/graphic-16.gif [13]: /embed/graphic-17.gif [14]: /embed/inline-graphic-2.gif [15]: /embed/graphic-18.gif [16]: /embed/inline-graphic-3.gif [17]: /embed/inline-graphic-4.gif