Modeling the dynamics of COVID-19 in the presence of Delta and Omicron variants with vaccination and non-pharmaceutical interventions

Since its inception in December 2019, many safe and effective vaccines have been invented and approved for use against COVID-19 along with various non-pharmaceutical interventions. But the emergence of numerous SARS-CoV-2 variants has put the effectiveness of these vaccines, and other intervention measures under threat. So it is important to understand the dynamics of COVID-19 in the presence of its variants of concern (VOC) in controlling the spread of the disease. To address these situations and to find a way out of this problem, a new mathematical model consisting of a system of non-linear differential equations considering the original COVID-19 strain with its two variants of concern (Delta and Omicron) has been proposed and formulated in this paper. We then analyzed the proposed model to study the transmission dynamics of this multi-strain model and to investigate the consequences of the emergence of multiple new SARS-CoV-2 variants which are more transmissible than the previous ones. The control reproduction number, an important threshold parameter, is then calculated using the next-generation matrix method. Further, we presented the discussion about the stability of the model equilibrium. It is shown that the disease-free equilibrium (DFE) of the model is locally asymptotic stable when the control reproduction is less than unity. It is also shown that the model has a unique endemic equilibrium (EEP) which is locally asymptotic stable when the control reproduction number is greater than unity. Using the Center Manifold theory it is shown that the model also exhibits the backward bifurcation phenomenon when the control reproduction number is less than unity. Again without considering the re-infection of the recovered individuals, it is proved that the disease-free equilibrium is globally asymptotically stable when the reproduction threshold is less than unity. Finally, numerical simulations are performed to verify the analytic results and to show the impact of multiple new SARS-CoV-2 variants in the population which are more contagious than the previous variants. Global uncertainty and sensitivity analysis has been done to identify which parameters have a greater impact on disease dynamics and control disease transmission. Numerical simulation suggests that the emergence of new variants of concern increases COVID-19 infection and related deaths. It also reveals that a combination of non-pharmaceutical interventions with vaccination programs of new more effective vaccines should be continued to control the disease outbreak. This study also suggests that more doses of vaccine should provide to combat new and deadly variants like Delta and Omicron.