Research Article | Open Access

Lin Wang, "Existence of Periodic Solutions of Seasonally Forced SEIR Models with Pulse Vaccination", *Discrete Dynamics in Nature and Society*, vol. 2020, Article ID 9381375, 11 pages, 2020. https://doi.org/10.1155/2020/9381375

# Existence of Periodic Solutions of Seasonally Forced SEIR Models with Pulse Vaccination

**Academic Editor:**Zuonong Zhu

#### Abstract

In this paper, we are interested in finding the periodic oscillation of seasonally forced SEIR models with pulse vaccination. Many infectious diseases show seasonal patterns of incidence. Pulse vaccination strategy is an effective tool to control the spread of these infectious diseases. Assuming that the seasonally dependent transmission rate is a T-periodic forcing, we obtain the existence of positive T-periodic solutions of seasonally forced SEIR models with pulse vaccination by Mawhin’s coincidence degree method. Some relevant numerical simulations are presented to illustrate the effectiveness of such pulse vaccination strategy.

#### 1. Introduction

It is a common phenomenon that the incidence of many infectious diseases often changes periodically with the seasonal cycle, such as measles, chickenpox, mumps, rubella, pertussis, and influenza [1–3]. In order to understand the mechanisms responsible for seasonal disease incidence and the epidemiological consequences of seasonality, a large number of mathematical models of infectious diseases with periodic transmission rates have been established [4–7]. Dietz [8] was the first to investigate the effects of one-year periodic contact rate in the classical SIR and SEIR models. Dietz considered a periodical contact rate given by

The periodically forced nonlinear effects in epidemic models have been studied extensively in the mathematical literature [9, 10].

Pulse vaccination strategy (PVS) is an effective tool to control the spread of epidemics, for example, control of poliomyelitis and measles in Central [11] and South American [12] and the UK vaccination campaigns against measles in 1994 [13]. The theoretical study on pulse vaccination strategy was firstly presented by Agur et al. [14]. Shulgin et al. [15, 16] incorporated pulse vaccination into the SIR epidemic model. Nokes and Swinton studied the control of childhood viral infections by pulse vaccination [17]. d’Onofrio applied the pulse vaccination method for SIR and SEIR epidemic models [18, 19]. PVS consists of periodical repetitions of impulsive multiage cohort vaccinations in a population [18, 19]. PVS proposes to vaccinate a constant fraction of the entire susceptible people in a single pulse, which can be formulated aswhere is left continuous satisfying . Pulse vaccination gives life-long immunity to susceptibles who are transferred to the “recovered” class of the population, which can be formulated as

This kind of vaccination is called impulsive since all the vaccine doses are applied in a time which is short considering the dynamics of the disease. PVS has been further developed, for example, in [20–23]. A comprehensive introduction on vaccination strategies can be found in [24].

Many infectious diseases do not die out, but become endemic. For autonomous epidemic models, the existence of the positive equilibrium plays an important role. A positive periodic solution in the periodic model will play the same role as a positive equilibrium in the autonomous model [25, 26]. Recently, Katriel [27] proved that the seasonally forced SIR model with a T-periodic forcing has a periodic solution with periodic *T* by Leray–Schauder degree theory provided . Jódar et al. [26] obtained that a -periodic solution exists for a more general system by the famous Mawhin’s coincidence degree method, whenever the condition holds. Using Leray–Schauder degree theory, Zu and the author [28] established new results on the existence of at least one positive periodic solution for a seasonally forced SIR model with impact of media coverage. The author [29] proved the existence of positive periodic solutions of seasonally forced SIR models with impulse vaccination at fixed time by Mawhin’s coincidence degree method if the basic reproductive number . Coincidence degree theory has been applied to prove the existence of multiple periodic solutions of the epidemic model with seasonal periodic rate [30, 31]. There are some research activities about the existence of periodic solutions of impulsive differential equation [32–36].

The aim of this paper is to study the existence of periodic solution of seasonally forced SEIR models with pulse vaccination. The paper is organized as follows. A seasonally forced SEIR model with pulse vaccination is formulated and a suitable region to our problem is chosen in Section 2. The existence of periodic solutions of our impulsive systems is established in Section 3. Some numerical simulations are demonstrated to verify the effectiveness of our pulse vaccination strategy in Section 4. The relevant conclusion will be stated in Section 5.

#### 2. The Model and Prelimimaries

##### 2.1. The Seasonally Forced SEIR Model with Pulse Vaccination

In this paper, we focus on the existence of periodic solution of seasonally forced SEIR models with pulse vaccination; we consider models of the formmodeling the spread of infectious diseases with PVS under the following hypotheses:(i)The population is divided into four classes; , , , and denote, respectively, the fractions of the susceptible, exposed/latent, infective, and recovered population, and is invariant with , for all .(ii), , and denote the birth (death) rate, the rate of latent individuals becoming infectious, and recovery rate, respectively, which are positive constants.(iii) is the seasonally dependent transmission rate, which is a positive continuous -periodic function.(iv)The susceptible population can be divided into many groups and all the groups cannot be vaccinated at the same time, and the susceptible population will be vaccinated for several times:

Our vaccination strategies concern the impact of infected population, which can be formulated aswhere and the sensitivity coefficient is sufficiently large.

Denote the basic reproduction number:

Obviously, we have . Since , , and are fractions of the population, we have for all . Because does not appear in the first three equations in (4), system (4) reduces to the following 3-dimensional system:with

Obviously, finding the periodic solution of (8) is equivalent to finding the solutions of the following periodic boundary value problem:with periodic boundary condition

##### 2.2. The Suitable Region to Our Problem

In order to prove the existence of periodic solutions of (8), we consider the following auxiliary problem:where and .

Let be an open bounded subset of (will be denoted in Section 2.4) satisfying

Proposition 1. * is an invariant region with respect to (12). The disease-free equilibrium is the unique periodic solution of (12) satisfying , .*

*Proof. *First, we will prove that is an invariant region.

In fact, it follows from model (12) thatSince there is no impulsive motion for , andit is easy to conclude that every possible solution will remain in the region ultimately.

Second, we will prove that the disease-free equilibrium is the unique periodic solution of (12) satisfying .

We assume that is a solution of (12); this means at least one of the following conditions holds:(i)There exists so that (ii)There exists so that (iii)There exists so that (iv)There exists so that We now consider each of these four cases:

In case of (i), we have and . If , and . Thus, the only possible periodic solution of is . On the other hand, if , . Thus, it is easy to obtain that for sufficiently close to , which contradicts that is an invariant region.

In case of (ii), we have and . Thus, it is easy to obtain that for sufficiently close to , which contradicts that is an invariant region.

In case of (iii), we have and . Because the cases and have been discussed above, we have which again contradicts that is an invariant region.

In case of (iv), we getBecause has been discussed, we only discuss , , which contradicts that is an invariant region.

*Remark 1. *If the impulsive motion is not influenced by , which means , the system , can have a nonconstant periodic solution on , which is hard to handle. In fact, if there are no infectious patients, it is often meaningless to vaccinate the susceptible people.

In order to use continuity method, it is necessary to choose an open bounded set , such that there is no solution of (12) satisfying for any . Motivated by the idea of Katriel [27], we choose to be the open subset of given bywhere is to be fixed.

Proposition 2. *Let , . Then, there exists no solution of (12) satisfying , for any .*

*Proof. *Suppose . Then, either or .

In the first case, and Proposition 1 imply that there is no solution of (12) on .

In the second case, we have and , . Since there is no impulsive motion for and , after integrating the third equation of (12) over , we obtainSince there is no impulsive motion for and , after integrating the second equation of (12) over , we obtainWith the help of (17)–(19), we conclude thatFor , we get immediatelywhich is a contradiction to the assumption .

##### 2.3. Outline of Mawhin’s Coincidence Degree Method

We introduce a few definitions and recall the continuation theorem which will help us to prove the existence of positive solutions of system (10).

Consider the operator equation:where is a linear bounded operator, is a continuous operator, and and are Banach spaces.

*Definition 1. *(see [37]). The linear mapping *L* is called a Fredholm mapping of index zero if(1)(2)*is closed in*If is a Fredholm mapping of index zero, there exist continuous projectors and such that(1)(2)(3)(4)It follows that is invertible. We denote the inverse of that map by .

*Definition 2. *(see [37]). Let be an open bounded subset of . The mapping is called -compact on , if(1)*is bounded*(2)*is compact*Since is isomorphic to , there exists an isomorphism .

Theorem 1 (see [37]) (Mawhin’s continuation theorem). *Let be an open bounded set. Let be a Fredholm mapping of index zero and be -compact on . Assume that*(1)* For each, , *(2)

*(3)*

*For each*,*Then, the equation has at least one solution in .*

Lemma 4 (Arzela–Ascoli) (see [38]). *Let be a compact subset of and be the linear space of continuous functions which take D into ; any uniformly bounded equicontinuous sequence of functions in has a subsequence which converges uniformly on .*

##### 2.4. Notations

We now prepare the setting to apply Mawhin’s continuation theorem. Define exists for and , exist with , where *j* is a non-negative integer and , ,

Define a Banach spacewith the norm

Define another Banach spacewith the norm

Letwhere

Letwhere

Obviously, system (10) can be written by . By a simple calculation, we get

It is easy to see that

Since is closed, is Fredholm mapping of index 0. Let be the projector given by

Obviously,

Let be the projector given by

Obviously,

Furthermore, the generalized inverse (to ) exists given by

Then, has the following form:

By a direct calculation, we have

#### 3. Results

The following theorem gives the main results of this paper.

Theorem 2. *Let ; there exists at least one -periodic solution of (1), all of whose components are positive.*

*Remark 2. *When , , system (4) is the usual seasonally forced SEIR model without pulse vaccination.

The proof of Theorem 2 will be divided into 4 steps.

*Proof. * Step 1: is -compact on . First, it is obvious that is bounded. For any , Then, it is easy to check that is uniform bounded and equicontinuous on each . Assume that . Using Lemma 4, there exists a uniformly convergent subsequence denoted by on . Using Arzela–Ascoli lemma again on , we have a uniformly convergent subsequence which is also uniformly convergent on . Repeat it again and again, and we can prove that is uniformly convergent on . In this way, is compact. Step 2: For each , , , which has been proved by Proposition 2. Step 3: For each , . If , we have Assume ; we know that is a constant vector in . Thus, (40) is equivalent to We claim that there are exactly two solutions: and in satisfying According to the definition of , we know that . Denote We have Since is a monotonous function with respect to . Thus, is the unique solution of in . Since we have . Consequently, for each , . Step 4. There exists an isomorphism such thatWe will prove that . On account of the discussion in step 3, we know that is the unique solution of in .

A direct calculation shows thatThus, we conclude from Theorem 1 that the equation has at least one solution on .

From Proposition 2, we deduce that for each , , . It follows that has at least one solution in .

#### 4. Simulation

In this section, some relevant numerical simulations about the T-periodic solution of the seasonally forced SEIR models are presented to show the effectiveness of PVS. Furthermore, we will compare the effects of different parameters on the solutions of the seasonally forced SEIR models with PVS.

With the period of the forcing representing one year, we take corresponding to a 2-week infectious period. We set . Assume that there are three impulsive points at fixed time , , and with , . Let be divided into intervals equally. Given the initial point , which is the endemic equilibrium of the SEIR model without periodic transmission rate and pulse vaccination. The periodic solutions of system (10) can be solved by the Newton iteration method in which we set at fixed time , , and .

*Simulation 1. *Set , , ; we get the approximate susceptible population, exposed population, infective population, and recovered population of system (10) by Newton iteration. In Figure 1, we compare the T-periodic solution of the seasonally forced SEIR models with PVS in 1 year and 10 years. Figure 1 shows periodicity of the positive solution of the seasonally forced SEIR models with PVS.

**(a)**

**(b)**

*Simulation 2. *Set , , and ; we make 8 steps of Newton iteration to get the approximate infective population of system (10) with both (the surface at the bottom) and (the surface at the top). Obviously, the infective population of system (10) with pulse is lower than the infective population of system (10) without pulse in Figure 2. The susceptible population of system (10) with has periodic and impulsive properties. Thus, it is very effective to lower the infective population by PVS. Furthermore, Figure 2 shows the stability of the periodic solution by Newton iteration.

**(a)**

**(b)**

*Simulation 3. *Set , , and ; we simulate system (10) with , , and by Newton iteration. Obviously, it is very effective to lower the exposed population and the infective population by PVS in Figure 3.

**(a)**

**(b)**

*Simulation 4. *Set , , and ; we simulate system (10) with , , and by Newton iteration. Figure 4 shows the impact on exposed population, infective population, recovered population, and susceptible population by different transmission rates. As increases, the infective population increases while susceptible population decreases.

**(a)**

**(b)**

*Simulation 5. *Set , , and ; we simulate system (10) with , , and by Newton iteration. Figure 5 shows the impact on exposed population, infective population, recovered population, and susceptible population by different birth rates. As increases, the infective population and susceptible population increase, respectively.

**(a)**

**(b)**

*Simulation 6. *Set , , and ; we simulate system (10) with , , and by Newton iteration. Figure 6 shows the impact on systems by the different rates of latent individuals becoming infectious. As increases, the infective population increases while susceptible population decreases.

**(a)**

**(b)**

#### 5. Conclusion

We obtain the existence of positive T-periodic solutions of seasonally forced SEIR models with pulse vaccination by Mawhin’s coincidence degree method. Some relevant numerical simulations are presented to show the T-periodic solution of the seasonally forced epidemiological models and to illustrate the effectiveness of such pulse vaccination strategy.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The author declares that there are no conflicts of interest.

#### Acknowledgments

This research was supported by NSFC (grant nos. 11901052 and 11626043) and the fund of Jilin Provincial Education Department (grant no. JJKH20170535KJ).

#### References

- Z. Grossman, I. Gumowski, and K. Dietz, “The incidence of infectious diseases under the influence of seasonal fluctuations analytical approach,” in
*Nonlinear Systems and Applications*, pp. 525–546, Academic Press, New York, NY, USA, 1977. View at: Google Scholar - S. Altizer, A. Dobson, P. Hosseini, P. Hudson, M. Pascual, and P. Rohani, “Seasonality and the dynamics of infectious diseases,”
*Ecology Letters*, vol. 9, no. 4, pp. 467–484, 2006. View at: Publisher Site | Google Scholar - N. C. Grassly and C. Fraser, “Seasonal infectious disease epidemiology,”
*Proceedings of the Royal Society B: Biological Sciences*, vol. 273, no. 1600, pp. 2541–2550, 2006. View at: Publisher Site | Google Scholar - H. W. Hetcote, “Qualitative analyses of communicable disease models,”
*Mathematical Biosciences*, vol. 28, p. 335, 1976. View at: Google Scholar - Y. Lou, X. Q. Zhao, and X.-Q. Zhao, “Threshold dynamics in a time-delayed periodic SIS epidemic model,”
*Discrete & Continuous Dynamical Systems-B*, vol. 12, no. 1, pp. 169–186, 2009. View at: Publisher Site | Google Scholar - Y. Lou and X.-Q. Zhao, “A climate-based malaria transmission model with structured vector population,”
*SIAM Journal on Applied Mathematics*, vol. 70, no. 6, pp. 2023–2044, 2010. View at: Publisher Site | Google Scholar - J. Ma, Z. E. Ma, and Z. Ma, “Epidemic threshold conditions for seasonally forced SEIR models,”
*Mathematical Biosciences and Engineering*, vol. 3, no. 1, pp. 161–172, 2006. View at: Publisher Site | Google Scholar - K. Dietz, “The incidence of infectious diseases under the influence of seasonal fluctuations,”
*Mathematical Models in Medicine*, Springer, Berlin, Germany, 1976. View at: Google Scholar - H. W. Hethcote and S. A. Levin, “Periodicity in epidemiological models,” in
*Applied Mathematical Ecology*, pp. 193–211, Springer, Trieste, Italy, 1986. View at: Google Scholar - B. Buonomo, N. Chitnis, and A. d’Onofrio, “Seasonality in epidemic models: a literature review,”
*Ricerche di Matematica*, vol. 67, no. 1, pp. 7–25, 2018. View at: Publisher Site | Google Scholar - C. A. de Quadros, J. K. Andrus, J.-M. Olivé et al., “Eradication of poliomyelitis,”
*The Pediatric Infectious Disease Journal*, vol. 10, no. 3, pp. 222–229, 1991. View at: Publisher Site | Google Scholar - A. B. Sabin, “Measles killer of millions in developing countries strategies of elimination and continuing control,”
*European Journal of Epidemiology*, vol. 7, p. 1, 1991. View at: Publisher Site | Google Scholar - M. Ramsay, N. Gay, E. Miller et al., “The epidemiology of measles in England and Wales Rationale for 1994 national vaccination campaign,”
*Communicable Disease Report*, vol. 4, no. 12, pp. R141–R146, 1994. View at: Google Scholar - Z. Agur, L. Cojocaru, G. Mazor, R. M. Anderson, and Y. L. Danon, “Pulse mass measles vaccination across age cohorts,”
*Proceedings of the National Academy of Sciences*, vol. 90, no. 24, pp. 11698–11702, 1993. View at: Publisher Site | Google Scholar - B. Shulgin, L. Stone, and Z. Agur, “Pulse vaccination strategy in the SIR epidemic model,”
*Bulletin of Mathematical Biology*, vol. 60, no. 6, pp. 1123–1148, 1998. View at: Publisher Site | Google Scholar - L. Stone, B. Shulgin, and Z. Agur, “Theoretical examination of the pulse vaccination policy in the SIR epidemic model,”
*Mathematical and Computer Modelling*, vol. 31, no. 4-5, pp. 207–215, 2000. View at: Publisher Site | Google Scholar - D. J. Nokes and J. Swinton, “The control of childhood viral infections by pulse vaccination,”
*Mathematical Medicine and Biology*, vol. 12, no. 1, pp. 29–53, 1995. View at: Publisher Site | Google Scholar - A. d’Onofrio, “Stability properties of pulse vaccination strategy in SEIR epidemic model,”
*Mathematical Biosciences*, vol. 179, no. 1, pp. 57–72, 2002. View at: Google Scholar - A. d’Onofrio, “Pulse vaccination strategy in the SIR epidemic model: global asymptotic stable eradication in presence of vaccine failures,”
*Mathematical and Computer Modelling*, vol. 36, no. 4-5, pp. 473–489, 2002. View at: Google Scholar - Y. Zhou and H. Liu, “Stability of periodic solutions for an SIS model with pulse vaccination,”
*Mathematical and Computer Modelling*, vol. 38, no. 3-4, pp. 299–308, 2003. View at: Publisher Site | Google Scholar - J. Hui, L. S. Chen, and L. Chen, “Impulsive vaccination of SIR epidemic models with nonlinear incidence rates,”
*Discrete & Continuous Dynamical Systems-B*, vol. 4, no. 3, pp. 595–605, 2004. View at: Publisher Site | Google Scholar - M. Khasin, M. I. Dykman, and B. Meerson, “Speeding up disease extinction with a limited amount of vaccine,”
*Physical Review E*, vol. 81, p. 11, 2010. View at: Publisher Site | Google Scholar - R. Shi, X. Jiang, and L. Chen, “The effect of impulsive vaccination on an SIR epidemic model,”
*Applied Mathematics and Computation*, vol. 212, no. 2, pp. 305–311, 2009. View at: Publisher Site | Google Scholar - M. J. Keeling and P. Rohani,
*Modeling Infectious Diseases in Humans and Animals*, Princeton University Press, Princeton, NJ, USA, 2007. - T. Zhang, J. Liu, and Z. Teng, “Existence of positive periodic solutions of an SEIR model with periodic coefficients,”
*Applications of Mathematics*, vol. 57, no. 6, pp. 601–616, 2012. View at: Publisher Site | Google Scholar - L. Jódar, R. J. Villanueva, and A. Arenas, “Modeling the spread of seasonal epidemiological diseases: theory and applications,”
*Mathematical and Computer Modelling*, vol. 48, no. 3-4, pp. 548–557, 2008. View at: Publisher Site | Google Scholar - G. Katriel, “Existence of periodic solutions for the periodically forced SIR model,”
*Nelı̄nı̄ĭnı̄ Koliv*, vol. 16, pp. 359–366, 2013. View at: Google Scholar - J. Zu and L. Wang, “Periodic solutions for a seasonally forced SIR model with impact of media coverage,”
*Advances in Difference Equations*, vol. 136, p. 10, 2015. View at: Google Scholar - L. Wang, “Existence of periodic solutions of seasonally forced SIR models with impulse vaccination,”
*Taiwanese Journal of Mathematics*, vol. 19, no. 6, pp. 1713–1729, 2015. View at: Publisher Site | Google Scholar - Z. Bai and Y. Zhou, “Existence of two periodic solutions for a non-autonomous SIR epidemic model,”
*Applied Mathematical Modelling*, vol. 35, no. 1, pp. 382–391, 2011. View at: Publisher Site | Google Scholar - Z. Bai, Y. Zhou, and T. Zhang, “Existence of multiple periodic solutions for an SIR model with seasonality,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 74, no. 11, pp. 3548–3555, 2011. View at: Publisher Site | Google Scholar - D. D. Baĭnov, S. G. Hristova, S. C. Hu, and V. Lakshmikantham, “Periodic boundary value problems for systems of first order impulsive differential equations,”
*Differential Integral Equations*, vol. 2, pp. 37–43, 1989. View at: Google Scholar - Y. Dong and E. Zhou, “An application of coincidence degree continuation theorem in existence of solutions of impulsive differential equations,”
*Journal of Mathematical Analysis and Applications*, vol. 197, no. 3, pp. 875–889, 1996. View at: Publisher Site | Google Scholar - S. G. Hristova and D. D. Bainov, “Existence of periodic solutions of nonlinear systems of differential equations with impulse effect,”
*Journal of Mathematical Analysis and Applications*, vol. 125, no. 1, pp. 192–202, 1987. View at: Publisher Site |