key: cord-0802579-m64brxkb authors: Naim, Mouhcine; Lahmidi, Fouad; Namir, Abdelwahed title: Threshold Dynamics of an [Formula: see text] Epidemic Model with Nonlinear Incidence Rates date: 2021-09-15 journal: Differ Equ Dyn Syst DOI: 10.1007/s12591-021-00581-9 sha: 0f32774e07c1e90a9205eab27196adbd27896e1c doc_id: 802579 cord_uid: m64brxkb In this paper, we consider an SEIS epidemic model with infectious force in latent and infected period, which incorporates by nonlinear incidence rates. The local stability of the equilibria is discussed. By means of Lyapunov functionals and LaSalle’s invariance principle, we proved the global asymptotic stability of the disease-free equilibrium and the endemic equilibrium. An application is given and numerical simulation results based on real data of COVID-19 in Morocco are performed to justify theoretical findings. In our life there are various infectious diseases that we face, such as coronavirus that was first identified in the Chinese city of Wuhan in 2019, and it is a new strain that has not been previously identified in humans. COVID-19 is a fatal illness caused by a severe acute respiratory syndrome that is transmitted from human-to-human [1] . In epidemiology, mathematical models describing the population dynamics of infectious diseases plays an important role in analyzing the spread and controlling these diseases. Some diseases, such as tuberculosis, measles, AIDS, COVID-19, have a latent or incubation period when a susceptible has become infected but not yet infectious. COVID-19, for example, has about 5-6 days latent period [2] . This latency can be modeled by introducing a new class, called the exposed class, in which the susceptible remains for a given length of time before moving into the infective class. The resulting models are of SEI, SEIS, SEIR, SEIRS, or SEIRI type (here, S, E, I and R denote the numbers of susceptible, exposed, infectious and recovered individuals, respectively), see, e.g., [3] [4] [5] [6] [7] [8] . In many diseases the infection is not only in the infected period but even in the latent period, such as COVID-19 at present [9] . Many epidemic models with the infectious force in the latent period have been performed. Yuang et al. in [10] considered the local stability of an epidemic model having infectious force in both latent period and infected period. Guihua and Zhen [11, 12] studied the global stability of an SEI epidemic model having infectious force in both latent period and infected period. Global stability of an SEIS epidemic model with standard incidence rate having infectious force in both latent period and infected period given by Mukhopadhyay and Bhattacharyya [13] . The method used in [11] [12] [13] for proving global asymptotic stability of epidemic equilibrium is based on the geometric approaches developed by Smith [14] , Li and Muldowney [15] . The incidence of a disease is the number of new cases per unit time and plays an important role in the modelling of epidemic dynamics. In many epidemic models, the bilinear incidence rate SI ( denote the transmission rate, which measures the infection force of the disease) [16, 17] and the standard incidence rate SI N , where N is the total population size [18] , are frequently used. However, the bilinear incidence rate may face many problems and challenges when it is used to describe the evolution of infectious diseases among human. Consequently, the bilinear incidence rate should be modified into a nonlinear incidence rate which play a vital role in ensuring that the model can give a reasonable qualitative description of the capture the nonlinear nature of disease propagation law. Initially, Capasso and Serio [19] in the study of the cholera epidemic spread in Bari in 1973, observed that the incidence rate may increase more slowly as I increases rather than linearly increasing, and found the saturated incidence rate SI 1+ I which gives good agreement with observations, where the parameter measures the inhibitory effect either due to the crowding caused by high numbers of cases or from the behavioral change of susceptible individuals as their number increases. From then on, a number of authors have studied epidemiological models with several nonlinear incidence rates, see, e.g., [20] [21] [22] [23] [24] [25] [26] . In this paper we will consider an SEIS epidemic model with general incidence rates that the diseases such as COVID-19 can be infected in the latent period and the infected period. Hence, our study will be carried out on the following epidemic model: where Λ is the recruitment rate of the susceptible, is the natural death rate of the population, is the transfer rate from the exposed group to the infectious group so that 1 is the mean latent period, is the death rate due to disease, is the transfer rate from the infectious group to the susceptible group. In this work, we assume that and are nonnegative constants, and Λ , , and are positive constants. To study the dynamics of the model (1), we suppose that the functions f i (i = 1, 2) are continuously differentiable on [0, +∞) satisfying the following two hypotheses. (H1) f i (x) is strictly monotone increasing on [0, +∞) with f i (0) = 0. (H2) x is monotone decreasing on (0, +∞). (1) The hypothesis (H1) is a basic requirement for any biologically feasible incidence rate, since it is obvious that the disease cannot spread when the number of exposed and infected individuals is zero. In addition, the incidence Sf 1 (I) (resp., Sf 2 (E) ) becomes faster as the number of the infectious (resp., exposed) individuals increases. Since the incidence Sf 1 (I) (resp., Sf 2 (E) ) can be rewritten as I SI (resp., . Then, function E ) refers to effective contact rate between an infective (resp., exposed) and a susceptible individuals. The hypothesis (H2) implies that the effective contact rate between infective (resp., exposed) and susceptible individuals is decreasing with increase of infectives since the individuals tend to reduce the number of contacts among them per unit time when they are under intervention policies. The idea of hypothesis (H2) has first been explored in [19] . By the hypotheses (H1)-(H2), it can be seen that f i satisfying and consequently, the effective contact rates are bounded. The functions f i (i = 1, 2) includes a number of especial incidence rates. For instance: On the positivity and boundedness of solutions, we have the following result. Since the proof is simple, we here omit it. The rest of this paper is organized as follows. In the next section, we discuss the existence and the local stability of the equilibria. In Section 3, by using Lyapunov functionals and LaSalle's invariance principle, we establish the global asymptotic stability of the equilibria. An application and illustrative numerical simulations have been given in Section 4, and the numerical results of real data of COVID-19 in Morocco is conducted in the same section. Finally, conclusion is given in Section 5. In this section, we discuss the existence and the local stability of the equilibria of system (1). It is easy to visualize that the system (1) always has a disease-free equilibrium E 0 = (S 0 , 0, 0) , where S 0 = Λ , that is, there is no infection present in the population and all individuals are susceptible. By applying the next generation matrix approach provided by van den Driessche and Watmough [30] , the basic reproduction number of model (1) is defined as follows , where 1 = + and 2 = + + . Biologically, R 0 represents the average number of secondary infectious cases produced by an exposed individual and an infectious individual during their effective infectious period when introduced in a completely susceptible population [31] . Now, we obtain the following lemma which gives the result on the unique existence of the endemic equilibrium E * for R 0 > 1. Proof Assume that R 0 > 1 . Suppose that (S, E, I) is any positive equilibrium of system (1), then From the third equation of system (3), we have and by summing all equations of system (3), we find that where = 2 + + . We have S ≥ 0 implies that I ∈ 0, Λ . Hence, no positive equilibrium exists if I ≥ Λ . On the other hand, by the second equation of (3) and equation (4), it follows that Then, by equations (5) and (6), I is a positive zero of the function F, where By the hypothesis (H2), the function F is strictly monotone decreasing on 0, Λ satisfying and F( Λ ) < 0 . Hence, there is a unique I * ∈ 0, Λ such that F(I * ) = 0 if R 0 > 1 . This shows that model (1) has a unique positive endemic equilibrium E * = (S * , E * , I * ) when where E * and S * are determined by (4) and (5), respectively. Second we assume R 0 ≤ 1 . Then it is obvious that system (1) has no endemic equilibria. Hence the proof is complete. ◻ Next, we focus on the local stability of the disease-free equilibrium E 0 . The Jacobian matrix of system (1) at the equilibrium E 0 is as follows The other two eigenvalues of J E 0 are determined by the following equation , and this suggests that J E 0 admits a positive real eigenvalue. Consequently, we have the following result. Now, we analyze the local stability of the endemic equilibrium of the system (1) by assuming that R 0 > 1 . The Jacobian matrix of system (1) at the equilibrium E * is as follows Therefore, its characteristic equation at E * is where By (H2), x is decreasing on (0, +∞) , so Then it is easy to show that C 2 > 0, C 1 > 0 , C 0 > 0 and C 2 C 1 > C 0 . Thus, by the Routh-Hurwitz criterion, all roots j (j = 1, 2, 3) of (7) have negative real part. Hence, we have the following theorem. Theorem 2 When R 0 > 1 , the endemic equilibrium E * is locally asymptotically stable. The aim of this section is to analyze the global stability of the equilibria. First, we consider the global stability of model (1) at the disease-free equilibrium E 0 . According to (2), we have Therefore, R 0 ≤ 1 ensures that dU dt ≤ 0 . Furthermore, it is easy to verify that the singleton E 0 is the largest compact invariant set in (S, E, I) ∈ ℝ 3 + ∶ dU dt = 0 . By LaSalle's invariance principle [32] , we get that equilibrium E 0 of system (1) is globally asymptotically stable if R 0 ≤ 1. ◻ In the following theorem, we study the global stability of the endemic equilibrium E * . Proof Let V be the Lyapunov functional defined as where and By means of the equalities the derivative of V 0 is The derivative of V 1 is Then the sum of these two functions gives By the hypotheses (H1) and (H2), we obtain and using the arithmetic-geometric inequality, we have Hence, the condition (8) ensures that dV dt ≤ 0 . Further, the largest invariant set of (S, E, I) ∈ ℝ 3 + ∶ dV dt = 0 is the singleton E * . By applying the LaSalle's invariance principle, we can obtain that the endemic equilibrium E * of model (1) is globally asymptotically stable if R 0 > 1 and S * − I * ≥ 0. ◻ Remark 1 When = 0 , the system (1) reduces to an SEI model, and the endemic equilibrium E * is globally asymptotically stable if R 0 > 1. For > 0 , the endemic equilibrium E * is globally asymptotically stable if where given in the proof of Lemma 1. Proof We prove that the condition (8) in Theorem 4 holds if the condition R 0 ≤ 1 + holds. We have I * satisfies the following Therefore, we obtain that which completes the proof. ◻ The main purpose of this section is to give a particular nonlinear incidence rate to illustrate our main results. Consider the system 1 < R 0 ≤ 1 + , subject to initial conditions where 1 and 2 are the rates of the efficient contact in infected period and latent period, respectively. The positive constants 1 and 2 measure the saturation or inhibitory effect. The remaining parameters have the same biological significance as those in system (1). Model (9) is a particular case of system (1) by choosing f 1 (I) = 1 I 1+ 1 I and f 2 (E) = 2 E 1+ 2 E . Obviously, the hypotheses (H1)-(H2) are satisfied. The basic reproduction number of model (9) is Now, we give some numerical simulations in order to validate our theoretical results. Firstly, we simulate model (9) with the parameter values as shown in Table 1 . By calculation, we have R 0 = 0.5919 ≤ 1 . Hence, by Lemma 1, model (9) has only an disease-free equilibrium E 0 = (1384.45, 0, 0) . Fig. 1 shows that the trajectories of S(t), E(t) and I(t), which are the solutions of (9), from different initial values to E 0 . This supports the analytical result obtained in Theorem 3, i.e., the disease-free equilibrium E 0 is globally asymptotically stable provided R 0 ≤ 1 . Thus, under the condition R 0 ≤ 1 , the disease will be eradicated. Now, we simulate the case when the basic reproduction number is bigger than one. For this end, we consider model (9) with the parameter values as shown in Table 2 . By computing, we get R 0 = 4.2787 > 1 and S * − I * = 97.0392 ≥ 0 . Fig. 2 shows that the trajectories of S(t), E(t) and I(t), which are the solutions of (9), from different initial values to E * = (1192.7, 29.0039, 18.6941) . This confirms the global stability of (9) . Fig. 3 shows that the trajectories of S(t), E(t) and I(t) in the model under investigation (9), for the parameter values as shown in Table 3 , converges towards the endemic equilibrium E * = (1.9007 × 10 7 , 1.8539 × 10 6 , 2.3607 × 10 6 ). In this case, we have the basic reproduction number is greater than unity R 0 = 1.5159 > 1 and S * − I * ≥ 0 , which means that the disease persists. In this situation, it will be important to eventually undertake some strategies like quarantine, isolation, wearing of masks, disinfection and if it will be possible vaccination. This paper presents a mathematical study on the dynamical behavior of an SEIS epidemic model with infectivity in incubation period and general incidence rates. From the system, we have found the following basic reproduction number Table 1 . This plot shows the stability of the disease-free equilibrium E 0 from different values It determines whether the disease will die out or persist in the population as time increases. We have shown that if R 0 ≤ 1 , the disease-free equilibrium E 0 is globally asymptotically, which means that the disease cannot persist in the population, and the situation is under control. Moreover, we have proved that if R 0 > 1 , the disease-free equilibrium becomes unstable and the model has an endemic equilibrium E * which is globally asymptotically stable provided that S * ≥ I * , so that the disease, if initially present, will be persistent at the unique endemic equilibrium level. An application and some numerical simulations are presented to verify the feasibility of the theoretical results. We would like to point out here that Theorem 4 leaves us an open problem: whether the endemic equilibrium E * is globally asymptotically stable when R 0 > 1 without the additional condition S * − I * ≥ 0 ?. We leave this for future work. Table 2 . This plot shows the stability of the endemic equilibrium E * from different values Table 3 . This plot shows the dynamics of COVID-19 in Morocco The species severe acute respiratory syndromerelated coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2 Early transmission dynamics in Wuhan, China, of novel coronavirusinfected pneumonia Global analysis of SEI epidemic model with the constant inflows of different compartments Qualitative analysis of a SEIS epidemic model with nonlinear incidence rate Global stability for delay SIR and SEIR epidemic models with saturated incidence rates. Electron Global dynamics of an SEIRI epidemiological model with time delay Stability of SEIR model of infectious diseases with human immunity Global asymptotic stability for the SEIRS models with varying total population size National health commission of the people's Republic of China A kind of epidemic model having infectious force in both latent period and infected period Global stability of an SEI epidemic model Global stability of an SEI epidemic model with general contact rate Analysis of a spatially extended nonlinear SEIS epidemic model with distinct incidence for exposed and infectives Systems of ordinary differential equations which generate an order preserving flow A geometric approach to global stability problems Infectious Diseases of Humans: Dynamics and Control The mathematics of infectious diseases On the global stability of SIS, SIR and SIRS epidemic models with standard incidence A generalization of the Kermack-Mckendrick deterministic epidemic model Global stability for delay SIR and SEIR epidemic models with nonlinear incidence rate Lyapunov functional techniques for the global stability analysis of a delayed SIRS epidemic model Global stability for an SEI model of infectious disease with immigration A class of Lyapunov functions and the global stability of some epidemic models with nonlinear incidence Threshold dynamics of an SIRS model with nonlinear incidence rate and transfer from infectious to susceptible An SEIR epidemic model with relapse and general nonlinear incidence rate with application to media impact Stability analysis of a delayed fractional order SIRS epidemic model with nonlinear incidence rate Complete global stability for an SIRS epidemic model with generalized non-linear incidence and vaccination Dynamics analysis of an epidemiological model with media impact and two delays Global dynamics for an age-structured epidemic model with media impact and incomplete vaccination Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission Mathematical Epidemiology of Infectious Diseases: Model Building. Analysis and Interpretation The Stability of Dynamical Systems The authors declare that they have no competing interests.