Assessing the potential impact of immunity waning on the dynamics of COVID-19: an endemic model of COVID-19 =========================================================================================================== * Musa Rabiu * Sarafa A. Iyaniwura ## 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. Keywords * COVID-19 * endemic model * backward bifurcation * non-pharmaceutical intervention * vaccination * immunity waning ## 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 *R**o*. 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 (*S**v*). We also divided the infected population into three: exposed (not infectious) (*E*1), pre-symptomatic (infectious without symptoms) (*E*2) 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. View this table: [Table 1:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/T1) Table 1: Model variables and descriptions. Our model assumes that individuals in the exposed compartment (E1) are not infectious, while those in the pre-symptomatic compartment (E2) are infectious but do not show symptoms. 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 ![Formula][1] where *N* (*t*) = *S*(*t*) + *S**v*(*t*) + *E*1(*t*) + *E*2(*t*) + *I*(*t*) + *R*(*t*) denotes the total population at any time *t*. The differential equations of our model are given by ![Formula][2] ![Formula][3] ![Formula][4] ![Formula][5] ![Formula][6] ![Formula][7] under the following initial conditions ![Formula][8] 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 *E*1 to *E*2 and *E*2 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 *δ*. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F1) Figure 1: Model schematic. Compartments are as follows: Unvaccinated susceptible (S); vaccinated susceptible (Sv); exposed (E1); pre-symptomatic infectious (E2); symptomatic infectious (I) and recovered (R). Black solid arrows show the flow of individuals through the compartments of the model from susceptible to recovered. The green arrows show the return of individuals to the susceptible compartment due to immunity waning. The blue arrow indicates birth and immigration into the susceptible compartment, while the red arrows show natural death and death due to COVID-19. The rates of transitioning from one compartment to the other are shown beside the arrows (see equations (2.2)-(2.7) for more details). 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. View this table: [Table 2:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/T2) Table 2: Model parameters, descriptions, and values. The birth and immigration rate Λ, and the death rates µ and δ are based on South African demography, and were obtained from the official Republic of South Africa government websites [1, 61]. ## 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. Lemma 1. *Let ϒ represents the closed set defined as* ![Formula][9] *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 ![Formula][10] Observe that equation (3.2) is a first-order linear differential equation. The integrating factor is chosen as ![Formula][11] where *λ*(*η*) ≡ *λ*(*S, S**v*, *E*1, *E*2, *I, R*) so that the solution to (3.2) with equality is ![Formula][12] Evaluating the equation above and using the fact that *S*(0) ≥ 0, we have ![Formula][13] which implies that ![Formula][14] 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 *S**v*, *E*1, *E*2, *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 ![Formula][15] whose solution is given by ![Formula][16] where exp(*µt*) is the integrating factor. Taking the limit of both sides as *t* → ∞ gives ![Formula][17] This means that the expression in (3.5) is bounded above by ![Graphic][18] 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 ![Formula][19] where ![Graphic][20]. 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 ![Formula][21] 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 ![Formula][22] The first term in the control reproduction number (3.7) accounts for the infections caused by the pre-symptomatic individuals in *E*2, while the second term account for those caused by symptomatic individuals in (*I*). Using (3.7), we establish the following result. Theorem 3.1. *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 (*E*1, *E*2, and *I*) tends to zero while the susceptible populations (*S* and *S**v*) 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 *S**v*) 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. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F2) Figure 2: Disease-free dynamics. The time dynamics of the susceptible (S), susceptible-vaccinated (Sv) and recovered (R) populations (left panel), and those for the exposed (E1), pre-symptomatic (E2), and symptomatic (I) populations (right panel), showing the disease-free equilibrium dynamics. Vaccination rate is σ = 0.85, with vaccine efficacy, f = 0.85. The vaccine and post-recovery immunity waning rates are ϕ1 = ϕ2 = 0.3, corresponding to a control reproduction number of ℛcon = 0.9896. Remaining parameters are as given in Table 2. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F3) Figure 3: Endemic dynamics. The time dynamics of the susceptible (S), susceptible-vaccinated (Sv) and recovered (R) populations (left panel), and those for the exposed (E1), pre-symptomatic (E2), and symptomatic (I) populations (right panel), for when the disease is endemic in the population. Vaccination rate is σ = 0.85, with vaccine efficacy, f = 0.85. The vaccine and post-recovery immunity waning rates are ϕ1 = ϕ2 = 0.5, corresponding to a control reproduction number of ℛcon = 1.2322. Remaining parameters are given in Table 2. 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 ![Formula][23] 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 ![Formula][24] The model equations (2.2)-(2.7) can be re-expressed as ![Formula][25] ![Formula][26] ![Formula][27] ![Formula][28] ![Formula][29] ![Formula][30] where the newly defined parameters *g*1, *g*2, …, *g*7 are ![Formula][31] Solving (3.9)-(3.14) simultaneously gives ![Formula][32] Again, solving the equations in (3.15) simultaneously gives ![Formula][33] Substituting the expression for *S**v* in (3.15) into (3.16) and factorizing, we have ![Formula][34] Similarly, we get ![Formula][35] Substituting (3.18) into (3.17) and factorizing gives ![Formula][36] Similarly, we have ![Formula][37] ![Formula][38] 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 ![Formula][39] where ![Formula][40] and ![Formula][41] We can easily establish the positiveness of the coefficient *x* by substituting the respective value of all the parameters in *x* 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, *x*4 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 *x*1, *x*2 and *x*3. 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. View this table: [Table 3:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/T3) Table 3: Possible positive real roots of equation (3.22) for ℛcon < 1 and ℛcon > 1. Theorem 3.2. *The COVID-19 model* (2.2)*-*(2.7) 1. *has a unique endemic equilibrium if ℛ**con* > 1 *and whenever cases 1, 2, 6 and 8 of Table 3 are satisfied;* 2. *could have more than one endemic equilibrium if ℛ**con* > 1 *and cases 3, 4, 5 and 7 of Table 3 are satisfied;* 3. *could have 2 or more endemic equilibria if R**con* *<* 1 *and cases 2-8 of Table 3 are satisfied*. Lemma 2. *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. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F4) Figure 4: Bifurcation Diagrams. The bifurcation diagrams for different values of vaccine efficacy f and different transmission rates β1 and β2. The top-left panel is the original bifurcation diagram for f = 0.55 and σ = ϕ1 = ϕ2 = 0.5 while other parameter values remain unchanged. The top-right panel is the bifurcation diagram for f = 0.15 and σ = ϕ1 = ϕ2 = 0.5 while other parameter values remain unchanged. The bottom-left panel shows the bifurcation diagram for f = 0.65 and σ = ϕ1 = ϕ2 = 0.5 while other parameter values remain unchanged. The bottom-right panel shows the bifurcation diagram f = 0.55, σ = ϕ1 = ϕ2 = 0.5 and β1 = 0.0219, β2 = 0.0352 (which is 10% of its baseline value) while other parameter values remain unchanged. 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 *R**con* > 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. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F5) Figure 5: Effect of vaccination and vaccine immunity waning on the control reproduction number. Left panel: vaccine efficacy, f = 0.5. Right panel: vaccine efficacy, f = 0.85. Remaining parameters are given in Table 2. 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 (*σ*). ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F6) Figure 6: Effect of vaccination rate and disease transmission rates on the control reproduction number. Left panel: Pre-symptomatic transmission rate β1, and right panel: Symptomatic transmission rate β2. In the first row, the vaccine immunity waning rate is ϕ1 = 0.1, and ϕ1 = 0.5 in the second row. Vaccine efficacy is set as f = 0.85. Remaining parameters are given in Table 2. 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. View this table: [Table 4:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/T4) Table 4: Model initial conditions. These initial conditions 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]. The unvaccinated susceptible population was computed using S(0) = N (0) − (Sv(0) + E1(0) + E2(0) + I(0) + R(0)). ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F7) Figure 7: Effect of immunity waning on disease dynamics. The dynamics of the symptomatic population over time for different vaccine immunity waning rates (left panel) and post-recovery immunity waning rates (right panel). In the left panel, for the waning rates ϕ1 = 0.3, 0.5, 0.75, and 0.9, the corresponding control reproduction numbers are ℛcon = 0.9896, 1.2322, 1.4504, and 1.5516, respectively, while in the right panel, ℛcon = 1.2322 for all values of ϕ2. Vaccination rate used is σ = 0.85, with vaccine efficacy, f = 0.85. Remaining parameters are given in Table 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. ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F8) Figure 8: Effect of vaccination rate and immunity waning on disease dynamics. The dynamics of symptomatic population over time for different vaccine immunity waning rates (left panel) and post-recovery immunity waning rates (right panel). In the left panel, for the waning rates ϕ1 = 0.3, 0.5, 0.75, and 0.9, the corresponding control reproduction numbers are ℛcon = 1.2474, 1.5230, 1.7442, and 1.8392, respectively, while in the right panel, ℛcon = 1.5230 for all values of ϕ2. Vaccination rate is reduced to σ = 0.5, with vaccine efficacy, f = 0.85. Remaining parameters are given in Table 2. 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. ![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F9) Figure 9: Effect of vaccine efficacy and immunity waning on disease dynamics. The dynamics of symptomatic population over time for different vaccine immunity waning rates (left panel) and post-recovery immunity waning rates (right panel). In the left panel, for the waning rates ϕ1 = 0.3, 0.5, 0.75, and 0.9, the corresponding control reproduction numbers are ℛcon = 1.3768, 1.5623, 1.7292, and 1.8065, respectively, while in the right panel, ℛcon = 1.5623 for all values of ϕ2. Vaccination rate is σ = 0.85, with vaccine efficacy reduced to f = 0.65. Remaining parameters are given in Table 2. 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. ![Figure 10:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/26/2021.10.23.21265421/F10.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2021/10/26/2021.10.23.21265421/F10) Figure 10: Effects of disease transmission rates and immunity waning on disease dynamics. The dynamics of symptomatic population over time for different vaccine immunity waning rates (left) and post-recovery immunity waning rates (right). In the top panel, the pre-symptomatic transmission rate is β1 = 0.3284 and the symptomatic transmission rate is β2 = 0.3521, while in the lower panel, the pre-symptomatic transmission rate is β1 = 0.2189 and the symptomatic transmission rate is β2 = 0.4049. In the top left panel, for the waning rates ϕ1 = 0.3, 0.5, 0.75, and 0.9, the corresponding control reproduction numbers are ℛcon = 1.0680, 1.3297, 1.5652, and 1.6744, respectively, while in the top right panel, Rcon = 1.5623 for all values of ϕ2. Similarly, in the bottom left panel, for the waning rates ϕ1 = 0.3, 0.5, 0.75, and 0.9, the corresponding control reproduction numbers are ℛcon = 1.1146, 1.3877, 1.6336, and 1.7475, respectively, while in the bottom right panel, ℛcon = 1.3877 for all values of ϕ2. Remaining parameters are given in Table 2. ## 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 (*S**v*), exposed (not infectious) (*E*1), 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 *S**v* 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. * Received October 23, 2021. * Revision received October 23, 2021. * Accepted October 25, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. [1].STATISTICS SOUTH AFRICA. STATISTICAL RELEASE, Mid-year population estimates, 2019, (accessed August 28, 2021). [http://www.statssa.gov.za/publications/P0302/P03022019.pdf](http://www.statssa.gov.za/publications/P0302/P03022019.pdf). 2. [2]. Sean C Anderson, Andrew M Edwards, Madi Yerlanov, Nicola Mulberry, Jessica E Stockdale, Sarafa A Iyaniwura, Rebeca C Falcao, Michael C Otterstatter, Michael A Irvine, Naveed Z Janjua, et al. Quantifying the impact of covid-19 control measures using a bayesian model of physical distancing. PLoScomputational biology, 16(12):e1008274, 2020. 3. [3]. Ally Yeketi Ayinla, Wan Ainun Mior Othman, and Musa Rabiu. A mathematical model of the tuber-culosis epidemic. Acta Biotheoretica, pages 1–31, 2021. 4. [4]. Sarita Bugalia, Vijay Pal Bajiya, Jai Prakash Tripathi, Ming-Tao Li, and Gui-Quan Sun. Mathematical modeling of covid-19 transmission: the roles of intervention strategies and lockdown. MathematicalBiosciences and Engineering, 17(5):5961–5986, 2020. 5. [5]. Carlos Castillo-Chavez and Baojun Song. Dynamical models of tuberculosis and their applications. Mathematical Biosciences & Engineering, 1(2):361, 2004. 6. [6]. Monika Chaudhary, PR Sodani, and Shankar Das. Effect of covid-19 on economy in india: Some reflections for policy and programme. Journal of Health Management, 22(2):169–180, 2020. 7. [7]. Zhangkai J Cheng and Jing Shan. 2019 novel coronavirus: where we are and what we know. Infection, 48(2):155–163, 2020. 8. [8]. Pyoeng Gyun Choe, Chang Kyung Kang, Hyeon Jeong Suh, Jongtak Jung, Kyoung-Ho Song, Ji Hwan Bang, Eu Suk Kim, Hong Bin Kim, Sang Won Park, Nam Joong Kim, et al. Waning antibody responses in asymptomatic and symptomatic sars-cov-2 infection. Emerging infectious diseases, 27(1):327, 2021. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 9. [9].Republic of South Africa Department of health. COVID-19 Statistics in SouthAfrica, (accessed August 29, 2021). [https://sacoronavirus.co.za/2021/08/29/update-on-covid-19-sunday-29-august-2021/](https://sacoronavirus.co.za/2021/08/29/update-on-covid-19-sunday-29-august-2021/). 10. [10].Republic of South Africa Department of health. COVID-19 Total Vaccination Summary, (accessed August 29, 2021). [https://sacoronavirus.co.za/latest-vaccine-statistics/](https://sacoronavirus.co.za/latest-vaccine-statistics/). 11. [11].Republic of South Africa Department of Statistics. Mid-year population estimates, 2021, (accessed August 29, 2021). [http://www.statssa.gov.za/publications/P0302/P03022021.pdf](http://www.statssa.gov.za/publications/P0302/P03022021.pdf). 12. [12]. Amelie Desvars-Larrive, Elma Dervic, Nils Haug, Thomas Niederkrotenthaler, Jiaying Chen, Anna Di Natale, Jana Lasser, Diana S Gliga, Alexandra Roux, Johannes Sorger, et al. A structured open dataset of government interventions in response to covid-19. Scientific Data, 7(1):1–9, 2020. 13. [13]. Odo Diekmann, Johan Andre Peter Heesterbeek, and Johan AJ Metz. On the definition and the computation of the basic reproduction ratio r 0 in models for infectious diseases in heterogeneous populations. Journal of mathematical biology, 28(4):365–382, 1990. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00178324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2117040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DG35000001&link_type=ISI) 14. [14]. Mathias Dzobo, Itai Chitungo, and Tafadzwa Dzinamarira. Covid-19: a perspective for lifting lockdown in zimbabwe. The Pan African Medical Journal, 35(Suppl 2), 2020. 15. [15]. Lilybeth Fontanesi, Daniela Marchetti, Cristina Mazza, Serena Di Giandomenico, Paolo Roma, and Maria Cristina Verrocchio. The effect of the covid-19 lockdown on parents: A call to adopt urgent measures. Psychological Trauma: Theory, Research, Practice, and Policy, 12(S1):S79, 2020. 16. [16]. Tapiwa Ganyani, Cecile Kremer, Dongxuan Chen, Andrea Torneri, Christel Faes, Jacco Wallinga, and Niel Hens. Estimating the generation interval for coronavirus disease (covid-19) based on symptom onset data, march 2020. Eurosurveillance, 25(17):2000257, 2020. 17. [17]. Mukoro Duke George. A novel corona virus (covid-19) pandemic, pathogenesis, clinical features and management options, public health measures. International Journal of TROPICAL DISEASE & Health, pages 11–21, 2020. 18. [18]. Yair Goldberg, Micha Mandel, Yinon M Bar-On, Omri Bodenheimer, Laurence S Freedman, Eric Haas, Ron Milo, Sharon Alroy-Preis, Nachman Ash, and Amit Huppert. Waning immunity of the bnt162b2 vaccine: A nationwide study from israel. medRxiv, 2021. 19. [19]. Brian Greenwood. The contribution of vaccination to global health: past, present and future. Philosophical Transactions of the Royal Society B: Biological Sciences, 369(1645):20130433, 2014. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2013.0433&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24821919&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 20. [20]. John Guckenheimer and Philip Holmes. Nonlinear oscillations, dynamical systems, and bifurcations ofvector fields, volume 42. Springer Science & Business Media, 2013. 21. [21]. Ali Hamady, JinJu Lee, and Zuzanna A Loboda. Waning antibody responses in covid-19: what can we learn from the analysis of other coronaviruses? Infection, pages 1–15, 2021. 22. [22]. Nils Haug, Lukas Geyrhofer, Alessandro Londei, Elma Dervic, Amélie Desvars-Larrive, Vittorio Loreto, Beate Pinior, Stefan Thurner, and Peter Klimek. Ranking the effectiveness of worldwide covid-19 government interventions. Nature human behaviour, 4(12):1303–1312, 2020. 23. [23]. Mackenzie Henderson, Amanda Bragg, Germin Fahim, Monica Shah, and Evelyn R Hermes-DeSantis. A review of the safety and efficacy of vaccines as prophylaxis for clostridium difficile infections. Vaccines, 5(3):25, 2017. 24. [24]. Enahoro A Iboi, Calistus N Ngonghala, and Abba B Gumel. Will an imperfect vaccine curtail the covid-19 pandemic in the us? Infectious Disease Modelling, 5:510–524, 2020. 25. [25]. Hiroyasu Inoue and Yasuyuki Todo. The propagation of economic impacts through supply chains: The case of a mega-city lockdown to prevent the spread of covid-19. PloS one, 15(9):e0239251, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 26. [26]. Benjamin Ivorra, Miriam R Ferrández, María Vela-Pérez, and Angel Manuel Ramos. Mathematical modeling of the spread of the coronavirus disease 2019 (covid-19) taking into account the undetected infections. the case of china. Communications in nonlinear science and numerical simulation, 88:105303, 2020. 27. [27]. Sarafa Adewale Iyaniwura, Musa Abolade Rabiu, Jummy David, and Jude Dzevela Kong. Assessing the impact of adherence to non-pharmaceutical interventions and indirect transmission on the dynamics of covid-19: a mathematical modelling study. medRxiv, 2021. 28. [28]. Waasila Jassat, Caroline Mudara, Lovelyn Ozougwu, Stefano Tempia, Lucille Blumberg, Mary-Ann Davies, Yogan Pillay, Terence Carter, Ramphelane Morewane, Milani Wolmarans, et al. Difference in mortality among individuals admitted to hospital with covid-19 during the first and second waves in south africa: a cohort study. The Lancet Global Health, 2021. 29. [29]. Matthew D Johnston and Bruce Pell. A dynamical framework for modeling fear of infection and frustration with social distancing in covid-19 spread. arXiv preprint arxiv:2008.06023, 2020. 30. [30]. Natalie Kash, Michael A Lee, Ramya Kollipara, Christopher Downing, Jacqueline Guidry, and Stephen K Tyring. Safety and efficacy data on vaccines and immunization to human papillomavirus. Journal of clinical medicine, 4(4):614–633, 2015. 31. [31]. Dennis M Klinman, Mitsuhiro Takeno, Motohide Ichino, Mili Gu, Galina Yamshchikov, Gil Mor, and Jacqueline Conover. Dna vaccines: safety and efficacy issues. In Springer seminars in immunopathology, volume 19, pages 245–256. Springer, 1997. 32. [32]. Philip R Krause and Marion F Gruber. Emergency use authorization of covid vaccines—safety and efficacy follow-up considerations. New England Journal of Medicine, 383(19):e107, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMp2031373&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 33. [33]. Adam J Kucharski, Timothy W Russell, Charlie Diamond, Yang Liu, John Edmunds, Sebastian Funk, Rosalind M Eggo, Fiona Sun, Mark Jit, James D Munday, et al. Early dynamics of transmission and control of covid-19: a mathematical modelling study. The lancet infectious diseases, 20(5):553–558, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1473-3099(20)30144-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 34. [34]. Stephen A Lauer, Kyra H Grantz, Qifang Bi, Forrest K Jones, Qulu Zheng, Hannah R Meredith, Andrew S Azman, Nicholas G Reich, and Justin Lessler. The incubation period of coronavirus disease 2019 (covid-19) from publicly reported confirmed cases: estimation and application. Annals of internalmedicine, 172(9):577–582, 2020. 35. [35]. Shaoqing Lei, Fang Jiang, Wating Su, Chang Chen, Jingli Chen, Wei Mei, Li-Ying Zhan, Yifan Jia, Liangqing Zhang, Danyong Liu, et al. Clinical characteristics and outcomes of patients undergoing surgeries during the incubation period of covid-19 infection. EClinicalMedicine, 21:100331, 2020. 36. [36].Long-quan Li, Tian Huang, Yong-qing Wang, Zheng-ping Wang, Yuan Liang, Tao-bi Huang, Hui-yun Zhang, Weiming Sun, and Yuping Wang. Covid-19 patients’ clinical characteristics, discharge rate, and fatality rate of meta-analysis. Journal of medical virology, 92(6):577–583, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.2575&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 37. [37]. Qun Li, Xuhua Guan, Peng Wu, Xiaoye Wang, Lei Zhou, Yeqing Tong, Ruiqi Ren, Kathy SM Leung, Eric HY Lau, Jessica Y Wong, et al. Early transmission dynamics in wuhan, china, of novel coronavirus–infected pneumonia. New England journal of medicine, 2020. 38. [38]. Qun Li, Xuhua Guan, Peng Wu, Xiaoye Wang, Lei Zhou, Yeqing Tong, Ruiqi Ren, Kathy SM Leung, Eric HY Lau, Jessica Y Wong, et al. Early transmission dynamics in wuhan, china, of novel coronavirus–infected pneumonia. New England journal of medicine, 2020. 39. [39]. Zhihua Liu, Pierre Magal, Ousmane Seydi, and Glenn Webb. A covid-19 epidemic model with latency period. Infectious Disease Modelling, 5:323–337, 2020. 40. [40]. Shabir A Madhi, Vicky Baillie, Clare L Cutland, Merryn Voysey, Anthonet L Koen, Lee Fairlie, Sherman D Padayachee, Keertan Dheda, Shaun L Barnabas, Qasim Ebrahim Bhorat, et al. Safety and efficacy of the chadox1 ncov-19 (azd1222) covid-19 vaccine against the b. 1.351 variant in south africa. MedRxiv, 2021. 41. [41]. Matete R Magamela, Tafadzwa Dzinamarira, and Mbuzeleni Hlongwa. Covid-19 consequences on mental health: An african perspective. South African Journal of Psychiatry, 27(1):1–2, 2021. 42. [42]. Elisabeth Mahase. Covid-19: South africa pauses use of oxford vaccine after study casts doubt on efficacy against variant, 2021. 43. [43]. Daily Maverick. South Africa to Give First COVID-19 Vaccine Doses to President,Health Workers., (accessed on 29 August 2021). [https://www.dailymaverick.co.za/article/2021-02-17-south-africa-to-give-first-covid-19-vaccine-doses-to-president\-health-workers/](https://www.dailymaverick.co.za/article/2021-02-17-south-africa-to-give-first-covid-19-vaccine-doses-to-president\-health-workers/). 44. [44]. Liam G McCoy, Jonathan Smith, Kavya Anchuri, Isha Berry, Joanna Pineda, Vinyas Harish, Andrew T Lam, Seung Eun Yi, Sophie Hu, Laura Rosella, et al. Characterizing early canadian federal, provincial, territorial and municipal nonpharmaceutical interventions in response to covid-19: a descriptive analysis. CMAJ open, 8(3):E545, 2020. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiY21ham8iO3M6NToicmVzaWQiO3M6ODoiOC8zL0U1NDUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMC8yNi8yMDIxLjEwLjIzLjIxMjY1NDIxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 45. [45]. Malizgani Mhango, Itai Chitungo, and Tafadzwa Dzinamarira. Covid-19 lockdowns: impact on facility-based hiv testing and the case for the scaling up of home-based testing services in sub-saharan africa. AIDS and Behavior, 24:3014–3016, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10461-020-02939-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32488551&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 46. [46]. Rabiu Musa, Absalom E Ezugwu, and Godwin CE Mbah. Assessment of the impacts of pharmaceutical and non-pharmaceutical intervention on covid-19 in south africa using mathematical model. medRxiv, 2020. 47. [47]. Steady Mushayabasa, Ethel T Ngarakana-Gwasira, and Josiah Mushanyu. On the role of governmental action and individual reaction on covid-19 dynamics in south africa: A mathematical modelling study. Informatics in Medicine Unlocked, 20:100387, 2020. 48. [48]. Neha Pathak. The great invader: How covid attacks every organ, 2020. 49. [49]. Liangrong Peng, Wuyue Yang, Dongyan Zhang, Changjing Zhuge, and Liu Hong. Epidemic analysis of covid-19 in china by dynamical modeling. arXiv preprint arxiv:2002.06563, 2020. 50. [50]. Gashema Pierre, Annette Uwineza, and Tafadzwa Dzinamarira. Attendance to hiv antiretroviral collection clinic appointments during covid-19 lockdown. a single center study in kigali, rwanda. AIDSand Behavior, 24(12):3299–3301, 2020. 51. [51]. Musa Rabiu, Robert Willie, and Nabendra Parumasur. Analysis of a virus-resistant hiv-1 model with behavior change in non-progressors. Biomath, 9(1):2006143, 2020. 52. [52]. Musa Rabiu, Robert Willie, and Nabendra Parumasur. Mathematical analysis of a disease-resistant model with imperfect vaccine, quarantine and treatment. Ricerche di Matematica, 69(2):603–627, 2020. 53. [53]. Manish Sadarangani, Bahaa Abu Raya, Jessica M Conway, Sarafa A Iyaniwura, Rebeca Cardim Falcao, Caroline Colijn, Daniel Coombs, and Soren Gantt. Importance of covid-19 vaccine efficacy in older age groups. Vaccine, 39(15):2020–2023, 2021. 54. [54]. Muhammad Arif Nadeem Saqib, Shajee Siddiqui, Muhammad Qasim, Muhammad Azhar Jamil, Ibrar Rafique, Usman Ayub Awan, Haroon Ahmad, and Muhammad Sohail Afzal. Effect of covid-19 lockdown on patients with chronic diseases. Diabetes & Metabolic Syndrome: Clinical Research & Reviews, 14(6):1621–1623, 2020. 55. [55]. Sadeeq Sheshe, Abdullahi Nazifi, AM Labbo, GM Khalid, Abubakar Yahya, Umar Muhammad, and Aliyu Musa Haruna. Mechanism of antiviral immune response and covid-19 infection. Asian Journalof Immunology, pages 1–8, 2020. 56. [56]. Madhumita Shrotri, Annalan MD Navaratnam, Vincent Nguyen, Thomas Byrne, Cyril Geismar, Ellen Fragaszy, Sarah Beale, Wing Lam Erica Fong, Parth Patel, Jana Kovar, et al. Spike-antibody waning after second dose of bnt162b2 or chadox1. The Lancet, 398(10298):385–387, 2021. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S01406736(21)01642-1&link_type=DOI) 57. [57]. Surya Singh, Mujaheed Shaikh, Katharina Hauck, and Marisa Miraldo. Impacts of introducing and lifting nonpharmaceutical interventions on covid-19 daily growth rate and compliance in the united states. Proceedings of the National Academy of Sciences, 118(12), 2021. 58. [58]. Lauren Tindale, Michelle Coombe, Jessica E Stockdale, Emma Garlock, Wing Yin Venus Lau, Manu Saraswat, Yen-Hsiang Brian Lee, Louxin Zhang, Dongxuan Chen, Jacco Wallinga, et al. Transmission interval estimates suggest pre-symptomatic spread of covid-19. MedRxiv, 2020. 59. [59]. Pauline Van den Driessche and James Watmough. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Mathematical biosciences, 180(1-2):29–48, 2002. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0025-5564(02)00108-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12387915&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179220600004&link_type=ISI) 60. [60].Relief Web. COVAX Expects to Start Sending Millions of COVID-19 Vaccines toAfrica in February, (accessed on 29 August 2021). [https://reliefweb.int/report/world/covax-expects-start-sending-millions-covid-19-vaccines-africa-february](https://reliefweb.int/report/world/covax-expects-start-sending-millions-covid-19-vaccines-africa-february). 61. [61].Webometer. Total Coronavirus Deaths in South Africa, (accessed August 28, 2021). [https://www.worldometers.info/coronavirus/country/south-africa/](https://www.worldometers.info/coronavirus/country/south-africa/). 62. [62].World Health Organization (WHO). First COVID-19 COVAX vaccine doses administeredin Africa, (accessed September 16, 2021). [https://www.who.int/news/item/01-03-2021-first-covid-19-covax-vaccine-doses-administered-in-africa](https://www.who.int/news/item/01-03-2021-first-covid-19-covax-vaccine-doses-administered-in-africa). 63. [63].World Health Organization (WHO) et al. Naming the coronavirus disease(COVID-19) and the virus that causes it. 2020, 2020. [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it#:~:text=The%20International%20Committee%20on%20Taxonomy,two%20viruses%20are%20different](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it#:~:text=The%20International%20Committee%20on%20Taxonomy,two%20viruses%20are%20different). 64. [64]. Uwe Wollina, Ayşe Serap Karadağ, Christopher Rowland-Payne, Anca Chiriac, and Torello Lotti. Cutaneous signs in covid-19 patients: a review. Dermatologic therapy, 33(5):e13549, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) 65. [65].Worldometer. Corona-virus data, (accessed August 29, 2021). [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). 66. [66]. Chuang Xu, Yongzhen Pei, Shengqiang Liu, and Jinzhi Lei. Effectiveness of non-pharmaceutical interventions against local transmission of covid-19: An individual-based modelling study. InfectiousDisease Modelling, 6:848–858, 2021. 67. [67]. Zifeng Yang, Zhiqi Zeng, Ke Wang, Sook-San Wong, Wenhua Liang, Mark Zanin, Peng Liu, Xudong Cao, Zhongqiang Gao, Zhitong Mai, et al. Modified seir and ai prediction of the epidemics trend of covid-19 in china under public health interventions. Journal of thoracic disease, 12(3):165, 2020. 68. [68].New york times. They Were Infected With the Coronavirus. They Never Showed Signs., (accessed August 29, 2021). [https://www.nytimes.com/2020/02/26/health/coronavirus-asymptomatic.html](https://www.nytimes.com/2020/02/26/health/coronavirus-asymptomatic.html). 69. [69]. Zebin Zhao, Xin Li, Feng Liu, Gaofeng Zhu, Chunfeng Ma, and Liangxu Wang. Prediction of the covid-19 spread in african countries and implications for prevention and control: A case study in south africa, egypt, algeria, nigeria, senegal and kenya. Science of the Total Environment, 729:138959, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.scitotenv.2020.138959&link_type=DOI) 70. [70]. Haonan Zhong and Wendi Wang. Mathematical analysis for covid-19 resurgence in the contaminated environment. Mathematical Biosciences and Engineering, 17(6):6909–6927, 2020. 71. [71]. Lirong Zou, Feng Ruan, Mingxing Huang, Lijun Liang, Huitao Huang, Zhongsi Hong, Jianxiang Yu, Min Kang, Yingchao Song, Jinyu Xia, et al. Sars-cov-2 viral load in upper respiratory specimens of infected patients. New England Journal of Medicine, 382(12):1177–1179, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2001737&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F26%2F2021.10.23.21265421.atom) [1]: /embed/graphic-2.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 [7]: /embed/graphic-8.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-12.gif [10]: /embed/graphic-13.gif [11]: /embed/graphic-14.gif [12]: /embed/graphic-15.gif [13]: /embed/graphic-16.gif [14]: /embed/graphic-17.gif [15]: /embed/graphic-18.gif [16]: /embed/graphic-19.gif [17]: /embed/graphic-20.gif [18]: /embed/inline-graphic-1.gif [19]: /embed/graphic-21.gif [20]: /embed/inline-graphic-2.gif [21]: /embed/graphic-22.gif [22]: /embed/graphic-23.gif [23]: /embed/graphic-26.gif [24]: /embed/graphic-27.gif [25]: /embed/graphic-28.gif [26]: /embed/graphic-29.gif [27]: /embed/graphic-30.gif [28]: /embed/graphic-31.gif [29]: /embed/graphic-32.gif [30]: /embed/graphic-33.gif [31]: /embed/graphic-34.gif [32]: /embed/graphic-35.gif [33]: /embed/graphic-36.gif [34]: /embed/graphic-37.gif [35]: /embed/graphic-38.gif [36]: /embed/graphic-39.gif [37]: /embed/graphic-40.gif [38]: /embed/graphic-41.gif [39]: /embed/graphic-42.gif [40]: /embed/graphic-43.gif [41]: /embed/graphic-44.gif