Summaries, Analysis and Simulations of Recent COVID-19 Epidemics in Mainland China

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 higher than those of the mainland COVID-19 infected individuals. The blocking rates to symptomatic and asymptomatic mainland infections are lower than those of previous epidemics. 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 172022. Conclusions: Need to implement more strict prevention and control strategies to prevent the spread of COVID-19 mainland infected individuals. Need to keep the present prevention and therapy measures to foreign input COVID-19 infected infections until infected individuals to be cleared.


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][4][5][6][7][8][9][10][11]).
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. Section 3.1 modelings and simulates the dynamics of COVID-19 Epidemics in Mainland. Section 3.2 modelings and simulates the dynamics of foreign input COVID-19 Epidemics in China. Section 3.3 discuses and compares the simulation results, virtual simulations are implemented in the same Section. Concluding remarks are given in Section 4.

Materials and Methods
The dataset of the China COVID-19 epidemics from December 31, 2021 to March 31, 2022 was collected and edited from the National Health Commission of the People's Republic of China official website [12]. 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 [13,14]. Using virtual simulations estimates outcomes of the spreads of COVID-19 epidemics in mainland China. Simulations and figure drawings were implemented via Matlab programs.

Analysis and Simulations of COVID-19 Epidemics in China
COVID-19 Epidemics in Mainland Figure 1 shows outcomes of the numbers of the current symptomatic individuals (CSI) and the current asymptomatic individuals (CAI). Figure 2 shows outcomes of the numbers of the cumulative recovered symptomatic individuals (CCSI) and the cumulative asymptomatic individuals (CCAI) discharged from observations. 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 rapidly.
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 [13], [14]) to simulate the dynamics of spread of infection disease. Because we modeling uniformly the epidemic situations appeared in different provinces and regions, the 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 lth time interval, the model has the form (similar to [13,14]): CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.
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 twelve time intervals (see solid points in Figs. 1 and 2). We need to determine the parameters of equation Denote I c (t i ) to be the number of the reported current symptomatic infected individuals, and I cr (t i ) be the number of the reported current asymptomatic individuals charged in medical observations. Denote I cr(t i ) to be 3 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 7, 2022. the number of the reported current cumulative recovered symptomatic infected individuals, and I cra (t i ) be the number of the reported current cumulative asymptomatic individuals discharged in medical observations. Using the minimization error square criterion: 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 (3.1) are shown in Fig. 1 and Fig. 2. Observe that the simulation results of equation (3.1) were in good agreement with the data of the COVID-19 epidemics in mainland China (see the solid blue lines, the red lines and the legends in Fig. 1 and Fig. 2). 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. The epidemics have reached the first turning point of current symptomatic infections on day 16. After day 40, the number of the current symptomatic individuals began to increase until on day 74 to reach the second turning 4 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 7, 2022. ; l Days    For the foreign input COVID-19 infected individuals, they were discovered immediately and no further transmissions generated. Therefore the model has simply the form 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 twelve time intervals (see solid points in Fig. 3 and Fig. 4). We need to determine the parameters for equation 6 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 7, 2022. ; https://doi.org/10.1101/2022.04.03.22273225 doi: medRxiv preprint Denote I c (t l ) to be the number of the reported current symptomatic infected individuals, and I cr (t l ) be the number of the reported current asymptomatic individuals charged in medical observations. Denote I cr(t l ) 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: (I(t l ) − I c (t l )) 2 + (I r (t l ) − I cr (t l )) 2 + (I r (t l ) − I cr (t l )) 2 + (I ra (t l ) − I cra (t l )) 2 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 (3.2) are shown in Fig. 3 and Fig. 4. Observe that the simulation results of equation (3.2) were in good agreement with the data of the foreign input COVID-19 epidemics (see the solid blue lines and the red lines and the legends in Fig. 3 and Fig. 4).

Results and Discussions
Mainland Epidemics Recent China COVID-19 epidemics with both Omicron and Delta variants make more difficult to prevent spread of the diseases.
(1) The transmission rates of the symptomatic infections caused by the symptomatic individuals were increasing from day 21 to day 70, and then seems to stop (see β 11 (l) ′ s in Table 1).
(2) The transmission rates of the asymptomatic infections caused by the symptomatic individuals have obviously increased after day 56 (see β 12 (l) ′ s in Table 1).
(4) 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).
(5) The transmission rates of the asymptomatic infections caused by the asymptomatic individuals are very low (see β 22 (l) ′ s in Table 1).
(6) The blocking rates Θ 1 (l) ′ s and Θ 2 (l) ′ s to symptomatic and asymptomatic infections were not hight. Even on day 90, the blocking rates only reach to 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 [13][14][15]. (7) The recovery rates κ(l) and κ a (l) of the symptomatic individuals and asymptomatic individuals waved. The recovery rates reached maximums during days 21-30 and days 31-43, respectively (see Table 1).

Foreign input epidemics
It seems that the foreign input COVID-19 infected individuals have been obtained good managements and therapies.
(1) The input transmission rates of the symptomatic infection individuals and the asymptomatic infection individuals waved. The maximal input transmission rates reached during days 50-53 and days 54-64, respectively (see Table 2).
(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).
(3) The recovery rates κ(l) ′ s of the foreign input symptomatic individuals and the foreign input asymptomatic individuals were much higher than those of mainland symptomatic and asymptomatic individuals (see Tables  1 and 2).
• The average recovery rates κ(l) ′ s for foreign input COVID-19 symptomatic infected individuals was about 0.060. The average recovery rates κ(l) ′ s for mainland COVID-19 infected asymptomatic individuals was about 0.040.
• The average recovery rates κ a (l) ′ s for foreign input COVID-19 asymptomatic infected individuals was about 0.0497. The average recovery rates κ(l) ′ s for mainland COVID-19 infected asymptomatic individuals was about 0.24.
• The last eight days' prevention measures implemented to 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 turning point of the number of current symptomatic individuals appear.
• During days 83 to 90, the ratio of the blocking rates of the asymptomatic infection and the symptomatic infection was about 44.6%. Therefor need to increase largely the blocking rates of the asymptomatic infection.

Mainland epidemic virtual simulations
Assume that after day 90 (March 31, 2022), it still keeps the transmission rates β ′ ij s, the blocking rates Θ 1 (11), Θ 2 (11)), the recovery rates κ(11), and κ a (11)) until day 100 (10 April, 2022). The simulation results of equation (3.1) are shown in Fig. 1 and Fig. 2 by cyan lines and magenta lines, respectively. Calculated 8 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 7, 2022. ; https://doi.org/10.1101/2022.04.03.22273225 doi: medRxiv preprint results show that the numbers of current symptomatic and the asymptomatic infected individuals reach about 38020 and 106871, respectively. The numbers of cumulative recovered symptomatic individuals and cumulative asymptomatic individuals discharged in medical observations reach about 35650 and 23030, respectively.
Furthermore assume that after day 90 (March 31, 2022), it still keeps the transmission rates β ′ ij s, and recovery rates (κ(11), κ a (11)) but increases the blocking rates to (Θ 1 , Θ 2 ) = (99%, 99%) until day 100. The simulation results of equation (3.1) are shown in Fig. 1 and Fig. 2 by green lines and yellow lines, respectively. Calculated results show that the numbers of current symptomatic and the asymptomatic infected individuals reach about 20691 and 48854 , respectively. The numbers of cumulative recovered symptomatic individuals and cumulative asymptomatic individuals discharged in medical observations reach about 32284 and 17814 respectively.

Foreign input epidemic virtual simulations
Assume that after day 90 (March 31, 2022), it still keeps the transmission rates β 11 (11), β 22 (11), the recovery rates κ(11), and κ a (11)) until day 100 (10 April, 2022). The simulation results of equation (3.2) are shown in Fig.  3 and Fig. 4 by cyan lines and magenta lines, respectively. Calculated results show that the numbers of current symptomatic and the asymptomatic infected individuals reduce to 318 and 879, respectively. The numbers of cumulative recovered symptomatic individuals and cumulative asymptomatic individuals discharged in medical observations reach about 7147 and 5393, respectively.

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 March 31, 2022 in mainland China. It shows a clear picture to prevent and control the spread of the COVID-19 China epidemics [12].
(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 [12], 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 higher that 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. 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.
• Keeping the measures implemented during days 83-90 (March 24-31,, 2002) can be cleared shortly the current symptomatic and asymptomatic foreign input individuals.
(6) Different combinations of the eight parameters of Equation (3.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 [14,15]. 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 [14,15].
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.

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.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 7, 2022. ; https://doi.org/10.1101/2022.04.03.22273225 doi: medRxiv preprint