 Research
 Open Access
 Published:
Discrete epidemic models with two time scales
Advances in Difference Equations volume 2021, Article number: 478 (2021)
Abstract
The main aim of the work is to present a general class of two time scales discretetime epidemic models. In the proposed framework the disease dynamics is considered to act on a slower time scale than a second different process that could represent movements between spatial locations, changes of individual activities or behaviors, or others.
To include a sufficiently general disease model, we first build up from first principles a discretetime susceptible–exposed–infectious–recovered–susceptible (SEIRS) model and characterize the eradication or endemicity of the disease with the help of its basic reproduction number \(\mathcal{R}_{0}\).
Then, we propose a general full model that includes sequentially the two processes at different time scales and proceed to its analysis through a reduced model. The basic reproduction number \(\overline{\mathcal{R}}_{0}\) of the reduced system gives a good approximation of \(\mathcal{R}_{0}\) of the full model since it serves at analyzing its asymptotic behavior.
As an illustration of the proposed general framework, it is shown that there exist conditions under which a locally endemic disease, considering isolated patches in a metapopulation, can be eradicated globally by establishing the appropriate movements between patches.
Introduction
Infectious diseases such as SARSCoV2, AIDS, Ebola, or COVID19 are becoming a part of usual life. They can catastrophically spread and cause a significant number of deaths. It is crucial to use the best management methods to curtail their harmful consequences. The only way to try to compare the effectiveness of these methods is to formulate appropriate mathematical models that help us on making predictions [4]. Mathematical models in epidemiology have a long history of more than two centuries. Most of these models are formulated in continuous time, possibly because of the wealth of analytical tools available for their study. Nevertheless, at least in the last twenty years, timediscrete mathematical epidemics models have been also used with a significant and increasing frequency.
Formulating epidemic models in discrete time has some advantages over the differential equation models, specially when these latter are untractable analytically and one must resort to numerical simulations. Discretetime models are better implemented in computer simulations when needed, and their parameters can be more easily related to data due to the natural fit of discrete time units to the periodic data collection used in the laboratory or the field [15]. The formulation of discretetime models should be done directly from first principles and not as a discretization of continuous time models [7, 19]. The mere discretization of continuous models could lead to unfeasible results and usually deviate the focus from the relevant disease dynamics analysis.
Some authors have proposed discretetime epidemic models without demographic dynamics [1, 6]. The discretetime epidemic models that we propose in this work follow the literature [13, 27, 29] that does allow for demographic effects. The epidemic and the demographic processes are sequentially included in the model in the following way: at each time interval, first the individuals change their disease status according to the disease flow and then reproduce into the susceptible class and survive following the demographic flow. This constitutes a difference with the case of continuoustime models where these processes are supposed to act instantaneously and simultaneously.
It is frequent in epidemic models that parameters vary by orders of magnitude. In models including disease and demographic dynamics, it is typical that infectious periods have lengths of the order of days, whereas the life spans have them of the order of years [16, 17]. In vectorborne epidemic models in which the vector is an insect, the time scale associated with the vector is usually much faster than the host’s one [5, 21, 26]. Individual mobility and behavior can also entail the existence of time scales in the disease dynamics [12, 24]. The models proposed in the previous references are continuoustime and expressed in the form of two time scales models. The analytical tools at disposal for the analysis of this kind of model, quasisteadystate hypothesis [5, 26] or geometric singular perturbation theory [16, 17, 21, 24, 28], serve at simplifying its analysis by previously reducing the dimension of the system of ordinary differential equations to be studied.
We present in this work a family of timediscrete epidemic models with two time scales. We consider a population divided into groups that we call patches as if they constituted a generalized metapopulation. The patches could represent real spatial locations but also different individual activities or different individual behaviors. The epidemic process acts locally in each patch and may differ from patch to patch. The key point in considering the patches is that slow and fast time scales can be associated with, respectively, the disease dynamics and the process of patch changes that henceforth we call movements.
In [8, 23] it is shown how to construct the kind of discretetime model with two time scales that we are using. A reduction method that helps to carry out the analytical study of the model is also developed. It is assumed that within a slow time unit the slow process, epidemic, is defined by a map \(\mathcal{S}\) and, analogously, within a fast time unit the fast process, movements, is defined by a map \(\mathcal{F}\). Thus, the effect of the fast process along a slow time unit can be described by the kth iterate of \(\mathcal{F}\), \(\mathcal{F}^{(k)}\), where k approximates the time scales ratio. The combined effect of both processes during a slow time unit can then be seen sequentially as the occurrence of k movement episodes followed by a disease dynamics one. In terms of maps \(\mathcal{F}\) and \(\mathcal{S}\), the associated discretetime system, denoting by X the population vector and by t the slow time variable, takes the form
The slow dynamics, represented by \(\mathcal{S}\), corresponds to the local disease model in each patch. As it is not possible to use a general disease model, we have proposed as a rather general case a discretetime SEIRS model based upon the same assumptions of the continuoustime one used in [3]. We have built the model from first principles following the sequential inclusion of the epidemic and the demographic processes as in [13, 27, 29]. The model is the same in all patches, which are distinguished by different parameter values.
In this work the fourdimensional discretetime system representing the SEIRS model is studied with the help of the basic reproduction number \(\mathcal{R}_{0}\) obtaining sufficient conditions that ensure either the global asymptotic stability of the diseasefree equilibrium (DFE) or the uniform persistence of the disease. Thus, disease eradication or endemicity can be straightforwardly characterized.
The fast dynamics, represented by map \(\mathcal{F}\), corresponds to the movements of individuals between patches. They are defined for each of the disease compartments (S, E, I, and R) by a probability matrix, which in general can depend on the total number of individuals in each compartment across patches.
The proposed full model is represented by the 4ndimensional discretetime system (1), where n is the number of patches. The reduction method [20] applied to the full model leads to a reduced fourdimensional system. If the movement rates are constant, this system is similar to the SEIRS model already studied. The asymptotic analysis of the full model can then be undertaken by studying the asymptotic behavior of this reduced system.
The main aim of this work is to develop a general framework where two timescale discretetime epidemic models can be included. This is achieved by proceeding as described above. To complete the presentation, we have also treated a relevant application in a particular setting. We have shown that there exist conditions under which a locally endemic disease, considering every patch in isolation, can be eradicated globally by establishing the appropriate movements.
The structure of the paper is as follows: In Sect. 2, we first propose a discretetime SEIRS epidemic model, we compute its \(\mathcal{R}_{0}\), and we find conditions for the DFE to be globally asymptotically stable (GAS) when \(\mathcal{R}_{0}<1\) and for the disease to persist when \(\mathcal{R}_{0}>1\). Section 3 is devoted to the presentation and reduction of a model combining the disease dynamics at a slow time scale with a fast process interpretable as movements between patches. The previous general model is analyzed in Sect. 4 in the particular case of constant recruitment, standard incidence, and constant movement rates in order to illustrate the advantages of the reduction procedure. We summarize our results and perspectives in the final Discussion section.
Discrete SEIRS model
In this section we present a discretetime SEIRS epidemic model which is a discrete version of the continuous SEIRS model included in [3].
The individuals of the population are divided into four epidemiological classes that, at each time \(t\in \{0,1,2,\dots \}\), are represented by \(S(t)\), susceptible, \(E(t)\), exposed, \(I(t)\), infectious, and \(R(t)\), recovered. So, the state vector of the population is \(X(t)=(S(t),E(t),I(t), R(t))\), and the total population \(N(t)=S(t)+E(t)+I(t)+R(t)\).
Following [27, 29], we consider that in each time interval there exist two distinct temporal phases. In the first one, the disease dynamics acts and individuals can change from one epidemiological class to the next one, as shown in Fig. 1. Reproduction and survival happen in the second temporal phase. We assume that the disease does not affect the birth process. Concerning survival, to take into account both the natural and the disease induced mortalities in full generality, we associate different survival rates with each of the epidemiological classes. The corresponding diagram is included in Fig. 2.
The disease transmission is a function \(\Phi (I(t)/N(t))\) of the fraction of susceptible individuals that become exposed, where
A common choice, which corresponds to the socalled proportional or standard incidence, is \(\Phi (x)=\beta x\), where \(\beta \in (0,1]\) is the transmission parameter. When infections are modeled as Poisson processes, then \(\Phi (x)=1e^{\beta x}\) for \(\beta >0\) [29].
Upon infection, the transitions between classes are defined by parameters \(\gamma ^{C}\in (0,1)\) for \(C\in \{E,I,R\}\) that represent the fraction of the individuals of class C that pass to the next class per time unit (see Fig. 1).
Henceforth, we use C as a generic letter for an unspecified epidemiological class and call \(\boldsymbol{\mathcal{C}}:=\{S,E,I,R\}\).
Demography is included in the model in a simple form. It is assumed that there is no vertical transmission of the disease, so that all births occur into the susceptible class. The recruitment of individuals to the susceptible class per time unit is a function
of the total population N. We will use a constant recruitment function as a particular case. Other common choices, proposed in [27, 29], are geometric, Beverton–Holt, and Ricker recruitment functions.
Individuals in all epidemiological classes are subject to natural death, possibly affected by the disease. Parameters \(\sigma ^{C}\in (0,1)\), for \(C\in \boldsymbol{\mathcal{C}}\), represent the fraction of the individuals of class C that survive per time unit. Thus, the fraction of individuals dying per time unit is \(1\sigma ^{C}\in (0,1)\) (see Fig. 2).
The transitions between epidemiological classes in a time unit are defined by the following map:
and the demographic changes by the next one
We now take into account both processes sequentially: first epidemiological transitions followed by demography. The discretetime SEIRS epidemic model that we propose can so be expressed in terms of maps T and D as
or in a detailed form
The assumptions on the parameters of the model allow us to straightforwardly prove that \(T (\mathbf{R}_{+}^{4} )\subset \mathbf{R}_{+}^{4}\) and also \(D (\mathbf{R}_{+}^{4} )\subset \mathbf{R}_{+}^{4}\). Therefore, \(\mathbf{R}_{+}^{4}\) is forward invariant under the semiflow defined by the discretetime system (4).
Proposition 1
If the function B is bounded, then system (4) is dissipative.
Proof of Proposition 1
Let \(\hat{B}>0\) be such that \(B(x)\leq \hat{B}\) for all \(x\in {}[ 0,\infty )\), and \(\hat{\sigma }=\max_{C\in \boldsymbol{\mathcal{C}}}\{ \sigma ^{C}\}\in (0,1)\). Then
We now present some basic properties for the solution of the linear scalar difference equation
with \(0< a<1\), \(b\geq 0\) that will be used here and later on in the manuscript. The solution to (6) is
from where it follows that \(x(t)\) converges monotonically to \(b/(1a)\). In particular, for all \(x(0)\geq 0\), we have
Therefore, from (5) and (7) it follows that any solution of system (4), with initial conditions \(X(0)\in \mathbf{R}_{+}^{4}\), satisfies that
This inequality proves that, fixed any \(M>\hat{B}/(1\hat{\sigma })\), for every solution of (4), there is a time t̂ such that \(N(t)\leq M\) for all \(t\geq \hat{t}\) and, thus, the dissipativity of system (4). □
Function B is bounded in the particular case of having a constant, Beverton–Holt, or Ricker recruitment function.
Notice that in the proof it is shown that all nonnegative solutions of system (4) are attracted by the compact set \(K=\{(S,E,I,R)\in \mathbf{R}_{+}^{4}:N\in [0,\hat{B}/(1\hat{\sigma })] \}\).
We continue the analysis by finding conditions for the eradication or endemicity of the disease.
To consider disease eradication, we try to find an equilibrium of (4) with \(I=0\). We immediately obtain \(E=0\), \(R=0\), and \(S=B(S)+\sigma ^{S} S\). To guarantee the existence of a unique diseasefree equilibrium (DFE), we make the following assumption on the scalar difference equation representing the demography of the population without disease.
Hypothesis 2.1
The equation \(S(t+1) = \sigma ^{S} S(t)+B(S(t))\) possesses a unique positive equilibrium \(S^{*}\) that is hyperbolic and globally asymptotically stable (GAS) in \((0,\infty )\).
If the recruitment function is constant, then Hypothesis 2.1 holds. This also happens, for certain values of the parameters, in the cases of the Beverton–Holt and the Ricker recruitment functions.
Note that if Hypothesis 2.1 is met, then the unique DFE of system (4) is \(X_{0}^{*}=(S^{*},0,0,0)\).
As we will prove, the basic reproduction number \(\mathcal{R}_{0}\) of system (4) determines whether the disease is eradicated or it becomes endemic. We use the nextgeneration method to calculate \(\mathcal{R}_{0}\) of a discretetime system as developed in [2]. We consider as infected states E and I, and as uninfected states S and R. The socalled diseasefree system is the one associated with the uninfected states setting \(E=0\) and \(I=0\):
Proposition 2
If function B is bounded and Hypothesis 2.1is met, then \((S^{*},0)\) is a GAS equilibrium of system (9) in \(\Omega = \{ (S,R)\in \mathbf{R}_{+}^{2}:S>0 \} \).
Proof of Proposition 2
It is straightforward to show that \((S^{\ast },0)\) is an equilibrium of system (9) and, by linearization, that it is hyperbolic and locally asymptotically stable (LAS).
To prove that \((S^{*},0)\) attracts all points in Ω, we apply Theorem 2.1 in [30].
For that, let \((S_{0},R_{0})\in \Omega \) and define, for \(t=0,1,2,\ldots \) ,
that is, a continuous map in \(\mathbf{R}_{+}\), and the discrete dynamical process defined by \(\tau _{0}:=I\), the identity map, and \(\tau _{t}:=\sigma _{t1}\circ \sigma _{t2}\circ \cdots \circ \sigma _{1}\circ \sigma _{0}\) for \(t\geq 1\).
If \(\{(S(t),R(t)):t\geq 0\}\) is the orbit associated with \((S_{0},R_{0})\) in system (9), we have that \(R(t)=(\sigma ^{R}(1\gamma ^{R}))^{t} R_{0}\) and \(S(t)=\tau _{t}(S_{0})\).
Since \(0<\sigma ^{R}(1\gamma ^{R})<1\), we have \(\lim_{t\rightarrow \infty }R(t)=0\), and so we need to show that \(\lim_{t\rightarrow \infty }S(t)=S^{\ast }\).
We now prove that the discrete dynamical process \(\tau _{t}\) (\(t\geq 0\)) is asymptotically autonomous (Definition 2.1 in [30]) and that its limit discrete semiflow, \(\Sigma _{t}=\Sigma \circ \overset{(t)}{\cdots }\circ \Sigma \) (\(t \geq 0\)), is that generated by the continuous map \(\Sigma (S)=B(S)+\sigma ^{S}S\) in \(\mathbf{R}_{+}\). In order to do so, let \(x\in \mathbf{R}_{+}\), and let \(\{ x_{n} \} _{n=1}^{\infty }\) and \(\{ t_{n} \} _{n=1}^{\infty }\) be sequences in \(\mathbf{R}_{+}\) such that \(\lim_{n\rightarrow \infty }x_{n}=x\) and \(\lim_{n\rightarrow \infty }t_{n}=\infty \). We need to show that \(\lim_{n\rightarrow \infty }\Sigma _{t_{n}}(x_{n})=\Sigma (x)\), which follows immediately taking into account that \(0<\sigma ^{R}(1\gamma ^{R})<1\) and that B is continuous.
As the only Σinvariant subset of \(\mathbf{R}_{+}\) is \(\{S^{\ast }\}\), to complete the proof by applying Theorem 2.1 in [30], we just need to prove that the set \(\{S(t):t\geq 0\}\) is bounded.
Let B̂ be un upperbound of function B. Then, for every \(t\geq 0\), we have \(\sigma _{t}(S)\leq \hat{B}+R_{0}+\sigma ^{S}S\), and so \(\sigma _{t}(S_{0})\leq x(t)\), where \(x(t)\) is the solution to
which corresponds to the linear scalar equation (6) with \(a:=\sigma ^{S}\) and \(b:=\hat{B}+R_{0}\). Now, using (8), we have \(\sigma _{t}(S_{0})\leq x(t)\leq \max \{ S_{0},(\hat{B}+R_{0})/(1 \sigma ^{S}) \} \), and so
as we wanted to show. □
To proceed with the application of the nextgeneration method, in the equations for the infected compartments we must separate the terms \(\mathcal{F}_{E}\) and \(\mathcal{F}_{I}\), representing new infections, from the terms \(\mathcal{T}_{E}\) and \(\mathcal{T}_{I}\) associated with transitions between compartments:
Now, we can calculate matrices
and
to obtain the nextgeneration matrix \(Q=F(IdT)^{1}\) [2], where Id is the identity matrix of order 2 whose spectral radius is \(\mathcal{R}_{0}\).
The next result gives sufficient conditions for the disease eradication locally and globally.
Theorem 3
Let system (4) satisfy Hypothesis 2.1. Then

(a)
If \(\mathcal{R}_{0}<1\), then DFE \(X_{0}^{*}\) is LAS.

(b)
If \(\mathcal{R}_{0}< 1\) and \(\Phi ''(x)\leq 0\) for \(x>0\), then DFE \(X_{0}^{*}\) is GAS in
\(\Omega = \{ (S,E,I,R)\in \mathbf{R}_{+}^{4}:S>0 \} \).

(c)
If \(\mathcal{R}_{0}>1\), then DFE \(X_{0}^{*}\) is unstable.
Proof of Theorem 3
(a) and (c) are direct consequences of Theorem 2.1 in [2].
(b) \(\Phi \in C^{2}([0,1])\), \(\Phi (0)=0\), and \(\Phi ''(x)\leq 0\) imply that \(\Phi (x)\leq \Phi '(0)x\) for \(x\in [0,1]\).
Therefore
and so, using the equation for I in (4),
Matrix \(F+T\) can be considered to be the projection matrix of a standard linear matrix model of population dynamics [18], with matrices F and T representing the fertility and transition matrices respectively. The net reproductive rate of the model coincides with \(\mathcal{R}_{0}\) and, since \(\rho (T)=\max \{\sigma ^{E}(1\gamma ^{E}),\sigma ^{I}(1\gamma ^{I}) \}<1\) and \(\mathcal{R}_{0}<1\), we can apply Theorem 3.3 in [18], which yields that \(\rho (F+T)<1\) and, therefore, \((F+T)^{t}\underset{t\to \infty }{\longrightarrow }0\). This proves that \(E(t)\), \(I(t)\underset{t\to \infty }{\longrightarrow }0\).
To prove that also \(R(t)\) tends to 0, we use from the previous arguments that we can find \(\alpha \in (0,1)\) and \(K>0\) such that \(I(t)\leq K\alpha ^{t}\). Thus, substituting in the equation of \(R(t)\), we obtain \(R(t+1)\leq K\alpha ^{t} +\sigma ^{R}(1\gamma ^{R}) R(t)\), which implies that there exist \(\bar{\alpha }\in (\max (\alpha ,\sigma ^{R}(1\gamma ^{R})),1 )\) and \(\bar{K}>0\) such that \(R(t)\leq \bar{K}\bar{\alpha }^{t} \underset{t\to \infty }{\longrightarrow }0\).
Finally, to prove that \(S(t)\underset{t\to \infty }{\longrightarrow }S^{*}\), we just need to follow the reasoning in the proof of Proposition 2. Defining, for \((S_{0},E_{0},I_{0},R_{0})\in \Omega \) and \(t=0,1,2,\ldots \) ,
The corresponding discrete dynamical process \(\tau _{t}\) (\(t\geq 0\)) is asymptotically autonomous with limit discrete semiflow \(\Sigma ^{t}\) (\(t \geq 0\)) generated by the continuous map \(\Sigma (S)=B(S)+\sigma ^{S}S\) in \(\mathbf{R}_{+}\), which coincides with the one in the proof of Proposition 2. □
Note that the assumption \(\Phi ''(x)\leq 0\) for \(x>0\) is met in the case of standard incidence, i.e., if \(\Phi (x)=\beta x\).
The endemicity of the disease is represented in mathematical terms by the concept of uniform persistence. We use the persistence function \(\rho (S,E,I,R)=E+I\). Thus, system (4) is uniformly persistent [25] if there exists \(\varepsilon >0\) such that \(\liminf_{t\to \infty }(E(t)+I(t))>\varepsilon \) for any solution with \(E(0)+I(0)>0\). If lim inf is substituted by lim sup in the definition, the system is said to be uniformly weakly persistent.
In the next theorem we prove the uniform persistence of system (4) when \(\mathcal{R}_{0}>1\) in the case of constant recruitment function \(B(N)=B\) and standard incidence \(\Phi (x)=\beta x\), \(\beta \in [0,1)\).
Theorem 4
Let \(\Phi (x)=\beta x\), \(0<\beta \leq 1\), and \(B(N)=B\) be constant in system (4). If \(\mathcal{R}_{0}>1\), then (4) is uniformly persistent.
Proof of Theorem 4
Since all nonnegative solutions of (4) are attracted by a compact set and \(\{X\in \mathbf{R}_{+}^{4}:E+I>0\}\) is forward invariant, Corollary 4.8 in [25] establishes that it is sufficient to prove that it is uniformly weakly persistent to obtain that it is also uniformly persistent.
So, let us prove that (4) is uniformly weakly persistent. We argue by contradiction. Suppose that it is not. Then, for any arbitrary \(\varepsilon >0\), there exists a solution \(X(t)\) with \(E(0)+I(0)>0\) and \(\limsup_{t\to \infty }(E(t)+I(t))<\varepsilon \). Thus, there exists some \(t_{0}>0\) such that \(E(t)+I(t)<\varepsilon \) for all \(t\geq t_{0}\).
Also, for \(t\geq t_{0}\), \(R(t+1)< \sigma ^{R} (1\gamma ^{R})R(t)+\sigma ^{R}\gamma ^{I} \varepsilon \) with \(0<\sigma ^{R} (1\gamma ^{R})<1\). This implies, iterating the righthand side from \(t_{0}\) on,
Then, as \((\sigma ^{R} (1\gamma ^{R}) )^{tt_{0}} \underset{t\to \infty }{\longrightarrow }0\), there exists \(t_{1}>t_{0}\) such that \(R(t)< (1+ \frac{\sigma ^{R}\gamma ^{I} }{\gamma ^{R} +d^{R}\gamma ^{R} d^{R}} )\varepsilon :=\bar{r}\varepsilon \) for \(t\geq t_{1}\).
Thus, for any \(\varepsilon >0\), there exists \(t_{1}>0\) such that
Let us now establish a lower bound for the total population \(N(t)\). Reproducing, with the appropriate changes, the calculations in the proof of Proposition 1, we obtain the following inequality:
where \(\bar{\sigma }=\min_{C\in \boldsymbol{\mathcal{C}}}\{ \sigma ^{C}\}\in (0,1)\) that implies, considering the difference equation \(x(t+1)=\bar{\sigma }x(t)+B\) and using (8),
This inequality together with (11) yields
With the help of (11) and (12), we can find the following lower bound for the coefficient of \(I(t)\) in the E equation:
Let us define matrix
that satisfies, for \(t\geq t_{1}\),
As \(\bar{P}_{\varepsilon }\) is a primitive matrix, if we can find ε such that \(\rho (\bar{P}_{\varepsilon })>1\), we obtain that \(E(t)\), \(I(t)\underset{t\to \infty }{\longrightarrow }\infty \) whenever \(E(0)+I(0)>0\), which is the contradiction we were looking for. As seen in the proof of Theorem 3, we can check it by means of the net reproductive rate \(\mathcal{R}_{0,\varepsilon }\) of the model associated with matrix \(\bar{P}_{\varepsilon }\) (Theorem 3.3. in [18]).
Now \(\mathcal{R}_{0,\varepsilon }\) and \(\mathcal{R}_{0}\) satisfy
So, if we choose \(\bar{\varepsilon }=\frac{1}{2} (1\frac{1}{\mathcal{R}_{0}} )\frac{\bar{n}}{1+\bar{r}}\), and take ito account that \(\mathcal{R}_{0}>1\), we obtain the required result to complete the proof
□
The model
In this section we present a model of disease dynamics with two time scales. The population can be considered divided into groups that we call patches as if they were forming a sort of generalized metapopulation. The patches could represent real spatial locations and involve explicit movements of the individuals between them [3], but also different individual daily activities (places of residency, work, or business) [14] or different individual behaviors [12]. The key point in considering the generalized patches is that the disease dynamics occurs locally and at a slow time scale compared to the fast time scale associated with the movements between patches.
We consider that individuals move between n patches. In each patch the disease dynamics follows SEIRS model (4) with appropriate local parameters. We assume that movements are almost instantaneous with respect to the disease dynamics. Thus, the model takes the form of a timediscrete two time scale system (1) with movements being the fast process and the disease dynamics the slow process [8].
We denote the densities of susceptible, exposed, infectious, and recovered individuals in patch \(j\in \{1,\ldots ,n\}\) at time \(t\in \{0,1,2,\dots \}\) by \(S_{j}(t)\), \(E_{j}(t)\), \(I_{j}(t)\), and \(R_{j}(t)\), respectively, and the total population \(N_{j}(t)=S_{j}(t)+E_{j}(t)+I_{j}(t)+R_{j}(t)\).
We denote by \(\bar{x}_{C}=\operatorname{col} (C_{1},\ldots ,C_{n} )\in \mathbb{R}^{n\times 1}\) for \(C\in \boldsymbol{\mathcal{C}}\) (\(\boldsymbol{\mathcal{C}}=\{S,E,I,R\}\)) the state vectors of individuals in each compartment (susceptible, exposed, infective, and recovered) across the n patches. The population, or rather metapopulation, state vector is called
The existence of two time scales in the complete model that we propose leads to a reduced model for some global variables. In this case the global variables correspond to the total number of individual in each compartment:
that we collect in the vector of global variables
whose sum yields the total number of individuals in the metapopulation
It is straightforward to see that we can obtain the global variables from the state variables with the help of matrix \(U=\operatorname{diag} (\bar{1},\bar{1},\bar{1},\bar{1} )\in \mathbb{R}^{4\times 4n}_{+}\), where \(\bar{1}=(1,\overset{(n)}{\ldots },1)\in \mathbb{R}^{n}_{+}\) is a row vector
Fast process: movements
We assume that individuals in each compartment move between patches according to movement rates that can generally depend on the global variables Y. In this way, for each disease compartment, movements are represented by a regular stochastic matrix depending on \(Y\in \mathbb{R}^{4}_{+}\)
The movements of the whole metapopulation are then defined through the following matrix:
The state X of the metapopulation after one movement episode is defined by the following map:
that represents the fast process in system (1).
Slow process: disease dynamics
The slow process, the disease dynamics, is defined locally, i.e., in each patch \(j\in \{1,\ldots ,n\}\), by SEIRS model (4):
The recruitment functions \(B_{j}\) verify assumption (3), and the transmission functions \(\Phi _{j}\) assumption (2). All the parameters \(\sigma _{j}^{C}\) and \(\gamma _{j}^{C}\) are in \((0,1)\), \(j\in \{1,\ldots ,n\}\) and \(C\in \boldsymbol{\mathcal{C}}\).
To obtain the map \(\mathcal{S}\) representing the slow process, we need to appropriately reorder equations (14) for all n patches. Following the order of variables in the population vector X, we must include first the equations for variables \(S_{1},\dots ,S_{n}\) and then follow consecutively with those corresponding to compartments E, I, and R.
Finally, the complete two time scale model takes the form of system (1)
Note that
and so map \(\mathcal{F}\) keeps invariant the values of the global variables Y and, therefore, its kth iterate can be expressed in terms of the kpower \(M(Y)^{k}\) of matrix \(M(Y)\)
so that the complete model reads as follows:
Proposition 5
If functions \(B_{j}\) are bounded, \(j\in \{1,\ldots ,n\}\), then system (15) is dissipative.
Proof of Proposition 5
Let us represent with hat \(\hat{C}_{j}(t+1)\) (\(C\in \boldsymbol{\mathcal{C}}\), \(j\in \{1,\ldots ,n\}\)) the state variables after the fast process and before the disease dynamics acts.
We know that \(\sum_{j=1}^{n}\hat{C}_{j}(t+1)=\sum_{j=1}^{n}C_{j}(t)\) for \(C\in \boldsymbol{\mathcal{C}}\).
Moreover, following the proof of Proposition 1, we have, for every \(j\in \{1,\ldots ,n\}\), that
where \(\hat{B}_{j}>0\) is such that \(B_{j}(x)\leq \hat{B}_{j}\) for all \(x\in [0,\infty )\), and \(\hat{\sigma }_{j}=\max_{C\in \boldsymbol{\mathcal{C}}} \{\sigma _{j}^{C}\}\in (0,1)\). Now, calling \(\hat{B}=\hat{B_{1}}+\cdots +\hat{B_{n}}\) and \(\hat{\sigma }=\max_{j\in \{1,\ldots ,n\}}\{ \hat{\sigma }_{j}\}\), we obtain the following recurrent inequality for the total population:
Thus, the proof can be completed as in Proposition 1. □
Notice that all nonnegative solutions of system (15) are attracted by the compact set \(K=\{X\in \mathbf{R}_{+}^{4n}:N\in [0,\bar{B}/(1\bar{\sigma })]\}\).
Reduced system
The fact that the probability matrix \(M^{C}(Y)\) is primitive for every \(C\in \boldsymbol{\mathcal{C}}\) and \(Y\in \mathbb{R}^{4}_{+}\) implies that 1 is its strictly dominant eigenvalue, \(\bar{1}\) is an associated row left eigenvector, and there exists a unique column right eigenvector \(\bar{m}^{C}(Y)\), representing the corresponding stable probability distribution, that satisfies \(\bar{1}\bar{m}^{C}(Y)=1\).
In the reduction procedure of model (15) we need to calculate the limit of the iterates of map \(\mathcal{F}\) that, in this case, is equivalent to calculating the limit of the powers of matrix \(M(Y)\). This latter follows from the Perron–Frobenius theorem:
where \(\bar{M}(Y)=\operatorname{diag} (\bar{m}^{S}(Y),\bar{m}^{E}(Y),\bar{m}^{I}(Y), \bar{m}^{R}(Y) )\). Finally,
Thus, the fourdimensional reduced model to be used to study the asymptotic behavior of the solutions of system (15) is
The next result states how the analysis of the stability of the equilibria of system (17) extends to system (15). It is a direct translation to our setting of Theorem 2 in [20]. Loosely speaking, for k large enough, from the equilibria of the reduced system we can obtain good approximations of the equilibria of system (15) and, in case of asymptotic stability, of their basins of attraction.
Theorem 6
Let \(Y^{*} \in \mathbb{R}_{+}^{4}\) be a hyperbolic equilibrium point of system (17). Then there exists an integer \(k_{0} \geq 0\) such that, for all \(k \geq k_{0}\), system (15) has an equilibrium point \(X_{k}^{*}\) which is hyperbolic and satisfies
Moreover, the following hold:

(i)
\(X_{k}^{*}\) is asymptotically stable (resp. unstable) if and only if \(Y^{*}\) is asymptotically stable (resp. unstable).

(ii)
Let \(Y^{*}\) be asymptotically stable, and let \(X_{0} \in \mathbb{R}_{+}^{4n}\) be such that the solution \(\{ Y(t) \} _{t=0, 1, \ldots } \) to (17) corresponding to the initial data \(Y_{0} := U X_{0}\) satisfies \(\lim_{t \to \infty } Y(t) = Y^{*}\). Then, for all \(k \geq k_{0}\), the solution to (15) \(\{ X_{k}(t) \} _{t=0, 1, \ldots }\) with \(X_{k} (0)= X_{0}\) satisfies \(\lim_{t \to \infty } X_{k} (t)= X^{*}_{k}\).
Metapopulation SEIRS epidemic model with constant recruitment, standard incidence, and constant movement rates
We consider, for each patch \(j\in \{1,\ldots ,n\}\), a local SEIRS model (14) assuming constant recruitment function and standard incidence for disease transmission.
Parameters satisfy \(B_{j}>0\), \(\beta _{j}\in (0,1]\) and \(\sigma _{j}^{C},\gamma _{j}^{C}\in (0,1)\), \(j\in \{1,\ldots ,n\}\) and \(C\in \boldsymbol{\mathcal{C}}\).
The associated local basic reproduction number,
rules the disease dynamics. If \(\mathcal{R}^{j}_{0}<1\), then (Theorem 3) the DFE \(X^{j}_{0}=(B_{j}/(1\sigma _{j}^{S}),0,0,0)\) is GAS (disease eradication), and if \(\mathcal{R}^{j}_{0}>1\), then (Theorem 4) the system is uniformly persistent (disease endemicity).
We assume that movement rates are constant. The constant regular stochastic matrices describing the movements in each compartment are \(M^{S},M^{E},M^{I},M^{R}\in \mathbb{R}^{n\times n}_{+}\). Let \(M=\operatorname{diag} (M^{S},M^{E},M^{I},M^{R} )\in \mathbb{R}^{4n \times 4n}_{+}\) be the matrix defining the whole fast process \(F(X)=MX\).
The complete two time scale model takes the form of system (15)
where \(\mathcal{S}\) is the map representing the disease dynamics.
Reduced system
We follow the procedure described in Sect. 3.1. For every \(C\in \boldsymbol{\mathcal{C}}\), we denote \(\bar{m}^{C}=(m^{C}_{j})_{j\in \{1,\ldots ,n\}}\) the column right eigenvector of matrix \(M^{C}\) associated with eigenvalue 1 that satisfies \(\bar{1}\bar{m}^{C}=1\), and \(\bar{M}=\operatorname{diag} (\bar{m}^{S},\bar{m}^{E},\bar{m}^{I},\bar{m}^{R} )\). Thus, the reduced system associated with system (20) can be expressed in the form of (17) as
that corresponds to the following SEIRS model:
where the parameters are weighted means of the corresponding local parameters, with the weights being the elements of the stable probability distributions associated with the movements process:
System (21) is similar to the local system (18), the main difference between them being the transmission terms \(\bar{\beta }_{S}\) and \(\bar{\beta }_{E}\). Its recruitment constant is the sum of the local ones. Coefficients \(\delta _{R}^{S}\), \(\delta _{S}^{S}\), \(\delta _{E}^{E}\), \(\delta _{E}^{I}\), \(\delta _{I}^{I}\), \(\delta _{I}^{R}\), and \(\delta _{R}^{R}\) are all in \((0,1)\). The same happens to \(\bar{\beta }_{S}(Y(t))I(t)\) and \(\bar{\beta }_{E}(Y(t))I(t)\). In fact,
and \(\bar{\beta }_{E}(Y(t))I(t)\leq \sum^{n}_{j=1}\sigma _{j}^{E} \beta _{j} m^{S}_{j}\leq \max_{j}\{\sigma _{j}^{E}\beta _{j}\}<1\).
We now proceed to briefly analyze system (21) following the steps carried out in Sect. 2 with system (4). As they are both very similar, we only fall into the details in case of significant difference.
All parameters are positive. Also, in the S equation, \(\delta _{S}^{S}\bar{\beta }_{S}(Y(t))I(t)\) is positive for any \(Y(t)\in \mathbf{R}_{+}^{4}\). Therefore, \(\mathbf{R}_{+}^{4}\) is forward invariant under the semiflow defined by the discretetime system (21).
Proposition 7
System (21) is dissipative.
Proof of Proposition 7
Let us define \(\hat{\sigma }=\max_{C\in \boldsymbol{\mathcal{C}},j \in \{1,\ldots ,n\}}\{\sigma ^{C}_{j}\}\in (0,1)\).
Summing up the four equations of system (21)
The coefficients of \(S(t)\), \(E(t)\), \(I(t)\), and \(R(t)\) are all bounded by σ̂. Let us show it in the case of \(S(t)\). The rest are analogous. Using the fact that the sum of a convex combination is less than or equal to the maximum of the summands, we easily obtain
Therefore, the total population verifies the following recurrent inequality:
and the proof can be completed as in Proposition 1. □
Note from the proof that, as in the case of system (4), all the nonnegative solutions of system (21) are attracted by the compact set \(\{Y\in \mathbf{R}_{+}^{4}:N\in [0,\bar{B}/(1\hat{\sigma })]\}\).
The associated diseasefree system
is linear and, thus, it is straightforward to prove that it possesses the GAS equilibrium \((\bar{B}/(1\delta _{S}^{S}),0)\).
The unique DFE of system (21) is \(Y^{*}_{0}=(\bar{B}/(1\delta _{S}^{S}),0,0,0)\).
The nextgeneration method [2] allows us to calculate its basic reproduction number
where \(\bar{\beta }_{I}=\sum^{n}_{j=1} \sigma _{j}^{E} \beta _{j} m^{I}_{j}\in (0,1)\).
The next result characterizes the eradication/endemicity of the disease in terms of \(\overline{\mathcal{R}}_{0}\).
Theorem 8
Consider system (21). Then

(a)
If \(\overline{\mathcal{R}}_{0}<1\), then DFE \(Y_{0}^{*}\) is GAS.

(b)
If \(\overline{\mathcal{R}}_{0}>1\), then DFE \(Y_{0}^{*}\) is unstable and the system is uniformly persistent.
Proof of Theorem 8

(a)
Follow the proof of Theorem 3 using that
$$ \begin{aligned} E(t+1) & = \delta _{E}^{E} E(t) + \Biggl(\sum^{n}_{j=1}\sigma _{j}^{E} \frac{\beta _{j}m^{I}_{j}m^{S}_{j}S(t)}{ \sum_{C\in \boldsymbol{\mathcal{C}}}m^{C}_{j} C(t)} \Biggr)I(t) \\ & \leq \delta _{E}^{E} E(t) + \Biggl(\sum ^{n}_{j=1}\sigma _{j}^{E} \beta _{j}m^{I}_{j} \Biggr)I(t)= \delta _{E}^{E} E(t) +\bar{\beta }_{I} I(t). \end{aligned} $$ 
(b)
As in the proof of Theorem 4, we prove, arguing by contradiction, that (21) is uniformly weakly persistent. We omit the details that are very similar in both proofs.
For any \(\varepsilon >0\), there exist \(t_{1}>0\) and \(\bar{r}>0\) such that
$$ E(t)+I(t)+R(t)\leq (1+\bar{r})\varepsilon\quad \text{for } t\geq t_{1}. $$(24)We have the inequality \(N(t+1)\geq \hat{\sigma }N(t)+\bar{B}\), where \(\hat{\sigma }=\min_{C\in \boldsymbol{\mathcal{C}},j \in \{1,\ldots ,n\}}\{\sigma ^{C}_{j}\}\), \(\hat{\sigma }\in (0,1)\), which implies that \(N(t)\geq \bar{n}:=\min \{N(0),B/(1\hat{\sigma })\}>0\).
Thus, together with (24), this yields
$$ S(t)\geq \bar{n}(1+\bar{r})\varepsilon \quad \text{for } t\geq t_{1}. $$(25)A lower bound for the coefficient of I(t) in the E equation, calling
\(\bar{m}=\max \{m^{E}_{j},m^{I}_{j},m^{R}_{j}\}\), is the following:
$$ \bar{\beta }_{E} \bigl(Y(t) \bigr)S(t) \geq G(\varepsilon ):= \sum ^{n}_{j=1} \frac{\sigma ^{E}_{j}\beta _{j}m^{I}_{j}m^{S}_{j} (\bar{n}(1+\bar{r})\varepsilon )}{m^{S}_{j}(\bar{n}(1+\bar{r})\varepsilon )+\bar{m}(1+\bar{r})\varepsilon }. $$Matrix
$${\overline{P}}_{\epsilon}=\left[\begin{array}{cc}{\delta}_{E}^{E}& G(\epsilon )\\ {\delta}_{E}^{I}& {\delta}_{I}^{I}\end{array}\right]$$satisfies, for \(t\geq t_{1}\),
$$\left[\begin{array}{c}E(t+1)\\ I(t+1)\end{array}\right]\ge {\overline{P}}_{\epsilon}\left[\begin{array}{c}E(t)\\ I(t)\end{array}\right].$$Therefore, to get the required contradiction and complete the proof, it is enough to find ε such that the net reproductive rate \(\overline{\mathcal{R}}_{0,\varepsilon }\) of \(\bar{P}_{\varepsilon }\) is larger than 1. The facts that \(\overline{\mathcal{R}}_{0,\varepsilon }= \frac{G(\varepsilon )}{\bar{\beta }_{I}}\overline{\mathcal{R}}_{0}\), \(\overline{\mathcal{R}}_{0}>1\), and \(G(\varepsilon )\) is a continuous decreasing function with \(\lim_{\varepsilon \to 0^{+}}G(\varepsilon )=\bar{\beta }_{I}\) justify the existence of the required ε.
□
A direct consequence of Theorems 6 and 8 is the following result, which we apply in the next section.
Corollary 9
If \(\overline{\mathcal{R}}_{0}<1\), there exists an integer \(k_{0} \geq 0\) such that for all \(k \geq k_{0}\) system (20) has an equilibrium point \(X_{k}^{*}\) which is GAS and satisfies
We note that \(X^{\ast }\) is a DFE. Indeed:
where \((\bar{x}^{*}_{S} )_{j}= \frac{\bar{B}m_{j}^{S}}{1\delta _{S}^{S}}\), \(j=1,2,\ldots ,n\). Therefore, Corollary 9 states that a sufficient condition for the global disease eradication in system (20) is that the basic reproduction number of the reduced system \(\overline{\mathcal{R}}_{0}\) is less than one.
Results
To illustrate the use of the developed framework, we explore the possibility of eradication of the disease through appropriate movements when it is endemic in every isolated patch.
We consider system (20). The disease is endemic in each patch \(j\in \{1,\ldots ,n\}\) if the associated local basic reproduction numbers satisfy
We want to find conditions such that the appropriate movements of exposed and infectious individuals drive the disease to eradication. We are looking for conditions on parameters \(\sigma _{j}^{E}\), \(\sigma _{j}^{I}\), \(\gamma _{j}^{E}\), \(\gamma _{j}^{I}\), and \(\beta _{j}\), for which there exist stable equilibrium proportions of the distribution of exposed and infectious individuals, \(m^{E}_{j}\) and \(m^{I}_{j}\), such that
We make some simplifying assumptions to deal with the large number of parameters. Among the five parameters involved in the expression of \(\mathcal{R}^{j}_{0}\), we consider three of them, \(\sigma ^{E}\), \(\gamma ^{E}\) and β, to be the same in all patches, that is,
Calling
we obtain the following simplified expressions for the basic reproduction numbers:
Considering \(\mathcal{R}^{j}_{0}\) as a function of any of the five parameters involved in its expression, it can be checked that it is monotone on \([0,1]\). Keeping four out of these five parameters constant across patches would make it impossible to find conditions to obtain \(\overline{\mathcal{R}}_{0}<1\), since in that case
Another simplification leading to the same negative conclusion is assuming coefficients \(m^{E}_{j}\) and \(m^{I}_{j}\) to be equal, i.e., \(m^{E}_{j}=m^{I}_{j}=m_{j}\), since in this case \(\mathcal{R}^{j}_{0}>1\), \(j\in \{1,\ldots ,n\}\), implies \(A\sigma _{j}^{I}>1\sigma _{j}^{I}(1\gamma _{j}^{I})\), and
which yields again \(\overline{\mathcal{R}}_{0}>1\). Therefore, to obtain \(\mathcal{R}^{j}_{0}>1\) and \(\overline{\mathcal{R}}_{0}<1\), we need to assume that movements lead exposed and infectious individuals to different distributions among patches.
We analyze this situation, local endemicity in isolated patches with global eradication through appropriate movements, in a simple twopatch metapopulation.
We have
and, calling \(x=m^{E}_{1}\) and \(y=m^{I}_{1}\), the expression for \(\overline{\mathcal{R}}_{0}\) becomes
where x and y, the fractions of exposed and infectious individuals in patch 1, take values on \((0,1)\).
Expression A, depending on the values of parameters \(\sigma ^{E},\gamma ^{E},\beta \in (0,1)\), can also take any value on \((0,1)\). On the other hand, the expression \(\sigma _{j}^{I}/ (1\sigma _{j}^{I}(1\gamma _{j}^{I} ))\), depending on \(\sigma ^{I}_{j},\gamma ^{I}_{j}\in (0,1)\), can take any positive value. Therefore, independently of the value of A, it is always feasible to have \(\mathcal{R}^{1}_{0}>1\) and \(\mathcal{R}^{2}_{0}>1\). Without loss of generality we can assume
or, equivalently,
Now, we look for the condition that ensures \(\overline{\mathcal{R}}_{0}<1\), that is,
where x and y can be chosen in \((0,1)\). This condition can be easily expressed as
This minimum is attained at the boundary of \([0,1]\times [0,1]\) because \(\partial g/\partial x\) never changes sign. In the four sides of the boundary g becomes a monotone function of either variable x or y, and therefore the minimum must be in one of the four corners. We have already assumed that \(g(0,0)=\mathcal{R}_{0}^{2}/A>1/A\) and \(g(1,1)=\mathcal{R}_{0}^{1}/A>1/A\). So, to get \(\overline{\mathcal{R}}_{0}<1\), one of the following two inequalities must hold: either \(g(1,0)<1/A\) or \(g(0,1)<1/A\). Both cannot hold, since the first one together with \(g(0,0)>1/A\) implies that \(\sigma _{1}^{I}<\sigma _{1}^{I}\), whereas the second one and \(g(1,1)>1/A\) lead to the opposite inequality.
Therefore if inequalities (29) hold and \(g(1,0)<1/A\), that is,
then the values \((m^{E}_{1},m^{I}_{1})\in (0,1)\times (0,1)\) under the line
correspond to movements leading to the global eradication of the disease (see Fig. 3).
Analogously, if inequalities (29) hold and \(g(0,1)<1/A\), that is,
then the values \((m^{E}_{1},m^{I}_{1})\in (0,1)\times (0,1)\) over the line \(\overline{\mathcal{R}}_{0}=1\) correspond to movements leading to the global eradication of the disease (see Fig. 4).
Discussion
Contrary to what happens in continuous time [5, 16, 17, 24, 28], in the literature there are almost no references on discretetime epidemic models with two time scales. The authors did present some of them [9–11] to analyze the influence of a parasite on the dynamics of a community. In them, the effect of the parasite is represented by an SIS model, and it is considered fast in comparison to the community dynamics.
This work aims to establish a framework applicable to discretetime models similar to the one existing in continuous time. However, it cannot be as flexible and powerful, for including two time scales in a discrete model has to be done sequentially and that necessarily imposes constrains. Moreover, the analysis of these models requires the reduction of their dimension, and the results at disposal to relate the behavior of the original and the reduced models are scarce [22, 23].
The proposed framework consists in a disease process acting on the slow time scale coupled with a second process whose dynamics is fast in comparison. To show the capabilities of this framework, we have tried to be fairly general in our choice of both processes. A new SEIRS discretetime model is presented and analyzed. This model is chosen because together with its subcases (SI, SIS, SEI, SEIS, SIR, SIRS,…) it encompasses a rather large number of commonly used models. Anyway, changing the local disease model in our framework is straightforward. The second process is presented under the form of generalized individual movements between patches of a generalized metapopulation. This is the way of admitting processes as general as possible in such a way that they lead to two time scale systems susceptible of being reduced by the methods mentioned above.
The analysis of SEIRS model (4) characterizes the eradication or endemicity of the disease. Its \(\mathcal{R}_{0}\) (10) is calculated by the nextgeneration method and, together with some other mild assumptions, if it is less than 1, then its DFE is GAS. On the other hand, if \(\mathcal{R}_{0}>1\), the disease becomes endemic, i.e., the model is uniformly persistent in its infected compartments. When analyzing the full model with two time scales (20), the associated reduced model (21) is an SEIRS model similar to (4). Its basic reproduction number \(\overline{\mathcal{R}}_{0}\) (23) serves to decide whether the EE is GAS or the system is persistent on infection. The main reduction result that is applied, Theorem 6, allows us to transfer the results on the asymptotic behavior of the system having to do with equilibria from the reduced system to the complete system for k large enough. Thus, \(\overline{\mathcal{R}}_{0}<1\) ensures the eradication of the disease in the complete model. This makes \(\overline{\mathcal{R}}_{0}\) play the role of the basic reproduction number of model (20) with a good approximation. Therefore, the presented framework provides a method to approximate in a simpler form the basic reproduction number of a complex system.
In the general case of the complete model (15), if the reduced model (17) were to admit an endemic equilibrium (EE), possibly GAS, indicating the persistence of the disease, then Theorem 6 would apply. Thus, we could state the following corollary: If \(Y_{1}^{*}\) is an EE of the reduced model (17), then, for k large enough, the complete model (15) possesses an EE close to \(X_{1}^{\ast }:= \mathcal{S}(\bar{M} (Y_{1}^{*} )Y_{1}^{*})\), and if \(Y_{1}^{*}\) is GAS so is \(X_{1}^{\ast }\). An issue that it is not covered by Theorem 6, and we will address it in a future work, has to do with the property of uniform persistence. If the reduced model (17) is uniformly persistent, Theorem 6 gives no information on the persistence of the complete model (15).
The general framework of Sect. 3 is treated in Sect. 4 for the case of constant recruitment, standard incidence, and constant movement rates. This is done to present a final pertinent illustration of the application of the process. We address the situation of a disease that is endemic in each isolated patch. Conditions on disease local parameters are found that enable driving the disease to extinction globally through a set of appropriate movements. Certainly, this application is only a minimal example of the large number of relevant aspects of disease dynamics that can be treated on similar terms using our approach.
Availability of data and materials
Not applicable.
Abbreviations
 DFE:

Diseasefree equilibrium
 EE:

Endemic equilibrium
 GAS:

Globally asymptotically stable
 LAS:

Locally asymptotically stable
 SEIRS:

Susceptible–Exposed–Infectious–RecoveredSusceptible
References
 1.
Allen, L.J.: Some discretetime SI, SIR, and SIS epidemic models. Math. Biosci. 124, 83–105 (1994)
 2.
Allen, L.J., van den Driessche, P.: The basic reproduction number in some discretetime epidemic models. J. Differ. Equ. Appl. 14, 1127–1147 (2008)
 3.
Arino, J.: Diseases in metapopulations. In: Ma, Z., Zhou, Y., Wu, J. (eds.) Modeling and Dynamics of Infectious Diseases, pp. 65–123. World Scientific, New Jersey (2009)
 4.
Brauer, F.: Mathematical epidemiology is not an oxymoron. BMC Public Health 9, S2 (2009)
 5.
Brauer, F.: A singular perturbation approach to epidemics of vectortransmitted diseases. Infect. Dis. Model. 4, 115–123 (2019)
 6.
Brauer, F., CastilloChávez, C., Feng, Z.: Discrete epidemic models. Math. Biosci. Eng. 7, 1–15 (2010)
 7.
Brauer, F., CastilloChávez, C., Feng, Z.: Mathematical Models in Epidemiology. Springer, New York (2019)
 8.
Bravo de la Parra, R., Marvá, M., Sánchez, E., Sanz, L.: Reduction of discrete dynamical systems with applications to dynamics population models. Math. Model. Nat. Phenom. 8, 107–129 (2013)
 9.
Bravo de la Parra, R., Marvá, M., Sánchez, E., Sanz, L.: A discrete predatorprey ecoepidemic model. Math. Model. Nat. Phenom. 12, 116–132 (2017)
 10.
Bravo de la Parra, R., Marvá, M., Sánchez, E., Sanz, L.: Discrete models of disease and competition. Discrete Dyn. Nat. Soc. 2017, Article ID 531083 (2017)
 11.
Bravo de la Parra, R., Sanz, L.: A discrete model of competing species sharing a parasite. Discrete Contin. Dyn. Syst., Ser. B 25, 2121–2142 (2020)
 12.
CastilloChávez, C., Bichara, B., Morin, B.R.: Perspectives on the role of mobility, behavior, and time scales in the spread of diseases. Proc. Natl. Acad. Sci. 113, 14582–14588 (2016)
 13.
CastilloChávez, C., Yakubu, A.A.: Dispersal, disease and lifehistory evolution. Math. Biosci. 173, 35–53 (2001)
 14.
FalcónLezama, J.A., MartínezVega, R.A., KuriMorales, P.A., RamosCastaneda, J., Adams, B.: Daytoday population movement and the management of Dengue epidemics. Bull. Math. Biol. 78, 2011–2033 (2016). https://doi.org/10.1007/s1153801602096
 15.
Getz, W.M., Salter, R., Muellerklein, O., Yoon, H.S., Tallam, K.: Modeling epidemics: a primer and numerus model builder implementation. Epidemics 25, 9–19 (2018)
 16.
JardónKojakhmetov, H., Kuehn, C., Pugliese, A., Sensi, M.: A geometric analysis of the SIR, SIRS and SIRWS epidemiological models. Nonlinear Anal., Real World Appl. 58, 103220 (2021)
 17.
JardónKojakhmetov, H., Kuehn, C., Pugliese, A., Sensi, M.: A geometric analysis of the SIRS epidemiological model on a homogeneous network. J. Math. Biol. 83, 37 (2021). https://doi.org/10.1007/s00285021016645
 18.
Li, C.K., Schneider, H.: Applications of Perron–Frobenius theory to population dynamics. J. Math. Biol. 44, 450–462 (2002)
 19.
Martcheva, M.: An Introduction to Mathematical Epidemiology. Springer, New York (2015)
 20.
Marvá, M., Sánchez, E., Bravo de la Parra, R., Sanz, L.: Reduction of slow–fast discrete models coupling migration and demography. J. Theor. Biol. 258, 371–379 (2009)
 21.
Rocha, F., Mateus, L., Skwara, U., Aguiar, M., Stollenwerk, N.: Understanding Dengue fever dynamics: a study of seasonality in vectorborne disease models. Int. J. Comput. Math. 93, 1405–1422 (2016)
 22.
Sanz, L., Bravo de la Parra, R., Marvá, M., Sánchez, E.: Nonlinear population discrete models with two time scales: rescaling of part of the slow process. Adv. Differ. Equ. 2019, 401 (2019)
 23.
Sanz, L., Bravo de la Parra, R., Sánchez, E.: Approximate reduction of nonlinear discrete models with two time scales. J. Differ. Equ. Appl. 14, 607–627 (2008)
 24.
Schecter, S.: Geometric singular perturbation theory analysis of an epidemic model with spontaneous human behavioral change. J. Math. Biol. 82, 54 (2021). https://doi.org/10.1007/s00285021016052
 25.
Smith, H.L., Thieme, H.R.: Dynamical Systems and Population Persistence. Am. Math. Soc., Providence (2011)
 26.
Souza, M.O.: Multiscale analysis for a vectorborne epidemic model. J. Math. Biol. 68, 1269–1293 (2014)
 27.
van den Driessche, P., Yakubu, A.A.: Disease extinction versus persistence in discretetime epidemic models. Bull. Math. Biol. 81, 4412–4446 (2019)
 28.
Wang, X., Wei, L., Zhang, J.: Dynamical analysis and perturbation solution of an SEIR epidemic model. Appl. Math. Comput. 232, 479–486 (2014). https://doi.org/10.1016/j.amc.2014.01.090
 29.
Yakubu, A.A.: Introduction to discretetime epidemic models. In: Gumel, A.B., Lenhart, S. (eds.) Modeling Paradigms and Analysis of Disease Transmission Models, pp. 83–109. Am. Math. Soc., Providence (2010)
 30.
Zhao, X.Q.: Asymptotic behavior for asymptotically periodic semiflows with applications. Commun. Appl. Nonlinear Anal. 3, 43–66 (1996)
Acknowledgements
The authors would like to thank an anonymous referee for their suggestions to improve this work.
Funding
Authors are supported by Ministerio de Ciencia e Innovación (Spain), Project PID2020114814GBI00.
Author information
Affiliations
Contributions
Both authors have the same contributions. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bravo de la Parra, R., SanzLorenzo, L. Discrete epidemic models with two time scales. Adv Differ Equ 2021, 478 (2021). https://doi.org/10.1186/s13662021036330
Received:
Accepted:
Published:
Keywords
 Discretetime epidemic model
 Time scales
 Disease eradication or persistence