Next Article in Journal
Experimental Research on the Settlement Feature of Two Ground Deformation Modes Induced by Tunnelling
Next Article in Special Issue
Coexistence and Replacement of Two Different Maturation Strategies Adopted by a Stage-Structured Population
Previous Article in Journal
Optimal Plastic Reliable Design of Reinforced Concrete Beams Considering Steel Bars Volume Probability
Previous Article in Special Issue
Mathematical Model to Understand the Dynamics of Cancer, Prevention Diagnosis and Therapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Periodic Behaviour of an Epidemic in a Seasonal Environment with Vaccination

by
Miled El Hajji
1,2,*,
Dalal M. Alshaikh
1 and
Nada A. Almuallem
1
1
Department of Mathematics, Faculty of Science, University of Jeddah, P.O. Box 80327, Jeddah 21589, Saudi Arabia
2
ENIT-LAMSIN, Tunis El Manar University, BP. 37, Tunis-Belvédère, Tunis 1002, Tunisia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(10), 2350; https://doi.org/10.3390/math11102350
Submission received: 10 April 2023 / Revised: 11 May 2023 / Accepted: 15 May 2023 / Published: 18 May 2023
(This article belongs to the Special Issue Applications of Differential Equations to Mathematical Biology)

Abstract

:
Infectious diseases include all diseases caused by the transmission of a pathogenic agent such as bacteria, viruses, parasites, prions, and fungi. They, therefore, cover a wide spectrum of benign pathologies such as colds or angina but also very serious ones such as AIDS, hepatitis, malaria, or tuberculosis. Many epidemic diseases exhibit seasonal peak periods. Studying the population behaviours due to seasonal environment becomes a necessity for predicting the risk of disease transmission and trying to control it. In this work, we considered a five-dimensional system for a fatal disease in a seasonal environment. We studied, in the first step, the autonomous system by investigating the global stability of the steady states. In a second step, we established the existence, uniqueness, positivity, and boundedness of a periodic orbit. We showed that the global dynamics are determined using the basic reproduction number denoted by R 0 and calculated using the spectral radius of an integral operator. The global stability of the disease-free periodic solution was satisfied if R 0 < 1 , and we show also the persistence of the disease once R 0 > 1 . Finally, we displayed some numerical investigations supporting the theoretical findings, where the trajectories converge to a limit cycle if R 0 > 1 .

1. Introduction

The flu is a seasonal disease. Every year, this viral and contagious disease reappears, and, at the same time, from the end of the year until the beginning of spring. It affects all categories of the population, young and old alike. This pathology has been monitored by the WHO (World Health Organization) since 1952 [1]. This monitoring makes it possible to identify viral strains but also to better understand and anticipate its consequences and its severity worldwide. Although the flu is most often a mild illness, it can lead to serious complications and can even be fatal. In Europe, the epidemic occurs every year, generally between November and April. There are three types of Influenza virus: viruses A and B are the cause of the “real” flu, while type C is responsible for infections similar to a common cold that go unnoticed. Influenza viruses are characterized by their significant genetic variability, with more or less frequent appearance depending on the virus of mutations [2]. Influenza viruses are easily transmitted through contact with respiratory secretions: by the projection of droplets, with direct contact and indirect contact (the virus can survive for up to 5 min on the skin, a few hours in dried secretions, and 48 h on objects). The incubation of the flu (period without symptoms during which the carrier subject is contagious) is about 3 days [2]. Every year, millions of people are affected by winter viruses. Influenza, gastroenteritis, and bronchiolitis are the main infections due to winter viruses, and they have a strong impact on healthcare structures during the winter. The most effective way to protect the fragile people around you is vaccination. Every year, a vaccination campaign is launched at the end of the year. The vaccine produced is different each year, in order to adapt to the genetic evolution of the virus. Thus, as the mutation of the virus is constant, it is necessary to be vaccinated every year to be protected and to protect those around you. If your child is less than 6 months old, it is recommended that you, the parents, and those around you get vaccinated to protect them. Indeed, before the age of 6 months, the vaccine is not suitable, but your child is vulnerable to the flu [2]. The seasonal flu vaccine is offered every year, especially to the elderly. Its formulation varies from year to year, as the strains in circulation change. The effectiveness of the vaccine varies according to the age of the patient but also to the formulation. In general, vaccination against seasonal flu is not 100% effective. A 2014 Cochrane review analysed 90 clinical trials including 70,000 people. The authors concluded that the efficacy of the vaccine was “modest”. In 2017, a Canadian study found the vaccine’s effectiveness to be 42%. In 2014–2015, it was not as good due to a mutation of the virus [3]. However, in people over 65, this efficiency decreases further, due to the ageing of their immune system. Finnish and Swedish data for 2016/2017 suggested protection in the order of 20–30% for people aged over 65. For the United States, the World Health Organization (WHO), or the European Union (EU), knowing the evolution of a human epidemic (H1N1 Flu, Ebola Virus, Coronavirus), animal (Bird Flu, Swine Fever, Rabies), or plant flu is essential.
The mathematical modelling in epidemiology is a way to study how a disease is spread, predict the future behaviour, and propose control strategies. Several works qualitatively proposed and studied some mathematical models describing the dynamical behaviour of infectious disease transmission (see, for example, [4,5,6,7,8,9,10]). In particular, the SVEIR epidemic models with constant coefficients have been analysed in several works (see, for example, [11,12,13,14,15,16,17]). However, seasonality is very repetitive in each of the ecological, biological, and human systems [18]. In particular, in the climate variation patterns repeated every year by the same way, bird migration is repeated according to the repeated season variation, schools open and close almost periodically each year, etc. Among other things, these seasonal factors affect the pathogens’ survival in the environment, host behaviour, and the abundance of vectors and non-human hosts. Therefore, several diseases show seasonal behaviours. Taking into account the seasonality in mathematical modelling becomes very important. Note that even the simplest mathematical models that take into account seasonality present many difficulties to study [6]. In [19], Bacaër and Gomes discussed the periodic S-I-R model, a simple generalization of the classical model of Kermack and McKendrick [20]. In [21], the authors studied a SEIRS epidemic model with periodic fluctuations. They calculated the basic reproduction number R 0 of the time-averaged system (autonomous). Then, they proved a sufficient but not necessary condition ( R 0 < 1 ) such that the disease could not persist in the population in a seasonal environment. In [22], Guerrero-Flores et al. considered a class of SIQRS models with periodic variations in the contact rate. They proved the existence of periodic orbits by using Leray–Schauder degree theory. Zhang and Teng [23] studied an alternative SEIRS epidemic model in a seasonal environment and established some sufficient equivalent conditions for the persistence and the extinction of the disease. These results were improved by Nakata and Kuniya in [4] by giving a threshold value between the uniform persistence and the extinction of the disease. In [7], Bacaër and Guernaoui gave the definition of the basic reproduction number in seasonal environments. In 2008, Wang and Zhao [24] defined R 0 for several compartmental epidemic models in seasonal environments. All these definitions were different, in several cases, from the basic reproduction number defined for the time-averaged system. By considering general compartmental epidemic models in seasonal environments, Wang and Zhao [24] showed that R 0 was the threshold value for proving or not the local stability of the disease-free periodic trajectory.
As seasonality is very repetitive in the environment, which affects several diseases that show seasonal behaviours, taking seasonality into account in mathematical modelling becomes a necessity. In this paper, we proposed an extension of the SVEIR model proposed in [8] by taking into account the seasonal environment. We studied, in a first step, the autonomous system by investigating the global stability of the steady states. In a second step, we showed that the disease-free periodic solution is globally asymptotically stable if R 0 is less than 1, and, if R 0 is greater than 1, the disease persists. The rest of the paper is organized as follows. In Section 2, we introduced the mathematical model. In Section 3, we studied the case of an autonomous system, where all parameters are supposed to be constants. In Section 4, we considered the non-autonomous system, gave some basic results, and gave the definition of R 0 . We showed that the value of R 0 around one was a threshold value between the disease’s extinction and the disease’s uniform persistence. We gave numerical examples that supported the theoretical findings in Section 5. Section 6 provided a brief conclusions of our obtained results.

2. Mathematical Model and Properties

This study makes it possible to predict the evolution of the disease over time, and its main purpose is to guide leaders in decision-making in terms of public health. It should now be specified that mathematical models are simplified models, with limits that are necessary to be aware about. However, they remain useful in trying to predict the evolution of an epidemic. It would be tempting to mathematically model the spread of seasonal flu. realistically while increasing the complexity of the model. We will content ourselves below with an epidemic model with direct and non-vectorial transmission in a periodic seasonal environment, which can be applied to the case of seasonal flu. In particular, we consider a generalisation of the model given in [8], taking into account the epidemics’ seasonal features, and then we study the behaviour of solutions of the new model designed in Figure 1 and presented by the system of ordinary differential Equation (1).
S ( t ) , V ( t ) , E ( t ) , I ( t ) , and R ( t ) describe susceptible, vaccinated, exposed, infectious, and recovered compartments, respectively (see Table 1 for the epidemiological meaning of the model parameters).
  • All compartments have instantaneous common death rates, d ( t ) .
  • Vaccination against seasonal flu is not 100% effective. Therefore, vaccinated individuals can be infected again.
  • Recovered individuals will not re-infect.
The proposed “SVEIR” mathematical model for the spread of the seasonal flu is:
S ˙ ( t ) = d ( t ) s i n ( t ) ( d ( t ) + p ( t ) ) S ( t ) β ( t ) μ ( I ( t ) ) S ( t ) , V ˙ ( t ) = p ( t ) S ( t ) d ( t ) V ( t ) θ ( t ) β ( t ) μ ( I ( t ) ) V ( t ) , E ˙ ( t ) = β ( t ) μ ( I ( t ) ) ( S ( t ) + θ ( t ) V ( t ) ) ( d ( t ) + δ ( t ) + ε ( t ) ) E ( t ) , I ˙ ( t ) = δ ( t ) E ( t ) ( d ( t ) + κ ( t ) ) I ( t ) , R ˙ ( t ) = ε ( t ) E ( t ) + κ ( t ) I ( t ) d ( t ) R ( t ) ,
with the positive initial condition ( S 0 , V 0 , E 0 , I 0 , R 0 ) R + 5 .
Assumption 1.
  • The function μ is non-negative C 1 ( R + ) , increasing bounded and concave with μ ( 0 ) = 0 .
  • The functions s i n ( · ) , p ( · ) , β ( · ) , θ ( · ) , d ( · ) , δ ( · ) , κ ( · ) , and ε ( · ) are continuous, positive, and T-periodic.
μ satisfies the following results.
Lemma 1.
The function μ satisfies μ ( x ) x μ ( x ) μ ( 0 ) x for all x > 0 .
Proof. 
See [25]. □

3. The Case of Constant Parameters

Consider System (1) for the case where all parameters are constants. The model becomes exactly the one considered in [8] but with a possible recovery of exposed individuals.
S ˙ = d s i n ( d + p ) S β μ ( I ) S , V ˙ = p S d V θ β μ ( I ) V , E ˙ = β μ ( I ) ( S + θ V ) ( d + δ + ε ) E , I ˙ = δ E ( d + κ ) I , R ˙ = ε E + κ I d R ,
with the positive initial condition ( S 0 , V 0 , E 0 , I 0 , R 0 ) R + 5 . It is necessary that the state variables S ( t ) , V ( t ) , E ( t ) , I ( t ) , and R ( t ) remain non-negative for all t 0 . We have the following result.
Proposition 1.
The compact set
Ω = { ( S , V , E , I , R ) R + 5 / S + V + E + I + R = s i n ; S d s i n d + p , V p s i n d + p }
is positively invariant for System (2).
Proof. 
See [8]. □
The basic reproduction number R 0 of System (2) is given [8] as follows:
R 0 = β δ s i n μ ( 0 ) ( d + θ p ) ( d + δ + ε ) ( d + κ ) ( d + p ) .
Proposition 2.
System (2) always has a trivial equilibrium point E 0 = d s i n d + p , p s i n d + p , 0 , 0 , 0 . Moreover, if R 0 > 1 , System (2) has a unique endemic steady state E * = ( S * , V * , E * , I * , R * ) , where I * ( 0 , s i n ) is the unique solution of f ( I ) = 0 , and f is an increasing function given by
f ( I ) = 1 β δ d s i n μ ( I ) I ( d + θ p + θ β μ ( I ) ) ( d + κ ) ( d + δ + ε ) ( d + p + β μ ( I ) ) ( d + θ β μ ( I ) ) .
S * , V * , E * , and R * are given by
S * = d s i n d + p + β μ ( I * ) , V * = p d s i n ( d + θ β μ ( I * ) ) ( d + p + β μ ( I * ) ) , E * = d β s i n μ ( I * ) ( d + θ p + θ β μ ( I * ) ) ( d + δ + ε ) ( d + p + β μ ( I * ) ) ( d + θ β μ ( I * ) ) , R * = β s i n μ ( I * ) ( ε d + ε κ + κ δ ) ( d + θ p + θ β μ ( I * ) ) ( d + κ ) ( d + δ + ε ) ( d + p + β μ ( I * ) ) ( d + θ β μ ( I * ) ) .
Hereafter, we discuss the global stability of the steady states, with respect to the value of R 0 . We started by giving the first main result for the autonomous System (2).
Theorem 1.
The equilibrium point E 0 is globally asymptotically stable once R 0 1 ; however, it is unstable once R 0 > 1 .
Proof. 
Let the function F 0 be given by: F 0 = δ E + ( d + δ + ε ) I . As μ ( I ) μ ( 0 ) I , the time-derivative of F 0 is given by:
F ˙ 0 = δ E ˙ + ( d + δ + ε ) I ˙ = δ ( β μ ( I ) ( S + θ V ) ( d + δ + ε ) E ) + ( d + δ + ε ) ( δ E ( d + κ ) I ) = δ β μ ( I ) ( S + θ V ) ( d + δ + ε ) ( d + κ ) I δ β μ ( 0 ) I ( S + θ V ) ( d + δ + ε ) ( d + κ ) I ( δ β μ ( 0 ) ( S + θ V ) ( d + δ + ε ) ( d + κ ) ) I ( d + δ + ε ) ( d + κ ) ( δ β μ ( 0 ) ( S + θ V ) ( d + δ + ε ) ( d + κ ) 1 ) I ( d + δ + ε ) ( d + κ ) δ β μ ( 0 ) ( d s i n d + p + θ p s i n d + p ) ( d + δ + ε ) ( d + κ ) 1 I = ( d + δ + ε ) ( d + κ ) R 0 1 I , ( S , V , E , I , R ) Ω .
If R 0 1 , then F ˙ 0 0 for all S , V , E , I , R > 0 . Let W 0 = { ( S , V , E , I , R ) : F ˙ 0 = 0 } = { E 0 } . Therefore, by applying the LaSalle’s invariance principle [26], we deduce that the steady state E 0 is globally asymptotically stable when R 0 1 . □
Give the second main result for the autonomous System (2), as, in the first phases of an epidemic, when the number of infections and exposed individuals remains small compared to their number more later, and as the number of susceptible and vaccinated individuals represents the most total size of the population. Thus, we consider the following set
Ω 0 = ( S , V , E , I , R ) Ω / S S * , V V * , E E * , I I * .
Lemma 2.
The function μ satisfies μ ( I ) μ ( I * ) I I * for all 0 I I * .
Proof. 
Let 0 I I * , and let us define the function g ( I ) = μ ( I ) I μ ( I * ) I * then g ( I ) = μ ( I ) I μ ( I ) I 2 0 by Lemma 1 and then g ( I ) g ( I * ) = 0 for all
0 I I * . Note that μ ( I ) μ ( I * ) I I * = I μ ( I * ) μ ( I ) I μ ( I * ) I * = I μ ( I * ) g ( I ) 0 for all 0 I I * . Therefore, μ ( I ) μ ( I * ) I I * for all 0 I I * . □
Theorem 2.
If R 0 > 1 , then the steady state E * is globally asymptotically stable in Ω 0 .
Proof. 
Define the Lyapunov function F * in Ω 0 as follows
F * = S S * ln ( S S * ) + V V * ln ( V V * ) + E E * ln ( E E * ) + d + δ + ε δ I I * ln ( I I * ) .
The function F * admits its minimum value F 1 m i n = S * + V * + E * + d + δ + ε δ I * once S = S * , V = V * , E = E * , I = I * . Let us calculate the derivative of F * along trajectories of System (2).
F ˙ * = ( 1 S * S ) S ˙ + ( 1 V * V ) V ˙ + ( 1 E * E ) E ˙ + d + δ + ε δ ( 1 I * I ) I ˙ = ( 1 S * S ) ( d s i n ( d + p ) S β μ ( I ) S ) + ( 1 V * V ) ( p S d V θ β μ ( I ) V ) + ( 1 E * E ) ( β μ ( I ) ( S + θ V ) ( d + δ + ε ) E ) + d + δ + ε δ ( 1 I * I ) ( δ E ( d + κ ) I ) = d s i n d S d s i n S * S + d S * + p S * + β μ ( I ) S * d V p S V * V + d V * + θ β μ ( I ) V * β μ ( I ) ( S + θ V ) E * E + ( d + δ + ε ) E * d + δ + ε δ ( d + κ ) I ( d + δ + ε ) E I * I + d + δ + ε δ ( d + κ ) I *
= d S * + d V * + β μ ( I * ) ( S * + θ V * ) d S d S * S * S d V * S * S β μ ( I * ) ( S * + θ V * ) S * S + d S * + d V * + θ β μ ( I * ) V * + β μ ( I ) S * d V d V * S S * V * V θ β μ ( I * ) V * S S * V * V + d V * + θ β μ ( I ) V * β μ ( I ) ( S + θ V ) E * E + β μ ( I * ) ( S * + θ V * ) β μ ( I * ) ( S * + θ V * ) I I * β μ ( I * ) ( S * + θ V * ) E E * I * I + β μ ( I * ) ( S * + θ V * ) = d S * ( 2 S S * S * S ) + d V * ( 3 S * S V V * S S * V * V ) + β μ ( I * ) S * + θ β μ ( I * ) V * β μ ( I * ) S * 2 S θ β μ ( I * ) S * V * S + θ β μ ( I * ) V * + β μ ( I ) S * θ β μ ( I * ) V * S S * V * V + θ β μ ( I ) V * β μ ( I ) E * S E θ β μ ( I ) E * V E + β μ ( I * ) S * + θ β μ ( I * ) V * β μ ( I * ) I S * I * θ β μ ( I * ) I V * I * β μ ( I * ) E E * I * I S * θ β μ ( I * ) E E * I * I V * + β μ ( I * ) S * + θ β μ ( I * ) V * = d S * ( 2 S S * S * S ) + d V * ( 3 S * S V V * S S * V * V ) + β μ ( I * ) S * 3 S * S E E * I * I I I * E * S E S * + θ β μ ( I * ) V * 4 S * S S S * V * V E E * I * I I I * E * V E V * + β μ ( I * ) S * μ ( I ) μ ( I * ) I I * 1 S S * E * E + θ β μ ( I * ) V * μ ( I ) μ ( I * ) I I * 1 V V * E * E .
As we have s 1 + s 2 + s 3 + + s m m s 1 . s 2 . s 3 s m m for all s 1 , s 2 , s 3 , , s m 0 , and m is an integer, then 2 S S * S * S 0 , 3 S * S V V * V * V S S * 0 ,
3 S * S I I * S S * E * E E E * I * I 0 , and 4 V * V S S * I I * V V * E * E E E * I * I S * S 0 . Therefore, by Lemma 2, F ˙ 2 0 in Ω 0 . Finally,
F ˙ * ( S , V , E , I , R ) = 0 ( S = S * , V = V * , E = E * , I = I * , R = R * ) .
Thus, { ( S , V , E , I , R ) | F ˙ * = 0 } = { E * } . Therefore, by applying LaSalle’s invariance principle [26], we deduce that the steady state E * is globally asymptotically stable in Ω 0 when it exists ( R 0 > 1 ). □

4. Seasonal Environment and Periodic Solution

Return now to the main Model (1), where all parameters are continuous and T-periodic positive functions. Let ( R m , R + m ) to be the ordered m-dimensional Euclidean space associated with the norm · . For X 1 , X 2 R m , we denote by X 1 X 2 if X 1 X 2 R + m . We denote by X 1 > X 2 if X 1 X 2 R + m \ { 0 } . We denote by X 1 X 2 if X 1 X 2 Int ( R + m ) . Consider a T-periodic m × m matrix function denoted by C ( t ) , which is continuous, irreducible, and cooperative. Let us denote by ϕ C ( t ) as the fundamental matrix, the solution of the following system
x ˙ ( t ) = C ( t ) x ( t ) .
Let us denote the spectral radius of the matrix ϕ C ( T ) by r ( ϕ C ( T ) ) . Therefore, all entries of ϕ C ( t ) are positive for each t > 0 . Let us apply the theorem of Perron–Frobenius [27] to deduce that r ( ϕ C ( T ) ) is the principal eigenvalue of ϕ C ( T ) (simple and admits an eigenvector y * 0 ). For the rest of the paper, the following lemma will be useful.
Lemma 3
([28]). There exists a positive T-periodic function y ( t ) such that x ( t ) = y ( t ) e k t will be a solution of system (4) where k = 1 T ln ( r ( ϕ C ( T ) ) ) .
Let us start by proving the existence (and uniqueness) of the disease-free periodic solution of System (1). Let us consider the following two-equations system
S ˙ ( t ) = d ( t ) s i n ( t ) ( d ( t ) + p ( t ) ) S ( t ) , V ˙ ( t ) = p ( t ) S ( t ) d ( t ) V ( t ) ,
with the initial condition ( S 0 , V 0 ) R + 2 . System (5) admits a unique T-periodic trajectory ( S * ( t ) , V * ( t ) ) , with S * ( t ) , V * ( t ) > 0 , and it is globally attractive in R + 2 . Therefore, System (1) admits a unique disease-free periodic trajectory ( S * ( t ) , V * ( t ) , 0 , 0 , 0 ) . For any continuous, positive T-periodic function h ( t ) , let h u = max t [ 0 , T ) h ( t ) and h l = min t [ 0 , T ) h ( t ) .
Let N ( t ) = S ( t ) + V ( t ) + E ( t ) + I ( t ) + R ( t ) be the population size at time t and N ¯ = d u d l s i n u . Then, we obtain the following lemma.
Lemma 4.
Ω u = { ( S , V , E , I , R ) R + 5 ; 0 S + V + E + I + R N ¯ } is a positively invariant attractor set for System (1). Furthermore, we have
lim t ( S ( t ) + V ( t ) + E ( t ) + I ( t ) + R ( t ) ( S * ( t ) + V * ( t ) ) ) = lim t ( N ( t ) ( S * ( t ) + V * ( t ) ) ) = 0 .
Proof. 
From System (1), we have
N ˙ ( t ) = d ( t ) s i n ( t ) d ( t ) N ( t ) d u s i n u d l N ( t ) 0   once d u d l s i n u N ( t ) .
This means that Ω u is a forward invariant compact absorbing set of all solutions of System (1). Next, we set y ( t ) = N ( t ) ( S * ( t ) + V * ( t ) ) for t 0 . Then, we obtain y ˙ ( t ) = d ( t ) y ( t ) , and this means that lim t y ( t ) = lim t ( N ( t ) ( S * ( t ) + V * ( t ) ) ) = 0 . □
Next, in Section 4.1, we will give the definition of the basic reproduction number R 0 , and we will prove that once R 0 is smaller than 1. Therefore, the disease-free periodic trajectory ( 0 , 0 , S * ( t ) , V * ( t ) , 0 ) is globally asymptotically stable, which means that the disease dies out. Then, in Section 4.2, we will prove that once R 0 > 1 , then System (1) is uniformly persistent. Therefore, we deduce that R 0 is the threshold parameter between the uniform persistence and the extinction of the disease.

4.1. Disease-Free Periodic Solution

We start by giving the definition of the basic reproduction number R 0 of System (1), using the theory given in [24], where F ( t , X ( t ) ) = β ( t ) μ ( I ( t ) ) ( S ( t ) + θ ( t ) V ( t ) ) δ ( t ) E ( t ) 0 0 0 , V ( t , X ( t ) ) = ( d ( t ) + δ ( t ) + ε ( t ) ) E ( t ) ( d ( t ) + κ ( t ) ) I ( t ) β ( t ) μ ( I ( t ) ) S ( t ) + ( d ( t ) + p ( t ) ) S ( t ) d ( t ) V ( t ) + θ ( t ) β ( t ) μ ( I ( t ) ) V ( t ) d ( t ) R ( t ) and V + ( t , X ( t ) ) = 0 0 d ( t ) s i n ( t ) p ( t ) S ( t ) ε ( t ) E ( t ) + κ ( t ) I ( t ) with X ( t ) = E ( t ) I ( t ) S ( t ) V ( t ) R ( t ) .
Our aim is to satisfy Conditions (A1)–(A7) in [24] (Section 1). Note that System (1) can have the following form
X ˙ ( t ) = F ( t , X ( t ) ) V ( t , X ( t ) ) = F ( t , X ( t ) ) V ( t , X ( t ) ) + V + ( t , X ( t ) ) .
The first five Conditions (A1)–(A5) are satisfied.
The System (8) admits a periodic solution (disease-free) X * ( t ) = 0 0 S * ( t ) V * ( t ) 0 . Let us define f ( t , X ( t ) ) = F ( t , X ( t ) ) V ( t , X ( t ) ) + V + ( t , X ( t ) ) , and M ( t ) = f i ( t , X * ( t ) ) X j 3 i , j 5 , where f i ( t , X ( t ) ) and I ( t ) are the i-th component of f ( t , X ( t ) ) and X ( t ) , respectively. By a simple calculus, we obtain M ( t ) = ( d ( t ) + p ( t ) ) 0 0 p ( t ) d ( t ) 0 0 0 d ( t ) ; thus, r ( ϕ M ( T ) ) < 1 . This means that X * ( t ) is linearly asymptotically stable in the subspace
Γ s = ( 0 , 0 , S , V , R ) R + 5 .
Therefore, Condition (A6) in [24] (Section 1) is satisfied.
Now, let us define F ( t ) and V ( t ) to be 2 by 2 matrices given by F ( t ) = F i ( t , X * ( t ) ) X j 1 i , j 2 , and V ( t ) = V i ( t , X * ( t ) ) X j 1 i , j 2 , where F j ( t , X ( t ) ) and V j ( t , X ( t ) ) are the j-th component of F ( t , X ( t ) ) and V ( t , X ( t ) ) , respectively. By an easy calculus, we obtain from System (8) that
F ( t ) = 0 β ( t ) μ ( 0 ) ( S * ( t ) + θ V * ( t ) ) δ ( t ) 0 and V ( t ) = d ( t ) + ε ( t ) + δ ( t ) 0 0 d ( t ) + κ ( t ) .
Consider Z ( t 1 , t 2 ) to be the two-by-two matrix solution of the system d d t Z ( t 1 , t 2 ) = V ( t 1 ) Z ( t 1 , t 2 ) for any t 1 t 2 , with Z ( t 1 , t 1 ) = I , the two by two identity matrix. Thus, condition (A7) was satisfied.
Let us define C T to be the ordered Banach space of T-periodic functions defined on R R 2 , associated with the maximum norm . and the positive cone C T + = { ψ C T : ψ ( s ) 0 ,   for   any   s R } . Consider the linear operator K : C T C T by
( K ψ ) ( s ) = 0 Z ( s , s w ) F ( s w ) ψ ( s w ) d w , s R , ψ C T
Let us now define the expression of the basic reproduction number, denoted by R 0 of the System (1), in terms of the spectral radius of the operator K as follows
R 0 = r ( K ) .
Therefore, we can conclude on the local stability of the disease-free periodic trajectory E 0 ( t ) = ( S * ( t ) , V * ( t ) , 0 , 0 , 0 ) for (1) as follows.
Theorem 3
([24] (Theorem 2.2)).
  • R 0 < 1 r ( ϕ F V ( T ) ) < 1 .
  • R 0 = 1 r ( ϕ F V ( T ) ) = 1 .
  • R 0 > 1 r ( ϕ F V ( T ) ) > 1 .
Therefore, E 0 ( t ) is unstable if R 0 > 1 , and it is asymptotically stable if R 0 < 1 .
Theorem 4.
The disease-free periodic solution E 0 ( t ) is globally asymptotically stable once R 0 < 1 . It is unstable if R 0 > 1 .
Proof. 
Using Theorem 3, E 0 ( t ) is locally stable once R 0 < 1 , and it is unstable once R 0 > 1 . Therefore, it remains to prove the global attractivity of E 0 ( t ) when R 0 < 1 . Consider the case where R 0 < 1 . Using the limit (6) in Lemma 4, for any δ 1 > 0 , there exists a positive constant T 1 > 0 satisfying S ( t ) + V ( t ) + E ( t ) + I ( t ) + R ( t ) S * ( t ) + V * ( t ) + δ 1 for t > T 1 . Then S ( t ) + θ V ( t ) S * ( t ) + θ V * ( t ) + δ 1 , and we deduce that
E ˙ ( t ) β ( t ) μ ( I ( t ) ) ( S * ( t ) + θ V * ( t ) + δ 1 ) ( d ( t ) + δ ( t ) + ε ( t ) ) E ( t ) , I ˙ ( t ) = δ ( t ) E ( t ) ( κ ( t ) + d ( t ) ) I ( t )
for t > T 1 . Let M 2 ( t ) to be the following two by two matrix function
M 2 ( t ) = 0 β ( t ) μ ( 0 ) 0 0 .
Again, using Theorem 3, we obtain r ( φ F V ( T ) ) < 1 . Let us choose δ 1 > 0 so that
r ( φ F V + δ 1 M 2 ( T ) ) < 1 . Let the two-equations model be hereafter
E ¯ ˙ ( t ) = β ( t ) μ ( I ¯ ( t ) ) ( S * ( t ) + θ V * ( t ) + δ 1 ) ( d ( t ) + δ ( t ) + ε ( t ) ) E ¯ ( t ) , I ¯ ˙ ( t ) = δ ( t ) E ¯ ( t ) ( κ ( t ) + d ( t ) ) I ¯ ( t ) .
Applying Lemma 3 and using the standard comparison principle, we deduce that there exists a positive T-periodic function y 1 ( t ) , satisfying x ( t ) y 1 ( t ) e k 1 t , where x ( t ) = E ( t ) I ( t ) and k 1 = 1 T ln r ( φ F V + δ 1 M 2 ( T ) < 0 . Thus, lim t E ( t ) = 0 and lim t I ( t ) = 0 . Therefore, we deduce that lim t R ( t ) = 0 . Furthermore, we have lim t ( S ( t ) + V ( t ) ) ( S * ( t ) + V * ( t ) ) = lim t N ( t ) E ( t ) I ( t ) R ( t ) ( S * ( t ) + V * ( t ) ) = 0 . Therefore, we deduce that the disease-free periodic trajectory E 0 ( t ) is globally attractive, which completes the proof. □
For the following subsection, we consider only the case where R 0 > 1 .

4.2. Endemic Periodic Solution

From Lemma 4, System (1) admits a positively invariant compact set Ω u . Now, since the recovery variable does not affect the other equations of System (1), then the Model (1) will be reduced as follows.
S ˙ ( t ) = d ( t ) s i n ( t ) ( d ( t ) + p ( t ) ) S ( t ) β ( t ) μ ( I ( t ) ) S ( t ) , V ˙ ( t ) = p ( t ) S ( t ) d ( t ) V ( t ) θ ( t ) β ( t ) μ ( I ( t ) ) V ( t ) , E ˙ ( t ) = β ( t ) μ ( I ( t ) ) ( S ( t ) + θ ( t ) V ( t ) ) ( d ( t ) + δ ( t ) + ε ( t ) ) E ( t ) , I ˙ ( t ) = δ ( t ) E ( t ) ( d ( t ) + κ ( t ) ) I ( t ) ,
with the initial condition ( S 0 , V 0 , E 0 , I 0 ) R + 4 such that S 0 > 0 , V 0 > 0 , E 0 > 0 and I 0 > 0 .
Let us define function P : R + 4 R + 4 to be the Poincaré map associated with System (13) with X 0 u ( T , X 0 ) , where u ( t , X 0 ) is the unique solution of the reduced Model (13) with the initial value u ( 0 , X 0 ) = X 0 R + 4 .
Let us define
Γ = ( S , V , E , I ) R + 4 , Γ 0 = I n t ( R + 4 )   and   Γ 0 = Γ \ Γ 0 .
Note that both Γ and Γ 0 are positively invariant from Lemma 4. P is point dissipative. Define
M = ( S 0 , V 0 , E 0 , I 0 ) Γ 0 : P n ( S 0 , V 0 , E 0 , I 0 ) Γ 0 ,   for   any   n 0 .
In order to use the uniform persistence theory detailed by Zhao [29] (or in [28] (Theorem 2.3)), we prove that
M = ( S , V , 0 , 0 ) , S 0 , V 0 .
Note that M ( S , V , 0 , 0 ) , S 0 , V 0 .
To show that M \ ( S , V , 0 , 0 ) , S 0 , V 0 = , let us consider ( S 0 , V 0 , E 0 , I 0 ) M \ ( S , V , 0 , 0 ) , S 0 , V 0 .
If I 0 = 0 and 0 < E 0 , then E ( t ) > 0 for any t > 0 . Then, it holds that I ˙ ( t ) | t = 0 = δ ( 0 ) E 0 > 0 . If I 0 > 0 and E 0 = 0 , then I ( t ) > 0 and S ( t ) > 0 for any t > 0 . Therefore, for any t > 0 , we have
E ( t ) = E 0 + 0 t β ( ω ) μ ( I ( ω ) ) ( S ( ω ) + θ ( ω ) V ( ω ) ) ) e 0 ω ( δ ( s ) + d ( s ) + ε ( s ) ) d s d ω × e 0 t ( δ ( s ) + d ( s ) + ε ( s ) ) d s > 0
for all t > 0 . This means that ( S ( t ) , V ( t ) , E ( t ) , I ( t ) ) Γ 0 for 0 < t 1 . Therefore, Γ 0 is positively invariant, from which we deduce Equation (14). Using the previous discussion, we deduce that there exists one fixed point ( S * ( 0 ) , V * ( 0 ) , 0 , 0 ) of P in M . We deduce, therefore, the uniform persistence of the disease as follows.
Theorem 5.
Consider the case R 0 > 1 . Equation (13) admits at least one positive periodic trajectory and γ > 0 , satisfying ( S 0 , V 0 , E 0 , I 0 ) R + 2 × I n t ( R + 2 ) ,
lim inf t I ( t ) γ > 0 .
Proof. 
Let us start by proving that P is uniformly persistent, respecting ( Γ 0 , Γ 0 ) , which will prove that the trajectory of the reduced Equation (13) is uniformly persistent, respecting ( Γ 0 , Γ 0 ) using [29] (Theorem 3.1.1). Recall that we obtain r ( φ F V ( T ) ) > 1 using Theorem 3. Therefore, η > 0 small enough and satisfies r ( φ F V η M 2 ( T ) ) > 1 . Let us consider the perturbed system
S ˙ α ( t ) = d ( t ) s i n ( t ) ( d ( t ) + p ( t ) ) S α ( t ) β ( t ) μ ( α ) S α ( t ) , V ˙ α ( t ) = p ( t ) S α ( t ) d ( t ) V α ( t ) θ ( t ) β ( t ) μ ( α ) V α ( t ) .
The function P associated with the perturbed System (15) has a unique positive fixed point ( S ¯ α 0 , V ¯ α 0 ) that it is globally attractive in R + 2 . Applying the implicit function theorem to deduce that ( S ¯ α 0 , V ¯ α 0 ) is continuously respecting α . Therefore, we can chose α > 0 to be small enough, satisfying S ¯ α ( t ) > S ¯ ( t ) η , and V ¯ α ( t ) > V ¯ ( t ) η , t > 0 . Let M 1 = ( S ¯ 0 , V ¯ 0 , 0 , 0 ) . As the trajectory is continuously respecting the initial condition, α * satisfies ( S 0 , V 0 , E 0 , I 0 ) Γ 0 with ( S 0 , V 0 , E 0 , I 0 ) u ( t , M 1 ) α * , and it holds that
u ( t , ( S 0 , V 0 , E 0 , I 0 ) ) u ( t , M 1 ) < α   for   0 t T .
We prove by contradiction that
lim sup n d ( P n ( S 0 , V 0 , E 0 , I 0 ) , M 1 ) α *   for   any   ( S 0 , V 0 , E 0 , I 0 ) Γ 0 .
Suppose that lim sup n d ( P n ( S 0 , V 0 , E 0 , I 0 ) , M 1 ) < α * for some ( S 0 , V 0 , E 0 , I 0 ) Γ 0 . In particular, we can suppose that d ( P n ( S 0 , V 0 , E 0 , I 0 ) , M 1 ) < α * , n > 0 . Therefore, we obtain u ( t , P n ( S 0 , V 0 , E 0 , I 0 ) ) u ( t , M 1 ) < α n > 0   and   0 t T . For all t 0 , let t = n T + t 1 , with t 1 [ 0 , T ) and n = [ t T ] (greatest integer t T ). Then, we get u ( t , ( S 0 , V 0 , E 0 , I 0 ) ) u ( t , M 1 ) = u ( t 1 , P n ( S 0 , V 0 , E 0 , I 0 ) ) u ( t 1 , M 1 ) < α , t 0 . Let ( S ( t ) , V ( t ) , E ( t ) , I ( t ) ) = u ( t , ( S 0 , V 0 , E 0 , I 0 ) ) . Therefore, 0 I ( t ) α for all t 0 and
S ˙ ( t ) d ( t ) s i n ( t ) ( d ( t ) + p ( t ) ) S ( t ) β ( t ) μ ( α ) S ( t ) , V ˙ ( t ) p ( t ) S ( t ) d ( t ) V ( t ) θ ( t ) β ( t ) μ ( α ) V ( t ) .
The fixed point S ¯ α 0 of the function P associated with the perturbed System (15) is globally attractive such that S ¯ α ( t ) > S ¯ ( t ) η , then T 2 > 0 large enough and satisfies
S ( t ) > S ¯ ( t ) η , t > T 2 .
Therefore, for t > T 2
E ˙ ( t ) β ( t ) μ ( I ( t ) ) ( S ¯ ( t ) η ) δ ( t ) E ( t ) d ( t ) E ( t ) , I ˙ ( t ) = δ ( t ) E ( t ) κ ( t ) I ( t ) d ( t ) I ( t ) .
Note that we have r ( φ F V η M 2 ( T ) ) > 1 . Applying Lemma 3 and the comparison principle, there exists a positive T-periodic trajectory y 2 ( t ) , satisfying J ( t ) e k 2 t y 2 ( t ) with k 2 = 1 T ln r φ F V η M 2 ( T ) > 0 , which implies that lim t I ( t ) = , which is impossible since the trajectories are bounded. Therefore, the inequality Equation (16) is satisfied, and P is weakly uniformly persistent respecting to ( Γ 0 , Γ 0 ) . By applying Lemma 4, P has a global attractor. We deduce that M 1 = ( S ¯ 0 , V ¯ 0 , 0 , 0 ) is an isolated invariant set inside X and W s ( M 1 ) Γ 0 = . All trajectory inside M converges to M 1 , which is acyclic in M . Applying [29] (Theorem 1.3.1 and Remark 1.3.1), we deduce that P is uniformly persistent, respecting ( Γ 0 , Γ 0 ) . Furthermore, using [29] (Theorem 1.3.6), P admits a fixed point ( S ˜ 0 , V ˜ 0 , E ˜ 0 , I ˜ 0 ) Γ 0 . Note that ( S ˜ 0 , V ˜ 0 , E ˜ 0 , I ˜ 0 ) R + 2 × I n t ( R + 2 ) . We prove also by contradiction that S ˜ 0 > 0 . Assume that S ˜ 0 = 0 . Using the first equation of the reduced System (13), S ˜ ( t ) verifies
S ˜ ˙ ( t ) d ( t ) s i n ( t ) β ( t ) μ ( I ˜ ( t ) ) S ˜ ( t ) d ( t ) S ˜ ( t ) ,
with S ˜ 0 = S ˜ ( m T ) = 0 , m = 1 , 2 , 3 , . Applying Lemma 4, δ 3 > 0 , there exists T 3 > 0 large enough and satisfying I ˜ ( t ) N ¯ + δ 3 , t > T 3 . Then, we have
S ˜ ˙ ( t ) d ( t ) s i n ( t ) β ( t ) μ ( N ¯ + δ 3 ) S ˜ ( t ) d ( t ) S ˜ ( t ) ,   for   t T 3
There exists a large enough m ¯ , satisfying m T > T 3 for all m > m ¯ . Applying the comparison principle, we deduce that
S ˜ ( m T ) = S ˜ 0 + 0 m T d ( ω ) s i n ( ω ) e 0 ω ( β ( u ) μ ( N ¯ + δ 3 ) + d ( u ) ) d u d ω × e 0 m T ( β ( u ) μ ( N ¯ + δ 3 ) + d ( u ) ) d u > 0
for any m > m ¯ , which is impossible. Therefore, S ˜ 0 > 0 and ( S ˜ 0 , V ˜ 0 , E ˜ 0 , I ˜ 0 ) is a positive T-periodic trajectory of the reduced System (13). □

5. More Examples, Numerical Results

We performed numerical simulations on the System (1) using the classical Monod function to express the transmission rate
μ ( I ) = μ m a x I k + I
where k and μ m a x are constants. Note that μ satisfies Assumption A1. The parameters of the model are T-periodic functions having the following forms:
s i n ( t ) = s i n 0 ( 1 + s i n 1 cos ( 2 π ( t + ϕ ) ) ) , d ( t ) = d 0 ( 1 + d 1 cos ( 2 π ( t + ϕ ) ) ) , p ( t ) = p 0 ( 1 + p 1 cos ( 2 π ( t + ϕ ) ) ) , β ( t ) = β 0 ( 1 + β 1 cos ( 2 π ( t + ϕ ) ) ) , θ ( t ) = θ 0 ( 1 + θ 1 cos ( 2 π ( t + ϕ ) ) ) , δ ( t ) = δ 0 ( 1 + δ 1 cos ( 2 π ( t + ϕ ) ) ) , ε ( t ) = ε 0 ( 1 + ε 1 cos ( 2 π ( t + ϕ ) ) ) , κ ( t ) = κ 0 ( 1 + κ 1 cos ( 2 π ( t + ϕ ) ) ) .
s i n 1 , d 1 , p 1 , β 1 , θ 1 , δ 1 , ε 1 , and κ 1 measure the amplitude of the seasonal variations in each of the parameters with | s i n 1 | < 1 , | d 1 | < 1 , | p 1 | < 1 , | β 1 | < 1 , | θ 1 | < 1 , | δ 1 | < 1 , | ε 1 | < 1 , and | κ 1 | < 1 . ϕ is the phase shift. Some fixed constants used for the numerical simulations are given in Table 2.
We will consider three cases. The first case is dedicated for the case of constant parameters (autonomous system) to validate the obtained theoretical results, concerning the local and global stability of the equilibrium points E 0 and E * . The second case deals only with a seasonal contact ( β ( t ) is a periodic function), where the other parameters are constants (partially non-autonomous system). The third case considers all parameters as periodic functions (non-autonomous system). For the numerical simulation, we used the standard solver “ode45” of MATLAB for ordinary differential equations, which implements an explicit Runge–Kutta (4,5) formula with a variable time step for efficient computation.

5.1. The Case of the Autonomous System

In a first step, we consider that all parameters of the System (1) are constants ( s i n 1 = d 1 = p 1 = β 1 = θ 1 = δ 1 = ε 1 = κ 1 = 0 ). Thus, the model is given by
S ˙ ( t ) = d 0 s i n 0 ( d 0 + p 0 ) S ( t ) β 0 μ ( I ( t ) ) S ( t ) , V ˙ ( t ) = p 0 S ( t ) d 0 V ( t ) θ 0 β 0 μ ( I ( t ) ) V ( t ) , E ˙ ( t ) = β 0 μ ( I ( t ) ) ( S ( t ) + θ 0 V ( t ) ) ( d 0 + δ 0 + ε 0 ) E ( t ) , I ˙ ( t ) = δ 0 E ( t ) ( d 0 + κ 0 ) I ( t ) , R ˙ ( t ) = ε 0 E ( t ) + κ 0 I ( t ) d 0 R ( t ) ,
with the positive initial condition ( S 0 , V 0 , E 0 , I 0 , R 0 ) = ( 10 , 10 , 6 , 8 , 10 ) R + 5 .
We give some numerical results that confirm the stability of the equilibrium points of Equation (22). In Figure 2, we give the results for the case where R 0 < 1 . The approximated solution of the given Model (22) approaches the equilibrium E 0 , which confirms that E 0 is globally asymptotically stable once R 0 1 . In Figure 3, we give the results for the case where R 0 > 1 . The approximated solution of the given Model (22) approaches asymptotically to E * . In Figure 4, we give the results for different initial conditions, The components of the approximated solutions converge asymptotically to the components of E * , which confirms that E * is globally asymptotically stable once R 0 > 1 .

5.2. The Case of the Partially Non-Autonomous System

In a second step, we perform numerical simulations on System (1), using a linear function to express the transmission rate, where only the seasonally forced T-periodic function β depends on time, t. The other parameters are constant ( d 1 = s i n 1 = δ 1 = κ 1 = 0 ). Thus, the model is given by
S ˙ ( t ) = d 0 s i n 0 ( d 0 + p 0 ) S ( t ) β 0 ( 1 + β 1 cos ( 2 π ( t + ϕ ) ) ) μ ( I ( t ) ) S ( t ) , V ˙ ( t ) = p 0 S ( t ) d 0 V ( t ) θ 0 β 0 ( 1 + β 1 cos ( 2 π ( t + ϕ ) ) ) μ ( I ( t ) ) V ( t ) , E ˙ ( t ) = β 0 ( 1 + β 1 cos ( 2 π ( t + ϕ ) ) ) μ ( I ( t ) ) ( S ( t ) + θ 0 V ( t ) ) ( d 0 + δ 0 + ε 0 ) E ( t ) , I ˙ ( t ) = δ 0 E ( t ) ( d 0 + κ 0 ) I ( t ) , R ˙ ( t ) = ε 0 E ( t ) + κ 0 I ( t ) d 0 R ( t ) ,
with positive initial condition ( S 0 , V 0 , E 0 , I 0 , R 0 ) = ( 10 , 10 , 6 , 8 , 10 ) R + 5 , where β 1 = 0.8 . The basic reproduction number R 0 is approximated using the time-averaged system. We give some numerical results that confirm the asymptotic behaviour of the solution of Equation (23). In Figure 5, we give the results for the case where R 0 > 1 . The approximated solution of the given Model (23) asymptotically approaches the periodic solution with the persistence of the disease. In Figure 6, we give the results for different initial conditions. The components of the approximated solutions asymptotically converge to the components of the periodic solution with the persistence of the disease. In Figure 7, we give the results for the case where R 0 < 1 . The approximated solution of the given Model (23) approaches the disease-free trajectory E 0 = d 0 s i n 0 d 0 + p 0 , p 0 s i n 0 d 0 + p 0 , 0 , 0 , 0 once R 0 1 .

5.3. The Case of Totally Non-Autonomous System

In a third step, we performed numerical simulations on System (1) using the classical Monod function to express the transmission rate, where all parameters were T-periodic functions. Thus, the model was given by
S ˙ ( t ) = d ( t ) s i n ( t ) ( d ( t ) + p ( t ) ) S ( t ) β ( t ) μ ( I ( t ) ) S ( t ) , V ˙ ( t ) = p ( t ) S ( t ) d ( t ) V ( t ) θ ( t ) β ( t ) μ ( I ( t ) ) V ( t ) , E ˙ ( t ) = β ( t ) μ ( I ( t ) ) ( S ( t ) + θ ( t ) V ( t ) ) ( d ( t ) + δ ( t ) + ε ( t ) ) E ( t ) , I ˙ ( t ) = δ ( t ) E ( t ) ( d ( t ) + κ ( t ) ) I ( t ) , R ˙ ( t ) = ε ( t ) E ( t ) + κ ( t ) I ( t ) d ( t ) R ( t ) ,
with positive initial condition ( S 0 , V 0 , E 0 , I 0 , R 0 ) = ( 10 , 10 , 6 , 8 , 10 ) R + 5 . Additional constants used for the numerical simulations in this step are given in Table 3.
We gave some numerical results that confirmed the asymptotic behaviour of the solution of Equation (24). The basic reproduction number R 0 was approximated using the time-averaged system.
In Figure 8, we give the results for the case where R 0 < 1 . The approximated solution of the given Model (24) approaches the disease-free periodic trajectory E 0 ( t ) = ( S * ( t ) , V * ( t ) , 0 , 0 , 0 ) once R 0 1 . In Figure 9, we give the results for the case where R 0 > 1 . The approximated solution of the given Model (24) approaches asymptotically to a periodic solution with the persistence of the disease. In Figure 10, we give the results for different initial conditions. The components of the approximated solutions asymptotically converge to the components of the periodic solution with the persistence of the disease.

6. Conclusions

In this work, we proposed an extension of the SVEIR epidemic model already considered in [8] by taking into account the seasonal environment. This model could describe seasonal flu behaviour. In the first step, we studied the case of the autonomous system, where all parameters were supposed to be constants. In the second step, we considered the non-autonomous system, gave some basic results, and defined the basic reproduction number R 0 . We showed the global asymptotic stability of the disease-free periodic solution once R 0 < 1 ; however, the disease persisted once R 0 was greater than 1. Finally, we gave some numerical examples that supported the theoretical findings, including the autonomous system, the partially non-autonomous system, and the full non-autonomous system. It was deduced that the trajectories converged to one of the equilibriums of the System (2), according to Theorems 1 and 2, if the system was autonomous. However, if at least one of the model parameters was periodic, the trajectories converged to a limit cycle according to Theorems 4 and 5. This study remained theoretical, and we need real data to compare, improve, and validate the model.

Author Contributions

Conceptualization, M.E.H., D.M.A. and N.A.A.; methodology, M.E.H., D.M.A. and N.A.A.; writing—original draft, M.E.H., D.M.A. and N.A.A.; writing—review and editing, M.E.H., D.M.A. and N.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the unknown referees for the many constructive suggestions, which helped to improve the presentation of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. Influenza (Seasonal). Available online: https://www.emro.who.int/health-topics/influenza/influenza-seasonal.html (accessed on 30 March 2023).
  2. Centers for Disease Control and Prevention; National Center for Immunization and Respiratory Diseases (NCIRD). Types of Influenza Viruses. Available online: https://www.cdc.gov/flu/about/index.html (accessed on 30 March 2023).
  3. Demicheli, V.; Jefferson, T.; Ferroni, E.; Rivetti, A.; Pietrantonj, C.D. Vaccines for preventing influenza in healthy adults. Cochrane Database Syst. Rev. 2018, 2018, CD001269. [Google Scholar] [CrossRef] [PubMed]
  4. Nakata, Y.; Kuniya, T. Global dynamics of a class of SEIRS epidemic models in a periodic environment. J. Math. Anal. Appl. 2010, 363, 230–237. [Google Scholar] [CrossRef]
  5. Bacaër, N.; Ouifki, R. Growth rate and basic reproduction number for population models with a simple periodic factor. Math. Biosci. 2007, 210, 647–658. [Google Scholar] [CrossRef]
  6. Bacaër, N. Approximation of the basic reproduction number R0 for vector-borne diseases with a periodic vector population. Bull. Math. Biol. 2007, 69, 1067–1091. [Google Scholar] [CrossRef] [PubMed]
  7. Bacaër, N.; Guernaoui, S. The epidemic threshold of vector-borne diseases with seasonality. J. Math. Biol. 2006, 53, 421–436. [Google Scholar] [CrossRef]
  8. El Hajji, M.; Albargi, A.H. A mathematical investigation of an “SVEIR” epidemic model for the measles transmission. Math. Biosci. Eng. 2022, 19, 2853–2875. [Google Scholar] [CrossRef]
  9. El Hajji, M. Modelling and optimal control for Chikungunya disease. Theory Biosci. 2021, 140, 27–44. [Google Scholar] [CrossRef]
  10. Alsolami, A.A.; El Hajji, M. Mathematical Analysis of a Bacterial Competition in a Continuous Reactor in the Presence of a Virus. Mathematics 2023, 11, 883. [Google Scholar] [CrossRef]
  11. Nkamba, L.; Ntaganda, J.; Abboubakar, H.; Kamgang, J.; Castelli, L. Global Stability of a SVEIR Epidemic Model: Application to Poliomyelitis Transmission Dynamics. Open J. Model. Simul. 2017, 5, 98–112. [Google Scholar] [CrossRef]
  12. Tang, Y.; Xiao, D.; Zhang, W.; Zhu, D. Dynamics of epidemic models with asymptomatic infection and seasonal succession. Math. Biosci. Eng. 2017, 14, 1407–1424. [Google Scholar] [CrossRef]
  13. Adda, P.; Nkague Nkamba, L.; Sallet, G.; Castelli, L. A SVEIR model with Imperfect Vaccine. In Proceedings of the CMPD 3 Conference on Computational and Mathematical Population Dynamics, Bordeaux, France, 31 May–4 June 2010. [Google Scholar]
  14. Wei, H.; Jiang, Y.; Song, X.; Su, G.; Qiu, S. Global attractivity and permanence of a SVEIR epidemic model with pulse vaccination and time delay. J. Comput. Appl. Math. 2009, 229, 302–312. [Google Scholar] [CrossRef]
  15. Gumel, A.B.; McCluskey, C.C.; Watmough, J. An SVEIR model for assessing potential impact of an imperfect anti-SARS vaccine. Math. Biosci. Eng. 2006, 3, 485–512. [Google Scholar] [PubMed]
  16. El Hajji, M.; Zaghdani, A.; Sayari, S. Mathematical analysis and optimal control for Chikungunya virus with two routes of infection with nonlinear incidence rate. Int. J. Biomath. 2022, 15, 2150088. [Google Scholar] [CrossRef]
  17. El Hajji, M.; Sayari, S.; Zaghdani, A. Mathematical analysis of an SIR epidemic model in a continuous reactor—Deterministic and probabilistic approaches. J. Korean Math. Soc. 2021, 58, 45–67. [Google Scholar]
  18. Xiao, D. Dynamics and bifurcations on a class of population model with seasonal constant-yield harvesting. Discret. Contin. Dyn. Syst. B 2016, 21, 699–719. [Google Scholar] [CrossRef]
  19. Bacaër, N.; Gomes, M. On the Final Size of Epidemics with Seasonality. Bull. Math. Biol. 2009, 71, 1954–1966. [Google Scholar] [CrossRef] [PubMed]
  20. Kermack, F.; McKendrick, D. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 1927, 115, 700–721. [Google Scholar]
  21. Ma, J.; Ma, Z. Epidemic threshold conditions for seasonally forced SEIR models. Math. Biosci. Eng. 2006, 3, 161–172. [Google Scholar] [CrossRef]
  22. Guerrero-Flores, S.; Osuna, O.; de Leon, C.V. Periodic solutions for seasonal SIQRS models with nonlinear infection terms. Electron. J. Differ. Equations 2019, 2019, 1–13. [Google Scholar]
  23. Tailei Zhang, Z.T. On a nonautonomous SEIRS model in epidemiology. Bull. Math. Biol. 2007, 69, 2537–2559. [Google Scholar] [CrossRef]
  24. Wang, W.; Zhao, X. Threshold dynamics for compartmental epidemic models in periodic environments. J. Dyn. Differ. Equ. 2008, 20, 699–717. [Google Scholar] [CrossRef]
  25. Alshehri, A.; El Hajji, M. Mathematical study for Zika virus transmission with general incidence rate. AIMS Math. 2022, 7, 7117–7142. [Google Scholar] [CrossRef]
  26. LaSalle, J. The Stability of Dynamical Systems; SIAM: Philadelphia, PA, USA, 1976. [Google Scholar]
  27. Frobenius, G. Uber Matrizen aus nicht negativen Elementen. Sitz. Preuss. Akad. Wiss. 1912, 26, 456–477. [Google Scholar]
  28. Zhang, F.; Zhao, X. A periodic epidemic model in a patchy environment. J. Math. Anal. Appl. 2007, 325, 496–516. [Google Scholar] [CrossRef]
  29. Zhao, X. Dynamical Systems in Population Biology; CMS Books in Mathematics; Springer: New York, NY, USA, 2003; Volume 16. [Google Scholar]
Figure 1. Epidemic compartmental model with seasonal environment (inspired from El Hajji and Albargi [8], where the authors consider a similar model with fixed parameters).
Figure 1. Epidemic compartmental model with seasonal environment (inspired from El Hajji and Albargi [8], where the authors consider a similar model with fixed parameters).
Mathematics 11 02350 g001
Figure 2. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.3 then R 0 0.7187 < 1 .
Figure 2. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.3 then R 0 0.7187 < 1 .
Mathematics 11 02350 g002
Figure 3. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.8 then R 0 1.9167 > 1 .
Figure 3. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.8 then R 0 1.9167 > 1 .
Mathematics 11 02350 g003
Figure 4. Behaviour of the solution components for different initial conditions, where k = 2 and μ m a x = 0.8 ( R 0 1.9167 > 1 ). The different graphs (colors) reflect different initial conditions.
Figure 4. Behaviour of the solution components for different initial conditions, where k = 2 and μ m a x = 0.8 ( R 0 1.9167 > 1 ). The different graphs (colors) reflect different initial conditions.
Mathematics 11 02350 g004
Figure 5. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.8 then R 0 1.9167 > 1 .
Figure 5. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.8 then R 0 1.9167 > 1 .
Mathematics 11 02350 g005
Figure 6. Behaviour of the solution components for different initial conditions where k = 2 and μ m a x = 0.8 ( R 0 1.9167 > 1 ). The different graphs (colors) reflect different initial conditions.
Figure 6. Behaviour of the solution components for different initial conditions where k = 2 and μ m a x = 0.8 ( R 0 1.9167 > 1 ). The different graphs (colors) reflect different initial conditions.
Mathematics 11 02350 g006
Figure 7. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.3 then R 0 0.7187 < 1 .
Figure 7. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.3 then R 0 0.7187 < 1 .
Mathematics 11 02350 g007
Figure 8. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.3 then R 0 0.7187 < 1 .
Figure 8. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.3 then R 0 0.7187 < 1 .
Mathematics 11 02350 g008
Figure 9. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.8 then R 0 1.9167 > 1 .
Figure 9. Behaviour of the solution of System (1) for k = 2 and μ m a x = 0.8 then R 0 1.9167 > 1 .
Mathematics 11 02350 g009
Figure 10. Behaviour of the solution components for different initial conditions where k = 2 and μ m a x = 0.8 ( R 0 1.9167 > 1 ). The different graphs (colors) reflect different initial conditions.
Figure 10. Behaviour of the solution components for different initial conditions where k = 2 and μ m a x = 0.8 ( R 0 1.9167 > 1 ). The different graphs (colors) reflect different initial conditions.
Mathematics 11 02350 g010
Table 1. Parameters and variables of System (1).
Table 1. Parameters and variables of System (1).
NotationDefinition
S ( t ) Susceptible individuals
V ( t ) Vaccinated individuals
E ( t ) Exposed individuals
I ( t ) Infected individuals
R ( t ) Recovered individuals
μ ( · ) Saturated incidence rate
β ( t ) Disease transmission coefficient
p ( t ) Vaccination rate
s i n ( t ) Instantaneous recruitment rate
θ ( t ) [ 0 , 1 ] Instantaneous effectiveness rate of the vaccination
d ( t ) Common instantaneous death rate
δ ( t ) 1/the instantaneous duration time spent in compartment E
ε ( t ) Exposed instantaneous recovery rate
1 / ε ( t ) The instantaneous duration time spent in compartment E
κ ( t ) Infected instantaneous recovery rate
1 / κ ( t ) The instantaneous duration time spent in compartment I
Table 2. Some fixed parameters for numerical simulations.
Table 2. Some fixed parameters for numerical simulations.
Parameter s i n 0 d 0 p 0 β 0 θ 0 δ 0 ε 0 κ 0 ϕ
Value10 0.8 0.7 2 0.5 1 0.2 0.8 4
Table 3. Additional parameters for numerical simulations of the totally non-autonomous system.
Table 3. Additional parameters for numerical simulations of the totally non-autonomous system.
Parameter s i n 1 d 1 p 1 β 1 θ 1 δ 1 ε 1 κ 1
Value 0.75 0.6 0.4 0.8 0.8 0.6 0.6 0.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

El Hajji, M.; Alshaikh, D.M.; Almuallem, N.A. Periodic Behaviour of an Epidemic in a Seasonal Environment with Vaccination. Mathematics 2023, 11, 2350. https://doi.org/10.3390/math11102350

AMA Style

El Hajji M, Alshaikh DM, Almuallem NA. Periodic Behaviour of an Epidemic in a Seasonal Environment with Vaccination. Mathematics. 2023; 11(10):2350. https://doi.org/10.3390/math11102350

Chicago/Turabian Style

El Hajji, Miled, Dalal M. Alshaikh, and Nada A. Almuallem. 2023. "Periodic Behaviour of an Epidemic in a Seasonal Environment with Vaccination" Mathematics 11, no. 10: 2350. https://doi.org/10.3390/math11102350

APA Style

El Hajji, M., Alshaikh, D. M., & Almuallem, N. A. (2023). Periodic Behaviour of an Epidemic in a Seasonal Environment with Vaccination. Mathematics, 11(10), 2350. https://doi.org/10.3390/math11102350

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop