Simulation and predictions of a new COVID-19 pandemic wave in Ukraine with the use of generalized SIR model =========================================================================================================== * Igor Nesteruk ## ABSTRACT A new wave of the COVID-19 pandemic in Ukraine, which began in the summer of 2021, was characterized by almost exponential growth of smoothed daily numbers of new cases. This is a matter of great concern and the need to immediately predict the dynamics of the epidemic in order to assess the possible maximum values of new cases, the risk of infection and the number of deaths. The generalized SIR-model and corresponding parameter identification procedure was used to simulate and predict the dynamics of this epidemic wave. The new COVID-19 epidemic wave in Ukraine will begin to subside in mid-October 2021, but its duration will be quite long. Unfortunately, new cases may appear by the summer of 2022. Keywords * COVID-19 pandemic * epidemic dynamics in Ukraine * mathematical modeling of infection diseases * SIR model * parameter identification * statistical methods ## Introduction The COVID-19 pandemic dynamics in Ukraine was discussed in [1-14]. To predict the first wave of the pandemic, the classical SIR model [15-17] and the statistics-based method of its parameter identification [18] were used. To simulate new epidemic waves, a numerical method of their detection [3, 19], a generalized SIR-model [20] and a corresponding parameter identification procedure [21] were developed. In particular, eleven epidemic waves in Ukraine were calculated [4, 7-10]. The calculations of the 11-th pandemic wave (based on the accumulated numbers cases reported by Ukrainian national statistics [22, 23] in the period May 23 – June 5, 2021) predicted the end of this wave on August 25, 2021 with the number of cases 2,226,797 (see [10]). As of August 25, 2021 the real number of cases accumulated in Ukraine was 2,278,171 (see Table 1). It means that the predicted saturation level was exceeded only 2.26% (after 81 days of observation). The obtained high accuracy of the method allows us to hope for a fairly accurate forecast for the next 12th epidemic wave in Ukraine, to which this study is devoted. View this table: [Table 1.](http://medrxiv.org/content/early/2021/10/18/2021.10.13.21264949/T1) Table 1. Cumulative numbers of laboratory-confirmed Covid-19 cases in Ukraine *V**j* in the summer and autumn of 2021 according to the national statistics, [22, 23]. ### Data We will use the data set regarding the accumulated numbers of laboratory-confirmed COVID-19 cases in Ukraine from national sources [22, 23]. The corresponding numbers *V**j* and moments of time *t**j* (measured in days) are shown in Table 1. The values for the period *T**c11*: May 23 – June 5, 2021 have been used in [10] for SIR simulations of the tenth epidemic wave in Ukraine. Here we use the dataset, corresponding to the period *T**c12*: September 29 – October 12, 2021 to simulate the 12-th wave. Other *V**j* and *t**j* values will be used to control the accuracy of predictions. To estimate the mortality rate (ratio of accumulated number of cases *V**j* to the accumulated number of deaths *D**j*), let us take figures *D**j* corresponding different days: 52,286 (June 26, 2021); 52,665 (July 13, 2021); 52,981 (August 2, 2021); 53,789 (August 30, 2021); 54,550 (September 15, 2021); 57,526 (October 5, 2021); 59,052 (October 11, 2021), [22, 23]. Taking corresponding *V**j* values from Table 1, we can calculate the mortality rates *m**i* =*D**j* *1000/*V**j* (per thousand of cases) for the listed days: 23.40; 23.49; 23.50; 23.52; 23.40; 23.17; 23.05. Thus the mortality rate is rather stable (its variation during June to October 2021 is only 0.47). We will use the average value *m*= 23.36 to predict the number of deaths in Ukraine during the new 12-th pandemic wave. ### Generalized SIR model and data smoothing procedure The generalized SIR-model relates the number of susceptible *S*, infectious *I* and removed persons *R* for a particular epidemic wave *i* [8, 20]. The exact solution of the set of non-linear differential equations uses the function ![Formula][1] corresponding to the number of victims or the cumulative confirmed number of cases [8, 20]. Its derivative: ![Formula][2] yields the estimation of the average daily number of new cases. When the registered number of victims *V**j* is a random realization of its theoretical dependence (1), the exact solution presented in [8, 20] depends on five parameters (*α**i* is one of them). The details of the optimization procedure for their identification can be found in [21]. Since daily numbers of new cases are random and characterized by some weekly periodicity, we will use the smoothed daily number of accumulated cases: ![Formula][3] and its numerical derivative: ![Formula][4] to estimate the smoothed number of new daily cases [3, 4, 9, 19]. ## Results and discussion The optimal values of the general SIR model and other characteristics of the 12-th pandemic wave in Ukraine are calculated and listed in Table 2 (middle column). The same characteristics of the 11-th wave calculated in [10] are presented in the last column of Table 2 for comparison. The corresponding SIR curves are shown in Figure by blue and black lines. It can be seen that the optimal values of SIR parameters are very different for 11-th and 12-th pandemic waves. Close values were obtained only for the average times of spreading the infection 1/ *ρ*11 =4.1 days (in comparison with 1/ *ρ*12 =3.3 days for the twelfth wave). The assessment of the 12-th epidemic wave duration corresponding the moment when the number of infectious persons becomes less that unit) is very pessimistic (June 16, 2022). A similar long epidemic wave was also predicted for India [24]. ![Figure.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/18/2021.10.13.21264949/F1.medium.gif) [Figure.](http://medrxiv.org/content/early/2021/10/18/2021.10.13.21264949/F1) Figure. 11-th and 12-th waves of the COVID-19 pandemic in Ukraine in the summer and autumn of 2021. The results of SIR simulations of the eleventh [10] and twelfth waves are shown by black and blue lines respectively. Numbers of victims *V(t)=I(t)+R(t)* – solid lines; numbers of infected and spreading *I(t)* multiplied by 5 – dashed; derivatives *dV/dt* (eq. (2)) multiplied by 100 – dotted. The magenta line represents the estimation of the accumulated number of deaths during the 12-th epidemic wave multiplied by 10. Black and blue “circles” correspond to the accumulated numbers of cases registered during the periods of time taken for SIR simulations (May 23 – June 5, 2021) and September 29 – October 12, 2021, respectively. “Stars” corresponds to *V**j* values beyond these time periods. The red “crosses” show the first derivative (4) multiplied by 100. View this table: [Table 2.](http://medrxiv.org/content/early/2021/10/18/2021.10.13.21264949/T2) Table 2. Optimal values of parameters and other characteristics of the eleventh and twelfth COVID-19 pandemic wave in Ukraine. The saturation level (final size) of the 12-th wave can approach 3 millions. As of the end of 2021 the accumulated numbers of laboratory-confirmed cases and deaths can exceed 2,943,700 and 68,764, respectively (see blue and magenta solid lines in Figure). The numbers of infectious persons and average daily new cases (see blue dashed and dotted lines) will stop to increase around 17 and 14 October 2021, respectively with the maximal values of 45,545 and 13,930, respectively (see blue dashed and dotted lines). As of the end of 2021 these values are expected to be around 505 and 113, respectively (see blue dashed and dotted lines). To illustrate the accuracy of simulations the accumulated number of cases (red “stars”) and averaged daily numbers of new cases (eq. (4), red “crosses”) are also presented in Figure. Comparisons with corresponding solid and dotted lines show that the theoretical curves are in good agreement with the observations of the pandemic dynamics (markers). To estimate the accumulated number of deaths during the 12-th wave, we have multiplied the theoretical *V(t)=I(t)+R(t)* values by the coefficient *m*= 23.36 calculated before. The magenta line represents expected accumulated numbers of deaths multiplied by 10. ## Conclusions The new COVID-19 epidemic wave in Ukraine will begin to subside in mid-October 2021, but its duration will be quite long. Unfortunately, new cases may appear by the summer of 2022. ## Data Availability data is listed in the text ## Acknowledgements The author is grateful to Oleksii Rodionov for his help in collecting and processing data. * Received October 13, 2021. * Revision received October 13, 2021. * Accepted October 18, 2021. * © 2021, 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. ## References 1. 1.Nesteruk I., Kydybyn I., Demelmair G. GLOBAL STABILIZATION TRENDS OF COVID-19 PANDEMIC. KPI Science News. 2020. No. 2. pp. 55–62. DOI: 10.20535/kpi-sn.2020.2.205124 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/kpi-sn.2020.2.205124&link_type=DOI) 2. 2.Nesteruk I. Simulations and predictions of COVID-19 pandemic with the use of SIR model. Innov Biosyst Bioeng. 2020. vol. 4. no. 2. 110–121. doi: 10.20535/ibb.2020.4.2.204274. [http://ibb.kpi.ua/article/view/204274](http://ibb.kpi.ua/article/view/204274) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/ibb.2020.4.2.204274&link_type=DOI) 3. 3.Nesteruk I. Coronasummer in Ukraine and Austria. [Preprint.] ResearchGate. 2020 June. DOI: 10.13140/RG.2.2.32738.56002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.13140/RG.2.2.32738.56002&link_type=DOI) 4. 4.Nesteruk I. COVID-19 pandemic dynamics. Springer Nature. 2021. DOI: 10.1007/978-981-33-6416-5. [https://link.springer.com/book/10.1007/978-981-33-6416-5](https://link.springer.com/book/10.1007/978-981-33-6416-5) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-981-33-6416-5&link_type=DOI) 5. 5.Pardhan S,. Drydakis N. Associating the Change in New COVID-19 Cases to GDP per Capita in 38 European Countries in the First Wave of the Pandemic. Front Public Health. 2021 Jan 20;8:582140. doi: 10.3389/fpubh.2020.582140. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpubh.2020.582140&link_type=DOI) 6. 6.Chintala S., Dutta R., Tadmor D. COVID-19 spatiotemporal research with workflow-based data analysis. Infect Genet Evol. 2021 Mar;88:104701. doi: 10.1016/j.meegid.2020.104701. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.meegid.2020.104701&link_type=DOI) 7. 7.Nesteruk I.. Benlagha N. PREDICTIONS OF COVID-19 PANDEMIC DYNAMICS IN UKRAINE AND QATAR BASED ON GENERALIZED SIR MODEL. Innov Biosyst Bioeng. 2021. vol. 5. no. 1. p. 37–46. DOI: 10.20535/ibb.2021.5.2.230487 [http://ibb.kpi.ua/article/view/230487](http://ibb.kpi.ua/article/view/230487) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/ibb.2021.5.1.228605&link_type=DOI) 8. 8.Nesteruk I. Visible and real sizes of new COVID-19 pandemic waves in Ukraine Innov Biosyst Bioeng. 2021. vol. 5. no. 2. pp. 85–96. DOI: 10.20535/ibb.2021.5.2.230487 [http://ibb.kpi.ua/article/view/230487](http://ibb.kpi.ua/article/view/230487) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/ibb.2021.5.2.230487&link_type=DOI) 9. 9.Nesteruk I. Detections and SIR simulations of the COVID-19 pandemic waves in Ukraine. Comput. Math. Biophys. 2021;9:46–65. [https://doi.org/10.1515/cmb-2020-0117](https://doi.org/10.1515/cmb-2020-0117) 10. 10.Nesteruk I. Will a natural collective immunity of Ukrainians restrain new COVID-19 waves? MedRxiv, July 2021. DOI: 10.1101/2021.07.20.21260840 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wNy4yMC4yMTI2MDg0MHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMTgvMjAyMS4xMC4xMy4yMTI2NDk0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. 11.Nesteruk I, Rodionov O. The impact of demographic factors on the accumulated number of COVID-19 cases per capita in Europe and the regions of Ukraine in the summer of 2021. medRxiv, July 2021. DOI: 10.1101/2021.07.04.21259980 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wNy4wNC4yMTI1OTk4MHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMTgvMjAyMS4xMC4xMy4yMTI2NDk0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 12. 12.Nesteruk I, Rodionov O, Nikitin AV. The impact of seasonal factors on the COVID-19 pandemic waves. medRxiv, August 2021. DOI:10.1101/2021.08.06.21261665 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wOC4wNi4yMTI2MTY2NXYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTAvMTgvMjAyMS4xMC4xMy4yMTI2NDk0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 13. 13.Nesteruk I. Comparison of the COVID-19 pandemics dynamics in Ukraine and Israel in the summer of 2021. Research Gate, September 2021. [https://www.researchgate.net/publication/354497030\_Comparison\_of\_the\_COVID-19\_pandemics\_dynamics\_in\_Ukraine\_and\_Israel\_in\_the\_summer\_of\_2021](https://www.researchgate.net/publication/354497030\_Comparison\_of\_the\_COVID-19\_pandemics\_dynamics\_in\_Ukraine\_and\_Israel\_in\_the_summer_of_2021) 14. 14.Nesteruk I. Effects of testing and vaccination levels on the dynamics of the COVID-19 pandemic and the prospects for its termination. MedRxiv, September 2021. DOI: [https://doi.org/10.1101/2021.09.20.21263823](https://doi.org/10.1101/2021.09.20.21263823) 15. 15.Kermack WO, McKendrick AG. A Contribution to the mathematical theory of epidemics. J Royal Stat Soc Ser A. 1927;115:700–21. 16. 16.Murray JD. Mathematical Biology I/II. New York: Springer; 2002. 17. 17.Langemann D, Nesteruk I, Prestin J. Comparison of mathematical models for the dynamics of the Chernivtsi children disease. Mathematics in Computers and Simulation. 2016;123:68–79. DOI: 10.1016/j.matcom.2016.01.003 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.matcom.2016.01.003&link_type=DOI) 18. 18.Nesteruk I. Statistics based models for the dynamics of Chernivtsi children disease. Naukovi Visti NTUU KPI. 2017;5:26–34. DOI: 10.20535/1810-0546.2017.5.108577 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/1810-0546.2017.5.108577&link_type=DOI) 19. 19.Nesteruk I. Identification of the New Waves of the COVID-19 Pandemic. In book: COVID-19 Pandemic Dynamics, Springer Nature, 2021. DOI: 10.1007/978-981-33-6416-5_8. [https://link.springer.com/chapter/10.1007/978-981-33-6416-5\_8](https://link.springer.com/chapter/10.1007/978-981-33-6416-5_8) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-981-33-6416-5_8&link_type=DOI) 20. 20.Nesteruk I. General SIR Model and Its Exact Solution. In book: COVID-19 Pandemic Dynamics, Springer Nature, 2021. DOI: 10.1007/978-981-33-6416-5_9. [https://link.springer.com/content/pdf/10.1007%2F978-981-33-6416-5\_9.pdf](https://link.springer.com/content/pdf/10.1007%2F978-981-33-6416-5_9.pdf) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-981-33-6416-5_9&link_type=DOI) 21. 21.Nesteruk I. Procedures of Parameter Identification for the Waves of Epidemics. In book: COVID-19 Pandemic Dynamics, Springer Nature, 2021. DOI: 10.1007/978-981-33-6416-5_10. [https://link.springer.com/chapter/10.1007%2F978-981-33-6416-5\_10](https://link.springer.com/chapter/10.1007%2F978-981-33-6416-5_10) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-981-33-6416-5_10&link_type=DOI) 22. 22.Coronavirus in Ukraine - Statistics - Map of infections, graphs [Internet]. Index.minfin.com.ua. 2021. [https://index.minfin.com.ua/ua/reference/coronavirus/ukraine/](https://index.minfin.com.ua/ua/reference/coronavirus/ukraine/) 23. 23.Cabinet of Ministers of Ukraine – Home [Internet]. [https://www.kmu.gov.ua/](https://www.kmu.gov.ua/) 24. 24.Nesteruk I. The COVID-19 pandemic storm in India. [Preprint] medRxiv Posted May 08, 2021. doi: [https://doi.org/10.1101/2021.05.06.21256523](https://doi.org/10.1101/2021.05.06.21256523) [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif