Differential and Time-Discrete SEIRS Models with Vaccination: Local Stability, Validation and Sensitivity Analysis Using Bulgarian COVID-19 Data

Bulgaria has the lowest COVID-19 vaccination rate in the European Union and the second-highest COVID-19 mortality rate in the world. That is why we think it is important better to understand the reason for this situation and to analyse the development of the disease over time. In this paper, an extended time-dependent SEIRS model SEIRS-VB is used to investigate the long-term behaviour of the COVID-19 epidemic. This model includes vaccination and vital dynamics. To apply the SEIRS-VB model some numerical simulation tools have been developed and for this reason a family of time-discrete variants are introduced. Suitable inverse problems for the identification of parameters in discrete models are solved. A methodology is proposed for selecting a discrete model from the constructed family, which has the closest parameter values to these in the differential SEIRS-VB model. To validate the studied models, Bulgarian COVID-19 data are used. To obtain all these results for the discrete models a mathematical analysis is carried out to illustrate some biological properties of the differential model SEIRS-VB, such as the non-negativity, boundedness, existence, and uniqueness. Using the next-generation method, the basic reproduction number associated with the model in the autonomous case is defined. The local stability of the disease-free equilibrium point is studied. Finally, a sensitivity analysis of the basic reproduction number is performed.