Abstract

Since the outbreak of the novel coronavirus disease (COVID-19) at the beginning of December 2019, there have been more than 28.69 million cumulative confirmed cases worldwide as of 12th September 2020, affecting over 200 countries and regions with more than 920,463 deaths. The COVID-19 pandemic has been sweeping worldwide with unexpected rapidity. In this paper, a hybrid modelling strategy based on tessellation structure- (TS-) configured SEIR model is adopted to estimate the scale of the pandemic spread. Building on the data pertaining to the global pandemic transmission over the last six months around the world, key impact factors in the transmission and control procedure have been analysed, including isolation rate, number of the infected cases before taking prevention measures, degree of contact scope, and medical level, so as to capture the fundamental factor influencing the pandemic. The quantitative evaluation allowed us to illustrate the magnitude of risks of pandemic and to recommend appropriate national health policy of prevention measures for effectively controlling both intra- and interregional pandemic spread. Our modelling results clearly indicate that the early-stage preventive measures are the most effective action to be taken to contain the pandemic spread of the highly contagious nature of the COVID-19.

1. Introduction to the COVID-19 Pandemic Development

Since the outbreak of novel coronavirus disease (COVID-19) at the beginning of December 2019, over 29.96 million people are infected with the disease, with 920,463 death cases around the world as of 12th September 2020. There are 87 countries whose confirmed cases exceed 10,000, affecting over 200 countries in the world at amazing width and rapidity against expectation. The risk level of COVID-19 pandemic was raised to “very high” on 28th February 2020 by the World Health Organization (WHO) since the global community is being rampaged by the disease [1]. The current distribution of COVID-19 pandemic cases around the world is shown in Figure 1 [2].

The COVID-19 virus has been proven as a pandemic at the global level [4, 5]. As it is known, the rapid declination in the number of the new infection cases in China was ascribed to the strict quarantine measures (i.e., including lockdown of Wuhan and suspension of all outdoor activities and public entertainment facilities) and timely construction of makeshift hospitals to confine infected patients, while the situation in other countries, i.e., South Korea, Japan, Iran, Italy, and France, was deteriorating rapidly at the same time. The preventive measures taken by the government in dealing with the virus prevention and control are beneficial for understanding the transmission feature of the COVID-19 for the purpose of effective containment [4, 6].

The pandemic control experiences learned from China would be invaluable to other countries [7, 8]. However, interventionist versus laissez-faire are the two different attitudes among the infected countries since the outbreak of COVID-19. The former is concerned with active inventions including radical confinement and quarantine measures, whereas the latter takes rather a sangfroid position, i.e., loose control or only providing suggestion in dealing with the virus transmission. For instance, Italy has taken similar strict prevention measures as in China, i.e., to close down the seriously infected regions and block the interpersonal spread at the early stage. Iranian government has also taken a series of intervention measures, e.g., closed public places, cancelled sports events, and suspended school nationwide. It should be noted that the pandemic prevention measures issued by the government would strictly be abided by the Iranian public. In contrast, only medical advice was provided by the government to the public in certain countries, without any further compulsory action taken to stop the community transmission at the early stage, allowing pandemic deteriorated rapidly.

Which measures are the appropriate ones to be taken by the governments have controversial debates, leading to various development directions of the highly contagious pandemic. In the meantime, the key factors influencing the COVID-19 transmission routes and infection rates should be investigated thoroughly to maximize the effectiveness of policy intervention with reasonable acceptable social cost to constraint the COVID-19 spread. Different modelling methods are applied to study the dynamics of the spread of infectious diseases such as susceptible-exposed-infectious-removed (SEIR) model [9], stochastic transmission model [10], and neural network model. However, the contact network effect, i.e., interpersonal contact times and scope, has seldom been taken into the modelling consideration. Therefore, to which certain level of the prevention measure should be adopted remains inconclusive without reaching a unanimous recognition although the control measures from certain countries have been proved to be effective.

In order to bridge this knowledge gap, an improved SEIR (susceptible-exposed-infectious-recovered) model embedded with tessellation structures (TSs) is developed to estimate the transmission tendency of the COVID-19 considering the key parameters in the procedure under different scenarios. It allows the random behaviour of interpersonal contact to be considered as probabilistic network parametric impact along with other factors, including isolation rate, contact number during gathering at regularity, and the size of the initial infected cases, from which the effective containment strategies, i.e., contact tracing or isolation measures, can be discussed and recommended accordingly. Based on the pandemic situation dealt in China and other early outbreak countries, the results presented in the paper clearly indicate that proactive prevention measures and effective medical resources providence are the most important actions to be adopted under high risk of COVID-19, especially in religious gathering areas and densely inhabited districts.

The remainder of the paper is organized as follows. Section 2 introduces the involved method with SEIR model combined with TS to form the virus transmission network. The key factors influencing the transmission speed and scope are analysed in Section 3. Discussion dealing with COVID-19 locally and globally to determine which measures should be adopted in different scenarios is provided in Section 4. Section 5 concludes the paper.

2. The COVID-19 Transmission Modelling Strategy

In this study, the standard SEIR (susceptible-exposed-infectious-recovered) model is coupled with a tessellation structure algorithm (TS) for describing each state in order to simulate the transmission of COVID-19 [11]. TSs are spatially and temporally discrete, abstract computational systems which can be used to describe nonlinear dynamics and composed of a set of homogenous and simple units. They evolve in parallel at discrete time steps considering the states of cells in their local neighbourhood, which can compute functions and emulate a universal computable problem, such as the coronavirus transmission pattern impacted by many stochastic factors [12]. In addition, the discrete random network characteristics of TS are particularly appropriate for investigating the effectiveness of the prevention and control strategy.

The developed TS-SEIR model has 7 states, such as S, E, I, R, Eq, H, and D, where the state setting of TS is designed based on four primary factors during COVID-19 control: (1) a 14-day quarantine period for an individual who comes from high-risk area or susceptible to be exposed by the infected; (2) individuals are highly suggested not going outside unless necessary; (3) individuals are restricted to participate in gatherings; and (4) patients should be admitted timely. In a typical mathematical model for an epidemic study, the number of the susceptible, latent, infectious, and removed individuals is represented as S, E, I, and R, respectively [13], while Eq, H, and D are the quarantined latent, quarantined infectious, and dead individuals, as illustrated in Figure 2.

There are total 14 rules to describe the state transmission: ① the susceptible individuals become exposed and enter the latent period without quarantine after contact with the infected; ② the susceptible individuals are quarantined and enter the latent period after contact with the infected; ③ the exposed individuals are quarantined; ④ the infected individuals are quarantined; ⑤ the exposed individuals without being in quarantine were infected after 14 days (average period of time); ⑥ the quarantined of the exposed individuals are confirmed infection and begin a quarantine of higher level; ⑦ the exposed individuals were infected and quarantined; ⑧ the infected individuals without quarantine recovered; ⑨ the quarantined of the exposed individuals are recovered without infection (small probability); ⑩ the quarantined of the infected individuals recovered; ⑪ the infected individuals without quarantine died; ⑫ the quarantined of the infected died; ⑬ the exposed individuals without quarantine infect the susceptible through contact; and ⑭ the infected individuals without quarantine infect the susceptible through contact.

It is assumed that the number of individuals in contact with other individuals in the unit time (one day) obeys a Poisson distribution, and then, the probability of m individuals in contact with j individuals per day is calculated aswhere is the node degree of individual m in the interpersonal contact network and the average number of each individual’s contact with other people in one day is denoted as . Thus, the probability of the ith individual being infected iswhere n is the total number of the individual nodes, and the infection probability can be calculated by the following equation:

It should be noted that there are small portions of people who have not been isolated in latent states but, with a larger node degree, would become more dangerous in an interpersonal contact network [14]. The developed TS-SEIR model can simulate the virus propagation, where each individual is treated as a cell evolving different states varied along with the time and space. The rules of the state variation are designed based on the characteristics of virus transmission. The implementation of programming the TS-SEIR model in our study follows the strategy as described below.

Formally, a TS is represented by a 4-tuple . represents the number of discrete regular grid of cells, i.e., the number of simulated population; represents the set of cell states; represents the neighbourhood; and represents the transition functions or transition rules. Here, it is set as , meaning one-dimensional TS. The unit for describing the transmission network having 7 states is written as

N is an adjacency matrix of the unit with 10000 × 10000 dimensions. Then, the transition rules (see Figure 2) of TS-SEIR model are described as Pseudocode 1.

  Parameters settings:
  q1 ⟵ 0.9; q2 ⟵ 0.99;
normalContactTimes ⟵ 6;
gatheringContactTimes ⟵ 80;
  initialInfectNum ⟵ 100;
  medicalLevel ⟵ 0.0413;
  
  //Initialize the adjacency matrix according to the value of normalContactTimes and gatheringContactTimes
  //Initialize the number of cell with states E and I according to the value of initialInfectNum
  for t ⟵ 1 to simulatedTimeStep
   do for i ⟵ 1 to populationNum
     if state(t, i) = S && anyNeighborOfCell(i) = E or I && infectRate ≥ rand
    then state(t + 1, i) ⟵ E;
    if state(t, i) = E && probTurnToInfected ≥ rand && q2 ≥ rand then state(t + 1, i) ⟵ H;
    if state(t, i) = E && probTurnToInfected ≥ rand %% q2 < r and then state(t + 1, i) ⟵ I;
   if state(t, i) = S or E && anyNeighborOfCell(i) = I or H && daysAfterContact(i, indexOfThisNeighbor) < 14 && q1 > rand then state(t + 1, i) ⟵ Eq;
    if state(t, i) = I && anyNeighborOfCell(i) = H then state(t + 1, i)⟵H;
    if state(t, i) = I && 1/(1 + exp((10 − timeLongInfected(i))/2)) < rand && q2 > rand
    then state(t + 1, i) ⟵ H;
    if state(t, i) = E && probTurnToInfected < rand && probESelfHealing ≥ rand
       then state(t + 1, i) ⟵ R;
   if state(t, i) = I && 1/(1 + exp((10 − timeLongInfected(i))/2))> = rand
    && 1 − medicalLevel < rand then state(t + 1, i) ⟵ R;
  if state(t, i) = I && 1/(1 + exp((10 − timeLongInfected(i))/2)) ≥ rand
    && 1 − medicalLevel ≥ rand then state(t+1, i) ⟵ D;
  if state(t, i) = Eq && probTurnToInfected ≥ rand then state(t + 1, i) ⟵ H;
  if state(t, i) = Eq && probTurnToInfected < rand && probEqSelfHealing ≥ rand
    then state(t + 1, i) ⟵ R;
    if state(t, i) = H && 1/(1 + exp((10 − timeLongInfected(i))/2)) ≥ rand.
    && 1 − medicalLevel < rand then state(t + 1, i) ⟵ R.
    if state(t, i) = H && 1/(1 + exp((10 − timeLongInfected(i))/2)) ≥ rand
    && 1 − medicalLevel ≥ rand then state(t + 1, i) ⟵ D
   //update daysAfterContact matrix and adjacency matrix

Several parameters for describing the transmission degree are involved in the TS-SEIR model, which are listed in Table 1. is the isolation rate, where q1 is the isolation rate in the latent period and q2 is the isolation rate in the infection period. The 1st group experiments are focused on the degree of the variation rate where the medical level index is set as 0.9455 based on the fact. For instance, is set due to the strict isolation measures taken in certain countries, and it is also taken in the following experiments. The 2nd group experiments deal with the comparison of normal contact times and gathering contact times during virus transmission, where the contact times are varied as (3, 6, 9) and varied gathering times. In the 3rd group, the initial infected number is varied as (40, 100, 160) to test the effect of the infected number during the pandemic transmission but set as 100 in the other experiments. The effect of medical level to the virus control and transmission is shown in the 4th group experiments [15].

The virus transmission is estimated via the typical SEIR models with the differential equation in which the parameters are fixed without considering the randomness of the virus transmission process. However, the comprehensive influencing factors are fully considered in the developed TS-SEIR model. Therefore, the parameters involved in the TS-SEIR model can be updated through iteration as a dynamic process for describing the gathering, medical conditions, and infected number in such varied states in the actual virus transmission situations.

3. Simulation Results and Analysis

3.1. Isolation Rate

It is known that isolation is the most effective measure to block the virus transmission in the absence of vaccine [16, 17]. TS-SEIR model is used to simulate the pandemic development trend with different isolation rates (see Figure 3), as set in the 1st group experiments, listed in Table 1. It can be seen from Figure 3 that the higher the isolation ratio, the lower the total number of infected individuals by quarantine during the COVID-19.

Besides, the recovery number with the highest isolation rate can also fast achieve balance compared with those of other isolation rates. For instance, the traffic restrictions were applied throughout China one week after the Wuhan City was closed down on 23rd January 2020, which would greatly decrease the contact possibility and scope. After a 17-day buffer period, the already spread pandemic situation from 14,380 (1st Feb 2020) to 72,436 (17th Feb 2020) has finally been stabilized and then gradually declined, which has proved a positive effect of the high-level quarantine on the spread controlling the pandemic after the outbreak stage (24th January–31st January 2020).

3.2. Infection Number prior to Policy Intervention

The number of infection cases before intervention, which depends largely on the reaction speed and intervention measures put in place by local authorities during the initial stage of the epidemic outbreak [18], plays a crucial role in shaping the consequent infections. The development trend of pandemic under different initial infected cases (40, 60, and 100, respectively) is simulated for illustrative purposes. Intuitively, the more the number of people infected before the intervention, the faster the rate of early cases increased, as illustrated in Figure 4. In other words, the larger the proportion of the population infected, the more difficult the epidemic development can be contained in the latter stage. This highlights the importance of timing of initial control, i.e., early intervention measures being put in place by the health authority. More specifically, the lower the number of initial infection cases, the more likely the peak of infections may be delayed (as well as smaller magnitude of the peak); this may contribute to winning precious time for local authorities to organize and coordinate cross-sectorial actions in fighting against a virulent pandemic such as COVID-19. In short, the initial situation of the pandemic outbreak at a given region has a profound effect on the effectiveness of public health crisis management.

To illustrate the timing effect, here we compare the pandemic data of South Korea and Germany during the early stage of COVID-19 development. It is found that South Korea had more than 200 confirmed patients when it took strict measures such as entry control on 21st February 2020; however, the domestic infections kept increasing rapidly since then. By contrast, Germany adopted strict immigration control measures at a similar date, on 27th February 2020 when only 40 patients were confirmed, and the consequent growth of confirmed cases in Germany was relatively slow for a certain period of time. Thus, a key takeaway is that the number of infected before the draconian pandemic prevention measures being implemented may determine the initial spread rate of the pandemic, leading to markedly different pathways of the epidemic spread and development.

3.3. Degree of Contact with Infected People

Close contact with a large number of susceptible people would highly increase the infection ratio, such as mass gathering in Church worships in Daegu of South Korea [19]. However, it is rather difficult to implement such mandatory measures in certain parts of the world.

Here, the TS-SEIR model is used to simulate the effect of the contact number on the pandemic control with a relatively high isolation rate of 0.9/0.99, where the results of the development trend with different contact number are depicted in Figure 5. The number of the susceptible, exposed, and infected people with/without quarantine would be greatly reduced with less contact number. Hence, the “mobile contact network,” as the mobile infection source, can greatly increase the risk/probability of being infected, especially in intensive social interconnections because of living closely to each other, i.e., India as the second largest population with confirmed infection on 8th September 2020, reported by the WHO [2].

Another example is South Korea; strict policy had been implemented for the importation of pandemic at the early stage of the outbreak. However, domestic mobility and mass gatherings are not prohibited, such as the congregants of the Shincheonji Church of Jesus without appropriate protection measures, resulting in more than 6,000 new confirmed COVID-19 cases in two weeks, as the largest number of infections outside China at 18th February 2020. The explosive growth of infected cases of COVID-19 in South Korea during the first week of nation’s outbreak evidently demonstrates the effectiveness of noninterpersonal contact on the pandemic spread control.

In order to eliminate the risk of fast spread rate in densely populated countries in the Asia Pacific, region strict quarantine measures, i.e., prohibiting unnecessary outside activities, were forced to be implemented in Europe, like Spain and Italy, to stop the sharply interpersonal contact; otherwise, the global situation would have been lost control entirely.

3.4. Medical Capacity

It is known that the medical level plays a crucial role in managing the epidemic emergency while preventing the virus from uncontrolled spreading [20], such as the number of medical personnel involved in the treatment, the length of the confirmed diagnosis, and the number of isolation wards with negative pressure. The influence of the medical level on pandemic development with the TS-SEIR model is simulated under different medical levels, as shown in Figure 6. The higher the medical level index is, the lower the death rate of the patients is expected to be.

At the initial stage of COVID-19 outbreak, the epicentre of the epidemic, the Chinese Hubei Province, was running severely short of sufficient medical resources, and the local medical capacity could not cope with the overwhelming coronavirus outbreak, leading to a relatively high mortality rate during the initial period of pandemic development. The Chinese government realised the serious medical challenge and immediately commanded building makeshift hospitals which began receiving patients on 8th February 2020, and the nationwide medical forces were rapidly mobilised and dispatched to assist with the saturated medical staff in Hubei Province so that the medical resources of Hubei Province were replenished and the treatment of patients soon returned to adequate level. From 8th February 2020, the number of deaths per day started to slowly increase and gradually declined, thanks to the improved medical level, which fully reflected the impact of medical conditions on the mortality rate associated with the pandemic.

The parameter-wise adjustment simulation experiment results clearly demonstrate that the restrained interpersonal contact or interactions are probably the most effective strategy to mitigate the spread of the virus along with sufficient medical treatment to cut off the transmission route and protect the population from the susceptible and the infected [16]. To summarize, a key lesson of the TS-SEIR model suggests that the proactive measures such as prevention and control prior to large-scale epidemic outbreak are more important than the medical countermeasures including diagnostics and clinical treatment for containing the spread of a virulent pandemic such as COVID-19 because the vaccines are still under development and specific cure is not yet available [15].

4. Discussion of the COVID-19 Situation

In the aspect of proactive pandemic management, emergency response should be launched in a timely manner in conjunction with preemptive alarm response and a priori risk assessment system. The primary objective of the preventive intervention measures is to cut off the route of transmission and to protect the susceptible population by eliminating the source of infections based on the integrated approach, such as reducing public gatherings and close interpersonal contact, good hygiene practice (wearing masks outdoors and washing hands frequently), suspending social and religious gathering activities, limiting traffic and mobility, strengthening quarantine at urban and intercity transportation hubs, closing temporarily access to educational institutions and public places, and carrying out thorough disinfection wherever necessary [21]. In addition, increasing the efficiency of isolation of infected population, early diagnosis, and treatment are highly recommended, including routine temperature check of vulnerable population, screening, and monitoring of symptom-susceptible patients and centralised isolation of the suspected cases and confirmed patients.

In the event of high risk of contamination, draconian quarantine measures such as home confinement of the entire infected areas with tight mobility restriction should be implemented (such as the lockdown policies adopted by a number of European countries to prohibit all nonessential outdoor activities) [22]. Both our modelling results and the experience of the Chinese city of Wuhan (one of the main reasons of the high mortality rate of the infected patients in Wuhan during the early stage of the outbreak of COVID-19) were related to the saturation in local medical resources including overwhelmed treatment centres and lack of caring personnel. The situation only started to improve after the Chinese central government sent more than 50,000 medical staff to the city and more than 15 makeshift hospitals started operation to admit tens of thousands of patients with mild symptoms of novel coronavirus for the purpose of isolation. Thus, the treatment capacity is a determining factor in reducing the mortality of infected patients. This invaluable experience gained from China’s pandemic combat may be shared with the countries facing similar challenges), suggesting that all nonnecessary outdoor activities, in particular, the social gatherings, and mobility should be suspended immediately to minimise the interpersonal transmission likelihood. Also, nationwide medical resources have to be mobilised timely to provide necessary medical support to the most infected areas to gain the struggle with the virus transmission. Furthermore, temporary medical isolation centres may be built to increase the patient intake capacity in the event of serious outbreak.

Another controversial issue is related to the prevention and control measure (P&CM) options with regard to the tradeoff between proactive and laissez-faire attitude. Indeed some European countries (e.g., UK and Sweden) had initially tempted the so-called nonpharmaceutical interventions (NPIs), or herd immunity strategy faced with rapidly increased infection cases and high pressure on local medical services (e.g., hospitalisation infrastructures to intake patients with critical situation), due to the fact that the local treatment capacities had been quickly saturated as a consequence of unpreparedness and lack of coordination of social resources while the optimal timing of strong intervention was missed out; however, the ultimate outcome of this passive tackling largely relying on the NPIs finally proved disastrous, and policymakers had to make a policy U-turn by re-enforcing proactive pandemic control measures to prevent the health crisis from deteriorating. This is in line with our modelling results which clearly stressed the importance of initial actions in retarding the peak of infections, which is crucial to managing the long-term epidemic spread as a whole.

The success of joint implementation of the pandemic control guideline depends upon the response of emergency plan such as cutting off the route of transmission and protecting the susceptible population with isolation and treatment of infectious diseases [22]. An integrated resources management and allocation system needs to be established to allow local medical and social workers to jointly complete the work of screening and controlling the sources of infection.

With the rapid spread of novel coronavirus worldwide, this pandemic was no longer a single country’s affair but has been developing into a global pandemic requiring cooperation and control in all countries [23]. The pandemic has so far spread across over 210 countries in the world. In many countries and regions, effective prevention and control are still “new topic” that has not been addressed in the previous pandemic outbreaks such as H1N1 or Ebola virus. The pandemic has not been fully brought under control even in many developed countries such as the United States, where the total infections exceeded 6 million with hundreds of thousands of infected people perished their lives due to ineffective control measures during the initial stage of the outbreak. In Europe, the situations in countries such as United Kingdom, Italy, Germany, France, and Sweden are also concerning, and each of them has recorded more than 250,000 confirmed cases. Based on our modelling analysis of the impact of aggregation factors on the development of the pandemic, it is likely to expect a second-wave outbreak across European and the Middle East countries, should no radical health policy be undertaken by the national governments in due course.

Another critical issue in terms of lifecycle management of pandemic is related to the precautionary measures to prevent a potential second-wave outbreak after the previously infected region or country has reached the stabilised state; e.g., all domestic infectious origins have been cleared off. In other words, how to minimise the risk of importing new cases for the disinfected origin of outbreak poses a serious challenge for local public health. This is particularly relevant to the current situations in East Asian countries with a large degree of increased mobility and high risk of exposure to international travels, which is facing a high risk of imported cases. Any mishandling or improper prevention of the imported cases, a secondary transmission of novel coronavirus, might lead to the outbreak anew, and it is noteworthy that a reinfection case has been recently confirmed in Hong Kong (on 1st March 2020, one imported case from UK has been found in Shenzhen; on 2nd March 2020, Qingtian County, Zhejiang Province, had 7 imported confirmed cases from Italy; before that, Beijing and Ningxia had 3 imported confirmed cases from Iran; these figures seem to be small compared with the overall scale of the existing confirmed cases in China, while they bring new challenges to the pandemic prevention and control in China) [24].

Last but not the least, the transparent and fluent information flow plays a critical role in coordinating pandemic control measures across different stakeholders, especially on social network platforms which may contribute to relaying and processing pandemic risk-related information through the whole process of pandemic control campaign for the well-being of international community. On the one hand, the timely release of information on the outbreak of highly contagious epidemic in an initially infected area is conducive to the prevention and control of infectious diseases in other countries and regions, which is instrumental to the global governance of a COVID-19-alike pandemic. Also, nonbiased and transparent disease-related intelligence information can help scientific understanding of new infectious diseases and pandemic situation of the public and stabilizing social mood, on the other hand. How to prevent the psychological panic from spreading across mass population is essential for health authorities when dealing with an emergent public health crisis with an unprecedented scale such as COVID-19. To gain legitimacy and public trust in local authority and health care institutions, the golden triangle of punctuality, accuracy, and transparency of pandemic-related information should be consistently maintained at all times.

5. Conclusion

The global experience of fighting against the COVID-19 pandemic across various countries clearly indicates that the efficient management of the person-to-person transmission risk of the epidemic requires taking timely preventive measures in global and national public health system, which is in line with our hybrid modelling results. In particular, serious consequences would be generated if the optimal window of intervention were forfeited at the early stage of outbreak, in particular, in the densely populated areas (i.e., the intranetwork transmission probability may be increased exponentially) with scarce medical resources.

The most effective P&CM globally of such infectious diseases as novel coronavirus mainly includes enhancing proactive pandemic management to maximising the efficiency of integrated public health intervention measures [25]. During the initial period of P&CM implementation, quantitative analysis can be carried out with a SEIR-type model introduced in the present study, i.e., TS-SEIR model, by focusing on key control parameters such as the estimated transmission rate, evaluating the possible impact of the pandemic according to the preliminary epidemiological investigation, and predicting the required resources, scale, and time to formulate and implement the intervention plan aided by medical experts and public health policymakers.

The interdisciplinary and intergovernmental coordination is also a prerequisite of the success in combatting the global public health crisis. The shared knowledge of pandemic control parameters may provide useful lessons for countries on the set of increasing risks of large infections. As such, it is an imperative for decision makers from international and national institutions to undertake efficient countermeasures based on robust modelling of the likelihood of pandemic spread pathways, whereby the frontier between policy arena and academia should be connected in an collaborative and open spirit for preserving the global public goods of human health and social welfare.

Data Availability

All data used to support the findings of the study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

JL and YMZ conceived and designed the study. ZGC, QSL, and LJY set up the experiments. XDW, HYN, and CNF collected the data. ZGC, HYN, and CNF analysed the data. JL, YMZ, ZGC, and LJY interpreted the results and wrote the manuscript. Yimin Zhou, Jun Li and Lingjian Ye are contributed equally.

Acknowledgments

This work was partially supported by the National Science Foundation of China (Grant no. 61973296).