key: cord-0705734-l6v5kijs
authors: Song, Haitao; Jia, Zhongwei; Jin, Zhen; Liu, Shengqiang
title: Estimation of COVID-19 outbreak size in Harbin, China
date: 2021-04-10
journal: Nonlinear Dyn
DOI: 10.1007/s11071-021-06406-2
sha: 665da5066a0ecb3fc9a3f43c4e521b4880a051fa
doc_id: 705734
cord_uid: l6v5kijs

Since the first level response to public health emergencies was launched on January 25, 2020, in Heilongjiang province, China, the outbreak of COVID-19 seems to be under control. However, an outbreak of COVID-19 caused by imported cases developed in Harbin during April 2020. A mathematical model is established to investigate the transmission of COVID-19 in Harbin. Based on the dynamical analysis and data fitting, the research investigates the outbreak of COVID-19 in Harbin and estimates the outbreak size of COVID-19 in Harbin. The outbreak size estimated of COVID-19 in Harbin reaches 174, where 54% of infected cases were identified while 46% of infected cases were not found out. We should maintain vigilance against unfound infected people. Our findings suggest that the effective reproduction number decreased drastically in contrast with the value of 3.6 on April 9; after that the effective interventions were implemented by the Heilongjiang province government. Finally, the effective reproduction number arrived at the value of 0.04 which is immensely below the threshold value 1, which means that the Heilongjiang province government got the outbreak of COVID-19 in Harbin under control.

Since December 2019, Coronavirus Disease 2019 spreads all provinces in China [1] . In March 2020, China government got the epidemic of COVID-19 under control by implementing the firstlevel response to public health emergencies (FLRPHE) and strict control measures [2] . However, the COVID-19 spreads rapidly all over the world, and 100,221,840 confirmed cases and 2,154,967 deaths were reported as of January 27, 2021.

An outbreak of COVID-19 in Harbin was caused by an imported case, the first related case was diagnosed and two asymptomatic cases were found on April 9 [3] [4] [5] . As of April 30, 68 confirmed cases and 23 asymptomatic cases were reported from the COVID-19 outbreak in Harbin. On March 19, an American student (H) carrying COVID-19 came back to Harbin from New York and then infected her neighbor (C) during the confinement period. Due to the infection within the family, G was infected with COVID-19 and then infected C by the dine together. Subsequently, C with cerebral apoplexy sought medical advice in the second hospital in Harbin and the first affiliated hospital of Harbin medical university, which caused that many patients, health care providers, doctors and nurses were infected with COVID-19 and then led to the local outbreak of COVID-19 in Harbin. The important events related to the outbreak of COVID-19 in Harbin are shown in Table 1 .

Mathematical modeling is used usually to explore the transmission of diseases such as COVID-19 and also predict the development trend of disease with the help of previous related information [6] [7] [8] [9] [10] [11] . According to the confirmed cases outside of mainland China by 18th January, Imai et al. [12] inferred the epidemic size of COVID-19 in Wuhan. Using public information and mathematical models, Wu et al. [13] estimated the clinical severity of COVID-19. Tang et al. [14] evaluated the epidemic size in China, and the effect of control measures was assessed in China. With the help of the mathematical model and public information, Song et al. [9] estimated the epidemic size of COVID-19 in China and predicted the potential second epidemic in China. Subsequently, Song et al. [11] established a mathematical model, computed the basic reproduction number, and estimated the epidemic size of COVID-19 in Wuhan as of January 23, 2020. In addition, the bilinear neural network method can be used to study the differential equations [15, 16] . The transmission of COVID-19 in China has been largely studied [9, [11] [12] [13] [14] , but there is no related study on mathematical modeling of COVID-19 transmission in Harbin by now.

To study the spread of COVID-19 in Harbin, the basic reproduction number and the effective reproduction number were computed. The basic reproduction number is defined as the expected number of secondary cases produced by a single infection in a completely susceptible population [17, 18] . While the effective reproduction number is the mean number of secondary cases an infected person can cause in a population where there is some immunity or some interventions in place [19] . It is useful to compute the basic reproduction number to estimated the transmission capacity of COVID-19 when the outbreak occurs. However, as a result of interventions such as FLRPHE and isolation measures, the effective reproduction number changes in time, and then, it is necessary to investigate the effective reproduction number to combat COVID-19 [20] . Most studies investigated the spread of COVID-19 using susceptible-infectious-recovered (SIR) or susceptibleexposed-infectious-recovered (SEIR) model [13, 14, [21] [22] [23] [24] [25] [26] [27] [28] . For example, Liu et al. [26] predicted the cumulative number of reported cases to a final size using the SIR model. Wu et al. [21] used the SEIR model to forecast the potential and international spread of COVID- 19 in Wuhan, where they assumed that the infected people in the incubation period were not infectious. Here, we established a susceptible-unfound infected-found infected-removed (SIFR) model, where the unfound infected people include infected people in the incubation period and unfound asymptomatic and symptomatic infected people.

To estimate the epidemic size of COVID-19 and assess the effectiveness of interventions in Harbin, a mathematical model explaining the transmission dynamics of COVID-19 is established. The basic reproduction number and the effective reproduction number are computed, and dynamics are analyzed by rigorous mathematical analysis. Using the mathematical model and public information, the epidemic size of COVID-19 in Harbin is estimated and the effect of interventions on the transmission of COVID-19 in Harbin is evaluated.

The rest of the paper is organized as follows. In Sect. 2, the mathematical model is established and dynamics are analyzed. Section 3 gives the parameter estimation. The estimation of the epidemic size of COVID-19 in Harbin is given and the effectiveness of interventions is assessed. Section 5 presents the discussion and conclusion.

Based on the found infected cases in Harbin in April 2020, a mathematical model is used to compute the basic reproduction number, estimate the epidemic size of COVID-19 in Harbin, and assess the effect of interventions on the transmission of COVID-19 in Harbin. The report of the WHO-China Joint Mission on COVID-19 [29] shows that the infected people in the incubation period have infectivity. In the model, several assumptions are given in the following. The natural birth rate and death rate are not considered since the period of the epidemic is short. The unfound infected people with COVID-19 have the same infectivity with infectious people with COVID-19 and found infected people are quarantined and could not infect healthy people. Removed people are not infected with COVID-19 again. The unfound infected people include infected people in the incubation period and unfound asymptomatic and symptomatic infected people. The population (N ) is divided into susceptible people (S), unfound infected people (I ), found infected people (F) and removed people (R), respectively. Here,

Susceptible people become infected people by the transmission rate β after contact with unfound infected people. Infected people are found out by the found rate p. Infected people are removed by the removed rates γ 1 and γ 2 , respectively. Here, the removed rate includes recovery rate and death rate due to disease. The flow diagram is shown in Fig. 1 .

Based on the epidemiological patterns of COVID-19 in Harbin and previous work [21, 28, 30] , the transmission dynamics of COVID-19 in Harbin is presented by the following differential equations

where the nonnegative initial values

Theorem 2.1 For system (1) with initial nonnegative initial values (2) , the solutions of system (1) are nonnegative and ultimately bounded.

Proof Through the Theorem 5.2.1 in [31] , the solutions S(t), I (t), F(t) and R(t) are nonnegative instantly.

From the first equation of model (1), we obtain

From the equations of model (1), we have

Then N (t) = N 0 , which means that I (t), F(t) and R(t) ultimately bounded. The proof is completed.

For system (1), the disease-free equilibrium E 0 = (S 0 , 0, 0, R 0 ) with S 0 > 0 and R 0 > 0. With the help of the next generation matrix theory [17, 18] , the basic reproduction number is calculated, where

Then the effective reproduction number R e (t) is

.

In what follows, the stability of the disease-free equilibrium E 0 is proved. (1),

Proof Since the stability of S(t) and R(t) is determined by the stability of I (t) and F(t), and N (t) = S(t) + I (t) + F(t) + R(t), then system (1) is simplified into the following equations

The Jacobian matrix of system (3) at E 0 is obtained, where

Thus, the characteristic equation at E 0 is

If R 0 < 1, then the roots β − γ 1 − p < 0 and −γ 2 < 0. Therefore, the disease-free equilibrium E 0 is locally asymptotically stable. Now we prove the global stability of the disease-free equilibrium E 0 . Define the Lyapunov function

Apparently, L(t) ≥ 0(∀t ≥ 0). Differentiating L(t) along the solutions of system (3) yields

When R 0 < 1, L (t) ≤ 0. Furthermore, the largest compact invariant set in {L (t) = 0} is the singleton {E 0 }. Using the LaSalle invariance principle [32] and the local stability of E 0 , the disease-free equilibrium point E 0 is globally asymptotically stable. The proof is completed.

In the following, it follows from [33] that the epidemic size and peak value of COVID-19 in Harbin could be obtained.

The sum of the first and second equations of system (1) is

which means that I ∞ = lim t→∞ I (t) = 0 and lim t→∞ (S(t) + I (t)) = S ∞ . Integration of the equation (4) from 0 to ∞ leads to

From the first equation of model (3), we have

It follows from equation (5) that the relation between the basic reproduction number and epidemic size of the COVID-19 is given. Furthermore, integration of the first equation of model (3) from 0 to t yields

and then

When the derivative of I is zero (that is, S = N 0 (γ 1 + p) β ), we obtain the maximum number of infectives which is

Simulations are carried out to verify the theoretical results. Setting the parameters β = 0.1403, p = 0.1128, γ 1 = 0.1, γ 2 = 0.1 and the initial value (S(0), I (0), F(0), R(0)) = (25, 11, 4, 0) . Figure 2 shows that the disease-free equilibrium E 0 is stable when R 0 < 1. Susceptible and removed people approach to S 0 and R 0 , respectively.

Setting the parameters β = 0.6403, p = 0.1128, γ 1 = 0.1, γ 2 = 0.1 and the initial value (S(0), I (0), F(0), R(0)) = (173, 1, 0, 0). Figure 3 gives the relation between S and I which describes the orbits of 

Data on found infected cases of COVID-19 from April 9 to April 30, 2020, in Harbin were obtained from the Health Commission of Heilongjiang Province [3] . The data set includes the cumulative number of found infected cases, newfound infected cases and cured cases. These data used are from publicly available data sources. The recovery rate of 0.1 -- [34] unfound infected people γ 2 The removed rate of 0.1 -- [34] found infected people 

As of April 30, 68 confirmed cases and 23 asymptomatic cases were reported from the COVID-19 outbreak in Harbin [3] . According to the study in [34] , the average time of treatment is 10 days (γ 1 = 0.1 and γ 2 = 0.1). As of April 9, one confirmed case and three asymptomatic cases were reported. Then the initial values on April 9 are F(0) = 4, R(0) = 0.

The initial values for two state variables and two unknown parameters were estimated using Bayesian methods. Multivariate Gaussian was chosen as the prior distribution of two unknown parameters. We chose the mean of the posterior distribution as the estimated value of the initial values and parameters, which were estimated by the Markov Chain Monte Carlo (MCMC) method. In view of the mathematical model and confirmed cases, using the MCMC method employing the adaptive Metropolis-Hasting algorithm with 20000 iterations and a 10000 iteration burn-in period [35] , parameter values β, p and initial values S(0), I (0) are estimated. In addition, the mean value, standard deviation (STD) and 95% confidence interval (95% CI) are given in Table 2 .

Concerning the uncertainty of estimated parameters and initial values, the MCMC method is used to evaluate the performance of our model (1) by the estimated parameters and initial values in Table 2 . Figure 4 shows the estimated cumulative infected cases and read data of COVID-19 in Harbin. Simulations are consistent with the reported cases, which validates the accuracy of our model.

Simulation results of the estimated cumulative infected cases are shown in Fig. 4 . Then we compute the basic reproduction number of 3.6 on April 9, 2020. Figure 4 shows that the cumulative number of infected people reached 174, the cumulative number of found infected people was 94 and the cumulative number of unfound infected people was 80. The effective reproduction number R e (t) for our model is shown in Fig. 5 . On April 9, the effective reproduction number was 3.6 which was the maximum of R e (t).

As time went on, R e (t) decreased quickly and was less than the threshold value 1 after April 15. Subsequently, R e (t) arrived at the minimum of 0.04 on April 30 which immensely below the threshold value 1. The outbreak of COVID-19 is under control when the effective reproduction number R e (t) is less than 1. This implies that the outbreak of COVID-19 in Harbin in April 2020 was under control when effective interventions were implemented.

Since the first case of COVID-19 was reported in Harbin on April 9, 2020, COVID-19 caused 68 confirmed cases and 23 asymptomatic cases in Harbin from March to April 2020. In the paper, using the public information and our mathematical model, we estimated the COVID-19 outbreak size in Harbin in April 2020.

Results show that the cumulative number of infected people reached 174, the cumulative number of found infected people was 94 and the cumulative number of unfound infected people was 80. To assess the effectiveness of interventions in Harbin, the effective reproduction number was estimated based on the public information and our mathematical model. This is the first study to estimate the transmission potential of the COVID-19 outbreak in Harbin in April 2020. The cumulative number of infected people finally reached 174, where 54% of infected people were found and 46% of infected people were not found out. Although all close contacts tracked were detected, some infected people were not found out. Indeed, the unfound infected people might be the infected people in the incubation period, unfound asymptomatic and symptomatic infected people. Therefore, it is dangerous for public health to ignore unfound infected people. We must maintain vigilance against unfound infected people.

Our findings indicate that the effective reproduction number on April 9 reaches 3.6 which is consistent with the estimated value in China [14, 21, 30, 36] . Besides, the estimated basic reproduction number was 6.8 in Hubei province [34] , and 3.6 in New York [37] . However, when the effective interventions were implemented by the Heilongjiang provincial government, the effective reproduction number R e (t) drastically dropped and finally reached 0.04 which is greatly below the threshold value 1, which suggest that the outbreak of COVID-19 in Harbin in April 2020 was under control and no subsequent outbreak in Harbin.

The mathematical modeling used in the study is analogous to the transmission dynamics model of COVID-19 in [12, 13, [22] [23] [24] [25] [26] [27] [28] [38] [39] [40] [41] [42] [43] [44] . The SIFR model helps us to estimate the cumulative number of infected cases of COVID-19 in Harbin in April and the effective reproduction number using found infected cases. Nonetheless, there are several limitations. First, we assumed that infected people with COVID-19 in the incubation period have the same infectivity as infec-tious people with COVID-19, which caused that the outbreak size of COVID-19 in Harbin was overestimated. Second, the detailed interventions were not incorporated into our model, which might lead to overestimating the outbreak size of COVID-19 in Harbin in April 2020. Third, very little is known about the effect of temperature and precipitation on the transmission of COVID-19. Our estimation could be untrustworthy if temperature and precipitation have a strong impact on the transmission of COVID-19.

Although the transmission of COVID-19 has been under control in China, the epidemic situation of COVID-19 all over the world is serious. Now, an increasing number of imported infected cases got into China and a growing number of asymptomatic infected people were found, which might increase the risk of a local outbreak of COVID-19 in China. Therefore, we should stay alert in case that unfound infected people might cause local outbreaks of COVID-19 in China such as the outbreak of COVID-19 in Harbin in April 2020 [9, 45] . The resurgence of COVID-19 in Beijing was likely caused by a polluted environment-to-human transmission from food via cold chain logistics [46] . The polluted environment-to-human transmission from food via cold chain logistics brings new challenges and its effect on the transmission of COVID-19 would be assessed, which we leave as our future work. Also, the differences of spread of COVID-19 in different hospitals and more modeling and verification for different cities in China will bring challenges from analysis as well as regression, and we will leave this as a future work.

A Novel Coronavirus from Patients with Pneumonia in China

Tibet activates highest-level public health alert

Chinese Center for Disease Control and Prevention

National Health Commission of the People's Republic of China

Global dynamics of two heterogeneous SIR models with nonlinear incidence and delays

Global dynamics of a Vector-Borne disease model with two delays and nonlinear transmission rate

Global dynamics of a multistage SIR model with distributed delays and nonlinear incidence rate

The Impact of Isolation on the Transmission of COVID-19 and Estimation of Potential Second Epidemic in China

Modeling the effect of temperature on dengue virus transmission with periodic delay differential equations

Using travellerderived cases in Henan Province to quantify the spread of COVID-19 in Wuhan

Estimating the potential total number of novel Coronavirus (2019-nCoV) cases in Wuhan City

Estimating clinical severity of COVID-19 from the transmission dynamics in Wuhan China

The effectiveness of quarantine and isolation determine the trend of the COVID-19 epidemics in the final phase of the current outbreak in China

Bilinear neural network method to obtain the exact analytical solutions of nonlinear partial differential equations and its application to p-gBKP equation

Rogue wave solutions and the bright and dark solitons of the (3+ 1)-dimensional Jimbo-Miwa equation

On the definition and the computation of the basic reproduction ratio R 0 in models for infectious diseases in heterogeneous populations

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

Epidemiology of transmissible diseases after elimination

Estimating the effective reproduction number of the 2019-nCoV in China

Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study

Modeling analysis of COVID-19 based on morbidity data in Anhui China

The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study

Current trends and future prediction of novel coronavirus disease (COVID-19) epidemic in China: a dynamical modeling analysis

A mathematical model for simulating the phase-based transmissibility of a novel coronavirus

Predicting the cumulative number of cases for the COVID-19 epidemic in China from early data

Assessment of monthly economic losses in Wuhan under the lockdown against COVID-19

An investigation of transmission control measures during the first 50 days of the COVID-19 epidemic in China

Using observational data to quantify bias of traveller-derived COVID-19 prevalence estimates in Wuhan China

Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems: An Introduction to the Theory of Competitive and Cooperative Systems

The stability of dynamical systems

Mathematical models for communicable diseases

Estimation of the timevarying reproduction number of COVID-19 outbreak in China

DRAM: efficient adaptive MCMC

Studies of the strategies for controlling the COVID-19 epidemic in China: Estimation of control efficacy and suggestions for policy makers

Epidemic model guided machine learning for COVID-19 forecasts in the United States

Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2)

Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions

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

A simple model for COVID-19

Modeling COVID-19 scenarios for the United States

Dynamics of COVID-19 under social distancing measures are driven by transmission network structure

Integrated vaccination and physical distancing interventions to prevent future COVID-19 waves in Chinese cities

Modelling COVID-19 transmission: from data to intervention

Cold-chain food contamination as the possible origin of COVID-19 resurgence in Beijing

The authors would like to thank the editor and anonymous referees for their valuable comments and suggestions. This study is supported by the National Natural Sci- 

Conflict of interest All authors declare no conflicts of interest in this paper.