First peak of COVID-19 outbreak in Japan might pass as of April 26, 2020 ======================================================================== * Junko Kurita * Tamie Sugawara * Yasushi Ohkusa ## Abstract **Background** The COVID-19 outbreak has two inconsistent phenomena: its reproduction number is almost two; and it shows earlier and lower peaks for new cases and the total number of patients. **Object** To resolve this difficulty, we constructed a mathematical model that can explain these phenomena. **Method** We applied a susceptible–infected–recovery model with the proportion of asymptomatic patients among the infected (q) as a key parameter for estimation as well as the basic reproduction number (R). **Results** We observed the first outbreak peak in Japan on April 3 for those infected on March 29. Their R and q were estimated respectively as 2.048 and 99.987%. Discussion and Conclusion: By introducing a very high proportion of asymptomatic cases, two inconsistent phenomena might be resolved. This hypothesis should be verified through additional study. Keywords * COVID-19 * voluntary event cancellation * outbreak in Japan * SIR model * asymptomatic case * peak * COVID-19 * SIR model * proportion of asymptomatic cases ## Introduction The initial case of COVID-19 in Japan was a patient who showed symptoms when returning from Wuhan, China on January 3, 2020. Subsequently, as of April 26, 2020, the Ministry of Labour, Health and Welfare (MLHW) in Japan announced that there were 7,741 cases in Japan, including asymptomatic people, but excluding those infected on a large cruise ship: the Diamond Princess [1]. In fact, the first peak was observed at the end of March. Subsequently, despite a declining epidemic curve, the Japanese government declared a state of emergency on April 7. The COVID-19 outbreak has two associated and inconsistent phenomena: its reproduction number is almost two; and it shows earlier and lower peaks of new cases and the total number of patients. That reproduction number of two indicates that the peak will be reached when half of the population is infected. However, the total number of patients in Japan until the peak has been approximately 7000. A similar phenomenon was confirmed in Wuhan, China. How can one reconcile these two inconsistent characteristics of the outbreak of COVID-19? To resolve the difficulty, we constructed a mathematical model to explain the two phenomena. ## Methods We applied a simple susceptible–infected–recovery (SIR) model [2–4] to the epidemic curve of Japan and its 120 million population. We assume an incubation period that conforms to the empirical distribution in the early stage of the outbreak. The proportion of asymptomatic people among those infected (q) is a crucially important parameter for estimation in the model, as is the basic reproduction number (R). For simplicity, we assumed no infectiousness among asymptomatic cases. As data, we used the number of symptomatic patients reported by MLHW for January 14 – April 24 published [1] as of April 26. We excluded some patients who were presumed to be persons infected abroad or as passengers on the Diamond Princess. They were presumed not to be community-acquired in Japan. For onset dates of some symptomatic patients that were unknown, we estimated their onset date from an empirical distribution with duration extending from the onset to the report date among patients for whom the onset date had been reported. Onset dates were sometimes unreported. Therefore, we estimated the onset date. To do so, we inferred a distribution extending from onset to reporting based on patients for whom onset dates were available. Then we applied this empirical distribution to patients for whom onset dates were not available. Letting *f*(*k*) represent this empirical distribution and letting *Nt* denote the number of patients for whom onset dates were not available published at date *t*, then the number of patients for whom the onset date was known is *t*-1. The number of patients for whom onset dates were not available was estimated as *f*(1)*Nt*. Similarly, the number of patients with onset date *t*-2 and whose onset dates were not available were estimated as *f*(2)*Nt*. Therefore, the total number of patients for whom the onset date was not available, given an onset date of *s*, was estimated as Σ*k*=1*f*(*k*)*Ns*+*k* for the long duration passing from *s*. Moreover, the reporting delay for published data from MHLW might be considerable. In other words, if *s*+*k* was larger than that in the current period *t*, then *s*+*k* represents the future for period *t*: therefore, *Ns*+*k* was not observable. Such a reporting delay leads to underestimation bias in the number of patients. Therefore, it must be adjusted as ![Graphic][1]. Similarly, patients for whom the onset dates were available are expected to be affected by the reporting delay. Therefore, we have ![Graphic][2], where *Ms*|*t* represents the reported number of patients for whom onset dates were within period *s*, extending until the current period *t*. We sought R and q to fit the data to minimize the sum of the absolute values of discrepancies among the bootstrapped epidemic curve and the fitted values. The estimated distributions of the three reproduction numbers were calculated using 10,000 fully replicated iterations of bootstrapping for the empirical epidemic curve of symptomatic patients. ## Results During January 14 – April 26 in Japan, 12,936 community-acquired cases were identified, excluding asymptomatic cases. Figure 1 depicts the empirical distribution of incubation period among 91 cases for which the exposure date and onset date were published by MHLW. The mode of the distribution was six days; the average was 6.6 days. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/01/2020.04.26.20081315/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/05/01/2020.04.26.20081315/F1) Figure 1: Epidemic curve and infection dates published by the Ministry of Labour, Health and Welfare (MLHW), Japan Figure 1 presents the epidemic curve and infection dates of onset patients inferred from the empirical distribution of the incubation period depicted in Figure 2. The curve implies that the first peak of the outbreak began in Japan on April 3, with those infected on March 29. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/01/2020.04.26.20081315/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/05/01/2020.04.26.20081315/F2) Figure 2: Empirical distribution of duration from onset to report by MLHW, Japan. Figure 3 depicts an empirical distribution of the duration of onset to report. The maximum delay was 30 days. The R was estimated as 2.048, with a range of [2.048, 2.024]. However, q was estimated as 99.987142%, with a range of [99.986819%, 99.987142%]. Figure 4 portrays a fitted line based on estimated R and q and their ranges. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/01/2020.04.26.20081315/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/05/01/2020.04.26.20081315/F3) Figure 3: Empirical distribution of incubation period published by MLHW, Japan. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/01/2020.04.26.20081315/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/05/01/2020.04.26.20081315/F4) Figure 4: Fitted line base on the model and epidemic curve. ## Discussion We observed the first peak of the outbreak as occurring with the April 3 onset date and March 29 infection date. Because Figure 3 depicts almost all cases reported up to 30 days, the first peak might not change over time. Therefore, we conclude that the first peak has already passed. It is noteworthy that no countermeasure was implemented when the peak was reached. In fact, during February 27 – March 20, schools were closed; events were cancelled voluntarily. Moreover, a state of emergency was declared on April 7. During March 21 through April 6, the only measure was a recommendation that large events be cancelled. Therefore, the peak infection date of the end of March was not regarded as a result of strong countermeasures. Nevertheless, school closure and voluntary event cancellations have apparently been effective because the epidemic curve is below the fitted line, as shown by the model in Figure 4. We applied a simple SIR model including the proportion of asymptomatic cases that had not been incorporated into the model to date. Earlier studies [5–7] have estimated R for COVID-19 as 2.24–3.58 in Wuhan. Our obtained R was similar but slightly smaller. Conversely, an earlier study [8] estimated R in Japan as 0.6. That figure might lead to misguided policies for countermeasures in Japan. If results of the present study are correct, then they would necessitate adherence to contact tracing to detect clusters when more than 60 million people would become infected. We estimated the proportion of asymptomatic cases among infected people as 99.985%, which means that the ratio of asymptomatic cases to symptomatic cases was 6482:1. In other words, for every symptomatic patient confirmed, 6482 asymptomatic cases are assumed to exist. An earlier study [9] revealed the proportion of asymptomatic cases among infected people as 3/23. That proportion represents a huge difference from our results presented herein. A possible reason might be laboratory test procedures. Earlier studies used PCR tests that detected infection at the time of specimen collection. Therefore, a negative PCR test result does not contradict a past infection. To confirm our results through laboratory testing, complete laboratory-based surveillance using an IgG antibody test, not PCR, in the community is expected to be necessary. Such a trail had started in New York City and found that about 15% of residents were positive by the antibody testing [10]. At that time, In New York State, the prevalence was reported as 0.88%. Therefore, q was 94% in New York state. It may be similar with our results of q. Recently, Keio University Hospital reported that about 6% of newly administrated and non-COVID-19 patients were infected asymptomatically during April 13–19, 2020 [11]. During that period, Figure 1 implies that 2095 patients were reported, representing 11% (=2095 × 6985/120 millions) of the total population. The figure was slightly higher than 6%. However, they were not healthy people. Therefore, their protection against infection might be more of a concern; their usual activity might be less than that of healthy people. Therefore, incidence among those patients might be lower. The present study has some limitations. First, these results obtained with a very high proportion of asymptomatic cases are merely hypothetical. This hypothesis should be verified through additional study. Secondly, the peak on April 2 might be only the first peak: a second and third wave might occur. Moreover, a second or third wave peak might be higher than the first peak. In fact, the peak of the entire outbreak might eventually be that second or third peak. One must particularly consider that April is the first month of the school year and the fiscal year in Japan. Consequently, the population in Japan reshuffles many of its activities at this time. Many new students and new employees move to Tokyo from outside Tokyo. Also, residents of Tokyo move away from the city. Therefore, the outbreak can be expected to increase again in middle or late April. Evaluation of the outbreak of COVID-19 in Japan in its entirety will have to be postponed until the outbreak ends. Thirdly, as described above, school closure and voluntary event cancellation are apparently effective. Therefore, their effects must be incorporated as effects influencing the model. Assessment of those effects constitutes our next challenge for future research efforts. Fourthly, for simplicity, we assumed no infectiousness of asymptomatic cases. However, another earlier study [9] found that such cases have equal infectiousness to that of symptomatic cases. We must modify the model for future studies. Fifthly, though we obtained very high proportion of asymptomatic cases, it may include some effect of under-ascertainment [12]. In other words, very mild cases probably were included as asymptomatic cases those who were not counted as infected in the present study. Isolation of those very mild cases and truly asymptomatic cases should be the next challenge. ## Conclusion We found that the first peak of COVID-19 outbreak was April 2. The central government of Japan declared an emergency on April 7. However, as shown in Figure 1, the number of symptomatic patients had already declined. Probably, the declaration served to mitigate the outbreak. By introducing a very high proportion of asymptomatic cases, two inconsistent phenomena might be resolved as a result of this study: the high reproduction number and low peak. The results also explain aspects of the outbreak in Wuhan. Nevertheless, it is currently only a hypothesis. Its validity must be verified using data of outbreaks in other countries including China and the United States, or prefectures in Japan. The present study is based on the authors’ opinions, but does not reflect any stance or policy of their professionally affiliated bodies. ## Data Availability Japan Ministry of Health, Labour and Welfare. Press Releases. [https://www.mhlw.go.jp/stf/newpage\_10723.html](https://www.mhlw.go.jp/stf/newpage_10723.html) ## Acknowledgments We acknowledge the great efforts of all staff at public health centers, medical institutions, and other facilities who are fighting the spread and destruction associated with COVID-19. * Received April 26, 2020. * Revision received April 26, 2020. * Accepted May 1, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## Reference 1. 1.Japan Ministry of Health, Labour and Welfare. Press Releases. [https://www.mhlw.go.jp/stf/newpage\_10723.html](https://www.mhlw.go.jp/stf/newpage_10723.html) (in Japanese) [accessed on 26 March, 2020] 2. 2.Kurita J, Sugawara T Ohkusa Y. Preliminary evaluation of voluntary event cancellation as a countermeasure against the COVID-19 outbreak in Japan as of 11 March, 2020. [https://www.medrxiv.org/content/10.1101/2020.03.12.20035220v1](https://www.medrxiv.org/content/10.1101/2020.03.12.20035220v1) 3. 3.Kurita J, Sugawara T Ohkusa Y. Forecast of the COVID-19 outbreak, collapse of medical facilities, and lockdown effects in Tokyo, Japan. [https://medrxiv.org/cgi/content/short/2020.04.02.20051490v1](https://medrxiv.org/cgi/content/short/2020.04.02.20051490v1) 4. 4.Ohkusa Y, Sugawara T, Taniguchi K, Okabe N. Real-time estimation and prediction for pandemic A/H1N1(2009) in Japan. J Infect Chemother. 2011;17:468–72. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21387184&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20081315.atom) 5. 5.Zhao S, Lin Q, Ran J, Musa SS, Yang G, Wang W, Lou Y, Gao D, Yang L, He D, Wang M. Preliminary Estimation of the Basic Reproduction Number of Novel Coronavirus (2019-nCoV) in China, From 2019 to 2020: A Data-Driven Analysis in the Early Phase of the Outbreak. Int J Infect Dis 2020[Online ahead of print] 6. 6.Liu Y, Gayle AA, Wilder-Smith A, Rockly J. The reproductive number of COVID-19 is higher compared to SARS coronavirus. J Travel Med. 2020.DOI: 10.1093/jtm/taaa021 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jtm/taaa021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32052846&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20081315.atom) 7. 7.Lai C, Shih T, Ko W, Tang H, Hsueh P. Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and Coronavirus disease-2019 (COVID-19): The Epidemic and the Challenges. Int J Antimicrob Agents. DOI:10.1016/j.ijantimicag.2020.105924 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijantimicag.2020.105924&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32081636&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20081315.atom) 8. 8. H Nishiura, H Oshitani, T Kobayashi, T Saito, T Sunagawa, Takaji Wakita, MHLW COVID-19 Response Team, Mi Suzuki. Closed environments facilitate secondary transmission of coronavirus disease 2019 (COVID-19). doi: https://doi.org/10.1101/2020.02.28.20029272 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMi4yOC4yMDAyOTI3MnYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMDEvMjAyMC4wNC4yNi4yMDA4MTMxNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 9. 9.Kimball A, Hatfield KM, Arons M, James A, Taylor J, Spicer K, Bardossy AC, Oakley LP, Tanwar S, Chisty Z, Bell JM, Methner M, Harney J, Jacobs JR, Carlson CM, McLaughlin HP, Stone N, Clark S, Brostrom-Smith C, Page LC, Kay M, Lewis J, Russell D, Hiatt B, Gant J, Duchin JS, Clark TA, Honein MA, Reddy SC, Jernigan JA; Public Health ? Seattle & King County; CDC COVID-19 Investigation Team. Asymptomatic and Presymptomatic SARS-CoV-2 Infections in Residents of a Long-Term Care Skilled Nursing Facility - King County, Washington, March 2020. MMWR Morb Mortal Wkly Rep. 2020 Apr 3;69(13):377–381. doi: 10.15585/mmwr.mm6913e1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/mmwr.mm6913e1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32240128&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F01%2F2020.04.26.20081315.atom) 10. 10.Reuter. New York survey suggests nearly 14% in state may have coronavirus antibodies. [https://www.reuters.com/article/us-health-coronavirus-usa-new-york/new-york-survey-suggests-nearly-14-in-state-may-have-coronavirus-antibodies-idUSKCN2252WN](https://www.reuters.com/article/us-health-coronavirus-usa-new-york/new-york-survey-suggests-nearly-14-in-state-may-have-coronavirus-antibodies-idUSKCN2252WN) [accessed on April 26, 2020] 11. 11.The Asahi Shimbun. Number of new patients testing positive for virus alarms hospital. [http://www.asahi.com/ajw/articles/13321354](http://www.asahi.com/ajw/articles/13321354) [accessed on April 26, 2020] 12. 12.Omori R, Mizumoto K, Ascertainment rate of novel coronavirus disease (COVID-19) in Japan. medRxiv 2020.03.09.20033183; doi: https://doi.org/10.1101/2020.03.09.20033183 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4wOS4yMDAzMzE4M3YxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMDEvMjAyMC4wNC4yNi4yMDA4MTMxNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif