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.
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.
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 Vj and moments of time tj (measured in days) are shown in Table 1. The values for the period Tc11: 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 Tc12: September 29 – October 12, 2021 to simulate the 12-th wave. Other Vj and tj values will be used to control the accuracy of predictions.
To estimate the mortality rate (ratio of accumulated number of cases Vj to the accumulated number of deaths Dj), let us take figures Dj 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 Vj values from Table 1, we can calculate the mortality rates mi =Dj *1000/Vj (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 corresponding to the number of victims or the cumulative confirmed number of cases [8, 20]. Its derivative: yields the estimation of the average daily number of new cases. When the registered number of victims Vj 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: and its numerical derivative: 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].
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.