ABSTRACT
The COVID-19 pandemic dynamics in Qatar in the second half of May and the first half of June 2021 was compared with the published results of SIR-simulations based on the data from the period April 25 - May 8, 2021. Forecast verification showed very good agreement with the real number of cases (which can exceed the laboratory-confirmed one more than 5 times). The positive effect of mass vaccination became visible in June 2021.
Introduction
The COVID-19 pandemic dynamics in Qatar was simulated with the use of SEIR-model (susceptible-exposed-infected-removed) [1] and SEIRD-model (susceptible-exposed-infected-removed-dead) [2]. Some recent SIR-simulations [3, 4] were based on the datasets about the number of cases in December 2020 and April 25 – May 8, 2021 available in COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU), [5]. In particular, the incompleteness of the statistical information has been taken into account in [4] with the use of method proposed in [6]. It was shown that the actual accumulated number of COVID-19 cases may exceed the laboratory-confirmed number of cases registered in Qatar more than 5 times. Comparison with the corresponding SIR-curves has showed that the effect of mass vaccination was not evident during 4 months after its onset in December 2020, [4]. In this study, we will follow the epidemic in Qatar until mid-June 2021 in order to verify the predictions and to investigate the impact of vaccination.
Data and monitoring changes in epidemic dynamics
We will use tree data sets regarding the accumulated numbers of confirmed COVID-19 cases Vj, number of vaccinated people Pj and number of vaccinations Qj in Qatar from JHU, [5]. These values and corresponding moments of time tj (measured in days) are shown in Table 1. For SIR-simulations, the values of Vj and tj corresponding to the time period: April 25 – May 8, 2021 were used in [4]. The values presented in Table 1 and datasets available in [4] were used only for monitoring the changes in epidemic dynamics, which can be caused by vaccination, changes in quarantine restrictions, emergence of new strains of the pathogen, etc. To control these changes and to identify the epidemic waves, we can use daily or weekly numbers of new cases and their derivatives (see, e.g., [7-9]). Since these values are random, we need some smoothing (especially for daily amounts, which are also characterized by some weekly periodicity). For example, we can use the smoothed daily number of accumulated cases:
Cumulative numbers for Qatar from [5]. The laboratory-confirmed Covid-19 cases Vj, number of vaccinations Qj, the number of vaccinated people Pj, and corresponding moments of time.
The first and second derivatives can be estimated with the use of following formulas:
Generalized SIR model and parameter identification procedure
In [3, 4] the generalized SIR-model and the exact solution of the set of non-linear differential equations relating the number of susceptible S, infectious I and removed persons R were used in order to simulate different epidemic waves (see, e.g., [6, 10]). This solution uses the function
corresponding to the number of victims or the cumulative confirmed number of cases. 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 (4), the exact solution presented in [6, 10] depends on five parameters (αi is one of them, i is the number of the epidemic wave). The details of the optimization procedure for their identification can be found in [6, 11].
If we assume, that data set Vj is incomplete and there is a constant coefficient βi ≥ 1, relating the registered and real number of cases during the i-th epidemic wave:
the number of unknown parameters increases by one. The values Vj, corresponding to the moments of time tj from the period April 25 – May 8, 2021 have been used in [4] to find the optimal values of these parameters corresponding to the third epidemic wave in Qatar. In particular, the optimal value of the visibility coefficient β3 = 5.308 was calculated.
Results and Discussion
The results of SIR simulations performed in [4] are shown in the figure. Red lines represent the estimations of the real sizes of the third COVID-19 epidemic wave in Qatar for the optimal value β3 = 5.308 in relationship (6): V(t)=I(t)+R(t) – solid; the dashed line represents the number of infectious persons multiplied by 10, i.e. I (t) × 10; the dotted line shows the derivative (dV / dt)×100 calculated with the use of formula (5). “Circles” and “stars” correspond to the accumulated numbers of cases registered during the period of time taken for SIR simulations (April 25 – May 8, 2021) and beyond this time period, respectively.
Numbers of victims V(t)=I(t)+R(t) – the red solid line; numbers of infected and spreading I(t) multiplied by 10 – dashed; derivative dV/dt (eq. (5)) multiplied by 100 – dotted, [4]. Markers show accumulated numbers of laboratory-confirmed cases Vj (from Table 1 and [4]) and derivatives. Red “circles” correspond to the accumulated numbers of cases taken in [4] for calculations of the third wave; blue “stars” – number of cases beyond this period. “Crosses” show the values of the first derivative (eq. (2)) multiplied by 100β3. The blue solid line represents smoothed accumulated number of laboratory-confirmed cases (1) multiplied by β3.
To illustrate the accuracy of calculations, the real accumulated number of laboratory confirmed cases Vj form Table 1 and Table 1 in [4] was smoothed with the use of (1) and multiplied by the “visibility” coefficient β = 5.308, calculated in [4]. The result is shown in the figure by the solid blue line and is in very good agreement with the theoretical estimation (4) obtained in [4] and shown by the solid red line. We have calculated also the first derivative (2) with the use of eq. (1). The corresponding values have been multiplied by 100β and shown by blue “crosses”. These values are in rather good agreement with the theoretical estimation of the real daily number of new cases (eq. (5), the dotted red line, [4]).
The epidemic dynamics in June 2021 shows that the averaged number of new cases (eqs. (1) and (2), blue “crosses”) started to deviate from the theoretical estimation (eq. (5), the dotted red line) obtained in [4] with the use of data set corresponding the time period April 25 – May 8, 2021. Probably it is a result of mass vaccination, since the number of vaccinated people exceeded the half of Qatar population in June 2021 (see green markers). It would be very interesting to follow the further dynamics of the pandemic in Qatar in order to assess the long-term impact of vaccination
Unfortunately, mathematical modeling cannot answer the question about possibility new epidemic waves after June 2021, because even strict quarantine measures does not prevent the emergence of new strains of coronavirus. The number of infectious persons in Qatar is still high. According to the estimation (shown by the red dashed line), the number of infectious persons was around 6,000 in mid-June 2021. This fact and the unfavorable course of the pandemic in India [12, 13] make the probability of the emergence of new coronavirus strains and new epidemic waves rather high. In the case of an optimistic scenario (without the emergence of new waves), the end of the epidemic in Qatar should be expected in March 2022 (see [4]). The actual number of cases may increase during this rather long period (from June 13, 2021) by 70,000 (the number of new registered cases will be about five times less). The government should closely monitor the average daily number of cases and its derivative (which can be calculated by formulas (1)-(3)) and immediately intensify quarantine measures in case of their increase.
Data Availability
data is in the text
Acknowledgements
The author is grateful to Oleksii Rodionov for his help in collecting and processing data.
