Forecast verifications for the real COVID-19 pandemic dynamics in Qatar ======================================================================= * Igor Nesteruk ## 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. Keywords * COVID-19 pandemic * vaccination efficiency * epidemic dynamics in Qatar * mathematical modeling of infection diseases * SIR model * parameter identification * statistical methods ## 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 *V**j*, number of vaccinated people *P**j* and number of vaccinations *Q**j* in Qatar from JHU, [5]. These values and corresponding moments of time *t**j* (measured in days) are shown in Table 1. For SIR-simulations, the values of *V**j* and *t**j* 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: ![Formula][1] View this table: [Table 1.](http://medrxiv.org/content/early/2021/06/22/2021.06.16.21259018/T1) Table 1. Cumulative numbers for Qatar from [5]. The laboratory-confirmed Covid-19 cases *V**j*, number of vaccinations *Q**j*, the number of vaccinated people *P**j*, and corresponding moments of time. The first and second derivatives can be estimated with the use of following formulas: ![Formula][2] ![Formula][3] ### 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 ![Formula][4] corresponding to the number of victims or the cumulative confirmed number of cases. Its derivative: ![Formula][5] 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 (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 *V**j* is incomplete and there is a constant coefficient *β**i* ≥ 1, relating the registered and real number of cases during the *i-th* epidemic wave: ![Formula][6] the number of unknown parameters increases by one. The values *V**j*, corresponding to the moments of time *t**j* 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. ![Figure.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.16.21259018/F1.medium.gif) [Figure.](http://medrxiv.org/content/early/2021/06/22/2021.06.16.21259018/F1) Figure. Real sizes of the third pandemic wave in Qatar at the calculated optimal value of the visibility coefficient *β*3 = 5.308. SIR simulations (red lines, [4]), numbers of vaccinations (yellow, Table 1) and vaccinated persons (green, Table 1). 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 *V**j* (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 *V**j* 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. * Received June 16, 2021. * Revision received June 16, 2021. * Accepted June 22, 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. Ahmed E. Fahmya, Mohammed M. El-desoukya, Ahmed S.A. Mohamed. Epidemic Analysis of COVID-19 in Egypt, Qatar and Saudi Arabia using the Generalized SEIR Model. MedRxiv. Posted August 22, 2020. Doi: [https://doi.org/10.1101/2020.08.19.20178129](https://doi.org/10.1101/2020.08.19.20178129) 2. 2. Ryad Ghanam, Edward L. Boone, Abdel-Salam G. Abdel-Salam. SEIRD MODEL FOR QATAR COVID-19 OUTBREAK: A CASE STUDY. Posted 26 May 2020. arXiv:2005.12777v1 3. 3.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, pp. 37–46. Doi: 10.20535/ibb.2021.5.1.228605 [http://ibb.kpi.ua/article/view/228605](http://ibb.kpi.ua/article/view/228605) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/ibb.2021.5.1.228605&link_type=DOI) 4. 4.Nesteruk I. Impact of vaccination and undetected cases on the COVID-19 pandemic dynamics in Qatar in 2021. [Preprint] MedRxiv, June 2021. Doi: 10.1101/2021.05.27.21257929 [https://medrxiv.org/cgi/content/short/2021.05.27.21257929v1](https://medrxiv.org/cgi/content/short/2021.05.27.21257929v1) [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wNS4yNy4yMTI1NzkyOXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMjIvMjAyMS4wNi4xNi4yMTI1OTAxOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. 5.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) 6. 6.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) 7. 7.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) 8. 8.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) 9. 9.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) 10. 10.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) 11. 11.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) 12. 12.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) 13. 13.Nesteruk I. The COVID-19 pandemic storm in India subsides, but the calm is still far away. [Preprint] medRxiv June, 2021. Doi: [https://doi.org/10.1101/2021.06.01.21258143](https://doi.org/10.1101/2021.06.01.21258143) [1]: /embed/graphic-1.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif