Waves of COVID-19 pandemic. Detection and SIR simulations

Data

The official information regarding the accumulated numbers of confirmed COVID-19 cases V_j in Ukraine and in the world from WHO daily situation reports (numbers 101-194), [1] is presented in Table 1. The corresponding moments of time t_j (measured in days, zero point corresponds to January 20, 2020) are also shown in this table. Some time periods were used for calculations. Other values were used only for verifications of their results.

Epidemic waves detection

Changes in quarantine conditions, human behavior, pathogen activity, weather etc. can cause changes in the epidemic dynamics, namely the so-called epidemic waves. The simplest way of their detection is to find some changes in the dependences of the number of registered cases on time. Since the number of cases is random, its time dependence needs some smoothing. We can use the method proposed in [16]: and the derivatives of the smoothed values in order to detect the changes in epidemics dynamics.

SIR model

The SIR model for an infectious disease [2-5] relates the number of susceptible persons S (persons who are sensitive to the pathogen and not protected); the number of infected is I (persons who are sick and spread the infection; please don’t confuse with the number of still ill persons, so known active cases) and the number of removed R (persons who no longer spread the infection; this number is the sum of isolated, recovered, dead, and infected people who left the region):

Parameters α_i and ρ_i are supposed to be constant for every wave of the epidemics, i.e. for the time periods t_i ≤ t ≤ t_i+1, i = 1, 2,3,….

Parameters α_i show how quick the susceptible persons become infected (see (4)). Large values of this parameter correspond to severe epidemics with many victims. These parameters accumulate many characteristics. First they shows how strong (virulent) is the pathogen and what is the way of its spreading. Parameters α_i accumulates also the frequency of contacts and the way of contacting. In order to decrease the values of α_i, we have to minimize the number of our contacts and change our contacting habits. For example, we have to avoid the public places and use masks there, minimize or cancel traveling. We have to change our contact habits: to avoid handshakes and kisses. First, all these simple things are very useful to protect yourself. In addition, if most people follow these recommendations, we have chance to diminish the values of parameter α_i and reduce the negative effects of the pandemic.

The parameters ρ_i characterize the patient removal rates, since eq. (6) demonstrates the increase rate for R. The inverse values 1/ ρ_i are the estimations for time of spreading the infection τ. So, we are interested in increasing the values of parameters ρ_i and decreasing 1/ ρ_i. People and public authorities should work on this and organize immediate isolation of suspicious cases.

Since the derivative d (S + I + R) / dt is equal to zero (it follows from summarizing Eqs. (4)-(6)), the sum must be constant for every wave and is not the volume of population (see also [15]).

Analytical solution of SIR equations

To determine the initial conditions for the set of equations (4)–(6), let us suppose that at the beginning of every epidemic wave t_i :

In particular, when the first wave of the epidemic starts with one infected person, the initial conditions (8) can be written as follows:

Equations (9) were used in [5-16] to simulate the first waves of epidemics in different countries.

It follows from (4) and (5) that

Integration of (10) with the initial conditions (8) yields:

It follows from (11) that function I has a maximum at S =ν_i and tends to zero at infinity. The corresponding number of susceptible persons at infinity S_i∞ > 0 can be calculated from a non-linear equation

Formula (12) follows from (11) at I=0.

As in [5] we solve (4)-(6) by introducing the function V (t) = I (t) + R(t), corresponding to the number of victims or cumulative confirmed number of cases. It follows from (5)-(7) and (11) that:

Integration of (13) yield an analytical solution for the set of Equations (4)–(6):

Thus, for every set of parameters N_i, I_i, R_i, ν_i, α_i, t_i and a fixed value of V, integral (15) can be calculated and a corresponding moment of time can be determined from (14). Then functions I(t) and R(t) can be easily calculated with the of formulas (11) and:

The final number of victims (final accumulated number of cases) can be calculated from:

To estimate the duration of an epidemic outbreak, we can use the condition: which means that at t > t _final less than one person still spreads the infection.

Parameter identification procedure

In the case of a new epidemic, the values of its independent six parameters are unknown and must be identified with the use of limited data sets. For the first wave of an epidemic starting with one infected person, the number of unknown parameters is only four, since I₁ = 1 and R₁ = 0. The corresponding statistical approach was used in [6-16] to estimate the values of four unknown parameters for the first waves of Covid-19 pandemic in different regions.

For the next epidemic waves (i > 1), the registered number of victims V_j corresponding to the moments of time t_j can be used in eq. (15) in order to calculate F_{i, j} = F_i (V_j, N_i, ν_i, I_i, R_i) for every fixed values of N_i, ν _i, I_i, R_i and then to check how the registered points fit the straight line (14). For this purpose the linear regression can be applied, e.g., [17], and the optimal straight line, minimizing the sum of squared distances between registered and theoretical points, can be defined. Thus we can find the optimal valuesα_i, t_i and calculate the correlation coefficient r_i for the of linear dependence (14).

Then the F-test may be applied to check how the null hypothesis that says that the proposed linear relationship (14) fits the data set. The experimental values of the Fisher function can be calculated with the use of the formula: where n_i is the number of observations for the i-th wave of epidemis, m = 2 is the number of parameters in the regression equation, [17]. The corresponding experimental value F_i has to be compared with the critical value F_C (k₁, k₂) of the Fisher function at a desired significance or confidence level (k₁ = m −1, k₂ = n_i − m), [18]. When the values n_i and m are fixed, the maximum of the Fisher function coincides with the maximum of the correlation coefficient. Therefore, to find the optimal values of parameters N_i, ν_i, I_i, R_i, we have to find the maximum of the correlation coefficient for the linear dependence (14). To compare the reliability of different predictions (with different values of n_i) it is useful to use the ratio F_i / F_C (1, n_i − 2) at fixed significance level. We will use the level 0.001; corresponding values of F_C (1, n_i − 2) can be taken from [18]. The most reliable prediction yields the highest F_i / F_C (1, n_i − 2) ratio.

In the case of sequential calculation of epidemic waves i = 1,2,3 …, it is possible to avoid determining the four optimal unknown parameters N_i, ν_i, I_i, R_i, thereby reducing the amount of calculations. For parameters I_i, R_i it is possible to use the numbers of I and R calculated for the previous wave of epidemic at the moment of time when the following wave began. In this study we will use this approximate approach.

Detection of Covid-19 pandemic waves in Ukraine and in the world

Applications of formulae (1)-(3) for the pandemic dynamics in Ukraine and in the world are shown in Fig. 1 and 2. The accumulated numbers of cases (blue “circles”) were smoothed with the use of eq. (1) and shown by blue lines. Red and black markers represent the results of differentiation (2) and (3) respectively. The make the results more visible, the first derivative (2) is multiplied by 100, the second one (formula (3)) – by 1000.

• Download figure
• Open in new tab

Fig 1. Epidemic dynamics in Ukraine versus time in days.

Accumulated number of cases (blue markers and line, eq. (1)). Red markers represent first derivative (eq. (2)), black one show the second derivative (eq. (3)).

• Download figure
• Open in new tab

Fig 2. Pandemic dynamics in the world versus time in days.

Accumulated number of cases (blue markers and line, eq. (1)). Red markers show first derivative (eq. (2)) multiplied by 100, black one - the second derivative (eq. (3)) multiplied by 1000.

Fig. 1 and 2 demonstrate that the second derivatives increase after epidemic outbreak, then become smaller and negative. Such behavior is typical for the first wave of the epidemics (before May 17, 2020 in Ukraine, see Fig. 1 and before May 12 in the world, see Fig. 2). The jumps in the values of the second derivative indicate changes in the conditions of the pandemic (for example, the weakening of quarantine) and the transitions to the next waves with other values of the parameters of mathematical models. These jumps occurred on May 16, May 29, June 8, July 3 and July 19-20 in Ukraine (see Fig. 1) and on May 12 and 28; June 13; July 5 and 14 in the world (see Fig. 2). Therefore these days can be treated as the beginning of the second, third etc. waves.

Thus four waves were calculated for Ukraine and presented in the next Section. The distance between some jumps of the second derivative, which took place in the second half of May, was too small to make statistical estimates of the parameters, so the individual waves of the epidemic during this period were not isolated. The same situation occurred for the world dynamics in July. The distance between jumps of second derivative was too short in order to isolate two different waves. Thus five waves were calculated for the world pandemic dynamics. These facts can reduce the accuracy of SIR simulations. At the end of June 2020, the number of days of observation of the fifth wave in Ukraine (after July 19) was still insufficient for statistical analysis, so this wave can be calculated later.

The second epidemic wave in Ukraine was caused by the cancelation of the national lockdown on May 10. After the incubation period (approximately after May 16), the number of new cases began to grow faster. Further waves of the epidemic in Ukraine are associated with further easing of quarantine and mass non-compliance with social distancing. The fifth epidemic wave in Ukraine can be explained by the consequences of the holiday season, which increased the number of trips and violations of social distancing. Analyze the causes of new waves in the world is difficult because of the large number of countries with different pandemic dynamics. Probably the fourth pandemic wave after June 13 is associated with mass protests in the United States.

Hidden periods of Covid-19 pandemic

It was already mentioned that during the initial stages of epidemic the registered number of cases is much lower than the real one, since there are a lot of infected persons without symptoms. In every country some hidden periods occurred, before the first case was confirmed [14, 15]. To make reliable estimations of hidden periods, we have to use the datasets corresponding to the period when the number of registered cases approaches the real one during the first wave of epidemic. In prediction 21 for Ukraine (see Table 2) and in prediction 2 for the world (see Table 4), the periods of time just before the start of the second wave were used for SIR calculations. The results are shown in Fig. 5.

• Download figure
• Open in new tab

Fig 5. Hidden periods in Ukraine and in the world. SIR curves for the first wave of pandemic (lines) and accumulated number of cases (markers) versus time.

Blue and magenta colors correspond to Ukraine (Table 1) and the world (Table 2) respectively. Numbers of infected and spreading I (dashed lines) and victims V=I+R (solid lines). “Stars” represent the accumulated number of cases V_j from [1]; “squares” show the number of laboratory confirmed cases in Wuhan, China calculated in [9] with the use of daily data presented in [19].

In particular the hidden period in the world according to the second simulation of the first wave is approximately 60 days longer in comparison this the first prediction for the first wave (see Table 4, [15]). Ukrainian statistics yields similar results. Estimation No. 8 available in [13-15] yielded February 16 as the day of the epidemic outbreak in Ukraine. The results of calculations, performed with the use of the more recent dataset (days immediately preceding the beginning of the second wave) and shown in Table 2, demonstrate that the first cases in Ukraine probably appeared in the beginning of January, 2020 (see solid blue line in Fig. 5).

Solid magenta line in Fig. 5 demonstrates that Covid-19 probably started to spread in August 2019. On December 31, 2019 – the day when China notified WHO about the situation in Wuhan - more than 2,000 persons could spread the infection (see dashed magenta line in Fig. 5). In the period from 18 to 27 October, 2019 (when the Military World Games held in Wuhan with the participation of 9,300 athletes from more than 140 countries) the number of infected persons can be estimated by 50-80 (see dashed magenta line in Fig. 5). May be some participants got the infection and passed it on to their families, [20].

Long-term predictions for Ukraine and the world

Table 6 illustrates a rather high accuracy of calculations. It can be seen that after one week of observations (after the last day of the periods taken for calculations) the real accumulated number of cases exceed the calculated ones by 0.54% and 1.35% for Ukraine and the world respectively. After two weeks of observations the discrepancies are higher (3.31% for Ukraine and 4.25% for the world).

Fig. 6 shows the long-term predictions for Ukraine and the world for the period August-December, 2020. Despite the relatively high accuracy of pandemic wave calculations (see Table 6), these long-term forecasts given should be considered as preliminary and overly optimistic. The new waves are expected due to the quarantine easing and changes in social distancing. Gradual weakening of quarantines may not cause pronounced new waves. In this case, the parameters of the pandemic dynamics will change continuously, which makes the proposed approach unsuitable. The expected emergence of vaccines could also change the course of the pandemic. If mass vaccination begins in 2021, the total number of cases will still exceed 20 million. In any case the global and Ukrainian dynamics must be updated with the use of new data sets.

• Download figure
• Open in new tab

Fig 6. Long-term predictions for Ukraine (fourth wave) and the world (fifth wave). SIR curves (lines) and accumulated number of cases (markers) versus time.

Magenta and brown colors correspond to Ukraine (Table 3) and the world (Table 5) respectively. Numbers of infected and spreading I (dashed lines) and victims (accumulated number of conformed cases) V=I+R (solid lines).

Very long duration of the pandemic requires correction of our behavior, we can not live as before it occurred. Decreased feelings of insecurity and non-compliance with social distancing may further increase the pandemic duration and the number of the coronavirus victims. Total closure of settlements or regions can be recommended only in the event of a sharp increase in the number of cases. There are many things that can be done without loss to the economy and our daily lives:

1. Minimize the number of contacts and trips, not visit crowded places. Work and study remotely where possible.

2. Refrain from shaking hands and kisses during meetings. Use masks in transport and crowded areas.

3. If you (or others) have any suspicious symptoms, do your best to avoid the spread of the infection.

It must be noted that the relatively small number of people who spread the coronavirus can cause a fairly long course of the epidemic. In particular, 20 such persons in Ukraine (estimation at the beginning of 2021, see Fig. 6) can cause new cases for another two and a half months. This explains the many new outbreaks in South Korea, Japan, China, Singapore and not all of them are caused by imported cases, [1, 21, 22]. To finally overcome the pandemic, it is very important to identify all the infected, and this is very difficult due to the large number of asymptomatic patients.

The dynamics of Covid-19 pandemic in different countries was simulated with the use of different mathematical models [23-37]. Nevertheless, in literature there are no long-time predictions for the final sizes and durations of the pandemic in Ukraine and in the world in order to compare with the presented results.

Probability of meeting an infected person

As long as there are infected people in the population, there is a chance to meet one of them. The corresponding probability can be estimated with the use of simple formula: where N_pop is the number of persons in a country or in a region (volume of population). With the use of corresponding SIR curves it is possible to estimate p(t) for every moment of time. For example, on July 19, 2020 these probabilities can be estimated as 1.04·10⁻⁴ and 9.28·10⁻⁵ for Ukraine and the world respectively. It means that Ukraine was 12% more dangerous in comparison with the average situation in the world.

Formula (20) can be useful for deciding which countries’ citizens are welcome guests after the resumption of international passenger traffic. If p_A (t) > p_B (t) for countries A and B, citizens of A are not welcome guests in country B at the moment of time t. Since such decisions are very important, we will propose some simple method of estimating p(t) without calculations of SIR parameters for epidemic waves in a country. According to (13):

Since function S(t) changes not very sufficient during a wave (compare initial S_i =N_i –I_i –R_i and final S_i∞ values of this function in Tables 2-5 and in corresponding Tables in [11, 15]), we can use the constant value S_i instead of S(t). To estimate dV/dt, it is possible to use the known values V_j or formula (2) for smoothed number of cases. Finally, with the use of (21), equation (20) can be written as follows:

The results of calculations for different countries and regions are presented in Table 7. It can be seen that coefficient k_p changes not very fast (compare values on different days for Ukraine and the world). It means that values of this coefficient can be used during rather long period of time before new exact estimations of I(t) will be available. According to eq. (23), the probability of meeting an infected person mostly depends on the daily increase in the number of cases. The values p were estimated for July 19, 2020 and presented in Table 7.

The p value for the fourth wave in Ukraine in 20% lower than the exact estimation 1.04·10⁻⁴ calculated before with the use of eq. (20). In the case of the world fifth wave, the difference between the exact estimation (20) and the approximate formula (23) is only 8%. It can be seen that the most infected country is USA. The average global figure is approximately 4 time less. The probabilities p in Kyiv, Moldova and Ukraine are around the world average value. Lower figures were calculated for Spain, Italy and UK. The safest counties are France, Austria and Germany.

Estimations of reproduction numbers

Effective reproduction number R_t(t) shows the average number of people infected by one person [38]. In terms of SIR model it can be estimated as:

Where τ is the average time of spreading the infection. Taking into account (5), τ ≈ 1/ ρ_I constant value S_i =N_i –I_i –R_i instead of S(t), equation (24) can be rewritten as follows: and using

Equation (25) was used for calculations of the effective reproduction number for different regions and moments of time. The results are shown in the last column of Table 7.

It can be seen that France, Italy, Germany and Austria have the lowest values of R_t which are close to the world figure (wave 5). Not much higher values were calculated for USA and Ukraine (wave 4). Spain, UK, Moldova and Kyiv demonstrate the highest figures. The large difference in Ukrainian and Kyiv values can be explained by very different days of estimation (July 3 and April 11 respectively). Earlier periods of pandemic have higher values of R_t. The same tendency demonstrates Table 7.

The R_t values corresponding to July 19, 2020, were estimated for Ukraine (wave 4, Table 3) and the world (wave 5, Table 5) with the use of exact dependences I(t) in formula (24) in order to compare with the results presented in Table 7. The values R_t = 0.0257 and R_t = 0.0059 for Ukraine and the world respectively demonstrate that the effective reproduction number can decrease very rapid over time.

The effective reproduction numbers of Covid-19 pandemic were estimated with the use of different mathematical models for different countries [37, 39, 40]. In particular, in [37] the R_t values were estimated as of May 10 for EU countries. Corresponding values for Austria, Spain, Germany, France and Italy vary from 0.45 to 0.74 and significantly exceed those shown in Table 7. Rapid changes in R_t values for the epidemic in China are reported in [39] (from 2.35 on January 16 to 1.05 on January 31).

The basic reproduction number R₀ for the initial periods of an epidemic outbreak can be estimated with the use of equation, [38]: where T_d is the doubling time (the period during which the number of cases duplicates; the exponential growth is assumed). In the first half of March 2020, the increase in cases in the United States and around the world was close to exponential, [10]. Taking into account the values T_d = 2.31 and T_d =3.65 days for USA and the world respectively (calculated in [10]) and estimation τ ≈ 1/ ρ₁, the values R₀ =1.55 (USA) and R₀ =1.99 (the world) can be obtained (we have used eq. (26) and values ρ₁ = 0.6877, [15]; ρ₁ = 0.275, Table 4).

The basic reproduction numbers of Covid-19 pandemic were estimated in [29, 41-44]. According to WHO calculations for China [41], values R₀ vary from 2 to 2.5. During the initial period of an epidemic outbreak many infected persons are undetected. This fact diminish the accuracy of R₀ estimations. For example, paper [44] reported the R₀ =7.97 for Indonesia and two different values 1.54 and 6.22 for Singapore.