Table of Contents Author Guidelines Submit a Manuscript
Abstract and Applied Analysis
Volume 2018, Article ID 8479638, 21 pages
Research Article

Global Dynamics of an SVEIR Model with Age-Dependent Vaccination, Infection, and Latency

1School of Mathematical and Statistical Sciences, North-West University, Mafikeng campus, Private Bag X2046, Mmabatho 2735, South Africa
2National Advanced School of Engineering of Maroua, University of Maroua, P.O. Box 46, Maroua, Cameroon

Correspondence should be addressed to Rodrigue Yves M’pika Massoukou;

Received 9 April 2018; Revised 26 June 2018; Accepted 18 July 2018; Published 14 August 2018

Academic Editor: Jaume Giné

Copyright © 2018 Rodrigue Yves M’pika Massoukou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Vaccine-induced protection is substantial to control, prevent, and reduce the spread of infectious diseases and to get rid of infectious diseases. In this paper, we propose an SVEIR epidemic model with age-dependent vaccination, latency, and infection. The model also considers that the waning vaccine-induced immunity depends on vaccination age and the vaccinated individuals fall back to the susceptible class after losing immunity. The model is a coupled system of (hyperbolic) partial differential equations with ordinary differential equations. The global dynamics of the model is established through construction of appropriate Lyapunov functionals and application of Lasalle’s invariance principle. As a result, the global stability of the infection-free equilibrium and endemic equilibrium is obtained and is fully determined by the basic reproduction number .

1. Introduction

Protection induced by vaccines plays a significant role in preventing and reducing the transmission of infectious diseases. One of the greatest successful events of vaccination is illustrated through the eradication of small-pox. It is reported in [1] that the case of small-pox was last recorded in 1977. Immunity conveyed by vaccination depends on different vaccines and vaccination policies—lifelong immunity occurs for certain vaccines while immunization period is induced by some vaccines. However, waning vaccine-induced immunity takes place (naturally) after immunization process. It is reported in [2] that a significant decay in the proportion of chicken pox took place in US in 1995 after conducting a universal vaccination campaign. But, surprisingly new cases of chicken pox appeared mainly in highly vaccinated school communities in US. Some studies were conducted and revealed waning vaccine-induced immunity in children under protection induced by vaccine. Moreover, this was also investigated in [35] and it was proved that such waning of immunity is attached to the time since vaccination and the age at vaccination. In this regard, it was published in [6, 7] that the time of vaccine-induced immunity depends on individuals features and vaccine age.

From the abovementioned statements and citations, we see that it is necessary to associate waning vaccine-induced immunity with vaccination in infectious disease models and it is interesting to investigate the impact of waning vaccine-induced immunity on the dynamics of infectious diseases. Many mathematical models on vaccination have been already investigated and, to cite a few of them, see [2, 816]. Some of the above-cited model considered either age-dependent vaccination, while some did not.

Despite vaccination age-structure being the main and appropriate feature required in the dynamics of infectious diseases with waning induced-vaccine immunity, most of epidemiological models with vaccination including waning induced-vaccine immunity were studied after assuming a constant rate of immunity loss (to name a few, see [9, 10, 15]). Age-dependent vaccination was considered in some epidemiological models studied recently in [1, 2, 1720]. However, some of these works considered either waning vaccine-induced immunity or not, either vaccine-age-dependent waning vaccine-induced immunity or not, either age-dependent latency or not, either age-dependent relapse or not, and either age-dependent infection or not. In [17] an SVIR epidemic model with continuous age-dependent vaccination was formulated to establish the global stability of equilibria. In [18] an SVIJS epidemic model with age-dependent vaccination was considered to study the asymptotical behavior of the equilibria, after assuming that age-dependent vaccine-induced immunity decays with time after vaccination. In [19] an SVIS epidemic model with age-dependent vaccination and vaccine-age-dependent waning vaccine-induced immunity and treatment was formulated to investigate backward bifurcations. In [2] an SVIR epidemic model with vaccination age was considered to establish global stability of equilibria, after assuming that vaccine-induced immunity decays with time after vaccination. In [20] an SVEIR epidemic model with age-dependence vaccination, latency, and relapse was formulated to establish the global stability of the equilibria. In [1] a multigroup SVEIR epidemic model with latent class and vaccination age was formulated to study global stability of equilibria, after assuming that vaccine-induced immunity decays with time after vaccination and likewise in [1719].

Recently, in [14] an SVEIR epidemic model with continuous age-structure in the latent and infectious classes and without continuous age-structure in the vaccinated class was formulated to prove the global stability of equilibria, while in [20] an SVEIR epidemic model with continuous age-structure in the latent, infectious, and recovered classes and with vaccine-age-dependent waning vaccine-induced immunity was formulated. Moreover, in [14] the latency and infection ages are denoted by the same variable . Similarly, in [20] the latency, relapse, and vaccination ages are denoted by the same variable . In spite of this, to the best of our knowledge, the global dynamics of an SVEIR epidemic model with continuous age-structure in latency, infection, vaccination, and vaccine-age-dependence waning vaccine-induced immunity has not yet been neither considered nor investigated using the approach of Lyapunov functionals. The aim of this work is to fill this gap by investigating the global dynamics of an SVEIR epidemic model as above defined. Motivated by [14, 20], we propose a new SVEIR epidemic model originated from an existing SVEIR formulated in [1], by considering continuous age-structure in latency and infection in addition to age-dependence vaccination and vaccine-age-dependence waning vaccine-induced immunity (which the authors took into account in [1]). However, the latency, infection, and vaccination ages are all denoted by , as in [14, 20]. Moreover, in this paper, we also consider a more significant incidence rate (taking into account transmission by both age-mates infective individuals and infective individuals of any age) of the form where and are defined below, instead of the classical incidence rate of the form where denotes the coefficient of disease transmission from infective individual, with infection age , to susceptible individual. The latter is considered in the references therein where continuous age-structure in infection is taken into account.

The model splits the total population into five epidemiological compartments, namely, the susceptible compartment, the vaccinated compartment, the latent compartment, infected compartment, and the removed compartment. Let and be the number of individuals in the susceptible and removed compartments at time , respectively. Let , , and be the density of vaccinated, (latently) infected, and (actively) infected individuals with vaccination, latency, and infection age at time , respectively. It follows that , , and defined by are the number of individuals in the vaccinated, latent, and infected compartments, respectively.

The model to be investigated consists of a hybrid system of nonlinear coupled ordinary differential equations and partial differential equations of the formwhere , , and , with boundary conditionsAnd initial conditionswhere and are initial size of susceptible and removed individuals, respectively, and , , and are initial age-density of vaccinated, latent, and infective individuals, respectively. Moreover, , , and are Lebesgue integrable functions, and it is assumed that the recruitment of newly vaccinated individuals in the vaccinated compartment is done at age zero.

The meaning of parameters in (4)-(5) is given below:: constant recruitment rate of susceptible individuals: rate of vaccination of susceptible individuals: natural mortality rate of individuals: age-specific rate of waning vaccine-induced immunity: age-specific natural mortality rate: age-specific removal rate: age-specific rate moving from latent to infective: age-specific infection rate of susceptible individuals by infective individuals (of the same age–intracohort contagion): probability that an infective individual of age will successfully infect a susceptible individual of age , after contact

In the sequel, we make following assumptions on parameters in (4)–(5):

A1(i) with .(ii) and with essential upper bounds and , respectively.(iii) are Lipschitz continuous on with coefficients , and , respectively.(iv)There exists such that .

A2. is a decreasing function of for any constant removal rate such that .

A3. , for every .

This paper is composed sections, in addition to the introduction, which are structured as follows. In Section 2, we present some preliminary result on compactness property of the semiflow generated by (4)-(5) and we discuss its asymptotic smoothness property. The uniform persistence property of (4)-(5) is addressed in Section 4. Section 3 deals with the existence of equilibria and the formulation of the threshold parameter (the basic reproduction number). Local stability of equilibria for (4) is established in Section 5, while global stability of equilibria for (4) is examined in Section 6.

2. Preliminaries

We consider the Banach space endowed with the norm where , for any . Let us denote by the positive cone of the Banach space such that For any initial value satisfying the conditions system (4) is well-posed, under assumption , due to [21]. Thus, a continuous semiflow is obtained and it is defined by (4) such thatNow, we introduce the functions By integrating the second, third, and fourth equations of (4) along the characteristic line , respectively, we getwhereTaking the norm of and using the positiveness of the components of , we getDifferentiating (15) with respect to leads toNext, we seek for the estimates of each time-derivative on the right hand side (16). First, we haveApplying the Leibniz Integral Rule to the first integral in (17) yieldsLikewise, we also haveandTherefore, we getUsing (iv) of , (21) yieldsthat is, Thus, we obtain where .

If we consider the state space of (4), defined bywe get for any , whenever . Moreover for any .

Then, we state the following result.

Lemma 1. The set is positively invariant for ; that is, Moreover, is point dissipative and the set attracts all points in .

As we aim to make use of the Lasalle’s Invariance Principle, we are required to establish the relative compactness of the orbit in due to the infinite dimensional Banach space . For this, we consider two operators and , (), such thatwhereThus, , for any ; and from the proof of [22, Proposition 3.13] and Lemma 1, we get to the following result.

Theorem 2. For , the orbit has a compact closure in if the following conditions are satisfied:(i)There exists a function such that, for any , , and if with , then for any .(ii)For any , maps any bounded sets of into a set with compact closure in .

For verifying conditions (i) and (ii) of Theorem 2, we need lemmas.

Lemma 3. For , let . Then . Then (i) of Theorem 2 holds.

Proof. Clearly, we see that . Making use of some equations in (13), we getTaking the initial condition such that and , we have

Lemma 4. For , maps any bounded set of into a set with a compact closure in .

Proof. Since and remain in the compact set by Lemma 1, it is sufficient to show that , , and remain in a precompact subset of , which does not depend on the initial data . To achieve this, the following conditions (see [23, Theorem B.2]) must be satisfied for , , and :(i)The supremum of with respect to is finite;(ii) uniformly with respect to ;(iii) uniformly with respect to ;(iv) uniformly with respect to , where . It follows from (13) and (30) thatand hence, using Lemma 1, we getand hence, (i), (ii), and (iv) follow. To establish (iii), we take a sufficiently small such that and we show thatIndeed, is a nondecreasing function of such that satisfying On the other hand, the Lipschitz continuity of is obtained from the first equation of (4), using the boundedness of the solution of (4) (see Lemma 1); that is, there exists such that for any .
Since does not depend on the initial data and as , it follows from (35) that (iii) is satisfied.

Therefore, From Lemma 1 and Theorem 2, the existence result of global attractors (see [24]) follows.

Theorem 5. The semiflow has a global attractor in , which attracts any bounded subset of .

3. Equilibria and Basic Reproduction Number

System (4) has a unique disease-free equilibrium , where

Apart from , system (4) could also have an endemic equilibrium. We suppose that there exists an endemic equilibrium for system (4) denoted by . Therefore, are satisfied. In addition, also satisfies (5), i.e.,The second equation of (38) and the first equation of (39) giveIt follows from the third equation of (38) thatEquations (5) and (41) together with the fourth equation of (38) yieldWe introduce parameters , , , and such thatBy substituting (42) into the second equation of (39), we haveAnd hence,Then, substituting (37), (40), (42), and (44) into the first equation of (38) yieldsTherefore, it easily follows that

A threshold condition is derived from the existence condition for the endemic equilibrium such that . Thus, the parameter , given by the first equation of (43), can be called the basic reproduction number of system (4). Moreover, can also be expressed aswhere

and can be understood, respectively, as the basic reproduction numbers for the corresponding model with purely intracohort infection mechanism (i.e., situation in which individuals can only be infected by their age-mates) and for the corresponding model with purely intercohort infection mechanism (i.e., situation in which individuals can be infected by those of any age).

Therefore, we state the following.

Theorem 6. If , then system (4) has only a disease-free equilibrium , while if , then system (4) also has an endemic equilibrium in addition to the disease-free equilibrium .

4. Uniform Persistence

This section is devoted to the uniform persistence of system (4) under the condition . For this, we introduce a function defined by where . Furthermore, we consider the set defined by such that as whenever .

Definition 7 (see [23, p. 61]). System (4) is uniformly weakly -persistent (respectively, uniformly strongly -persistent) if there exists a positive , independent of initial conditions, such that for .

Theorem 8. If , then system (4) is uniformly weakly -persistent.

Proof. We assume that, for any , we can find such that Since , then we can find a small enough such thatIn particular, we can find such that We can assume that, for any , . It follows from the first equation of (4) that We apply the Laplace transform to the above inequality and we obtain It follows that And hence, This yields . We can assume that, for any , Now, we consider the boundary condition defined by the second equation of (5) and we get We apply the Laplace Transform to the above inequality so that Dividing the above inequality by yields First, taking the limit inferior as on both sides of the above inequality, we obtain Next, take the limit as on both sides of the above inequality, we getwhich contradicts the inequality given in (54).

Combining the results from Theorems 5 and 8 with [25, Theorem 3.2] leads to the uniform (strong) -persistence as follows.

Theorem 9. If , then the semiflow is uniformly (strongly) -persistent.

Definition 10. A total trajectory of a continuous semiflow , defined by (11), is a function such that for any .

Note that a global attractor will only contain points with total trajectories through them as it needs to be invariant. So, the -limit point of a total trajectory , passing through , is given by

A total trajectory satisfies

Corollary 11. Let and be, respectively, a global attractor of in and a total trajectory of in . If , then there exists such that

Proof. We consider the boundary condition given by the second equation of (5). Using (25) and (ii) of , we getFrom the first equation of (4), we have that is, This yieldsTaking the limit inferior as in (72) leads to Therefore,