Next Article in Journal
Continuous Performance Improvement Framework for Sustainable Wastewater Treatment Facilities in Arid Regions: Case of Wadi Rumah in Qassim, Saudi Arabia
Next Article in Special Issue
Amikacin Therapy in Japanese Pediatric Patients: Narrative Review
Previous Article in Journal
Modeling Community Health with Areal Data: Bayesian Inference with Survey Standard Errors and Spatial Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling Analysis of COVID-19 Transmission and the State of Emergency in Japan

1
College of Engineering and Design, Hunan Normal University, Changsha 410081, China
2
School of Engineering and Technology, University of Washington, Tacoma, WA 98402, USA
3
State Key Laboratory of Developmental Biology of Freshwater Fish, Hunan Normal University, Changsha 410081, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(13), 6858; https://doi.org/10.3390/ijerph18136858
Submission received: 6 April 2021 / Revised: 14 June 2021 / Accepted: 21 June 2021 / Published: 26 June 2021
(This article belongs to the Special Issue State-of-the-Art Infectious Disease Epidemiology in Japan)

Abstract

:
To assess the effectiveness of the containment strategies proposed in Japan, an SEIAQR (susceptible-exposed-infected-asymptomatic-quarantined-recovered) model was established to simulate the transmission of COVID-19. We divided the spread of COVID-19 in Japan into different stages based on policies. The effective reproduction number R e and the transmission parameters were determined to evaluate the measures conducted by the Japanese Government during these periods. On 7 April 2020, the Japanese authority declared a state of emergency to control the rapid development of the pandemic. Based on the simulation results, the spread of COVID-19 in Japan can be inhibited by containment actions during the state of emergency. The effective reproduction number R e reduced from 1.99 (before the state of emergency) to 0.92 (after the state of emergency). The transmission parameters were fitted and characterized with quantifiable variables including the ratio of untracked cases, the PCR test index and the proportion of COCOA app users (official contact confirming application). The impact of these variables on the control of COVID-19 was investigated in the modelling analysis. On 8 January 2021, the Japanese Government declared another state of emergency. The simulated results demonstrated that the spread could be controlled in May by keeping the same strategies. A higher intensity of PCR testing was suggested, and a larger proportion of COCOA app users should reduce the final number of infections and the time needed to control the spread of COVID-19.

1. Introduction

The coronavirus disease (COVID-19) has continued to spread since the first recorded case at the end of 2019 [1,2]. The global cumulative number of reported confirmed cases has risen to over 168 million with 3.6 million deaths [3].
A number of containment strategies have been enforced to control the pandemic by many governments including restrictions on travel and public gatherings, school delays or closings, social distancing, and even lockdown measures [4,5,6,7,8]. Some extreme actions can have a significant impact on the socioeconomic systems [9], which may be continued by some authorities [10]. Thus, it is of importance to optimize the containment policies to stop the spread.
The COVID-19 infection in Japan started in the middle of January 2020 [11]. The Japanese government put forward basic control policies in February and March to curb the early transmission [12]. As the number of confirmed cases rapidly grew, and a large portion of infections could not be tracked, a state of emergency was declared by the central government on 7 April 2020 [12]. After the termination of the first state of emergency in Japan, the spread of COVID-19 was not well controlled. The number of newly confirmed cases quickly rose at the end of 2020. The number of cumulative confirmed cases in Japan hit 401,355 as of 6 February 2021 [13]. Another state of emergency was declared in early 2021 [14].
To assess the containment policies and evaluate the effectiveness of the state of emergency, appropriate modelling studies could offer useful instruction and guidance. A study adopted the susceptible–infected–removed (SIR) transmission model to analyze the effectiveness of the early strategies proposed by the Japan Government [6]. It was concluded that it was necessary to minimize the length of time people stayed in crowded places. However, the study did not highlight testing or isolation measures. In addition, the spread of COVID-19 was affected by the time-varying policies. Thus, the parameters associated with the model should be dynamic, and should be estimated accurately. Another study utilized a susceptible–exposed–infected–recovered (SEIR) model to forecast the peak time of transmission in Japan. However, the rate of transmission from the exposed population was not taken into consideration [15]. Moreover, the infectivity of the large portion of asymptomatic patients should not be neglected.
The SEIR-like models were also applied in some other studies for the analysis of the transmission of COVID-19 and of the control policies [16,17,18,19]. Jiao considered that patients were contagious during the incubation period [17], which was underestimated in many models. It was demonstrated that isolation should be strictly carried out to curb the pandemic. A modified SEIRV model examined the effect of environmental conditions, in which V represented the external concentration of the virus. The conclusion was drawn that long-term strategies were necessary since COVID-19 would, otherwise, remain endemic according to the simulation.
In this study, we established a susceptible–exposed–infected–asymptomatic–qurantine–recovered/removed (SEIAQR) model. Asymptomatic infection [20,21,22] was especially emphasized in the model. The policies of isolation measures and the Polymerase Chain Reaction (PCR) test [23] were quantitatively characterized and analyzed. The spreading period of COVID-19 was partitioned into several stages according to the containment policy announcements, including the declaration of the state of emergency. The parameters in the model associated with the dynamics of COVID-19 transmission were determined by optimizing and fitting the simulated results with the actual number of infected cases reported.
Based on the simulation results, the declaration of the state of emergency was effective, as shown by the reductions in the coefficients of transmission rates and the basic reproduction number. Then, we used the established model to analyze the second declaration of a state of emergency. Finally, the strategies, with different intensities of isolation, PCR test index and the percentage of users of tracking applications, were optimized to offer some instructions to control the COVID-19 pandemic.

2. Methods

2.1. Model Description

The compartmental model established for this study is illustrated in Figure 1. The compartments included susceptible (denoted by S), exposed (denoted by E), infected (denoted by I), asymptomatic (denoted by A), quarantined (denoted by Q) and recovered/removed (denoted by R) populations. The dynamics of the spread of COVID-19 and the interactions between compartments were characterized by the system of Equation (1) in this model. Specifically, asymptomatic infection and infectivity during the incubation period (exposed compartment) were taken into consideration, which improved the accuracy for simulating COVID-19 transmission.
The structure and dynamics of COVID-19 transmission are shown in Figure 1, and the discrete dynamic system is
S ( t + 1 ) = S ( t ) + Λ β e S ( t ) E ( t ) β i S ( t ) I ( t ) β a S ( t ) A ( t ) d S ( t ) E ( t + 1 ) = E ( t ) + β e S ( t ) E ( t ) + β i S ( t ) I ( t ) + β a S ( t ) A ( t ) μ E ( t ) d E ( t ) I ( t + 1 ) = I ( t ) + α μ E ( t ) γ I ( t ) d i I ( t ) A ( t + 1 ) = A ( t ) + ( 1 α ) μ E ( t ) δ A ( t ) ε A ( t ) d A ( t ) Q ( t + 1 ) = Q ( t ) + δ A ( t ) + γ I ( t ) λ Q ( t ) d i Q ( t ) R ( t + 1 ) = R ( t ) + ε A ( t ) + λ Q ( t ) R ( t ) ,
where S, E, I, A, Q and R represent the proportion of individuals that are susceptible, exposed, asymptomatic, documented and recovered/removed at time t (discrete time). The transmission rates from the compartments E, I and A are denoted by β e , β i and β a (per day), respectively. Λ (the natural birth rate) is about 2.0 × 10 5 [24]. d and d i denote the natural death rate and the COVID-19 death rate in Japan, respectively. d was set to be 3.0 × 10 5 [24] and d i was determined as 1.4 × 10 3 by fitting the ratio of deaths to the number of confirmed cases (Figure 2a). The incubation period 1 μ was set to be 5 days [15] ( μ is the rate of transmission from the exposed population to the infected individuals with symptoms). α and 1 α represent the proportion of infected patients with symptoms or asymptomatic patients. γ and δ denote the confirmation rate of the symptomatic infections and asymptomatic infections, respectively. ε (per day) is the recovery rate of asymptomatic infections and λ (per day) is the recovery rate of confirmed infections. The average recovery times 1 ε ( ε is the recovery rate for symptomatic infected patients) and 1 λ ( λ is the recovery rate for asymptomatic infected patients) were determined to be 11.5 days by data fitting (Figure 2b).
The confirmed asymptomatic or symptomatic cases were presumed to be isolated and thus could not spread COVID-19. The reproduction number R e , the expected value of secondary infections produced by one typical infected case, was characterized as the following [25]:
R e = β e Λ d ( μ + d ) + β i α μ Λ d ( γ + d i ) ( μ + d ) + β a ( 1 α ) μ Λ d ( δ + ϵ + d ) ( μ + d ) .

2.2. Determination of the Transmission Parameters

The parameters β e , β i , β a , α , γ , δ needed to be determined based on the real confirmed cases. Usually, the least-squares method is adopted to obtain the best fitting values based on the number of daily confirmed cases. However, the data on the daily reported confirmed cases for COVID-19 fluctuated greatly as can be seen in Figure 1. On the contrary, the change of cumulative confirmed cases was smooth. Therefore, the cumulative confirmed cases were utilized to calculate these parameters.
Set Y r ( t ) as the actual cumulative confirmed cases at time t, where t Z + and
Y r ( t ) = τ = 0 t D ( τ ) ,
where D ( τ ) is the daily confirmed cases on τ Z + . Set Y h ( t ) as the numerical cumulative confirmed cases at time t by solving the SEIAQR model (1), thus
Y r ( t ) = Y h ( t ) + ξ ( t ) ,
where ξ ( t ) is the error at time t and Y h ( t ) is expressed as
Y h ( t ) = τ = 0 t γ I ( τ ) + δ A ( τ ) .
Then, the problem of solving the parameters β e , β i , β a , α , γ , δ was transformed into a nonlinear square optimization problem as
min τ = 0 t ξ 2 ( t ) = min τ = 0 t ( Y r ( t ) Y h ( t ) ) 2 , for 0 < β e , β i , β a , α , γ , δ < 1 .

2.3. COVID-19 Data and Stages for Analysis

In this modelling study, the data were collected from the official website of the Japanese Ministry of Health, Labour and Welfare, including a reported confirmed number of symptomatic or asymptomatic cases, PCR testing, and so forth [26]. The COVID-19 pandemic timeline in Japan was divided into 16 stages, primarily based on the policies proposed by the government (Table 1). The timeline with 16 stages is shown in Figure 3. The dates segregating the stages are marked on the horizontal axis. The transmission parameters were determined dynamically in each stage. In particular, the emphasis in this study was on the analysis of the stages that involved the state of emergency.

3. Results

3.1. Simulated Results and Determined Model Parameters

The simulation of the spread of COVID-19 in Japan was performed using the established model. With the transmission parameters obtained (Table 2), the development of the number of confirmed cases over time fitted well with the reported data (Figure 4a). The relative error between the simulated and real number of confirmed cases did not exceed 5% in all the stages (Figure 4b). The effective reproduction numbers ( R e ) during the stages were calculated with a 95% confidence interval (CI) as shown in Table 3.
Upon the outbreak of COVID-19 in Japan, the transmission parameters β e , β i and β a were 0.0858, 0.2183 and 0.3217, respectively. In stage 2, when the basic policies were announced, the transmission parameters were β e = 0.1635 , β i = 0.4668 and β a = 0.5562 , respectively. The parameters were not lowered, illustrating that early containment measures did not control the disease spread. The confirmation rates increased in stage 2, suggesting a higher proportion of cases were tracked by testing. The effective reproduction number R e increased from 1.52 (95% CI 1.47–2.12) to 1.99 (95% CI 1.90–2.27). After declaring the first state of emergency (stages 3 and 4), the transmission parameters were lowered to β e = 0.0849 , β i = 0.1933 and β a = 0.2357 (stage 3); β e = 0.0534 , β i = 0.1929 and β a = 0.1800 , respectively (stage 4). The reproduction number R e declined to 0.92 (95% CI 0.83–1.00) in stage 3 and 0.59 (95% CI 0.53–0.61) in stage 4. Although a previous study suggested that following the early basic policies proposed in stage 2 could control the transmission of COVID-19 [40], the rapid spread of COVID-19 was not contained until taking effective measures during the state of emergency.
After the first state of emergency was lifted, the spread of COVID-19 continued to develop with reproduction numbers higher than 1 except for in stages 8 and 9 (Table 3). In stage 8, the PCR testing intensity was reinforced [32]. This may have resulted in a more effective control of the spread of COVID-19, compared to the previous stages. In stage 6, an application used for tracking the cases was adopted. The gradual increase of users of the app should play a positive role in controlling the disease.
By the end of 2020, the second wave of infections occurred in Japan. Another state of emergency was declared by the Japanese government on 7 January 2021. After declaring the second state of emergency (stages 15 and 16), the transmission parameters decreased from β e = 0.0633 , β i = 0.2795 , β a = 0.1932 (stage 14) to β e = 0.0489 , β i = 0.1721 and β a = 0.1897 (stage 15), respectively; β e = 0.0444 , β e = 0.1613 and β e = 0.1585 (stage 16), respectively. The reproduction number Re reduced to 0.91 (95% CI 0.89–1.90) in stage 15 and 0.67 (95% CI 0.38–0.71) in stage 16. Compared to the first state of emergency in stage 4, the R e was slightly larger and the confirmation rate was smaller, which may possibly be due to the different intensities of the distancing strategies or variations in the virus’ infectivity [41]. The proportion of asymptomatic cases ranged from less than 10% to over 80% based on different reports [42]. In this model, the proportion of asymptomatic cases ranged from 67% to 90% in various stages (Table 2).

3.2. Fitting of the Transmission Parameters

To investigate the impact of various strategies on the control of the COVID-19 spread, the transmission rates β e , β i and β a were fitted by the defined variables related to the containment policies. The variables used in the fitting included the ratio of untracked to tracked cases at stage p, the PCR test index χ and the proportion of COCOA (contact-confirming application) users in the population η . The official contact-confirming application was put forward in stage 6 and the percentage of users gradually increased. The values of the variables across the stages are shown in Figure 5.
By fitting the transmission parameters with the ratio of untracked to tracked cases p, the PCR test index χ , the proportion of COCOA users in the population η and the confirmation rate of symptomatic infected cases γ , the following equations were determined:
β e = e 2.79 + 0.72 γ 0.92 γ 2 0.57 p + 0.23 p 2 + 3.72 χ 3.69 χ 2 + 0.96 χ 3 3.44 η ,
β i = e 0.63 + 0.37 γ 0.91 γ 2 1.78 p + 0.61 p 2 + 3.25 χ 3.08 χ 2 + 0.75 χ 3 1.72 η ,
β a = e 1.70 + 2.44 γ 2.71 γ 2 0.96 p + 0.39 p 2 + 3.24 χ 3.39 χ 2 + 0.93 χ 3 2.45 η .
The R 2 (goodness-of-fit measure) for the multivariable regressions were 0.926, 0.917, 0.970 and 0.939, respectively. The confirmation parameters δ and γ are associated with the testing. The infected cases with symptoms could be more easily identified but the asymptomatic patients could hardly be identified other than with PCR testing, which should result in difficulties for disease control. The confirmation rate δ characterized the proportion of determined asymptomatic cases by testing the estimated total asymptomatic cases. The confirmation rate δ was fitted for the 16 stages in the following:
δ = 1 1 + e 1.06 1.87 γ + 2.65 γ 2 + 3.65 p 1.16 p 2 2.07 χ + 1.74 χ 2 0.43 χ 3 + 3.81 η .
The R 2 for the multivariable regressions was 0.938. It was found that the confirmation rate could be increased by a higher proportion of COCOA app users. Furthermore, the effects of p and χ were investigated in the following section. The strategies used in stage 16 were optimized by this model.
Given the previous analysis, in the current stage with enhanced testing parameters but where the basic production number has not been reduced to less than 1, we may consider using a combination of increased isolation level and maintaining a moderate testing intensity.

3.3. Prediction of Cumulative Confirmed Cases and Effect of Combined Strategies

In this section, the SEIAQR model was used to predict the spread of COVID-19 the current containment strategies in stage 16 are maintained or modified strategies are adopted. On the basis that the simulation of the previous stages fitted well with the trend of real reported confirmed cases (Figure 4), the impact of the ratio of untracked to tracked cases p, the PCR test index χ and the proportion of COCOA users in the population η on the development of predicted cumulative cases was analyzed by numerical simulation (Figure 6), which provided a sensitivity analysis of the parameters for the model and instructions to optimize the combination of control strategies.
By maintaining the same current containment strategy in stage 16 (the second state of emergency), the cumulative number of infections would reach a stable peak in early May. The predicted cumulative number of confirmed cases would be around 438,000 (red line in Figure 6). If the containment strategies were modified with the same ratio of untracked cases p = 0.9 , the proportion of COCOA users as η = 0.19 and with an increase in the intensity of the PCR testing index χ = 1.5 by 15% or 30% (dashed or solid cyan lines in Figure 6), the epidemic could be effectively controlled by mid-March or early March. The final number of cumulative infections would be around 428,900 and 421,700, respectively.
At the current stage, the proportion of COCOA users is 19.3%. If the proportion η was increased to 22.2% or 25.1% while keeping the same p = 0.9 and χ = 1.5 (dashed or solid blue lines), the epidemic would be effectively controlled in mid-April or late March. The final number of cumulative infections would be around 428,900 and 421,700, respectively. The dissemination of disease information could interact with the spread of the disease [43,44]. It was found that promoting the usage of the COCOA application regarding disease information could help to control the spread of COVID-19. If p was decreased by 15% or 30% (dashed or solid pink lines), the effect on controlling the disease spread was not as significant as strengthening the PCR testing or promoting the application. The number of cumulative infections at the peak would be around 428,900 and 421,700, respectively.
The contour of the effective reproduction number R e was plotted as a function of the proportion of COCOA users in the population η and the PCR test index χ with the other parameters kept the same as in stage 16 (Figure 7). These contour lines demonstrated that Re decreased as the parameter values η and χ increased. In order to control the development of the epidemic ( R e < 1 ), the parameters η and χ must be greater than 30 % and 1.8 , respectively. To achieve better control than the current state of emergency (stage 16), η must be increased by more than 0.65 and χ should be greater than 3.5. The combined control strategy in the following stages should be evaluated conveniently based on the contour.

3.4. Verification of the Model

The data from 15 January 2020 to 5 February 2021 were used for the previous fitting and analysis. To validate the model, the number of confirmed cases from 6 February 2021 to 18 February 2021 was predicted by using the same configurations as in stage 16. As can be seen in Figure 8, under the assumption that the containment strategies should not change during the second state of emergency, the real data on the number of cumulative infections fell in the range of 95% CI of the cumulative confirmed cases predicted by the model. Thus, in addition to the fitting of the data (Figure 3), this model was further validated. The strategies during the second state of emergency could be optimized as guided by the model.

4. Discussion

An SEIAQR model was established to analyze the spread of COVID-19 in Japan and to evaluate the policies to control the pandemic. Generally, SIR or SEIR-like models [6,17,20] are composed of three or four compartments. In contrast, this study adopted six compartments, emphasizing the impact of asymptomatic infection. Infectivity during the incubation period, the rates of natural deaths and deaths due to COVID-19 were incorporated into the system to improve the reliability and authenticity of the model, although this resulted in more difficulties for the determination of the parameters. The proportion of asymptomatic patients was reported as ranging from less than 10% to over 80% based on different reports [42]. The estimation of the asymptomatic percentage ranges from 67% to 90% in this modelling study, which agrees with the previous experimental studies.
The transmission of COVID-19 was divided into 16 stages primarily determined by the policies put forward by the Japanese Government. The declaration of the state of emergency (stages 3, 4, 15, 16) could inhibit the rapid development of the disease spread since the effective reproduction numbers R e were less than 1. The strategies in the state of emergency were effective, which could be shown by the decrease of the transmission parameters. After the lifting of the state of emergency, the reproduction numbers R e became larger than 1. However, in stages 8 and 9, COVID-19 was better controlled compared to the other stages where there was no state of emergency declared, which was probably due to a higher PCR test intensity.
Many variables should affect the spread of COVID-19 but are rarely measured. We selected quantifiable indexes including the ratio of the untracked to tracked cases, the PCR test index and the proportion of COCOA app users. The transmission rates and confirmation rates of asymptomatic patients were fitted using these quantifiable variables. Then, for the second declaration of the state of emergency, the number of cumulative cases was predicted. It was found that increasing the PCR test index by 30% could decrease the number of overall infected cases by around 3 × 10 4 . The state of emergency could also be lifted earlier. Although social distancing measures are critical [45], it was found that increasing the PCR test index should have a more significant effect, followed by the impact of promoting the COCOA app. As suggested in another study, a higher intensity of testing should be implemented during the early spread period [46].
It was inferred that herd immunity in Japan could be achieved without sacrificing the current healthcare system if the basic reproduction number is less than 2 [47]. Our simulation revealed that, under the state of emergency, the reproduction number could be reduced to less than 1. After the lifting of the state of emergency, the reproduction number should be controlled at less than 1 by keeping a relatively high PCR test index and increasing the proportion of COCOA app users.
This study could be improved in the future. First, it was assumed that permanent immunity was achieved after recovery and the impact of the vaccine was also not taken into consideration. The compartment of the immune population could be included in the future. In addition, the stages were segregated based on policies. The stages could be divided with multiple approaches, such as using additionally defined indicators [48] to better evaluate the status of the spread. The transmission parameters could be estimated more accurately since the impact of policies may have some delay. Lastly, we excluded the impact of policies announced by local governments [49], which may result in some deviations in the fitting of parameters.

Author Contributions

Conceptualization, Z.C. and J.P.; methodology, Z.C. and J.P.; software, Z.C. and J.P.; validation, Z.S., X.H. and K.P.; formal analysis, Z.C. and J.P.; investigation, Z.S.; resources, Z.C. and J.P.; data curation, X.H.; writing—original draft preparation, Z.C. and J.P.; writing—review and editing, Z.C., J.P., Z.S., X.H. and K.P.; visualization, Z.C. and J.P.; supervision, Z.C. and J.P.; project administration, Z.C. and J.P.; funding acquisition, K.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National College Student Innovation and Entrepreneurship Training Program (Grant No. S202010542031), the Doctoral Start-Up Foundation of Hunan Normal University (Grant No. 0531120-3827) and the Hunan Provincial Education Department (Grant No. HNKCSZ-2020-0813).

Institutional Review Board Statement

This study did not deal with animal or human subjects. Existing data were used for analysis. Therefore, the IRB statement is not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data, such as the number of positive cases, the number of people performing a PCR test, the number of people requiring inpatient treatment, the number of persons discharged or recovered, the number of the dead and the number of PCR tests performed can be accessed online at https://www.mhlw.go.jp/stf/covid-19/open-data.html. The data on new untraceable cases are available online at https://stopcovid19.metro.tokyo.lg.jp/en/cards/untracked-rate. The data on COCOA users are available online at https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/cocoa_00138.html.

Acknowledgments

We thank all the providers of the open access data for our analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X. Clinical features of patients infected with 2019 novel coronavirus in wuhan, china. Lancet 2020, 395, 497–506. [Google Scholar] [CrossRef] [Green Version]
  2. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R. A novel coronavirus from patients with pneumonia in china, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef] [PubMed]
  3. Global Cumulative Confirmed Covid-19 Deaths and Cases, Who. Available online: https://covid19.who.int/ (accessed on 18 May 2021).
  4. Djaoue, S.; Kolaye, G.G.; Abboubakar, H.; Ari, A.A.A.; Damakoa, I. Mathematical modeling, analysis and numerical simulation of the COVID-19 transmission with mitigation of control strategies used in Cameroon. Chaos Solitons Fractals 2020, 139, 110281. [Google Scholar] [CrossRef] [PubMed]
  5. González-Bustamante, B. Evolution and early government responses to covid-19 in south america. World Dev. 2021, 137, 105180. [Google Scholar] [CrossRef] [PubMed]
  6. Karako, K.; Song, P.; Chen, Y.; Tang, W. Analysis of covid-19 infection spread in japan based on stochastic transition model. Biosci. Trends 2020, 14, 134–138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ng, Y.; Li, Z.; Chua, Y.X.; Chaw, W.L.; Zhao, Z.; Er, B.; Pung, R.; Chiew, C.J.; Lye, D.C.; Heng, D. Evaluation of the effectiveness of surveillance and containment measures for the first 100 patients with covid-19 in singapore—2 January–29 February 2020. Morb. Mortal. Wkly. Rep. 2020, 69, 307–311. [Google Scholar] [CrossRef]
  8. Panovska-Griffiths, J.; Kerr, C.C.; Stuart, R.M.; Mistry, D.; Klein, D.J.; Viner, R.M.; Bonell, C. Determining the optimal strategy for reopening schools, the impact of test and trace interventions, and the risk of occurrence of a second covid-19 epidemic wave in the uk: A modelling study. Lancet Child Adolesc. Health 2020, 4, 817–827. [Google Scholar] [CrossRef]
  9. Wells, C.R.; Townsend, J.P.; Pandey, A.; Moghadas, S.M.; Krieger, G.; Singer, B.; McDonald, R.H.; Fitzpatrick, M.C.; Galvani, A.P. Optimal covid-19 quarantine and testing strategies. Nat. Commun. 2021, 12, 356. [Google Scholar] [CrossRef]
  10. Wilder-Smith, A.; Chiew, C.J.; Lee, V.J. Can we contain the covid-19 outbreak with the same measures as for SARS? Lancet Infect. Dis. 2020, 20, e102–e107. [Google Scholar] [CrossRef] [Green Version]
  11. The First Reported Confirmed Case in Japan Announced by the Minstry of Health, Labour and Welfare. Available online: https://www.mhlw.go.jp/stf/newpage_08906.html (accessed on 18 May 2021).
  12. The Basic Polices for Novel Coronavirus Disease Control by the Government of Japan. Available online: https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/newpage_00032.html (accessed on 18 May 2021).
  13. The Current Covid-19 Situation in Japan, Who. Available online: https://www.who.int/countries/jpn/ (accessed on 18 May 2021).
  14. Declaration of a State of Emergency in Response to the Novel Coronavirus Disease (January 13). Available online: https://japan.kantei.go.jp/ongoingtopics/_00039.html (accessed on 18 May 2021).
  15. Kuniya, T. Prediction of the epidemic peak of coronavirus disease in japan, 2020. J. Clin. Med. 2020, 9, 789. [Google Scholar] [CrossRef] [Green Version]
  16. Carcione, J.M.; Santos, J.E.; Bagaini, C.; Ba, J. A simulation of a covid-19 epidemic based on a deterministic SEIR model. Front. Public Health 2020, 8, 230. [Google Scholar] [CrossRef]
  17. Jiao, J.; Liu, Z.; Cai, S. Dynamics of an SEIR model with infectivity in incubation period and homestead-isolation on the susceptible. Appl. Math. Lett. 2020, 107, 106442. [Google Scholar] [CrossRef]
  18. Prem, K.; Liu, Y.; Russell, T.W.; Kucharski, A.J.; Eggo, R.M.; Davies, N.; Flasche, S.; Clifford, S.; Pearson, C.A.; Munday, J.D. The effect of control strategies to reduce social mixing on outcomes of the covid-19 epidemic in wuhan, china: A modelling study. Lancet Public Health 2020, 5, e261–e270. [Google Scholar] [CrossRef] [Green Version]
  19. Yang, Z.; Zeng, Z.; Wang, K.; Wong, S.-S.; Liang, W.; Zanin, M.; Liu, P.; Cao, X.; Gao, Z.; Mai, Z. Modified SEIR and ai prediction of the epidemics trend of covid-19 in china under public health interventions. J. Thorac. Dis. 2020, 12, 165. [Google Scholar] [CrossRef]
  20. Mizumoto, K.; Kagaya, K.; Zarebski, A.; Chowell, G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance 2020, 25, 2000180. [Google Scholar] [CrossRef] [Green Version]
  21. Rothe, C.; Schunk, M.; Sothmann, P.; Bretzel, G.; Froeschl, G.; Wallrauch, C.; Zimmer, T.; Thiel, V.; Janke, C.; Guggemos, W. Transmission of 2019-nCoV infection from an asymptomatic contact in Germany. N. Engl. J. Med. 2020, 382, 970–971. [Google Scholar] [CrossRef] [Green Version]
  22. Sutton, D.; Fuchs, K.; D’alton, M.; Goffman, D. Universal screening for SARS-CoV-2 in women admitted for delivery. N. Engl. J. Med. 2020, 382, 2163–2164. [Google Scholar] [CrossRef]
  23. Peto, J. Covid-19 mass testing facilities could end the epidemic rapidly. BMJ 2020, 368, m1163. [Google Scholar] [CrossRef] [Green Version]
  24. Death Rate of Japan. Available online: https://www.macrotrends.net/countries/JPN/japan/death-rate (accessed on 18 May 2021).
  25. Allen, L.J.; Van den Driessche, P. The basic reproduction number in some discrete-time epidemic models. J. Differ. Equ. Appl. 2008, 14, 1127–1147. [Google Scholar] [CrossRef]
  26. Official Statistics of Covid-19 from the Japanese Government. Available online: https://www.mhlw.go.jp/stf/covid-19/open-data.html (accessed on 18 May 2021).
  27. Basic Policies to Contain the Covid-19 Spread in Japan. Available online: https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/newpage_00032.html (accessed on 18 May 2021).
  28. Declaration of the State of Emergency in Several Prefectures. Available online: https://japan.kantei.go.jp/ongoingtopics/_00018.html (accessed on 18 May 2021).
  29. The Declaration of the State of Emergency Nationwide. Available online: https://japan.kantei.go.jp/ongoingtopics/_00020.html (accessed on 18 May 2021).
  30. The Promotion of Contact Confirming Application. Available online: https://www.mhlw.go.jp/stf/newpage_11954.html (accessed on 18 May 2021).
  31. Salivapcr Testing for Asymptomatic Infection. Available online: https://www.mhlw.go.jp/stf/newpage_12488.html (accessed on 18 May 2021).
  32. Strengthened the pcr Testing. Available online: https://www.mhlw.go.jp/stf/newpage_12892.html (accessed on 18 May 2021).
  33. Announced the Employment Subsidy. Available online: https://www.mhlw.go.jp/stf/houdou/houdou_list_202008.html (accessed on 18 May 2021).
  34. Extended Applications for Subsidies for Business Suspension. Available online: https://www.mhlw.go.jp/stf/newpage_13766.html (accessed on 18 May 2021).
  35. Ensured the Vacation and Welfare of the Patients. Available online: https://www.mhlw.go.jp/stf/newpage_13997.html (accessed on 18 May 2021).
  36. Abandoned the Policy of the Two-Week Quarantine Upon Entry. Available online: https://www.sohu.com/a/423581473_120707383 (accessed on 18 May 2021).
  37. Signed an Agreement on Provision of Information Sharing of Cluster Countermeasures for Covid-19. Available online: https://www.mhlw.go.jp/stf/newpage_14984.html (accessed on 18 May 2021).
  38. Announcement to Send Additional Medical Staff from the Self-Defense Forces. Available online: https://www.asahi.com/articles/ASND73G6NND7IIPE007.html (accessed on 18 May 2021).
  39. The Declaration of the Second State of Emergency. Available online: https://www.mhlw.go.jp/stf/newpage_15933.html (accessed on 18 May 2021).
  40. Declared the Second Nationwide State of Emergency. Available online: https://japan.kantei.go.jp/ongoingtopics/_00039.html (accessed on 18 May 2021).
  41. Korber, B.; Fischer, W.M.; Gnanakaran, S.; Yoon, H.; Theiler, J.; Abfalterer, W.; Hengartner, N.; Giorgi, E.E.; Bhattacharya, T.; Foley, B.; et al. Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus. Cell 2020, 182, 812–827.e19. [Google Scholar] [CrossRef]
  42. Meyerowitz, E.A.; Richterman, A.; Bogoch, I.I.; Low, N.; Cevik, M. Towards an accurate and systematic characterisation of persistently asymptomatic infection with SARS-CoV-2. Lancet Infect Dis. 2021, 21, e163–e169. [Google Scholar] [CrossRef]
  43. González-Rubio, J.; Navarro-López, C.; López-Nájera, E.; López-Nájera, A.; Jiménez-Díaz, L.; Navarro-López, J.D.; Nájera, A. A systematic review and meta-analysis of hospitalised current smokers and covid-19. Int. J. Environ. Res. Public Health 2020, 17, 7394. [Google Scholar] [CrossRef] [PubMed]
  44. Zhan, X.X.; Liu, C.; Zhou, G.; Zhang, Z.K.; Sun, G.Q.; Zhu, J.J.; Jin, Z. Coupling dynamics of epidemic spreading and information diffusion on complex networks. Appl. Math. Comput. 2018, 332, 437–448. [Google Scholar] [CrossRef] [PubMed]
  45. Inoue, H. Japanese strategy to COVID-19: How does it work? Glob. Health Med. 2020, 2, 131–132. [Google Scholar] [CrossRef] [PubMed]
  46. Chen, Y.T.; Yen, Y.F.; Yu, S.H.; Su, E.C. An Examination on the Transmission of COVID-19 and the Effect of Response Strategies: A Comparative Analysis. Int. J. Environ. Res. Public Health 2020, 17, 5687. [Google Scholar] [CrossRef]
  47. Akamatsu, T.; Nagae, T.; Osawa, M.; Satsukawa, K.; Sakai, T.; Mizutani, D. Model-based analysis on social acceptability and feasibility of a focused protection strategy against the COVID-19 pandemic. Sci. Rep. 2021, 11, 2003. [Google Scholar] [CrossRef]
  48. Nakano, T.; Ikeda, Y. Novel Indicator to Ascertain the Status and Trend of COVID-19 Spread: Modeling Study. J. Med. Internet. Res. 2020, 22, e20144. [Google Scholar] [CrossRef]
  49. Tanaka, T.; Yamaguchi, T.; Sakamoto, Y. Estimation of the percentages of undiagnosed patients of the novel coronavirus (SARS-CoV-2) infection in Hokkaido, Japan by using birth-death process with recursive full tracing. PLoS ONE 2020, 15, e0241170. [Google Scholar] [CrossRef]
Figure 1. The graphical illustration of the SEIAQR model.
Figure 1. The graphical illustration of the SEIAQR model.
Ijerph 18 06858 g001
Figure 2. Data fitting for the recovery rate and death rate of COVID-19 in Japan. (a) Data fitting for the recovery rate. (b) Data fitting for the death rate.
Figure 2. Data fitting for the recovery rate and death rate of COVID-19 in Japan. (a) Data fitting for the recovery rate. (b) Data fitting for the death rate.
Ijerph 18 06858 g002
Figure 3. Segregation of the spread of COVID-19 in Japan for the modelling study.
Figure 3. Segregation of the spread of COVID-19 in Japan for the modelling study.
Ijerph 18 06858 g003
Figure 4. The simulation of confirmed cases over time by the model. (a) The simulated number of cases fitted well with the real accumulated confirmed cases. (b) The relative error between the simulated and real data.
Figure 4. The simulation of confirmed cases over time by the model. (a) The simulated number of cases fitted well with the real accumulated confirmed cases. (b) The relative error between the simulated and real data.
Ijerph 18 06858 g004
Figure 5. The variables for the fitting of the model parameters and characterization of the strategies. (a) The ratios of untracked to tracked cases in the stages. (b) PCR test indexes in the stages. (c) The proportions of COCOA app users in the stages.
Figure 5. The variables for the fitting of the model parameters and characterization of the strategies. (a) The ratios of untracked to tracked cases in the stages. (b) PCR test indexes in the stages. (c) The proportions of COCOA app users in the stages.
Ijerph 18 06858 g005
Figure 6. The prediction of the development of cumulative confirmed cases when changing the intensity of the strategies. Red line represents using the current strategy in the state of emergency; Pink lines represent using the modified strategies and lowering the ratio of untracked cases; cyan lines represent using the modified strategies with strengthened PCR tests; blue lines represent using the modified strategies alongside promoting the COCOA application.
Figure 6. The prediction of the development of cumulative confirmed cases when changing the intensity of the strategies. Red line represents using the current strategy in the state of emergency; Pink lines represent using the modified strategies and lowering the ratio of untracked cases; cyan lines represent using the modified strategies with strengthened PCR tests; blue lines represent using the modified strategies alongside promoting the COCOA application.
Ijerph 18 06858 g006
Figure 7. The contour of reproduction number R e with different η and χ .
Figure 7. The contour of reproduction number R e with different η and χ .
Ijerph 18 06858 g007
Figure 8. Verification of the model by comparing the predicted and real number of infections.
Figure 8. Verification of the model by comparing the predicted and real number of infections.
Ijerph 18 06858 g008
Table 1. The stages of the spread of COVID-19 in Japan.
Table 1. The stages of the spread of COVID-19 in Japan.
StagePeriodNote
12020.01.06–03.27The first case appeared in Japan [11].
22020.03.28–04.06A series of basic policies were announced [27].
32020.04.07–04.15Declared the state of emergency in several prefectures [28].
42020.04.16–05.25Declared the first nationwide state of emergency [29].
52020.05.26–06.18Lifted the state of emergency [12].
62020.06.19–07.16Promoted cellphone app for contact information [30].
72020.07.17–08.07Started saliva PCR testing to detect asymptomatic infection [31].
82020.08.08–08.24Strengthened the testing intensity [32].
92020.08.25–09.25Announced employment subsidy [33].
102020.09.26–10.08Extended applications for subsidies for business suspension [34].
112020.10.09–11.01Ensured the vacation and welfare of the patients [35].
122020.11.02–11.20Abandoned two-week quarantine policy upon entry [36].
132020.11.21–12.07Signed an agreement on provision of information sharing of cluster countermeasures for COVID-19 [37].
142020.12.08–2021.01.07Announced to send additional medical staffs [38].
152021.01.08–01.15Declared another state of emergency in several prefectures [39].
162021.01.16–02.18Declared the second nationwide state of emergency [40].
Table 2. Parameters in different stages fitted by the model.
Table 2. Parameters in different stages fitted by the model.
Stage β e (per day) β i (per day) β a (per day) α γ δ
10.08580.21830.32170.10.52440.2057
20.16350.46680.55620.17120.63740.3701
30.08490.19330.23570.20740.47710.3829
40.05340.19290.18000.22740.46030.5485
50.09710.34040.31030.12320.31880.3695
60.12130.42350.38640.12950.32380.3304
70.10150.36250.31750.12510.38390.2930
80.03350.13940.13710.10020.65360.1275
90.04780.13170.11650.10190.91070.1020
100.04910.13050.19760.11280.13740.1877
110.07420.21870.21130.11210.74530.1655
120.08670.27270.29290.22150.42330.1700
130.07240.25740.16230.10060.93780.1273
140.06330.27950.19320.10420.05880.1625
150.04890.17210.18970.10010.56870.1969
160.04440.16130.15850.10000.61520.2652
Table 3. Effective reproduction number R e determined by the model.
Table 3. Effective reproduction number R e determined by the model.
Stage R e 95% CIStage R e 95% CI
11.52(1.47, 2.12)90.86(0.83, 0.88)
21.99(1.90, 2.27)101.03(0.96, 1.12)
30.92(0.83, 1.00)111.2(1.16, 1.26)
40.59(0.53, 0.61)121.52(1.46, 1.63)
51.23(1.19, 1.27)131.12(1.01, 1.15)
61.61(1.56, 1.67)141.44(1.31, 1.92)
71.39(1.37, 1.44)150.91(0.89, 1.90)
80.81(0.79, 0.89)160.67(0.38, 0.71)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Z.; Shu, Z.; Huang, X.; Peng, K.; Pan, J. Modelling Analysis of COVID-19 Transmission and the State of Emergency in Japan. Int. J. Environ. Res. Public Health 2021, 18, 6858. https://doi.org/10.3390/ijerph18136858

AMA Style

Chen Z, Shu Z, Huang X, Peng K, Pan J. Modelling Analysis of COVID-19 Transmission and the State of Emergency in Japan. International Journal of Environmental Research and Public Health. 2021; 18(13):6858. https://doi.org/10.3390/ijerph18136858

Chicago/Turabian Style

Chen, Zhongxiang, Zhiquan Shu, Xiuxiang Huang, Ke Peng, and Jiaji Pan. 2021. "Modelling Analysis of COVID-19 Transmission and the State of Emergency in Japan" International Journal of Environmental Research and Public Health 18, no. 13: 6858. https://doi.org/10.3390/ijerph18136858

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop