key: cord-0045063-mnjnra38 authors: Baniya, Vinod; Keval, Ram title: The influence of vaccination on the control of JE with a standard incidence rate of mosquitoes, pigs and humans date: 2020-06-09 journal: J Appl Math Comput DOI: 10.1007/s12190-020-01367-y sha: a9fa3139ff6ddee1ab6cba40a7ceeb55ef9606c9 doc_id: 45063 cord_uid: mnjnra38 In this article, a nonlinear mathematical model used for the impact of vaccination on the control of infectious disease, Japanese encephalitis with a standard incidence rate of mosquitoes, pigs and humans has been planned and analyzed. During the modeling process, it is expected that the disease spreads only due to get in touch with the susceptible and infected class only. It is also assumed that due to the effect of vaccination, the total human population forms a separate class and avoids contact with the infection. The dynamical behaviors of the system have been explored by using the stability theory of differential equations and numerical simulations. The local and global stability of the system for both equilibrium states under certain conditions has been studied. We have set up a threshold condition in the language of the vaccine-induced reproduction number [Formula: see text] , which is fewer than unity, the disease dies in the absence of the infected population, otherwise, the infection remains in the population. Furthermore, it is found that vaccine coverage has a substantial effect on the basic reproduction number. Also, by continuous efforts and effectiveness of vaccine coverage, the disease can be eradicated. It is also found a more sensitive parameter for the transmission of Japanese encephalitis virus by using sensitivity analysis. In addition, numerical results are used to investigate the effect of some parameters happening the control of JE infection, for justification of analytical results. Japanese encephalitis (JE) is an illness of the brain caused through the Japanese encephalitis virus (JEV) and is transmitted from the pig population to the human population through bloodsucking insects Culex mosquitoes, Centers for Disease Control, (CDC). The JEV is maintained in a natural enzootic cycle between reservoir population (especially pigs), mosquitoes and humans. The transmission of JEV, for the most part, happens in the country in horticultural zones, where floodwater system and rice development are related [1] [2] [3] [4] . In the areas where JEV is endemic the annual incident range from 1 to 10 cases per 10,000 person per year affected children are under 16 years of age [5] . The time duration of infection to illness in infected humans is about 5-15 days and initial symptoms are fever, headache and vomiting [3, [6] [7] [8] . Most humans who are infected with only a few symptoms or have no symptoms and JE can be more deadly in advanced stages, National Vector Borne Disease Control Program, (NVBDCP). The first human case of JEV was registered in Japan in 1935 and in India in 1955 [3, 9] . Unfortunately, there is no specific treatment available that can cure Japanese encephalitis, including the management of some basic supportive medical care and complications in medicine. Those who experience a more severe form of the disease will have to undergo treatment at the hospital. The finding of JE disease ought to be performed by utilizing a JEV-explicit immunoglobulin M (IgM) counteracting agent catch ELISA for serum and cerebrospinal liquid (CSF). JEV-explicit IgM can be estimated in the CSF of most patients by four days after the commencement of side effects and in serum by seven days after inception. The prevention of JE is mainly based on the control of mosquitoes and the vaccination of the human population. To control the mosquito population, some methods like using insecticides, reducing the number of breeding sites on your property, by killing the mosquito larvae where the insect's breed and bed nets are used in endemic regions. Mathematical modeling of infectious disease has a long history for understanding the effects of vaccination. Through the speedy rate of vaccine development, medical science has saved millions of lives. Smallpox, influenza, polio and lots of supplementary diseases are conquering examples in this regard [10] . In the present life, it is bizarre and uncommon for a youngster not to get an antibody. An individual can get accessible immunizations for a malady common in his general vicinity. Some sickness immunizations are as of now created and anybody can take the antibody if the specific malady is in danger for it. Currently, several groups of vaccines available in the world like inactivated Vero cell culture-derived JE vaccine, SA 14-14-2 live attenuated JE vaccine, Ixiaro JE vaccine, etc. But Ixiaro is the only licensed JE vaccine available in the United States [11] . It was avowed in march 2009 for the use of the human population aged more than 17 years. Later on, May 2013, avowed for the use of children whose aged lies between 2 months to 16 years [1, 6, 12] . The essential inoculation plans for Ixiaro is two portions intramuscularly administrated on 0-28 days. For youngsters from 2 months to 2 years, each dose is 0.25 ml, for grownups and kids for more than 3 years, each dose is 0.5 ml. Prior, In India, the mouse brain-derived JE antibody was utilized, yet the productivity was insufficient to meet the national prerequisite [13, 14] . The SA 14-14-2 live attenuated JE vaccine has been utilized in China for over 20 years with an effectiveness of 80% for single-dose and 97.5% for a two-doses routine with 1 year interim among youngsters under 15 years old [11, 14] . Thus, the effectiveness of the JE vaccine is more than 90%, but the time duration of protection with the JE vaccine is not clear, but its effectiveness appears to decrease over time [1, 6] . Japanese encephalitis spreading fast in all states of India except Meghalaya, Himachal Pradesh, Jammu and Kashmir, Dadra, Daman, Punjab, Diu, Nagar Haveli, Gujarat, Lakshadweep and Sikkim [1] . This disease is seasonal, most human cases reported in north India in May-October. The highest rate of this disease reported from the states of Uttar Pradesh, Bihar, Andhra Pradesh, Assam, Kerala, Goa, Haryana, Karnataka, Tamil Nadu and West Bengal [5, 9] . In modeling of such infectious disease, many authors have already worked on epidemic models to minimize the infection of JE and maximize the rate of recovered population [15] . In recent years, mathematical models have been used to understand the interplay between various traditional measures such as the mobility of epidemic outbreaks [16, 17] and vaccination [5, [18] [19] [20] quarantine, behavior change [21] . When the immunization is accessible, a general inquiry emerges; How to allocate and implement this immunization. Unquestionably, we can't inoculate all people to control the sickness from the human population on account of social and moral issues, significant expenses, and so forth. In this manner, the legitimate conveyance and usage systems are critical to controlling JEV using vaccines. The adequacy of such a vaccination assignment procedure can be resolved through the investigation of a mathematical model. Right now, we need to illuminate this significant issue and would like to give a valuable rule to the arrangement producer. To appropriately actualize vaccination battles in the human population, the recognized and insightful idea may be to inoculate individuals having a spot with explicit social affairs or territories vulnerable to contamination. On account of mosquito dependent infectious diseases, the vaccination can be an important tool to help control the spread of infection in such infectious disease and encephalitis is one of them. It is pointed out that in the above model, vaccination has been studied without considering the effective role of the mosquito population which depends on the human populationrelated factor. In this work, we construct and analyze a modified mathematical model of De et al. [9] , by incorporating the effect of vaccination on the spread of JE, which is the most effective control strategy against JE. In addition, we use a realistic standard incidence rate for direct contact between susceptible and infectious. We might want to express that we don't consider any structure in the model analyzed down right now, these references can be utilized to improve the present work. Kar et al. [19] have been studied the SIR epidemic mathematical model to minimize the infection rate and to maximize the recovered rate of the human population by using vaccination program and optimal control strategy to SIR model. Singh et al. [22] have found that the development rate of the mosquito population relies upon the increasing density of ecological discharge by the human population. Also, the growth rate of the mosquito populations can be controlled by controlling the increasing density of ecological discharge. Iannelli et al. [23] build an SIS mathematical model to contemplate the effect of inoculation by expecting the immunization age, super-infection, immaculate immunization and furthermore talked about the asymptotic behavior of the equilibriums state. Raul Prealta et al. [24] , developed an SVIR mathematical model to study the effect of the vaccine by considering the immune loss rate, which depends on the vaccine-age under the condition that during vaccination individuals have the same immunity. In the present investigation, modifications to the SVIR model for JEV-affected humans and SI model for JEV-affected pigs [9] are considered by reducing the vaccinated compartment of a pig. The main objective is to analyze the effect of vaccination on the human population for the prevention of Japanese encephalitis with a standard incidence rate of mosquitoes and pigs. In the case of mosquitoes there is nothing to explain because it is obvious that vaccination, medicine is not possible for that compartment. Hence we consider susceptible, vaccinated, infected and recovered (SVIR) compartment for humans and susceptible, infected (SI) compartment of mosquitoes and pigs (cf. the Fig. 1 ). Disease transmission of the susceptible human population depends on the standard incidence rate of infected mosquitoes. Similarly, JEV spreads among the susceptible mosquito population when it comes in contact with infected human and pig populations. In this way, JE transmission for the susceptible mosquito population should be dependent on the rate of standard incidence on infected pig and human populations. The organization of the manuscript is as follows: In the next Sect. 2, notations and assumptions along with full explanation are provided. In Sect. 3, the mathematical formulation of the JE endemic model with a schematic diagram of the transmission of JEV puts in place and derive the model equations. We provide the region of biologically feasible solutions in Sect. 4, together with their positivity and boundedness. Dynamical behavior (existence of equilibrium states, basic reproduction number and stability) of the framework has been discussed in Sect. 5. In Sect. 6, vaccine-induced reproduction number and model parameters sensitivity investigation has been discussed to study the effects of controllable parameters of the JE model on the basic reproduction number and infected human population in Sect. 7 while the simulations and discussions are presented in Sect. 8 and conclusions in Sect. 9. In our proposed model, it is assumed that the JEV starts its life cycle from the pig population and then spreads into the human population through Culex mosquitoes. It is also assumed that all the newborn's pigs and mosquitoes are susceptible, but newborn humans are not susceptible. In addition, JEV transmitted only in susceptible individuals. Many variables and constants(parameters) symbols have utilized in our proposed mathematical epidemiological model, including four compartments for humans, two compartments for pigs and mosquitoes, which are given in Tables 1 and 2. In our present investigation, we consider susceptible, vaccinated, infected and recovered (SVIR) model for human population, susceptible, infected (SI) models for mosquito and pig population. The model sub-divides the total human population The total mosquito population N m is divided into subpopulation of susceptible mosquitoes (S m (t)) and infected mosquitoes (I m (t)), so that N m = S m (t) + I m (t). Also, the total pig population N p is sub-divided into susceptible pigs S p (t) and infected pigs (I p (t)) such that N p = S p (t) + I p (t). In SVIR model, humans are recruited at vaccinated birth rate α 1 λ h and unvaccinated birth rate (1 − α 1 )λ h . It is assumed that the JE spreads among the S h (t) due to contribution of the I m (t) [5] . Hence, the transmission of JE disease to susceptible humans should be a standard incidence rate of an infected mosquito population and it is in the form of Bβ mh Similarly, the JE transmitted among susceptible mosquitoes (S m (t)) due to contact with infected pigs and human populations. Hence, the transmission of JE disease to susceptible mosquitoes should be standardized incidence rate of infected pig and human population and it is in the form of Cβ pm The interaction between I m (t) and S h (t), I p (t) and S m (t) are assumed to be rated B, C respectively. Ultimately, the JE spreads among S p (t) due to contact with I m (t). Hence, the JE transmission to susceptible pigs should be a standard incidence rate of infected mosquitoes and it is in the form of Cβ mp I m N p S p . Also, we consider the proportion of the infected human population become recovered naturally and with medicine control at a rate of τ 1 and τ 2 respectively. After that, some of the recovered individuals go to susceptible class at rate f 1 . Humans in each compartment are assumed to die with a natural mortality μ h and die with the disease at a rate of λ. In the SI model, the recruitment rate of mosquitoes and pigs is λ m and λ p respectively. The natural mortality of the mosquito population in each compartment is μ m . Also, the pig populations in each compartment are assumed to die with natural mortality of μ p . Further, we consider that in the whole modeling process the disease transmission is possible only in susceptible individual class, for other classes JE spreads among them after they become susceptible class only. So, we do not consider direct JE transmission from vaccinated Table 2 Model parameters, its description and their values used for numerical calculation The rate at which vaccinated human goes to susceptible class of human 0.001 [9] B Average biting rate on susceptible human by infected mosquitoes per day 0.7 [ 9] C Average biting rate on susceptible pig by infected mosquitoes per day 0.9 [ 9] λ Disease-related death rate of infected human population 0.0222 [9] τ 1 , τ 2 The rate at which infected human becomes susceptible without and with medicine respectively 0.5 [ 9] 1/η h Mean viremic period of JEV (in days) δ Effectiveness of medicine on the human population 0.10 [9] β mh JEV transmission potential rate from infected mosquitoes to susceptible humans 0.2 [ 9] β pm JEV transmission potential rate from an infected pigs to susceptible mosquitoes 0.3 [ 9] β mp JEV transmission potential rate from an infected mosquitoes to the susceptible pigs 0.3 [ 9] α 2 The rate at which infected pigs go to the susceptible class of pigs by medicine 0.5 The mathematical model is described by following system of ODEs: For our simplicity we can re-write Eqs. (1)- (8) in notations form as- where, Now, we present the positivity and boundedness of solutions of the system of Eqs. (1)- (8) . For the system of Eqs. (1)-(8) to be biologically significant and well-proposed, essentially we have to show that all the solutions with non-negative initial conditions will remain positive for all time t > 0. In the theoretical study of disease transmission, the boundedness of the system measures that the system is biologically legitimate and well-behaved [16] . Biologically, validation of the model will be appeared by the accompanying theorem. The proof is excluded. The feasible region for the system of Eqs. (1)-(8) is given by Clearly, the system (9)-(16) has following two equilibrium states, namely The existence of equilibrium E 0 is obvious. It demonstrates that if the model has an infected population, then the equilibrium magnitude of susceptible human, mosquito and pig population will reach the valueS h = a 11 a 22 +a 21 β 1 respectively and vaccinated human population will remain at its equilibriumV h = a 21 a 22 . Let us define threshold parameter R 01 , as Biologically, R 01 is the average number of secondary Culex mosquitoes produced by single Culex mosquitoes [25] . To maintain the Culex mosquitoes in the environment, R 01 must be greater than one. To measure the disease transmission potential in such a population, the basic reproduction number R 0 can be established by using the concept of the next-generation matrix approach method. Since the system of Eqs. (9)-(16) are consistent, therefore the conditions for new infection and transition terms are given by the matrices where, For infection-free equilibrium state, the next generation matrix (F V −1 ) is calculated as: The potent characteristic root of matrix F V −1 is the basic reproduction number R 0 [26] . It is calculated as Biologically, R 0 is the average number of secondary cases generated by one primary JEV infected person, which has been introduced in the susceptible population, in which some individuals have been vaccinated. From Eq. (17), it is found that as the increase with the parameters β 1 , μ h , μ m , μ p and α 2 , the value of R 0 decreases. Also the parameters B, C, β mh and β pm affects the value of R 0 significantly. From the expression of R 0 , it is also found that vaccination rate α 1 affects the value of R 0 . It is imperative to take note that, the above parameters assume a significant job in the investigation of global threshold dynamics with regards to R 0 . Pronouncedly, if R 0 < 1, the disease will stop after some time and R 0 > 1 disease will stay endemic and irresistible population will arrive at positive continued levels [2, 4, 27, 28] . Clearly, from the perspective of controlling the JE infection, one would normally need to decrease the R 0 . In this way, it is advantageous to research how we can successfully diminish R 0 by sensible usage of vaccines and different parameters. The variation of R 0 w.r.t. some μ m , B, μ h , μ p are shown in the Figs. 2, 3, 4 and 5 (Table 3) . The local stability of infection-free equilibrium state (E 0 ) is determined after the characteristic roots of the Jacobian matrix J (E 0 ) of system (9)-(16) by following theorem: See "Appendix A" for proof of this theorem. and . The infection-free equilibrium is locally asymptotically stable if R 0 < 1 and R 01 < 1. Biologically, R 01 < 1 implies the absence of Culex mosquito population. In this manner, the infection-free equilibrium must be steady without the Culex mosquito population, in any case, disease despite everything perseveres, regardless of whether the susceptible have no immediate contact with the infected population. In this section, we have studied the global behavior of infection-free equilibrium point for the system of Eqs. (9)-(16). If R 0 ≤ 1, the the infection-free equilibrium state E 0 is globally asymptotically stable (GAS) on D. See "Appendix B" for proof of this theorem. (2) Endemic or positive equilibrium Using the value of Eqs. (19) and (21) in the expression for I * m in the Eq. (20) . This shows that the endemic equilibrium of this model satisfy the Eq. where, The root I * m = 0 of Eq. (22) correspond to the infection-free equilibrium (E 0 ) and the positive roots of the equation correspond to endemic equilibria point. From Eq. (23) we have If R 0 > R 01 > 1, K 3 < 0 and hence (K 2 2 − 4K 1 K 3 ) > 0. Thus, the positive equilibrium point is as But if R 0 < R 01 < 1, K 3 > 0 and hence (K 2 2 − 4K 1 K 3 ) < 0, this gives imaginary value of I * m . To verify above numerically, we use the set of parametric values from See "Appendix C" for proof of this theorem. To verify Theorem 5.3 numerically, we use the set of parametric values from Table 2 . Thus, the conditions are as Theorem 5. 4 The endemic steady state (E 1 ) is GAS in the region D, if the following inequalities are satisfied: See "Appendix D" for proof of this theorem. To verify Theorem 5.4 numerically, we use the set of parametric values from Table 2 . The simulations are done in the following seven cases for different initial starts and graphically displayed in Figs. 6 and 7. In this figure, the variation of Table 4 The components of endemic equilibrium points for non-linear global asymptotically stable 275 207 20 30 35 40 17 20 250 207 20 30 25 40 17 20 301 207 20 30 30 40 17 20 301 195 23 30 35 40 17 20 301 200 15 30 35 40 17 20 301 203 17 30 35 40 17 (Table 4) . Also for following set of parameter values, the conditions for the global asymptotical stability of the endemic equilibrium E 1 are satisfied. From Eq. (17), we define the vaccine-induced reproduction number (R(α 1 )) as Now, From above it is clear that R(α 1 ) is a decreasing function of α 1 if α 1 ≥ 0. This signalize that the effect of vaccination in reducing vaccine-induced reproduction number. Also, when vaccination is zero i.e. (α 1 = 0) then Eq. (25) can be written as From Eqs. (25) and (26), we can write the vaccination as If R * c < 1 then R(α 1 ) < 1, therefore E 0 is locally asymptotically stable. Also if R * c = 1, then R(α 1 ) ≤ 1 this is also the case in which E 0 locally asymptotically stable. From Eqs. (25) and (26) we can write, . From this we find R(α 1 ) > 1, therefore no amount of vaccination can bring R(α 1 ) less then one. Hence, σ c defines the stationary values for the vaccine related reduction rate of infection. Therefore, if the vaccination level α 1 > σ c , then possibility of disease eradication from population can be maximized (Fig. 8) . Solving Eq. (17) for α 1 , we get a critical vaccination level α 1c 1 , is given by Here, α 1c 1 is always positive, since R * 2 c > R 2 0 . If α 1 > α 1c 1 , R(α 1 ) < 1 as R(α 1 ) is a decreasing function for α 1 ≥ 0. In this manner, in the event that the vaccination coverage α 1 surpasses the condition α 1c 2 , at that point disease can be eradicated. From Eq. (19) , we also find that This critical vaccination level α 1c 2 guarantees that the disease annihilation when inoculated people stay at its harmony. Along these lines, if the vaccination rate α 1 is with the end goal that α 1 > max{α 1c 1 , α 1c 2 }, at that point disease annihilation is conceivable in the human population. In this section, we want to examine the relationship between input parameters and the results of the proposed model using the reproduction number R 0 . To assist us with discovering the parameters that highly affect R 0 , we dissected the R 0 to decide if immunization, treatment of infectives and mortality of mosquitoes can prompt the successful disposal or control of the JE in the population. The parameters chosen to perform the sensitivity analysis were B, μ m , δ and α 1 . We evaluate the sensitivity indices for B = 0.7, μ m = 0.06, δ = 0.1, α 1 = 0.5 and all other parameters are same as in Table 2 . For sensitivity analysis of the model, we define called normalized forward sensitivity index for a variable x and it depends on a parameter s [29] . The sensitivity indices of R 0 w. r. t. the parameters B, μ m , δ and α 1 are calculated as The fact is that Z B R 0 = 0.01743 signify that 1% increase in B, will produce 0.01743% increase in R 0 , keeping all other parameters fixed. When the parameters μ m , δ and α 1 increase by 1%, while keeping all others fixed then the values of R 0 decrease by 0.01587%, 0.001987% and 0.00495% respectively. Overall a lower value of R 0 is better because it increases the chance of disease eradication in the areas. Therefore, all the above prevention practices should be considered to reduce the parameter B, while increases the parameters μ m , u 1 and α 1 , instead it should be favored. It is seen that the impact of vaccination on people assumes an overwhelming job in the control of the disease. On the off chance that the contact rate (B) between susceptible people and infected mosquitoes can be managed then Japanese encephalitis can be killed totally from the population. It is additionally seen that the infected human population consistently expanding as for B and it is progressively delicate for the spread of JE disease, (cf. Figs. 9, 10). That's why we have chosen the parameter B only for the model sensitivity investigation. Figure 9 shows the sensitivity plot of the basic reproduction number R 0 for α 1 and β 1 . From this, we found that a higher vaccination rate α 1 and a low contact rate β 1 will keep R 0 < 1. Also, in Fig. 10 , the sensitivity plot of the R 0 regarding B and α 1 has appeared. This implies a higher estimation of α 1 and a low contact pace of B will keep R 0 < 1. In this paragraph, we provide some simulations on the dynamics of the system (1)- (8) . Moreover, to check how the human, mosquito and pig populations are affected by various parameters, we perform simulations by using the parameter values as shown Fig. 9 The surface plot shows the variation of R 0 with respect to β 1 and α 1 Fig. 10 The surface plot shows the variation of R 0 with respect to B and α 1 in Table 2 , and numerical results are manifested graphically in Figs. 11, 12, 13, 14 , 15, 16, 17 and 18 by using MATLAB(R2014a) software. In Fig. 11 , the variation of the human population with vaccination is shown and in Fig. 12 , the variation of the human population without vaccination is shown. From these two figures, we can study the effect of vaccination on the human population. In the case of vaccination, the susceptible population is below 50 and the vaccinated population is above 50, whereas in case of without vaccination, the vaccinated level goes to zero and susceptible level increases. We also see that vaccination does not affect the infected and recovered human population. This signifies that vaccination only affects susceptible individuals 13, 14) show the stability of the endemic equilibrium when R 0 is larger than one and the threshold condition T 0 of the JE virus larger than one. This reflects that the virus can be dominant when the T 0 and R 0 are larger than one. We observed in Fig. 15 , that due to the control strategy, the number of susceptible human decreases in society. This signifies that the spread of JE infection in society can also be controlled through the proper implementation of vaccines, medicine and some effective prevention. This strategy also shows that the effect of vaccine and medicine control on susceptible humans is almost the same. This reflects that some effective medical treatment and medicine can also control the spread of JE infection. Whereas in Fig. 16 , the number of vaccinated human increases in the community. This shows that the community can make infection-free through the effective treatment and vaccination control of a human. From Fig. 17 , we observed that due to the control strategy, initially decreases the infected human but after some time it stabilized. This shows that the vaccination strategy has no significant impact on the infected human. Finally, in Fig. 18 we observed that due to controlled strategy, the number of recovered human increases in the community but after some time it stabilized. This shows that infected humans can be immune (but not completely) through some basic supportive medical care and treatment. In this article, perspective on human, mosquito and pig population, a nonlinear scientific model for JEV has been proposed and examined to contemplate the impact of vaccination on control of infectious disease, Japanese encephalitis. The model displays two equilibria, to be a specific disease-free and endemic equilibrium. The first equilibrium exists without any conditions but the second equilibrium exists if the threshold parameter R 01 > 1. We have discovered another threshold parameter in terms of vaccine, vaccine-induced reproduction number R(α 1 ) which is, if not as much as unity, the disease ceases to exist, without infected population otherwise the infection is kept up in the population and disease turns out to be progressively endemic. The system of the differential equations has a unique infection-free steady state which is locally asymptotically stable if R 01 < 1 and unstable if R 01 > 1. Furthermore if R 01 approaches to 1, E 0 becomes unstable but E 1 exist. If R 01 > 1 then E 1 exists and it is locally asymptotically stable under certain condition. Also for both the equilibrium point global stability have to show by using Liapunov's function in the region of attraction under certain conditions. Model analysis shows, if the vaccination program has adequate capacity to mobilize a high degree of population, there is an abundance of mosquitoes in the environment, then the disease can be controlled. Moreover, it is found that vaccine coverage has a substantial effect on the basic reproduction number, R 0 . Also, by continuous efforts and effectiveness of vaccine coverage, the disease can be eradicated. Also, found by sensitivity investigation, if 1% increase in the parameters μ m , δ and α 1 keeping all others parameters fixed, then the value of R 0 decreases by 0.01587%, 0.001987% and 0.00495% respectively. Also, if we increase 1% in B then R 0 increases by 0.01743%. This signifies that the parameter B is more sensitive to the transmission of JEV. Then again, in the event that the vaccination rate is over a critical level α 1c 1 , at that point the JE doesn't take a type of pestilence gave the density of mosquito population stays at its equilibrium level. During the entire investigation of the JEV model, it is discovered that vaccination is the main long haul control measure against JE in the population. The mosquitoes reduction strategies may incorporate the end of mosquito breeding sites, keeping environmental factors perfect and clean. On the other hand, vaccinating the susceptible individuals in the human population with the end goal that vaccination level is kept up over the critical vaccination level α 1c 2 and the vaccine efficacy is sufficiently high so that vaccinated individuals are not infected by infected mosquitoes. Along these lines, our present paper has been investigated various sorts of control procedures and their adequacy in controlling the JE disease by utilizing the steadiness hypothesis of differential conditions and the basic reproduction number. From the clinical perspective, we accept that our model examination is gainful for persons and control of the JE disease. Proof The Jacobian matrix (J (E 0 )) of the system (9)-(16) corresponding to equilibrium state E 0 is given by The characteristic roots of Jacobian matrix (J (E 0 )) are In this article, all parameters are considered positive, therefore all φ i (i = 1, 2, 3, 4, 5, 6, 7) < 0. Then, Routh-Hurwitz criterion all the characteristic roots of Jacobian matrix (J (E 0 )) have negative real part if φ 8 < 0 [27, 30] . Thus, the infectionfree steady state (E 0 ) is LAS if R 0 < R 01 < 1, and unstable if R 0 > R 01 > 1. ThusV (t) ≤ 0 whenever a 13Sh +a 42Sp μ m ≤ 1 i.e. R 0 ≤ 1. Therefore, LaSalle's invariance Principle, [18, 31] implies that E 0 is GAS in D. The proof is complete. Proof The Jacobian matrix for the model system (9)-(16) corresponding to equilibrium E 1 is given by Using the region of attraction D (as explain in Sect. 4), we conclude thatV will be negative definite if it satisfy the following inequalities: Stability and bifurcation analysis of Japanese encephalitis model with/without effect of some control parameters Modelling and analysis of the spread of Japanese encephalitis with environmental effects The impact of media on the spreading and control of Japanese encephalitis Analysis of a vaccination model for carrier dependent infectious disease with environmental effects World Health Organization, Japanese encephalitis fact sheets An SIRS epidemic model of Japanese encephalitis Japanese encephalitis: a review of the Indian perspective Application of various control strategies to Japanese encephalitis: a mathematical study with human, pig and mosquito Modeling the impact of vaccination on infectious disease dynamical Safety and immunogenicity of an inactivated Vero cell derived Japanese encephalitis vaccine in a pediatric population in JE non-endemic countries: An uncontrolled, open-label phase 3 study New initiatives for the control of Japanese encephalitis by vaccination: minutes of a WHO/CVI meeting Effectiveness of Japanese encephalitis SA 14-14-2 live attenuated vaccine among Indian children: retrospective 1:4 matched case-control study Acquired immunity dependent upon exposure in a SIRS epidemic model Mathematical Modeling in Population Biology and Epidemiology Modelling the spread of carrier dependent infectious diseases with environmental effect Modeling monocyte-derived dendritic cells as a therapeutic vaccine against HIV Stability analysis and optical control of an SIR epidemic model with vaccination Global stability of an SIR epidemic model with information dependent vaccination Effect of awareness program in controlling the prevalience of an epidemic with time delay Modelling and analysis of the spread of malaria: environmental and ecological effects Strain replacement in an epidemic model with super-infection and perfect vaccination Global stability results in a SVIR epidemic model with immunity loss rate depending on the vaccine-age A network model for control of dengue epidemic using sterile insect technique Dynamical behaviors of an influenza epidemic model with virus mutation Global dynamics and bifurcation in delayed SIR epidemic model Comparison of deterministic and stochastic SIS and SIR models discretetime Effect of active case finding on dengue control: implications from a mathematical model Effect of constant immigration in plant-pathogen-herbivore interactions The stability of dynamical system Modelling the dynamics of Hepatitis C virus with combined antiviral drug therapy: interferon and ribavirin Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Proof To establish the global stability of the infection-free steady state (E 0 ), consider the following Lyapunov function, Using the condition (18) and feasible region D (as explain in Sect. 4), Eq. (27) can be rewritten as Proof To establish the global stability of the endemic equilibrium state (E 1 ), consider the following positive definite function [29] .Diff. V w.r.t. time t along the solution of system of Eqs. (9)-(16), we geṫOn rearranging the terms,V reduced to, Hence, we shows that V is a Liapunov's function for the system of equations (9)-(16), provided the conditions (39)-(43) hold.