Abstract

Introduction

At the beginning of 2020, the COVID-19 pandemic was able to spread quickly in Wuhan and in the province of Hubei due to a lack of experience with this novel virus. For quite a bit of time, the current number of people infected was unknown. In fact, e.g., authorities in China found severe uncertainties regarding the dynamics and spread of the SARS-CoV-2 virus which causes COVID-19. In this paper, we give a more accurate estimate of the number of infected people in Wuhan and discuss the effects and the need of protective measures against this epidemic, e.g., such as the quarantines imposed by the Chinese government.

Several important studies about COVID-19 have already been published. Wu et al. [1] used the SEID metapopulation model applied to international travel data from Wuhuan and estimated that the basic reproductive number for COVID-19 was 2.68 (95% CrI 2.47-2.86); that means 75 815 individuals (95% CrI 37 304—130 330) were infected in Wuhan as of January 25th, 2020, see [1]. Li et al. [2] used observations of reported infections in China based on Chinese mobility data. By utilizing a network-based, dynamical metapopulation model using Bayesian inference, they obtained that 86% of all infections were simply not documented (95% CI: [82%-90%]) before the travel restrictions in Wuhan were arranged by January 23rd, 2020, see [2]. Lai et al. [3] investigated the temporal original rate of viral evolution and population dynamics of SARS-CoV-2 using 52 full genomes of viral strains sampled in different countries, and estimated the R value equal to 2.6 (range, 2.1-5.1) [3]. Qun Li et al. collected information on demographic characteristics, exposure history, and illness timelines of the first 425 laboratory-confirmed COVID-19 cases in China and estimated that the mean incubation period was 5.2 days (95% confidence interval [CI], 4.1 to 7.0) and the basic reproductive number was 2.2 (95% CI, 1.4 to 3.9) [4]. Sanche et al. combined case studies with model calculations and estimated that the incubation period is 4.2 days (95% CI 3.5-5.1 days) and the growth rate r is 0.29/day (95% CI 0.21-0.37/day), the estimated number of infected persons was ≈ 18 700 (95% CI 7 147-38 663) on January 23, 2020 [5]. Pan et al. [6] performed a cohort study using 32 583 laboratory-confirmed COVID-19 cases reported between December 8, 2019, and March 8, 2020 to determine the efficiency of public health interventions to control the COVID-19 outbreak in Wuhan [6].

In this paper, we use statistics of confirmed cases in China excluding Hubei and some other Chinese provinces as well as statistics from Beijing and Shanghai. Also, we use daily data generated by travel activity within China to analyze the daily number of infected persons in Wuhan City and Hubei Province from December 31, 2019 to January 23, 2020. In order to pursue, we develop a system dynamics (see, e.g., [7]) model that takes the infection and spread dynamics of the virus into account as well as the impact of protection measures against the epidemic such as quarantine. We employ multivariate Monte Carlo simulation (see, e.g., [8]) to parameterize the model using official statistics. In particular, we use the model to calculate the number of infections which already broke out but were not diagnosed due to various reasons.

Methods

A system dynamics model

System dynamics is an approach to understand the nonlinear behavior of complex systems over time using stocks, flows, feedback loops and non-linear functions, see, e.g., [7]. Based on the well-known SEIR model [12], we here develop a SEMIR model. This model takes the fact into account that in an epidemic situation there exist numerous infected people whose disease broke out but cannot be confirmed for various reasons. We put the emphasis on the city of Wuhan (W), the province Hubei excluding the city Wuhan (H) and China excluding the province Hubei (C) as three interconnected regions.

Our system dynamics model as shown in Fig 1 is a compact visualization of an integral equation system. Each node of the graph corresponds to an equation. For each of the so-called stocks, represented by a rectangle, an integral equation applies.

Open in a separate window

Fig 1

A system dynamics model for infection and spread dynamics of the SARS-CoV-2 virus.

Now we will explain the model step by step. An infection starts with an exposure. The exposure may take place locally or the infected persons can enter the region under consideration. The number of exposed cases is given by

$$\begin{array}{c}\hfill E_i\left(t\right)=E_i\left(0\right)+{\int }_0^t(f_i\left(\tau \right)+i_i\left(\tau \right)-o_i\left(\tau \right))d\tau \end{array}$$

(1)

where

$$\begin{array}{c}\hfill i=\{W,H,C\}.\end{array}$$

(2)

W stands for Wuhan, H for Hubei excluding Wuhan and C for China excluding Hubei, as mentioned above. The inflow of exposed cases is given by

$$\begin{array}{c}\hfill f_i\left(t\right)=f^{conv}[m_i^W\left(t\right)\frac{E_W\left(t\right)}{P_W}+m_i^H\left(t\right)\frac{E_H\left(t\right)}{P_H}-m_i\left(t\right)\frac{E_i\left(t\right)}{P_i}],\end{array}$$

(3)

where $m_i^W\left(t\right)$ and $m_i^H\left(t\right)$ are the daily immigrations into the region under consideration coming from Wuhan und Hubei, respectively. m_i(t) is the daily exmigration from the region under consideration, P_i the population of the region i. f^conv is the converting factor from the points to the number of persons.

The rate of local exposure is given by

$$\begin{array}{c}\hfill i_i\left(t\right)=i_i^E\left(t\right)+i_i^T\left(t\right)=p\left[q_i\left(t\right)\right(E_i\left(t\right)+M_i\left(t\right))+q_i^T{T}_i\left(t\right)]\end{array}$$

(4)

where $i_i^E\left(t\right)$ and $i_i^T\left(t\right)$ are the exposures caused by exposed cases, E_i(t), missed cases M_i(t), and the ones in treatment denoted by T_i(t). p is the probability of transmission per personal contact and q_i(t) and $q_i^T\left(t\right)$ denote the contact frequency of such cases. In order to simplify, we introduce the standard contact here, where the transmission probability p equals 1%. Note that all non-pharmaceutical interventions, e.g., quarantine, aim to lower the frequency q_i(t).

After a certain presymptomatic or incubation period t^E, an exposed case becomes infected and will be diagnosed:

$$\begin{array}{c}\hfill o_i\left(t\right)=\text{delay}3(f_i(t-\lambda \mu t^E)+i_i(t-\lambda t^E),(1-\lambda )t^E)\end{array}$$

(5)

where the function delay3 and the parameters λ are used to approximately depict the statisitical distribution of the incubation time of the exposed cases aorund the average value t^E [13]

Now we analyze the so-called missed cases. If a new epidemic breaks out, it’s likely that a certain amount of infected cases are not diagosed and registered due to a lack of experience, medical capacity, or other reasons. Our model takes such possible missed cases into account:

$$\begin{array}{c}\hfill M_i\left(t\right)=M_i\left(0\right)+{\int }_0^t(o_i\left(\tau \right)-l_i\left(\tau \right)-c_i\left(\tau \right))d\tau \end{array}$$

(6)

where

$$\begin{array}{c}\hfill c_i\left(t\right)=\{\begin{array}{cc}{o}_i\left(t\right)& \text{if}{m}^C=0\\ c_i^{data}\left(t\right)& \text{else}\end{array}\end{array}$$

(7)

In case of

$$\begin{array}{c}\hfill m^C=0,\end{array}$$

(8)

all infected cases are diagnosed and registered. Otherwise only some of the cases given by the data of new cases, $c_i^{data}\left(t\right)$, are found and some other cases end without being diagnosed and registered. That means:

$$\begin{array}{c}\hfill l_i\left(t\right)=\frac{M_i\left(t\right)}{t^T},\end{array}$$

(9)

where t^T is the average duration of the critical phase in which a patient may die.

The diagnosed cases are treated medically or at least observed. The following applies for the cases during treatment:

$$\begin{array}{c}\hfill T_i\left(t\right)=T_i\left(0\right)+{\int }_0^t(c_i\left(\tau \right)-d_i\left(\tau \right)-s_i\left(\tau \right))d\tau .\end{array}$$

(10)

The following applies to deceased and rescued cases:

$$\begin{array}{c}\hfill d_i\left(t\right)=r_i^D\frac{T_i\left(t\right)}{t^T},\end{array}$$

(11)

$$\begin{array}{c}\hfill s_i\left(t\right)=(1-r_i^D)\frac{T_i\left(t\right)}{t^T},\end{array}$$

(12)

where $r_i^D$ is the mortality of the disease under consideration.

Finally we discuss the features of the model when analyzing cases of death, recovered and cured cases. To be able using the official statistics for our calculation, we consider the cases of death and cured cases as well as the recovered cases but have not yet been released from medical treatment:

$$\begin{array}{c}\hfill D_i\left(t\right)=D_i\left(0\right)+{\int }_0^t{d}_i\left(\tau \right)d\tau ,\end{array}$$

(13)

$$\begin{array}{c}\hfill C_i\left(t\right)=C_i\left(0\right)+{\int }_0^t{r}_i\left(\tau \right)d\tau ,\end{array}$$

(14)

$$\begin{array}{c}\hfill R_i\left(t\right)=R_i\left(0\right)+{\int }_0^t(s_i\left(\tau \right)-r_i\left(\tau \right))d\tau ,\end{array}$$

(15)

where

$$\begin{array}{c}\hfill r_i\left(t\right)=s_i(t-t^R+t^T).\end{array}$$

(16)

t^R is the total time from the diagnosis of a case to its full recovery.

Parameterization

The model shown in Fig 1 and in the previous section contains a number of parameters. In the first step we neglect the exposure due to diagnosed cases, i.e. set $q_i^T\left(t\right)=0$, to focus on those parameters such that the model can reflect the actual infection and spread of COVID-19. These parameters include:

• E_W(0) is the number of exposed cases in Wuhan on Dec. 31th 2019

• f^conv is the converting factor from the points to number of persons

• q_W(0) and q_H(0) are the contact frequencies in Wuhan and Hubei before Jan. 23rd 2020

• t^E is the time between the exposure and the diagnosis

• λ is used to describe the the statisitical distribution of t^E

• μ is used to descibe the proportion of t^E which is spent in the target region

• q_C(t) is the time function of contact frequency in China excluding Hubei C which is given by four keypoints representing the contact frequency on Jan. 21st, Jan 25th, Feb. 10th and Feb. 27th, 2020.

By using multivariate Monte Carlo optimization, a large number (e.g., 10 000) of simulation runs are carried out with randomly generated values for the above mentioned parameters. The sensitivity results from these simulation runs are compared with the real data. (Fig 2). Hence, we obtain the appropriate values thereof and also confidence intervals for these parameters. However, since the daily departures from Wuhan and Hubei are considered by using random samples of the local population, the values for the confidence interval of the parameters that apply to Wuhan and Hubei increase accordingly. To keep the confidence interval within 15%, we focus on the period of time between January 17th through 26th where both the travel activity and the proportion of infected cases were at a relatively high level. Incidentally, a repetition of the simulation runs with other reasonable values of $q_i^T$ shows that the changes in the results compared to the confidence interval specified here are negligible.

Open in a separate window

Fig 2

Monte Carlo simulations for finding the parametrization.

Findings and interpretations

Using the method described in Section ‘Methods’, we find that the infection rate before the quarantine in Wuhan was 22.3% ± 3.3% (95% confidence interval, CI) per day. We determine the average time between exposure and confirmed diagnosis by 10.2 ± 1.5 days (95% CI). This yields to the reproduction number R equal to 2.27 ± 0.34 (95% CI) before the quarantine in Wuhan was arranged. Also, we estimate that the number of exposed cases based on the SEMIR model was 49 800 ± 7500 (95% CI) on January 23rd, 2020; this is the day the quarantine was arranged in Wuhan. More than 5 000 of these cases were active, although only 495 cases had been registered by that time. The number of infected cases which were not diagnosed for various reasons peaked at 28 200 ± 4 200 on February 04th 2020.

As shown by Fig 3, the model simulation can reproduce the actual course of the pandemic COVID-19 in China apart from Hubei. In addition, it provides information on the daily numbers of cases flowed in from Wuhan and Hubei and locally exposed cases. It is remarkable that non-pharmaceutical interventions actually work out in order to reduce the country-wide infection rate from 19% per day to ≤ 3% (or contact frequency from 19 per day to 3, a reduction of more than 80%) so that the spread of the pandemic was largely stopped within China.

Open in a separate window

Fig 3

Statistical data and simulated course of the pandemic COVID-19 in China excluding Hubei.

The situation in Wuhan until shortly after quarantine can also be determined with the present method. On February 3, there were almost 30 000 missed cases. However, we cannot determine how Wuhan then developed in this regard with the actual method. Fig 4 shows three different scenarios in which the contact frequency is 3.5, 4.0 and 4.5 per day, respectively.

Open in a separate window

Fig 4

Statistical data and simulated course of the pandemic COVID-19 in Wuhan.

Conclusions

In this paper, we came up with an estimate of the number of infected people in Wuhan to obtain a consistent, quantitative picture of the overall situation of a pandemic for connected areas. With our System Dynamics model, we took into account not only the infection dynamics, but also the spread through the travel activities and the anti-epidemic measures conducted. Our study revealed that the model was able to capture significant information to analyze the real number of infected people in Wuhan. Compared to foregoing studies discussed in the Section ‘Introduction’, our estimates are more accurate and the system dynamics model incorporates relevant factors such as quarantines imposed by the Chinese government.

Based on the model introduced in this article, we are closely monitoring the development of the COVID-19 pandemic in various countries and regions, observing the effect of different isolation methods and detection strategies in real time. Through data analysis, parameters such as infection rate, test coverage and immunity rate are obtained. Our aim is to contribute to contain the pandemic.

Funding Statement

Matthias Dehmer received funding from the Austrian Science Fund (P 30031). The funder did not play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

The data can be found at http://inform.bwi.unibw-muenchen.de/infoe/info.aspx?f=single&INFO=CBC9BC8F938999EE2CD6F0B193B. All data sources the authors used can be found through reference [6].

References