Final sizes and durations of new COVID-19 pandemic waves in Poland and Germany predicted by generalized SIR model ================================================================================================================= * Igor Nesteruk ## ABSTRACT New waves of the COVID-19 pandemic in Europe, which began in the autumn of 2021, are a matter of great concern and the need to immediately predict the epidemic dynamics 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 new epidemic waves in Poland and Germany. Results of calculations show that new cases in these countries will not stop to appear in 2022. Keywords * COVID-19 pandemic * epidemic waves * epidemic dynamics in Poland * epidemic dynamics in Germany * mathematical modeling of infection diseases * SIR model * parameter identification * statistical methods ## Introduction The COVID-19 pandemic dynamics in Poland was discussed in [1-9]. The pandemic spreading in Germany was investigated in [10-22]. In particular, to predict the first pandemic wave in Germany, the classical SIR model [23-25] and the statistics-based method of its parameter identification [26] were successfully used in [10, 14]. In the second SIR simulation for Germany [10, 14], the data set about the accumulated numbers of COVID-19 cases for the period: April 9–29, 2020 (according to WHO daily reports, [27]) was used. It was predicted that the first epidemic wave in Germany could stop on August 4, 2020 with the number of cases 176,952. The real accumulated number of cases was 212,022 (see [27]). It means that the accuracy of prediction was 16.5% after 95 days of pandemic monitoring (paper [10] was published in May 2020). Changes in quarantine restrictions, social behavior, testing and vaccination levels, appearance of new strains, etc. cause the changes in the epidemic dynamics, i.e., so known new waves. To simulate new epidemic waves, a numerical method of their detection [28] was developed. In particular, thirteen epidemic waves were simulated for Ukraine and six pandemic waves for the whole world [14, 31, 34] with the use of a generalized SIR-model [29], and corresponding parameter identification procedures [30]. The correct selection of all the COVID-19 pandemic waves in Poland and Germany is still ahead. Here we will consider the new waves which occurred in October-December 2021 and predict their durations and final sizes with the use of generalized SIR-model. ### Data We will use the data sets regarding the accumulated numbers of laboratory-confirmed COVID-19 cases in Poland and Germany reported by COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU), [32]. The corresponding numbers *V**j j* and moments of time *t**j* (measured in days) are shown in Tables 1 and 2 for the period of July to December 2021. Here we use the datasets, corresponding to the period of November 22 to December 5, 2021 to simulate the 4th wave in Poland and the 5-th wave in Germany. Other *V**j* and *t**j* values will be used to control the accuracy of predictions and the pandemic dynamics. View this table: [Table 1.](http://medrxiv.org/content/early/2021/12/15/2021.12.14.21267771/T1) Table 1. Cumulative numbers of laboratory-confirmed Covid-19 cases in Poland in the summer and autumn of 2021 according to the JHU datasets, [32]. View this table: [Table 2.](http://medrxiv.org/content/early/2021/12/15/2021.12.14.21267771/T2) Table 2. Cumulative numbers of laboratory-confirmed Covid-19 cases in Germany in the summer and autumn of 2021 according to the JHU datasets, [32]. ### 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*, [14, 29]. 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 laboratory-confirmed number of cases [14, 29]. 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 [14, 29] depends on five parameters (*α**i* is one of them). The details of the optimization procedure for their identification can be found in [30]. 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 [10, 14, 28]. ## Results and discussion The optimal values of the general SIR model and other characteristics of the 4th epidemic wave for Poland and the 5th one in Germany are calculated and listed in Table 3. For comparison, the results for the 13th pandemic wave in Ukraine (calculated in [31]) are also presented in the last column. The corresponding SIR curves are shown in the figure by red lines for Poland and black lines for Germany. The pandemic dynamics in both countries is rather irregular. In particular, the smoothed numbers of new daily cases (4) have multiple minima after August 15, 2021 (see red and black “crosses” in the figure). It means that pandemic dynamics changed very rapid and many short epidemic waves occurred. In particular, a distinct epidemic wave in Germany in August – September 2021 is visible (see black “crosses”). That is why we use different numbers for autumn epidemic waves in Poland and Germany. View this table: [Table 3.](http://medrxiv.org/content/early/2021/12/15/2021.12.14.21267771/T3) Table 3. Optimal values of parameters and other characteristics the COVID-19 pandemic waves in Poland, Germany and Ukraine For correct SIR simulations, we need to separate periods with more or less stable dynamics first (see [14, 34]). Unfortunately, for the period taken for SIR simulations of pandemic waves in Poland and Germany (November 22 – December 5, 2021), some local minima of the smoothed numbers of new daily cases (4) are visible (see red and black “crosses” in the figure). This fact reduces the accuracy of SIR simulations and predictions. In particular, some problems with the identifications of the optimal values of model parameters occurred. Table 3 demonstrates that the optimal values of SIR parameters are very different. Close values were obtained only for the average times of spreading the infection 1/ *ρ**i* for Germany and Ukraine. The assessments of the pandemic wave durations (corresponding the moment when the number of infectious persons becomes less that unit) are very pessimistic (April, 2023 for Germany and September 2025 for Poland). Similar long epidemic waves were predicted for India [33] and the whole world [31] with the use of the same generalized SIR model. In particular, if the global situation with vaccination, testing and treatment will not change, the COVID-19 pandemic could continue for another ten years [31]. ![Figure.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/15/2021.12.14.21267771/F1.medium.gif) [Figure.](http://medrxiv.org/content/early/2021/12/15/2021.12.14.21267771/F1) Figure. The COVID-19 pandemic waves in Poland (red) and Germany (black) in the autumn of 2021. The results of SIR simulations of the 4th wave in Poland and the 5th wave in Germany are shown by red and black lines, respectively. Numbers of victims *V(t)=I(t)+R(t)* – solid lines; numbers of infected and spreading *I(t)* multiplied by 10 – dashed; derivatives *dV/dt* (eq. (2), multiplied by 100) – dotted. “Circles” correspond to the accumulated numbers of cases registered during the period of time taken for SIR simulations. “Stars” correspond to *V**j* values beyond this time period. “Crosses” show the first derivative (4) multiplied by 100. “Stars” and “crosses” in the figure illustrate the accuracy of simulations for the accumulated number of cases and the averaged daily numbers of new cases (eq. (4)). Comparisons with corresponding solid and dotted lines show that the theoretical estimations are consistent with observations for Poland. In the case of Germany the theoretical simulation are close to observations after Nobember 1, 2021. The theory predicts the decrease of the daily numbers of new cases after December 1, 2021 in both countries (see black and red dotted lines). The rate of this decreasing will be higher in Germany and by the end of 2021, the daily number of new cases in Germany is expected to be lower than in Poland. A significant increase in contacts during the New Year and Christmas holidays or/and the wide spread of new coronavirus strains (e.g., omicron) could disrupt these positive trends. SIR simulations allows us to estimate the numbers of infectious persons *I(t)* and the probabilities to meet such a person in a population. Dashed lines show that the numbers of infectious in Poland are higher than in Germany and are expected to be much higher by the end of 2021. Taking into account the difference in populations of these countries, we can conclude, that contacts during the holidays in Poland are much more dangerous. Unfortunately, the general SIR model cannot predict the emergence of new epidemic waves. It simulates the dynamics for only the period with constant epidemic conditions. Therefore, permanent monitoring of the number of new cases is needed to determine changes in the epidemic dynamics. After that it will be possible to perform new simulations by means of the generalized SIR model with calculation and use of new values of its parameters. ## Conclusions The generalized SIR-model and corresponding parameter identification procedure was used to simulate and predict the dynamics of two new epidemic waves in Poland and Germany. If the situation with vaccination, testing and treatment will not change, the new COVID-19 will not stop to appear in 2022. ## Data Availability All data produced in the present work are contained in the manuscript ## Acknowledgements The author is grateful to Oleksii Rodionov for his help in collecting and processing data. * Received December 14, 2021. * Revision received December 14, 2021. * Accepted December 15, 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.Raciborski F, Pinkas J, Jankowski M, et al. Dynamics of the coronavirus disease 2019 outbreak in Poland: an epidemiological analysis of the first 2 months of theepidemic. Pol Arch Intern Med. 2020; 130: 615–621. doi:10.20452/pamw.15430 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20452/pamw.15430&link_type=DOI) 2. 2. M Wielechowski, K Czech, L Grzeda. Decline in Mobility: Public Transport in Poland in the time of the COVID-19 Pandemic. Economies 2020, 8(4), 78; [https://doi.org/390/economies804007810.3](https://doi.org/390/economies804007810.3) 3. 3.Gujski M. et al. Epidemiological Analysis of the First 1389 Cases of COVID-19 in Poland: A Preliminary Report J. Clin. Med. 2020, 9(10), 3245; [https://doi.org/10.3390/9103245jcm](https://doi.org/10.3390/9103245jcm) 4. 4.Raciborski F. et al. Dynamics of the coronavirus disease 2019 outbreak in Poland: an epidemiological analysis of the first 2 months of the epidemic. Pol Arch InternMed. 2020 Aug 27;130(7-8):615-621. doi: 10.20452/pamw.15430. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20452/pamw.15430&link_type=DOI) 5. 5.Pinkas J. et al. Public Health Interventions to Mitigate Early Spread of SARS-CoV-2 in Poland. Med Sci Monit. 2020 Apr 13; 26:e924730. doi: 10.12659/MSM.924730. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.12659/MSM.924730&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F15%2F2021.12.14.21267771.atom) 6. 6.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) 7. 7.Piasecki T, Mucha PB, Rosinska M. A new SEIR type model including quarantine effects and its application to analysis of Covid-19 pandemia in Poland in March-April 2020. Research gate, May 2020. [https://www.researchgate.net/publication/341785478\_A\_new\_SEIR\_type\_model\_including\_quarantine\_effects\_and\_its\_application\_to\_analysis\_of\_Covid-19\_pandemia\_in\_Poland\_in\_March-April\_2020](https://www.researchgate.net/publication/341785478\_A\_new\_SEIR\_type\_model\_including\_quarantine\_effects\_and\_its\_application\_to\_analysis\_of\_Covid-19\_pandemia\_in\_Poland_in_March-April_2020) 8. 8.Arti M.K., Antoni Wilinski. Mathematical modeling and estimation for next wave of COVID-19 in Poland. November 2021.Stochastic Environmental Research and Risk Assessment. DOI: 10.1007/s00477-021-02119-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00477-021-02119-5&link_type=DOI) 9. 9.Borucka A, Kozlowski E, Piotr Oleszczuk P, Mazurkiewicz D. The Use of Mathematical Models Describing the Spread of Covid-19 in Strategic State Security Management. European Research Studies Journal 2020(special issue 3):82–98. DOI: 10.35808/ersj/1855 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.35808/ersj/1855&link_type=DOI) 10. 10.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) 11. 11. Maria Vittoria Barbarossa et al. Modeling the spread of COVID-19 in Germany: Early assessment and possible scenarios. September 2020. PLoS ONE 15(9):e0238559 DOI:10.1371/journal.pone.0238559 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0238559&link_type=DOI) 12. 12. Kai Wang et al. Modelling the initial epidemic trends of COVID-19 in Italy, Spain, Germany, and France. November 2020. PLoS ONE 15(11):e0241743. DOI: 10.1371/journal.pone.0241743 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0241743&link_type=DOI) 13. 13. Goetz Th, Heidrich P. Early stage COVID-19 disease dynamics in Germany: models and parameter identification. December 2020. Journal of Mathematics in Industry 10(1). DOI: 10.1186/s13362-020-00088-y [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13362-020-00088-y&link_type=DOI) 14. 14.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) 15. 15.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) 16. 16. Andreas Schuppert et al. Effectiveness of extended shutdown measures during the ‘Bundesnotbrems’ introduced in the third SARS-CoV-2 wave in Germany. December 2021. Infection 49(6):1–5. DOI: 10.1007/s15010-021-01713-7 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s15010-021-01713-7&link_type=DOI) 17. 17.Schüler L, Calabrese JM, Attinger S. Data driven high resolution modeling and spatial analyses of the COVID-19 pandemic in Germany. August 2021. PLoS ONE 16(8):e0254660. DOI: 10.1371/journal.pone.0254660 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0254660&link_type=DOI) 18. 18.Radev ST et al. OutbreakFlow: Model-based Bayesian inference of disease outbreak dynamics with invertible neural networks and its application to the COVID-19 pandemics in Germany. October 2021. PLoS Computational Biology 17(10):e1009472. DOI: 10.1371/journal.pcbi.1009472 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1009472&link_type=DOI) 19. 19. Berit Lange. Assessment of the infection dynamics and the disease burden of COVID-19 in Germany. March 2021.Public Health Forum 29(1):57–60. DOI: 10.1515/pubhef-2020-0122 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1515/pubhef-2020-0122&link_type=DOI) 20. 20.Tomori V. et al. Individual social contact data and population mobility data as early markers of SARS-CoV-2 transmission dynamics during the first wave in Germany—an analysis based on the COVIMOD study. October 2021.BMC Medicine 19(1). DOI: 10.1186/s12916-021-02139-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12916-021-02139-6&link_type=DOI) 21. 21.Nakajima E. et al. Estimating Parameters of a SEIRD Model Applied to SARS-CoV-2 Infections in Germany based on the Particle Swarm Optimization Method. July 2021. Applied Mathematics & Information Sciences 15(4):423–428. DOI: 10.18576/amis/150403 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18576/amis/150403&link_type=DOI) 22. 22.Maria Teresa Signes Pont MTS, Cortes JJ, Mora HM. An Epidemic Grid Model to Address the Spread of Covid-19: A Comparison Between Italy, Germany and France. February 2021. Mathematical and Computational Applications 26(1):14. DOI: 10.3390/mca26010014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/mca26010014&link_type=DOI) 23. 23.Kermack WO, McKendrick AG. A Contribution to the mathematical theory of epidemics. J Royal Stat Soc Ser A. 1927;115:700–21. 24. 24.Murray JD. Mathematical Biology I/II. New York: Springer; 2002. 25. 25.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) 26. 26.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) 27. 27.Coronavirus Disease (COVID-19) Situation Reports [Internet]. Who.int. 2020 [cited 2020 May 9]. Available from: [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports/](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports/) 28. 28.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) 29. 29.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) 30. 30.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) 31. 31.Nesteruk I. Final sizes and durations of new COVID-19 pandemic waves in Ukraine and around the world predicted by generalized SIR model. MedRxiv. November 2021. DOI: 10.1101/2021.11.22.21266683 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4xMS4yMi4yMTI2NjY4M3YxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMTUvMjAyMS4xMi4xNC4yMTI2Nzc3MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 32. 32.COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU). [https://github.com/owid/covid-19-data/tree/master/public/data](https://github.com/owid/covid-19-data/tree/master/public/data) 33. 33.Nesteruk I. The COVID-19 pandemic dynamic in India in the spring and summer of 2021. J Bio Med Open Access. 2021;2(2):123. 34. 34.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) [1]: /embed/graphic-3.gif [2]: /embed/graphic-4.gif [3]: /embed/graphic-5.gif [4]: /embed/graphic-6.gif