The Dawn is Coming —— the Description and Prediction of Omicron SARS CoV-2 Epidemic Outbreak in Shanghai by Mathematical Modeling ===================================================================================================================================== * Minghao Jiang * Guoyu Meng * Geng Wu ## Abstract The Omicron outbreak in Shanghai has been going on for more than a month although the lockdown of the city has been enforced. In this work, we present a simple and easy-to-use model to evaluate and forecast the spreading trend. After data fitting, we concluded that the maximum overall number infected in Shanghai will exceed 300,000 and the turning point will be in the coming days around April 13-15, 2022. Our model would provide us important parameters such as the average quarantine ratio, the detection interval from being infected to being tested positive, and the spreading coefficient. We found that it would take 1 to 2 days before the infected patients could be tested positive, which was responsible for the rapid Omicron spread in Shanghai. Furthermore, the quarantine ratio was greater than 1, indicating that the lockdown policy was effective, making quelling the epidemic possible. More importantly, our model could not only enable us to predict the spreading precisely, but also lay the theoretical foundation for better policy making in the prevention of Omicron spread. KEYWORDS * SARS-CoV-2 Omicron * Shanghai outbreak * Prediction model ## Introduction The SARS-CoV-2 Omicron variant, reported as a new variant of COVID-19 in November, 2021[1], has broken out in Shanghai recently. This variant carries mutations to help the virus to resist or escape immunity provided by COVID-19 vaccines[1-3]. To prevent the spread of infection, the local government of Shanghai municipal decided to enforce a lockdown on March 28. However, in the past week, there were still over 20,000 Omicron cases confirmed per day, including both diagnosed and asymptomatic. The situation was quite severe when compared to other places across mainland China during the same period or that in Shanghai several months ago. Here, we propose a new mathematical model to evaluate the epidemic trend of Omicron spreading in Shanghai and make forecasts for the near future. ## Materials and methods ### Data used in modeling The case numbers were collected manually from officially released reports (**Table 1**). View this table: [Table 1.](http://medrxiv.org/content/early/2022/04/19/2022.04.13.22273788/T1) Table 1. The number of COVID-19 cases in Shanghai ### Mathematical modeling First, a differential equation is deducted to describe the spreading of COVID-19 as following. Let us denote *s*(*t*) as the number of people whose nucleic acid tests are positive, then ![Formula][1] in which *p*(*t*) is the total number of patients infected at time *t*, and *τ* is the time from being infected to being tested positive. Since the patients infected during the time *τ* will not be recorded as nucleic acid positive, the diagnosed or recorded patients (*s*(*t*) above) at time *t* is equal to the number of infected populationat time *t* − *τ*. During the time *dt*, the number of newly infected patients *dp*(*t*) is ![Formula][2] Where *k* is the spreading coefficient representing the average number of people one unquarantined patient can infect in a unit interval, and *f*(*t*) is the number of unquarantined patients at time *t. f*(*t*) is equal to *p*(*t*) subtracted by the number of quarantined patients *q*(*t*): ![Formula][3] We further assume that *q*(*t*) (including both nucleic acid positive patients and their close contacts) is proportional to *s*(*t*), and let *u* be the average quarantine ratio: ![Formula][4] Therefore, we get the final differential equation: ![Formula][5] In this equation, there are three unknown parameters: the spreading coefficient *k*, the average quarantine ration *u*, and the detection interval *τ*. These parameters may vary during different periods, *e*.*g*., before or after the city lockdown. Using the number of reported cases during a particular time in Shanghai, these three parameters can be estimated by data fitting. As a next step, we could forecast the trend of Omicron spreading in Shanghai with these estimated parameters. ## Results The Omicron epidemic in Shanghai has begun since early March, 2022. Until March 12, the number of confirmed cases per day did not exceed 100 (**Figure 1**). Due to the surge in COVID-19 positives since mid-March, the Shanghai municipal government made the lockdown decision on March 28 to control the spread. So, we fitted the data of total number of reported cases before and after 28 March to figure out to what extent the more stringent policy could affect, and we found that the fitting between the reported data and the curve calculated from our equation was quite close (**Figure 2**). We found that, in order to fit the data better, it was needed to alter the parameter *τ* in the two scenarios: i) *τ* = 2 (days) before the city lockdown and ii) *τ* = 1 (day) after the city lockdown. This can be construed as that the more stringent quarantine policy after the city lockdown, with nucleic acid tests for all citizens once or twice per day, reduced the detection interval *τ*. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/19/2022.04.13.22273788/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/04/19/2022.04.13.22273788/F1) Figure 1. The total number of confirmed cases in Shanghai since March1. The blue points and curve represent the data before the official lockdown, and the red ones represent the data after the lockdown. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/19/2022.04.13.22273788/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/04/19/2022.04.13.22273788/F2) Figure 2. Data fitting. The data were fitted separately before and after the lockdown. The blue points and curve represented the reported number of total confirmed cases, and the yellow ones represented the number calculated from our differential equation. The parameters were rounded to 3 significant digits. In addition, we found that the spreading coefficient *k* decreased after the city lockdown, from 1.63 to 1.49, approximately (**Figure 2**). Furthermore, in the early stage of the epidemic, the quarantine ratio *u* was 1.32 > 1, which means that almost all the patients tested positive as well as their close contacts are quarantined. In the later stage of the epidemic, the quarantine ratio*u* got smaller, ∼1.04. The reduced *u* may serve as a sign of the arriving turning point. Since *p*(*t*), the cumulative number of confirmed cases, is a non-decreasing function, its derivative, ![Graphic][6], has to be greater than or equal to 0. Therefore, *u* must be smaller than or equal to ![Graphic][7]. At the beginning of the epidemic, Omicron was spreading quickly, *p*(*t*)was much larger than *p*(*t* − *τ*), so *u* was relatively bigger. As a result of reaching the turning point, Omicron spreading has slowed down, thus *p*(*t*) was not much different from ![Graphic][8] was close to 1, which made the quarantine ratio *u* smaller than that in the early stage of the epidemic. With the three parameters estimated, the trend of Omicron spreading in the near future was predicted. We concluded that there will be more than 300,000 people infected in Shanghai during this Omicron outbreak (**Figure 3A**). And the number of confirmed cases per day will start to decrease around April 13 (**Figure 3B**), which is the so-called turning point. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/19/2022.04.13.22273788/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2022/04/19/2022.04.13.22273788/F3) Figure 3. Epidemic trend forecast. (A) The forecasted total number of confirmed cases, with the predicted data points colored in blue. (B) The reported daily increased number of cases (red) and those for upcoming days predicted by our model (blue). We found that the average detection interval *τ* was 1 or 2 days when fitting the known data, which means that, there would be 1 to 2 days of delay from infection to positive nucleic acid test, during which time the unnoticed patients were still infectious. We believe that this is one of the reasons why Omicron spread rapidly. Furthermore, the quarantine ratio *u* was found to be greater than 1, which means that most of the infected and close contacts had been quarantined, indicating the lockdown and other stringent policies were effective, making appearance of the turning point and quelling the epidemic possible. Moreover, the apparent R0 value, which is the number of people a patient could infect, could be obtained by multiplying the spreading coefficient *k* and the detection interval *τ*.After the lockdown, the value of R0 dropped from 3.27 to 1.49, providing further supporting evidence for the positive effects of the lockdown. ## Discussion The sudden Omicron epidemic outbreak in Shanghai has brought panic to the public, and great loss to business. In order to evaluate and forecast the spreading trend based on available confirmed data, we proposed a novel mathematical model, which took the government policies into account, which could provide us important parameters describing the Omicron epidemic such as the spreading coefficient *k*, the quarantine ratio *u*, and the detection interval *τ*. The predicted number of overall cases and the expected time of turning point may help the government to make judgment on the spreading and to revise the policies accordingly. Our model is easy to be applied to describe other epidemic diseases, and can help the society to avoid panic and build up confidence to fight with COVID-19. For example, the Guangzhou city was the latest metropolitan city in China that declared city-wise lockdown on Apr 11. We believe that our mathematical model presented in this report would help the public to have a better grasp of the current pandemic spread, and would undoubtedly instill confidence and calm that are urgently needed in the caught-up fight of COVID-19 in China. ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## Author contributions Conceived and designed the experiments: GM and GW. Performed the experiments: MJ. Analyzed the data, preparation of figures manuscripts and wrote the paper: MJ, GM, GW. Project supervision: GM. All authors read and approved the final manuscript. ## Compliance with ethics guidelines ### Conflict of interest This article does not contain any studies with human or animal subjects. The authors declare no competing interests. ## Acknowledgements This work was supported by research grants 81970132, 81770142, 81800144, 31800642, 82100153 and 82104582 from National Natural Science Foundation of China, a research grant 20JC1410600 from Shanghai Science and Technology Committee, Shanghai Guangci Translational Medical Research Development Foundation, a research grant 20152504 from “Shanghai Municipal Education Commission—Gaofeng Clinical Medicine Grant Support”, “The Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institute of Higher Learning”, Samuel Waxman Cancer Research Foundation. We thanks Professor Zhu Chen, Shanghai Institute of Hematology, for his kind support of this work. * Received April 13, 2022. * Revision received April 13, 2022. * Accepted April 19, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Team, C.C.-R., SARS-CoV-2 B.1.1.529 (Omicron) Variant -United States, December 1-8, 2021. MMWR Morb Mortal Wkly Rep, 2021. 70(50): p. 1731–1734. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/mmwr.mm7050e1&link_type=DOI) 2. 2.Karim, S.S.A. and Q.A. Karim, Omicron SARS-CoV-2 variant: a new chapter in the COVID-19 pandemic. Lancet, 2021. 398(10317): p. 2126–2128. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(21)02758-6/ATTACHMENT/FB9D1AF8-6927-48BD-B7B7-A5055C3743D5/MMC1.PDF&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F19%2F2022.04.13.22273788.atom) 3. 3.Del Rio, C., S.B. Omer, and P.N. Malani, Winter of Omicron-The Evolving COVID-19 Pandemic. JAMA, 2022. 327(4): p319–320. [1]: /embed/graphic-3.gif [2]: /embed/graphic-4.gif [3]: /embed/graphic-5.gif [4]: /embed/graphic-6.gif [5]: /embed/graphic-7.gif [6]: /embed/inline-graphic-1.gif [7]: /embed/inline-graphic-2.gif [8]: /embed/inline-graphic-3.gif