Abstract
We developed an endemic model of COVID-19 to assess the impact of vaccination and immunity waning on the dynamics of the disease. Our model exhibits the phenomenon of back-ward bifurcation and bi-stability, where a stable disease-free equilibrium co-exists with a stable endemic equilibrium. The epidemiological implication of this is that the control reproduction number being less than unity is no longer sufficient to guarantee disease eradication. We showed that this phenomenon could be eliminated by either increasing the vaccine efficacy or by reducing the disease transmission rate (adhering to non-pharmaceutical interventions). Furthermore, we numerically investigated the impacts of vaccination and waning of both vaccine-induced immunity and post-recovery immunity on the disease dynamics. Our simulation results show that the waning of vaccine-induced immunity has more effect on the disease dynamics relative to post-recovery immunity waning, and suggests that more emphasis should be on reducing the waning of vaccine-induced immunity to eradicate COVID-19.
1 Introduction
The Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) that infected millions of people around the world in 2002 and the Middle East Respiratory Syndrome Coronavirus (MERS-CoV) in 2012 are responsible for many cases of the common cold (Flu) in recent years. In mid-December 2019, twenty individuals suffering from pneumonia and seven other critically ill patients were admitted to the Wuhan Municipal Health Commission (WHC). The majority of them had prolonged exposure to animals such as snakes, poultry, bats etc., at the Huanan Seafood Market in Wuhan [7]. After extensive research by the World Health Organization (WHO) [37], these individuals were believed to be infected by a strange virus that was later named the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) by the WHO on February 11, 2020 [63]. The disease caused by the virus was also named COVID-19 on the same day [63]. As at August 29, 2021, there has been over 216,789,388 confirmed cases of COVID-19 globally, with 4,508,379 confirmed deaths and over 193,721,846 recovered individuals [65].
After infection, the corona-virus attacks the body cells and replicate quickly resulting in a severe cases of COVID-19 [48, 55]. According to the medical practitioners, COVID-19 victims experience mild or moderate respiratory problems and difficulty in breathing. Other symptoms include but not limited to tiredness, flu, cough, fever and shortness of breath [36, 64]. Infected individuals go through the incubation period of 2-14 days for developing symptoms, while some do not show symptoms at all during their infection (asymptomatic individuals) [34, 35]. It is believed that asymptomatic individuals account for more than 80% of the community transmissions [68]. The incubation period is generally referred to as asymptomatic stage since the infected person can still successfully transmit the disease around the community through effective contact [34, 35, 68]. The individuals in this category are generally encouraged by WHO to selfisolate for a minimum of 14 days and be in constant contact with the COVID-19 regulatory line in their community to forestall further transmission. People with severe illness, difficulty in breathing and those with underlying medical conditions are encouraged to visit the hospital with immediate alacrity [65]. According to the virologists, corona-virus is a single-stranded RNA virus while HIV is a retrovirus (ie. it has the ability to carry single-stranded RNA as its genetic material instead of double-stranded DNA carried by human cells) [16]. The implication of this is that the treatment used in treating diseases with RNA virus can also be used for COVID-19 treatment [17].
In an effort to curb the transmission of COVID-19, many countries and regions around the world implemented the non-pharmaceutical interventions (NPIs) such as physical distancing, wearing of face mask, quarantine, closure of schools/businesses, travel restriction, among others [12, 14, 22, 44, 45, 57]. The implementation of some of these NPIs have both physical and mental implications on the society. Some of these include loss of jobs, loneliness, lack of access to adequate medical treatment, economic meltdown, etc. [6, 15, 25, 41, 45, 50, 54]. The availability of COVID-19 vaccines has represented a unique opportunity in the fight against COVID-19. In addition to reducing disease transmission and burden, widespread vaccination will also allow countries to safely lift NPIs and promote the gradual recovery from the pandemic. Vaccination has been reported as one of the best development in the public health sector since around 1900 [23, 30, 31, 32]. It has helped in the eradication of several infectious diseases such as polio, measles, rubella, smallpox etc [19]. Despite the public health benefits derived from vaccination, there are still a lot of concerns about the need, efficacy and safety of many available vaccines [23, 30, 31, 32].
South Africa is one of the first African countries to receive the COVID-19 vaccine. They received their first shipment which contains a million doses of the AstraZeneca/Oxford COVID-19 vaccine, on February 1, 2020 [60, 62]. A few days later, the vaccination process was put on hold over concerns on the efficacy of the vaccine against the B.1.351 variant of SARS-CoV-2, which was the dominant strain of the virus in the country [40, 42]. The Johnson and Johnson COVID-19 vaccine, with an initial 80,000 doses [43], was rolled out on the February 17, 2021. This vaccine is observed to be more effective against the B.1.351 strain [40, 42].
Several mathematical models have been developed and analyzed to study the transmission dynamics of COVID-19 [2, 4, 26, 27, 33, 53, 39, 46, 66, 70]. Bugalia et al. [4] developed a COVID-19 model using an endemic framework considering the roles of intervention strategies such as lockdown, hospitalization and quarantine. They computed the reproduction number and calculated both the disease-free and endemic equilibriums of their model. They also established the epidemiological significance of the equilibrium stability, transcritical bifurcation and boundedness of solution. In order to validate their system of differential equations representing the COVID-19 dynamics, they fitted the cumulative and new daily cases in India after which they estimated the model parameters and predicted the near future scenario of the disease. The global sensitivity analysis was also established to observe the impact of different parameters of Ro. They also investigated the dynamics of disease in respect to partial lock-down, complete lock-down and no lock-down. Their analyses conclude that if there is partial or no lockdown case, then endemic level would be unbelievably high and that India may experience more that six million infections if care is not taken.
The effect of an imperfect vaccine in reducing COVID-19 transmission and eventually eradicating the disease in the United State of America was examined in [24]. Their mathematical model provides an estimate for the minimum number of susceptible individuals needed to be vaccinated to obtain vaccine-induced herd immunity. In addition, they investigated the epidemiological consequence of the herd immunity threshold and showed that the disease can be eliminated if the herd immunity threshold is achieved. They simulated their model using COVID-19 parameters related to states like Florida and New York, and estimated the herd immunity threshold for these states. The impact of government actions and individual reactions on the corona-virus dynamics in South Africa was assessed using a mathematical model by [47]. Their model examines the effects of individual behavioral reaction and non-pharmaceutical interventions such as quarantine and lockdown on the dynamics of COVID-19 in South Africa. Based on the South African officially published COVID-19 data for the period of March 2020 to early May 2020. They discovered that a detection rate of at least 0.5 per day may lead to a significant reduction in the total number of infected cases. They also observed that in the absence of intervention strategies, the peak number of cases could be attained earlier than expected which can also be delayed if the interventions are implemented. The importance of having a COVID-19 vaccine that is effect in older individuals was emphasized in [53].
Johnson and Bruce [29] proposed a framework for modeling fear of infection and frustration with social distancing during COVID-19 epidemic. They established that the SEIR behavior-perception model has three principal modes of qualitative behavior: no outbreak, controlled outbreak, and uncontrolled outbreak. They fitted their model to the cumulative cases of COVID-19 and mortality data for different regions, and showed that their model can produce and sustain waves of infections. Their results show that the region with significant decline after wave of infection may likely survive the second wave, while those with moderate success in controlling their initial COVID-19 outbreak are most likely to experience substantial second waves or devastating outbreaks.
In this work, we develop and analyze a novel endemic model of COVID-19 that incorporates vaccination as a pharmaceutical intervention strategy (PIS). This model considers the waning of both vaccine-induced immunity and post-recovery immunity. We derive the control reproduction number for our model and study the disease-free and endemic dynamics of the model with a focus on the South African COVID-19 scenario. Our goal is to use this model to study the effect of immunity waning on both the control reproduction number and the disease dynamics.
The rest of the paper is structured as follows. In Section 2, we present the model formulation and state the assumptions on the model. Section 3 entails the model analyses. In the first part, we discuss the non-negativity of solution and the invariant region, while in the second part we establish the bifurcation analysis and discuss its biological interpretations. In Section 4, we present numerical computations of the control reproduction number derived in Section 3, in terms of the vaccination and immunity waning parameters. Also in this section, we study the effects of vaccination and immunity waning of the dynamics of the disease. A brief discussion and future work conclude the paper in Section 5.
2 Mathematical model
We develop a six compartment endemic model of COVID-19 to study the effect of immunity waning on the dynamics of the disease. Our model incorporates vaccination as a pharmaceutical intervention strategy and assumes that the immunity provided by the vaccine (vaccine-induced immunity) and immunity due to infection (post-recovery immunity) wane over time [18, 8, 56, 21]. The implementation of vaccination in our model is such that it reduces the chances of susceptible individuals acquiring the infection. The susceptible compartment is divided into two: unvaccinated susceptible population (S) and vaccinated susceptible population (Sv). We also divided the infected population into three: exposed (not infectious) (E1), pre-symptomatic (infectious without symptoms) (E2) and symptomatic infectious (I). The recovered population is denoted by R. We assume that recovered individuals do not have permanent immunity to the disease, they become susceptible again at some rate (post-recovery immunity waning). The vaccine and post-recovery immunity waning assumptions in the model can also be interpreted as vaccinated and recovered individuals becoming susceptible to another variant of the virus. Table 1 gives the descriptions of the variables in our model.
Figure 1 shows a schematic diagram of the model, where the black and green arrows indicate the direction of the flow of individuals between the compartments at the rates indicated beside the arrows. The green arrows specifically show the flow of individuals due to immunity waning. The red arrows show natural death and death due to COVID-19, while the blue arrow shows birth and immigration into the susceptible population. Our model assumes that immigration into the population only happens in the susceptible population and individuals do not leave the population through emigration. The transmission rates for the pre-symptomatic and symptomatic individuals are respectively denoted by β1 and β2. Using these two rates, we define our force of infection, λ as where N (t) = S(t) + Sv(t) + E1(t) + E2(t) + I(t) + R(t) denotes the total population at any time t. The differential equations of our model are given by under the following initial conditions where σ is the vaccination rate and 0 ≤ f ≤ 1 is the vaccine efficacy. Here, f = 0 implies that the vaccine does not provide any protection from acquiring the disease, while f = 1 implies that it is 100% effective in preventing new infection. The vaccine and post-recovery immunity waning rates are ϕ1 and ϕ2, respectively. The natural death rate is denoted by µ, while the rates of transitioning from E1 to E2 and E2 to I are denoted by α1 and α2, respectively. Symptomatic infectious individuals recover from the disease at the rate ψ and die as a result of COVID-19 related complications at the rate δ.
In Table 2, we present the model parameters, their descriptions and values. The birth and immigration rate Λ, and the death rates µ and δ are related to South African demography. These parameters were obtained from the official Republic of South Africa websites. The vaccine and post-recovery immunity waning parameters are varied throughout this paper. The exact values used in each scenario will be specified in the figure caption.
3 Model analysis
In this section, we analyze our ODE model (2.2)-(2.7) and show the non-negativity of its solution and that of the invariant region that describes the solution. We also analyze the linear stability of disease-free equilibrium and the endemic equilibrium.
3.1 Non-negativity of Solution and the Invariant Region
Here, we establish the non-negativity of the solutions of the ODE model (2.2)-(2.7) as well as the invariant region that describes these solutions.
Let ϒ represents the closed set defined as then equation (3.1) is positively invariant and attracting with regards to (2.2)-(2.7).
Proof. Using (2.8) as the initial conditions of the ODE system (2.2)-(2.7), from (2.2), we have
Observe that equation (3.2) is a first-order linear differential equation. The integrating factor is chosen as where λ(η) ≡ λ(S, Sv, E1, E2, I, R) so that the solution to (3.2) with equality is
Evaluating the equation above and using the fact that S(0) ≥ 0, we have which implies that
Since S(0) is a non-negative constant (see (2.8)), (3.3) ensures that the solution S(t) is non-negative.
Following the same approach, the non-negativity of Sv, E1, E2, I and R can be established provided that the initial conditions exist for all time t > 0.
More so, summing up equations (2.2)-(2.7) gives whose solution is given by where exp(µt) is the integrating factor. Taking the limit of both sides as t → ∞ gives
This means that the expression in (3.5) is bounded above by in the domain described in Lemma 3.1. Hence, the ODE system (2.2)-(2.7) is well-posed epidemically and meaningful mathematically since all the state variables are non-negative for all t > 0. Thus, we proceed to analyze the model in ϒ. This completes the proof. □
3.2 Local Stability of Disease-free equilibrium (DFE)
The disease-free equilibrium of (2.2)-(2.7) is given by where . We calculate the control reproduction number using the next generation matrix operator [13, 59] on (2.2)-(2.7) as used by [46, 51]. The matrix for the rates of new infections ℱ and that of transfer of infected individuals, 𝒱 at the DFE are given by
The next generation matrix is defined by ℱ𝒱-1. Computing the spectral radius (dominant eigenvalue) of the next generation matrix, we obtain the control reproduction number
The first term in the control reproduction number (3.7) accounts for the infections caused by the pre-symptomatic individuals in E2, while the second term account for those caused by symptomatic individuals in (I). Using (3.7), we establish the following result.
The DFE of the COVID-19 model in (2.2)-(2.7) is locally asymptotically stable (LAS) if ℛcon < 1 and unstable if ℛcon > 1.
The proof of this theorem is standard and can be established using theorem 2 of [59]. Recall that the parameter ℛcon is used to measure the average number of new COVID-19 cases generated by a single infected individual introduced into a susceptible population where a certain fraction of individuals are vaccinated. The biological interpretation of Theorem 3.1 is that if the control reproduction number ℛcon < 1, the influx of a new COVID-19 case in the community will not overwhelm the community but the whole community will be endemic if ℛcon > 1.
In Figure 2, we present numerical simulations of the ODE system (2.2)-(2.7) showing the disease-free dynamics. The initial conditions used for this simulation is given in Table 4 and the parameters are given in Table 2. The vaccination rate is set as σ = 0.85 with vaccine efficacy f = 0.85, while the waning rate for both vaccine-induced and post-recovery immunity are ϕ1 = ϕ2 = 0.3. For these set of parameters, the corresponding control reproduction number is ℛcon = 0.9896, which is less than 1. As observed from the results in the right panel of this figure, the dynamics of the infected compartments (E1, E2, and I) tends to zero while the susceptible populations (S and Sv) initially decrease and then stabilize at some positive numbers (see left panel), which leads to the disease-free equilibrium (DFE), as stated in Theorem 3.1. Using the same set of parameters but increasing the immunity waning rates to ϕ1 = ϕ2 = 0.5, the corresponding control reproduction number increases to ℛcon = 1.2322. The results for this scenario are shown in Figure 3.As stated by Theorem 3.1, the DFE is unstable for this scenario since ℛcon > 1. From the right panel of Figure 3, we observe that unlike in Figure 2, the dynamics of the infected compartments do not tend to zero, rather they saturate as some positive numbers. Similar behaviours can be seen in the inset of the plot in the left panel of Figure 3 for the susceptible populations (S and Sv) as they reduced drastically but don’t tend to zero. This result shows that the DFE is unable when ℛcon > 1, and as a result, the disease cannot be eradicated from the population.
In the next subsection, we analyze the endemic equilibrium and show that our model undergoes the phe-nomenon of backward bifurcation.
3.3 The Backward Bifurcation and Endemic Equilibrium (EE)
Bifurcation is the change in the qualitative structure of the flow of a differential equation with respect to parameter values. The point where bifurcation occurs is called the bifurcation point. Examples of bifurcations in dynamical systems include transcritical, pitchfork, Hopf, saddle-node, forward, backward, e.t.c. [3, 20, 52]. However, the bifurcation type that is related to our study is the backward bifurcation, which is explained in [5] based on the provision of the center manifold theory.
The COVID-19 model presented in (2.2)-(2.7) has a positive and unique endemic equilibrium point, which can be interpreted as the point where at least one infected population exists and is non-zero. Let be the endemic equilibrium point. We seek to solve the ODE system (2.2)-(2.7) in terms of the parameter λc as follows: The force of infection at steady state can be written as
The model equations (2.2)-(2.7) can be re-expressed as where the newly defined parameters g1, g2, …, g7 are
Solving (3.9)-(3.14) simultaneously gives
Again, solving the equations in (3.15) simultaneously gives
Substituting the expression for Sv in (3.15) into (3.16) and factorizing, we have
Similarly, we get
Substituting (3.18) into (3.17) and factorizing gives
Similarly, we have
Substituting all the expressions in (3.18), (3.19), (3.20) and (3.21) into (3.8), it can be proved that the non-zero equilibria of the model satisfies where and
We can easily establish the positiveness of the coefficient x0 by substituting the respective value of all the parameters in x0 and subject the resultant expression into full expansion. The negative signs will disappear and the expression will become positive. The complete expansion is cumbersome and thus ignored. More-so, x4 is positive if and only if ℛcon < 1 and negative if and only if ℛcon > 1. This shows that the number of positive real roots of (3.22) is dependent on the signs of x1, x2 and x3. These signs can be analyzed using the Decartes Rule of signs on (3.22). The possible roots of equation (3.22) are given in Table 3. Hence, we have the following theorem.
The COVID-19 model (2.2)-(2.7)
has a unique endemic equilibrium if ℛcon > 1 and whenever cases 1, 2, 6 and 8 of Table 3 are satisfied;
could have more than one endemic equilibrium if ℛcon > 1 and cases 3, 4, 5 and 7 of Table 3 are satisfied;
could have 2 or more endemic equilibria if Rcon < 1 and cases 2-8 of Table 3 are satisfied.
The COVID-19 model has at least one positive endemic equilibrium whenever ℛcon > 1 but could have zero, two or four positive endemic equilibrium whenever ℛcon < 1.
It is worth noting that the exact expression for the endemic equilibria can be obtained by solving for λc in (3.22) to obtain the roots and then substitute its positive value into the expressions in (3.18), (3.19), (3.20) and (3.21). The existence of multiple endemic equilibria when the control reproduction number is less than unity established using Theorem 3.2, Lemma 2 and Table 3 suggests the possibility of backward bifurcation [3, 52]. As previously mentioned, backward bifurcation is a phenomenon where a stable DFE coexists with a stable endemic equilibrium when the control reproduction number is less than 1. We present the bifurcation diagrams in Figure 4.
The bifurcation diagram on the top left panel of Figure 4 shows the behaviour of the endemic model (2.2)-(2.7) when the control reproduction number ℛcon is less than, equal and greater than 1. The parameter values used are given in Table 1 with σ = ϕ1 = ϕ2 = 0.5 and vaccine efficacy set at f = 0.55. The red-circled line on the diagram shows that the endemic equilibrium is unstable for ℛcon < 1 while the blue-circled part of the diagram shows the stability of the endemic equilibrium for ℛcon > 1. The blue-dashed line shows the stability of the disease-free equilibrium for ℛcon < 1 while the red-dashed line shows that the disease-free equilibrium is unstable for Rcon > 1. We observe from this result that there are two steady-states of the system (DFE and endemic equilibrium) when ℛcon > 1. As ℛcon decreases, we have the emergence of an endemic equilibrium (red-circled line) at ℛcon = 1, which is unstable, that co-exist with a stable endemic equilibrium (blue-circled line) and a stable disease-free equilibrium (blue-dashed line). The backward movement of the equilibrium that emerge at ℛcon = 1 clearly shows that our model exhibits the phenomenon of backward bifurcation, where a stable disease-free equilibrium point co-exists with a stable endemic equilibrium point. The epidemiological implication of this is that the control reproduction number being less than unity is no longer sufficient to guarantee the disease eradication. More control measures would be required to eradicate the disease. For this reason, we considered the effects of increase in vaccine efficacy and reduction in transmission rates on the bifurcation diagram.
The top right panel shows the behavior of the bifurcation diagram when the vaccine efficacy f is reduced to 0.15 while other parameter values remain the same. We can see from this result that bi-stability persist, and now exist for all ℛcon < 1. The bottom left panel shows the behavior of the model when the vaccine efficacy f is increased to 0.65 while other parameter values remain the same. We can easily notice the disappearance of bi-stability. Also, we observe that the endemic equilibrium is unstable for some values of ℛcon > 1. For these values of ℛcon, the DFE equilibrium is stable (not shown). This clearly underlines the need for a highly efficacious vaccine in controlling the spread of the disease and its eventual eradication. The bottom right panel shows the behaviour of the model when the transmission rates β1 and β2 are reduced to 0.02189 and 0.03521, respectively (which is 10% of its baseline value) while other parameter values remain the same. We observe from this result that there is also no bi-stability in the system for this scenario. This result emphasizes the importance of reduction in community transmission in the fight against the disease. The lower/higher the transmission rates the lower/higher the disease prevalence. On the other hand the higher/lower the vaccine efficacy the higher/lower the protection against the disease. These two parameters are important in controlling the disease spread and its eventual eradication.
4 Numerical simulations
We present the numerical simulations of the endemic model (2.2)-(2.7) for different scenarios to study the effect of immunity waning on the dynamics of COVID-19. The parameters used are related to the demography of South Africa. These parameters are given in Table 2. We studied the effect of immunity waning on the control reproduction number and the entire diseases dynamics.
In Figure 5, we present the contour plots of the control reproduction number ℛcon in terms of the vaccine immunity waning parameter ϕ1 and the vaccination rate σ for vaccine efficacy f = 0.5 (left panel) and f = 0.85 (right panel). We observe from the results in this figure that the control reproduction number, ℛcon decreases as the vaccination rate increases and increases as the vaccine immunity waning rates increase. When the vaccination rate is relatively small, vaccine immunity waning has little effect on the reproduction number, compare to when a large fraction of the population is vaccinated. In this case, an increase in the waning rate increases the control reproduction number. In addition, comparing the results in the left panel of this figure for vaccine efficacy f = 0.5 to that of the right panel for f = 0.85, we notice that the vaccine immunity waning rate has more effect on the control reproduction number when the vaccine is more effective.
Similar results are shown in Figure 6 for the transmission and vaccination rates. In the top left panel of this figure, we show the control reproduction number as a function of the pre-symptomatic transmission rate (β1) and the vaccination rate (σ) for ϕ1 = 0.1. For the same vaccine waning rate, we present the reproduction number as a function of the symptomatic transmission rate (β2) and the vaccination rate (σ).
The remaining parameters are given in Table 2. As expected, the control reproduction number increases as the transmission rates β1 and β2 increase, although, the symptomatic transmission rate has more effect on the control reproduction number compared to the pre-symptomatic transmission rate. The results in the lower panel of Figure 6 are similar to those in the top panel but for a higher waning rate of vaccine-induced immunity (ϕ1 = 0.5). Comparing these results with those in the top panel, we notice that an increase in ϕ1 increases the dependence of the control reproduction number on the vaccination rate. Similarly, when ϕ1 is small (say ϕ1 = 0.1 as shown in the top panel of Figure 6), the transmission rates have more effect on the reproduction number when a small fraction of the population is vaccinated compared to when the vaccination rate is high. However, the effect of the transmission rates on the control reproduction number when vaccination rate is higher increases as the immunity waning rate increases. This can easily be seen by comparing the results in the top panel to those in the bottom panel for high vaccine efficacy.
Next, we study the effect of immunity waning on the dynamics of COVID-19 using South Africa as a case study. The results of this study are presented in terms of the time dynamics of the symptomatic infectious population. The initial conditions used are presented in Table 4 and are based on the South African COVID-19 statistics as of August 29, 2021, and were obtained from the Department of health COVID-19 report, Republic of South Africa [11]. In addition, the parameters used are also related to the demography of South Africa and are presented in Table 2. Figure 7 shows the dynamics of the symptomatic population over time for different vaccine immunity waning rate, ϕ1 (left panel) and post-recovery immunity waning rate, ϕ2 (right panel). We observe from these results that vaccine immunity waning has much more effect on the disease dynamics than the post-recovery immunity waning. For this scenario, we used ϕ1 = 0.3, 0.5, 0.75, and 0.9, which give the corresponding control reproduction numbers ℛcon = 0.9896, 1.2322, 1.4504, and 1.5516, respectively. By varying the vaccine waning rate from ϕ1 = 0.3 through ϕ1 = 0.9, the control reproduction number goes from ℛcon = 0.9896 (ℛcon < 1) to ℛcon = 1.5516 (ℛcon > 1). Changing the disease dynamics from disease-free to endemic (See the left panel of Figure 7). However, this is not the case for the post-recovery immunity waning rate ϕ2. This parameter has little effect on the dynamics of the disease (right panel of Figure 7). In addition, since this control reproduction number is independent of this parameter, ℛcon = 1.2322 for all values of ϕ2.
In Figure 8, we study the effect of the vaccination rate, σ and the immunity waning rates ϕ1 and ϕ2 on the disease dynamics. In this figure, we reduce the vaccination rate, from σ = 0.85 used in Figure 7 to σ = 0.5. Our goal is to study how a reduction in vaccination rate would affect the disease dynamics. Comparing the results in Figure 7 to those in Figure 8, we observe that a decrease in vaccination rate lead to more infections in the population. Unlike the case where vaccination rate is high, when vaccination rate is σ = 0.5, there are no disease-free dynamics for all values of the waning parameter ϕ1 used. This results suggests that we need to keep the vaccination rate high in order to eradicate the disease.
The result in Figure 9 is for when the vaccine efficacy is reduced from f = 0.85 used in Figure 7 to f = 0.65, while the vaccination rate remains σ = 0.85. Similar to the results in Figure 8, the disease is endemic in the population for all values of the waning parameter ϕ1. The increase in infections observed in Figure 9 relative to Figure 7 is due to a decrease in the efficacy of the vaccine, since the same vaccination rate is used for the two scenarios. This suggests that in addition to high vaccination rate, a high vaccine efficacy is also essential to eradicate the disease.
Lastly, we study the effect of the transmission rates on the disease dynamics. In the top panel of Figure 10, we present the results for when the pre-symptomatic transmission rate is increased to β1 = 0.3284, while the symptomatic transmission rate remains the same as in the remaining figures (β2 = 0.3521). On the other hands, in the lower panel we increased the symptomatic transmission rate to β2 = 0.4049, while the pre-symptomatic transmission rate remains the same as in the other figures (β1 = 0.2189). As expected, there are more infections are the immunity waning rates increase. Comparing the results in this figure to those in Figure 7, we observe that an increase in the transmission rates will increase the disease spread in the community making it difficult to eradicate the disease. In addition, an increase in the symptomatic transmission rate significantly affect the disease dynamics relative to the pre-symptomatic transmission rate. Note that for all the examples considered, the post-recovery immunity waning rate has little effect on the disease dynamics compared to the waning rate of vaccine-induced immunity.
5 Discussion
We have developed and analyzed an endemic model for COVID-19 to study the effect of immunity waning on the dynamics of thes disease. Our model has six compartment: unvaccinated susceptible (S), vaccinated susceptible (Sv), exposed (not infectious) (E1), pre-symptomatic (infectious but no symptoms), symptomatic (I), and recoversed (R). We considered two type of immunity waning in our model: vaccine immunity waning and post-recovery immunity waning. Vaccine immunity waning occurs when a vaccinated individual in the Sv compartment loses their immunity and return to the unvaccinated susceptible compartment S. On the other hand, post-recovery immunity waning occurs when recovered individuals lose their immunity and return to the unvaccinated susceptible compartment. Both immunity waning described above can also be interpreted as being susceptible to another strain of COVID-19. Our model implements a sterilizing vaccine that reduces the chances of getting infected. We apply our model to the South African COVID-19 scenario.
We analyzed our model rigorously by first computing the control reproduction number for the disease-free equilibrium. After which we analyzed the endemic equilibrium and show that our model exhibit the backward bifurcation phenomenon. We examined the condition for the existence of bifurcation and endemic equilibrium. These conditions are stated in Theorem 3.2. Furthermore, Lemma 2 provides the condition for the existence of multiple endemic equilibria with respect to the conditions specified on Table 3. This leads to the phenomenon where stable DFE co-exists with a stable endemic equilibrium when the control reproduction number is less than unity. Figure 4 shows the bifurcation diagrams for different scenarios. It can be understood from this diagrams that if the vaccine efficacy is reduced, bi-stability persist but when increased, the bi-stability disappears and the disease may die out. A better result is achieved when the transmission rates are reduced, thereby eliminating the bi-stability totally. This clearly underlines the need for a highly efficacious vaccine and reduction in disease transmission to achieve disease eradication.
We numerically studied the effects of vaccination and the waning of vaccine immunity on the control reproduction number. These results were presented in form of contour plots and show that the control reproduction number decreases as the vaccination rate increases. However, it increases as the vaccine immunity waning rate increases. As the vaccine efficacy increases, the effects these parameters (vaccination rate and vaccine immunity waning rate) become more pronounced. There is also a decrease in the control reproduction number, although, the decrease is more when the vaccination rate is high. When a small fraction of the population is vaccinated, vaccine immunity waning has little effect on the control reproduction number, compared to when vaccination rate is high. This is because when a large fraction of the population is vaccinated, we expect a reduction in infection. But when these people loose their immunity at a faster rate, the expected reduction in infection will not be realized, leading to a significant spread of the disease. This is not the case when only a small fraction of the population are vaccinated from the onset. In this case, there will be no much difference whether they loose their immunity slowly or at a faster rate. We also studied the effect of the transmission rates together with vaccination on the reproduction number. Our results show that the control reproduction number increases as the transmission rates increase, although, the symptomatic transmission rate has more effect on the control reproduction compared to the pre-symptomatic transmission rate. This may be because infectious individuals spend fewer days (2 days) in the pre-symptomatic compartment compared to the symptomatic compartment, where they spend 7 days. In addition, the transmission rates have little effect on the reproduction number when vaccination rate is high, and this effect increases as the vaccination rate drops. In other words, when we vaccinate a large fraction of the population with a relatively effective vaccine, the transmission rate has a little effect on the reproduction number. However, if we vaccinate only a small fraction of the population with the same vaccine when the transmission rate is higher, the control reproduction number increases too.
Furthermore, we investigated the effects of different parameters of the model on the disease dynamics. Our main focus is studying how vaccination and the waning of immunity affect the disease dynamics. The results of these investigations were presented in terms of the dynamics of the symptomatic population over time and show that vaccine immunity waning has more effect on the disease dynamics compared to the post-recovery immunity waning. By varying the vaccine immunity waning parameter, we realized that when this parameter is low, the diseased-free equilibrium (DFE) of the model is linearly asymptotically stable. This means that the disease can be eradicated when people retain their vaccine immunity for a longer time. However, as this waning rate is increased, the DFE becomes unstable, there by making the disease endemic in the population. On the other hand, the DFE remains unstable for all values of the post-recovery immunity waning parameter and a moderate vaccine waning parameter. By reducing the vaccination rate by ≈ 40%, we noticed that there is no value of vaccine immunity waning rate for which the DFE is linearly stable. As a result of this, the disease cannot be eradicated when the vaccination rate is low. In summary, keeping the vaccination rate high and preventing lose of immunity is required to eradicate COVID-19.
An interest future direction for this work is to fit our model to the reported cases of COVID-19 for a specific location/country. In this way, the parameters would be more related to that location and the predictions would be more representative of the location. Another interesting future direction is to extend the model to have explicit cycles of infection or loss of immunity. For this case, there will be all six compartments of our current model in each cycle such that individuals that lose their immunity in one cycle will go to the susceptible compartment of the next cycle rather than the susceptible compartment of their current cycle as it is in our current model, and the disease will progress from there to the other compartments of the next cycle. An advantage of structuring the model this way is that one can easily track the number of people getting infected at each infection cycle. In addition, different cycles of infection or loss of immunity can be interpreted as being susceptible to different variants of COVID-19, where different vaccine efficacy could also be used. This will increase the flexibility of the model. It would also be worthwhile to stratify the population by age and extend the model to have age-structure. This is particularly important since individuals younger than 12 years have not been approved for vaccination yet. It would also make the model more suitable for answering age related questions.
Some limitations of the modeling framework presented in this paper includes modeling vaccination as a single dose despite the fact that majority of the COVID-19 vaccines available today require more than one does to get a good level of protection. In addition, our model assumes that vaccine efficacy is the same, either one is getting the vaccine for the first time or one had already experienced waning of immunity. This may not be the case in reality as getting the vaccine the second time may lead to a better/stronger protection. Also the efficacy of the vaccines for different variants of COVID-19 may be different. Another limitation of our model is assuming that the population is well-mixed. This may not be so in reality as contact rates and mixing patterns may vary between individuals depending on their age and activity level. We would also like to include lack of stochasticity as a limitation to our modeling framework. Despite all these limitations, we believe that our modeling framework has provided some insights into the endemic dynamics of COVID-19 and has also provided some reasonable suggestions on the important parameters to be considered to effective eradicate the disease.
Data Availability
All data produced in the present work are contained in the manuscript
Competing Interest
We declare that there is no competing interest as far as this work is concerned.
Data Availability
All used data are properly referenced in the body of the work.
Acknowledgement
No funding is received for this project.
References
- [1].↵
- [2].↵
- [3].↵
- [4].↵
- [5].↵
- [6].↵
- [7].↵
- [8].↵
- [9].
- [10].
- [11].↵
- [12].↵
- [13].↵
- [14].↵
- [15].↵
- [16].↵
- [17].↵
- [18].↵
- [19].↵
- [20].↵
- [21].↵
- [22].↵
- [23].↵
- [24].↵
- [25].↵
- [26].↵
- [27].↵
- [28].
- [29].↵
- [30].↵
- [31].↵
- [32].↵
- [33].↵
- [34].↵
- [35].↵
- [36].↵
- [37].↵
- [38].
- [39].↵
- [40].↵
- [41].↵
- [42].↵
- [43].↵
- [44].↵
- [45].↵
- [46].↵
- [47].↵
- [48].↵
- [49].
- [50].↵
- [51].↵
- [52].↵
- [53].↵
- [54].↵
- [55].↵
- [56].↵
- [57].↵
- [58].
- [59].↵
- [60].↵
- [61].↵
- [62].↵
- [63].↵
- [64].↵
- [65].↵
- [66].↵
- [67].
- [68].↵
- [69].
- [70].↵
- [71].