key: cord-0686981-wjwjzphv authors: Cui, Ting; Liu, Peijiang; Din, Anwarud title: Fractal–fractional and stochastic analysis of norovirus transmission epidemic model with vaccination effects date: 2021-12-21 journal: Sci Rep DOI: 10.1038/s41598-021-03732-8 sha: 12f9611cd3ef9035feb844e114781b26efe30084 doc_id: 686981 cord_uid: wjwjzphv In this paper, we investigate an norovirus (NoV) epidemic model with stochastic perturbation and the new definition of a nonlocal fractal–fractional derivative in the Atangana–Baleanu–Caputo (ABC) sense. First we present some basic properties including equilibria and the basic reproduction number of the model. Further, we analyze that the proposed stochastic system has a unique global positive solution. Next, the sufficient conditions of the extinction and the existence of a stationary probability measure for the disease are established. Furthermore, the fractal–fractional dynamics of the proposed model under Atangana–Baleanu–Caputo (ABC) derivative of fractional order “[Formula: see text] ” and fractal dimension “[Formula: see text] ” have also been addressed. Besides, coupling the non-linear functional analysis with fixed point theory, the qualitative analysis of the proposed model has been performed. The numerical simulations are carried out to demonstrate the analytical results. It is believed that this study will comprehensively strengthen the theoretical basis for comprehending the dynamics of the worldwide contagious diseases. The NoV is disseminated by numerous factors that can promptly surge the transmission and henceforth the subsequent disease. The virus is known to be relatively sporadic, with more than half of the infections occurring during the cold season 1 . Seasonal fluctuation is attributed to ecological factors as well as demographic behaviour. Norovirus, for instance, is more easily feasting in cooler temperatures and might be aided by greater rainfall 2, 3 . Other factors in the population may influence the intensity of norovirus epidemics. The virus infects people of all ages, however it is most common in kids below the age of five 4 . According to epidemiological studies, that the very earliest norovirus disease arises in childhood. Serious problems and fatality are often more common among the aging and handicapped, according to 5, 6 . As a result, while disease is self-limiting in healthy people, the consequences for certain high-risk groups can be devastating. Such populations have also been observed to have much longer shedding intervals, which could lengthen an epidemic [7] [8] [9] [10] . Scholars from different fields have been working to prevent or minimize the rate of infection in certain groups. Mathematicians too have worked on the evaluation of relevant nonlinear dynamics of problems connected to infection, such as epidemics (see, for example [11] [12] [13] ). Since 1990, mathematicians and biologists have been working hard to learn more about how epidemics travel and how to prevent them from growing in the community. In the area of preventing diseases, mathematicians are also playing an important contribution by employing mathematical modeling approaches as well as optimal solutions [14] [15] [16] [17] [18] [19] [20] . Infection modelling and study have become increasingly popular in current years as a means of better comprehending the mechanisms of developing epidemics. The approaches of mathematical modeling and optimality may then be employed to come up with a decent control strategy or technique for managing infectious diseases. Most academics have highlighted the stability and optimality for viral models with a non-linear incident operator using deterministic model techniques 21, 22 . To come up with an effective control strategy or treatment for infectious disease, mathematical modelling in various dynamical systems, such as fractional and stochastic modelling, might be applied. Here W 1 (t), W 2 (t), W 3 (t), W 4 (t), W 5 (t) are independently standard Brownian motions, and ζ 1 , ζ 2 , ζ 3 , ζ 4 , ζ 5 are the intensities of standard Gaussain white noises, correspondingly. The parameters description given in Table 2 . Here, vaccinated people also become infected via contact with symptomatic people. Note that, 0 < τ < 1 means τ = 1 perfect vaccine, while τ = 0 represents a vaccine that offers no protection at all. The rest of the article layout has been established; In "Deterministic state stability" section, we determined the deterministic stability of system (1) . We prove the global positive solution for system (1) in "The existence and uniqueness of positive solution" section. We employed a stochastic threshold strategy to decrease the epidemic in "Extinction of the disease" section. And in "Stationary distribution and ergodicity" section we evaluate a stationary distribution of the stochastic model. In "A fractal-fractional NoV model with Mittag-Leffler kernel" section, we analyze the proposed model through Fractal-fractional Atangana-Baleanu operator. In "Parameter estimation" section, the application of the parameter estimation is presented for the proposed model. In "Numerical simulations" section, we give simulation results to support our theoretical results. In the "Conclusion" section, we are give some conclusions. The deterministic form of system (1) is governed by the following set of equations with initial conditions (1) (2) . www.nature.com/scientificreports/ The model equilibria provide useful information regarding the model trajectory over time. The model (2) posses two type of equilibrium points, the disease-free (DFE) and endemic equilibrium (EE) points. Disease-free equilibrium. The disease-free equilibrium (G 0 ) can be obtained by equating the right side of equations in system (2) to zero as follows: The DFE G 0 is locally and globally asymptotically stable if R D 0 < 1 , and unstable if R D 0 > 1. The Proof of the Lemma is simple so we omit it here. Using the next-generation matrix technique, we compute the basic reproductive ratio R D 0 using only the two equations corresponding to compartments U and A classes from system (2). The basic reproductive ratio, defined as the spectral radius of the matrix FV −1 , is obtained as Endemic equilibrium. The system (2) is also have EE equilibrium point which denoted by where Lemma 2 The EE G 1 is locally and globally asymptotically stable if R D 0 < 1 and unstable if R D 0 > 1. The Proof of the Lemma is simple so we omit it here. In order to addresses the existence and uniqueness of the stochastic model (1), we present the following theorem. The solution of the constructed stochastic epidemiological system (1) (H(t), V(t), U(t), A(t), C(t)) is unique for t ≥ 0 with initial condition (H(0), V (0), U(0), A(0), C(0)) ∈ C 5 + . Additionally, the solution will almost certainly remain in C 5 + with the unit probability, that is, (H(0), V (0), U(0), A(0), C(0)) ∈ C 5 + ∀ t ≥ 0 almost surely (a.s). , 0, 0, 0 . . , Proof The coefficients used in equations for the initial value of the state variables (H(t), V (t), U(t), A(t), C(t)) ∈ C 5 + are continuous and locally lipschitz. As a result, there should be a local unique solution of the system (H(t), V(t), U(t), A(t), C(t)) throughout t ∈ [0, τ e ) . The citations 22, 26 provide a detailed examination of the explosion duration τ e . To demonstrate the solution's global character, we must establish that τ e = ∞ a.s. Suppose we have a comparatively big nonnegative number k 0 such that each of the state's starting conditions are contained within the range [ 1 k 0 , k 0 ] . Let the final time be specified as k ≥ k 0 to every positive integer. We will use inf φ = ∞ throughout this article, whereas φ stands for the null set. The concept of τ k forces us to state that it rises as k reaches ∞ . Replacing τ ∞ = lim k→∞ for τ e ≥ τ ∞ a.s. After demonstrating that τ ∞ = ∞ a.s., we will assert that τ e = ∞ , and therefore (H(t), V(t), U(t), A(t), C(t)) will be found in C 5 + a.s. ∀t ≥ 0 . Hence, proving that τ e = ∞ a.s. suffices. Instead, two positive constants values ǫ from (0, 1) and T must exist, such that As a result, the integer k 1 ≥ k 0 exists in the following form After that, we'll look at how to interpret a C 2 -function H : C 5 + → C + in such a way that It is to be noted that the H is a nonnegative function, and it can be verified from the fact that 0 ≤ y − logy − 1, ∀ 0 < y . Assume that k 0 ≤ K and 0 < T are arbitrary. Upon applying Atô ′ s formula to Eq. (6) gives us In Eq. (7), LH : C 5 + → C + is defined by the following equation Setting � k = {τ k ≤ T} for k ≥ k 1 and by Eq. (5), P(� k ) ≥ ǫ . Note that for each ω from k there must exist one or more than one H(τ k , ω) , V (τ k , ω) , U(τ k , ω) , A(τ k , ω) and C(τ k , ω) which equals 1 k or k. As a result H(H(τ k ), V (τ k ), U(τ k ), A(τ k ), C(τ k )) is no less then 1 k − 1 + logk or k − 1 − logk . Therefore, By using Eqs. (5) and (9), we can write www.nature.com/scientificreports/ The indicator function of is represented as 1 �(ω) . As we get closer to ∞ , the contradiction ∞ > H H(0), V (0), U(0), A(0), C(0)) + MT = ∞ emerges, indicating that τ ∞ = ∞ a.s. This section focuses on the criteria for disease's extinction in system (1) . Prior to prove the major findings, let's look at an important lemmas. Let , then the root of system (1) fulfilling, so given as: It is simple to deduce below from the fourth equation of system (1). Furthermore, on both hand sides of the first equation in system (1), integrating from 0 to t and dividing with t obtains and considering (22) , and Lemma 3, it then follows that Similarly, we also can get The proof for Theorem 2 is finished. When it comes to stochastic systems, there are no endemic equilibria. As a result, the stability analysis could be utilised to investigate the disease's persistence. As a consequence, one should focus on the existence and uniqueness theory for the stationary distribution, which, in some ways, will help with disease persistence. We shall use Hasminskii's renowned finding 29 for this task. Let X(t) be a regular Markov process (time-homogeneous) in C n + for which the dynamics is as below: The diffusion matrix is of the form Describe a parameter is ergodic, and having a unique stationary distribution. Since R s 0 > 1 , π(.) is used. Proof To check the condition (2) of the Lemma 5, we should define a non-negative C 2 −function V : C 5 + → C + . For which we need to define here the positive constants c 1 , c 2 and c 3 must be calculated afterwards. We get the following results by utilising Itô's formula and the suggested system (1). Therefore, we have In addition, we obtain here c 4 > 0 is a constant that will be decided afterward. It's useful to illustrate that. It's not difficult to establish that V 2 has a unique stagnation point. 0)) contained in C 5 + . Following that, we'll define a non-negative C 2 −function V : C 5 + → C + as follows Considering Ito ′ s calculation and the proposed system, we reach at as a consequence of which the preceding assumption can be formed: where The subsequent step is to establish the set where ǫ i > 0 is a negligibly minor constant to be found afterwards for (i = 1, 2) . We'll split the entire C 5 + \D into the preceding regions for clarity's reason. www.nature.com/scientificreports/ For the smallest velue of ǫ 2 > 0 , so we can obtain LV < 0 for any arbitrary (H, V , U, A, C) ∈ D 1 0. As an outcome, we establish that a constant W > 0 is such that which assures Hence Assume that (H(0), V (0), U(0), A(0), C(0)) = (x 1 , x 2 , x 3 , x 4 , x 5 ) = x ∈ C 5 + \D , and τ x is that time for which a path start from x reach to the set D, One can get the following by taking integration of the both hand sides of the inequality (35) from zero to τ (n) (t) , considering expectation, and utilizing Dynkin's calculation. As V(x) is a non-negative, therefore We have P{τ e = ∞} = 1 as a result of the proof of Theorem 3. Conversely, the system (1) can be defined as regular. As a result, if we take t → ∞ and n → ∞ , we almost certainly obtain τ (n)(t) → τ x almost certainly. As a result, using Fatou's lemma, we arrive at Clearly, sup x∈K Uτ x < ∞, here K is a compact subset from C 5 + . It confirms Lemma 5 condition 2. Additionally, the diffusion matrix of the system (1) is where ρ = (ρ 1 , ρ 2 , ρ 3 , ρ 4 , ρ 5 ) ∈ C 5 + . This implies that Lemma 5 condition 1 is likewise true. Following to the preceding analysis, Lemma 5 indicates that the system (1) is ergodic and has just one stationary distribution. Fractional calculus (FC) has gained much interest from the researcher and scientists, because of its uses in different fields of real-world problems than that of integer order. FC has taken advantages and popularity of modelling with memory effects 30, 31 . Inspired from the work of Atangana 25, [31] [32] [33] , in the field of FC, we convert the proposed norovirus model to investigate its parameter with the available data in system (1) . The current section, deals with the approach of Atangana-Baleanu fractal-fractional (FF) derivative operator having fractional order p and fractal dimension q . We consider the presented model (2) to fractal-fractional order in sense of ABC operator. It is due to the fractional order derivatives have extra degree of freedom and some other characteristics of heredity, memory, and description of the past as well as present and future. This operator has non-singular kernel and is non-local operator. The NoV model can be shown through the following fractal-fractional differential system: Further, the system (36) can be write in the ABC fractal-fractional into the following way (37) www.nature.com/scientificreports/ The approximate function be K 1 on the interval [t ς , t ς +1 ] through the interpolation polynomial as follows which implies that Calculating I ς −1,p and I ς ,p we get and put t ς = ς�, we get and And similarly for the other classes V, U, A and C we may find the same scheme as In reference 34 , they study the data about a NoV infectious diarrhea incident reported in a middle school in a city. To guarantee the correctness and effectiveness of the methods, the data of a 2007 NoV outbreak in a middle school in one city is used as the real data to solve the inverse problem of the parameter estimation. For the daily reports and the data set in the Eq. (48), and the related relative error is used in the goodness of fit. (43) ��� . (44) ��� . (45) ��� . (46) ��� . (47) Table 2 (Set A), which ensure the conditions of Theorem 2, therefore, the infected population exponentially tends to zero with probability 1. The epidemic free equilibrium point shows global asymptotic stability in the related deterministic approach. As seen in Fig. 2 , the disease can die. We use parameter value from Table 2 (Set B) for the stochastic system (1) and compute R s 0 > 1 , which ensure the condition for the NoV persistent. As shown in Fig. 3 , the infection of system (1) will reside in the average, confirming the conclusions of Theorem 3. Observe 9000 attempts at t = 500 , then compute the average value. Theorem 3 suggests that the system (1) has an ergodic stationary distribution, as shown in Fig. 4 . Example 1 (Stochastic disease-free dynamical behavior) The parameter values are taken from Table 2 (Set A) . As a result, we obtain the reproduction number R s 0 < 1 , and the root of the model (1) may be satisfied by Theorem 2. and As a result, the pandemic will disappear from the population, as seen in Fig. 2 demonstrates that the numericalsimulation validates our strategy. Example 2 (Stochastic endemic dynamical behavior) We get the parameter values from Table 2 (Set B). We show that R s 0 > 1 , and that the illness will lie or stabilised according to Theorem 3, and we illustrate our findings in Fig. 3 . Theorem 3 states that the system (1) has just one stationary distribution, as shown in Fig. 3 . Numerical simulation for fractal-fractional system. In this section, we applied the novel numerical approach obtained above in (43)-(47) to simulate the proposed FF NoV model (36) having different arbitrary orders " p " and various dimensions " q " respectively given in Figs. 5 and 6, the initial and parameter value taken www.nature.com/scientificreports/ from Table 2 (Set C). We simulate the ABC fractal-fractional NoV model (36) when both the fractional order " p " and fractal order " q " are different. Figure 5a represents the susceptible class H(t) at different fractional orders and fractal dimensions. The class declines with the passage of time as the virus enters the society will transfer to the other classes of the system. This class converges quickly at low order and slowly at high orders. Figure 5b shows the vaccinated class V(t) growing at the beginning and then became stable at different fractal-fractional orders. As the disease control, the vaccination also controls and goes to their equilibrium point. Figure 5c is the representation of asymptotic or exposed class U(t) showing decrease in its behavior quickly with the passage of time like the behavior of susceptible class at different fractal-fractional orders. In Fig. 5 done can see the dynamical behavior of symptomatic or infectious class which also shows declines as controlled by vaccination and converging to their equilibrium point or became vanishes. Figure 5e shows the class of recovered population from the said epidemic at different fractional orders and fractal dimensions for the independent variable t. The cases of recovery increases at the beginning as more people have been vaccinated and the become stable. In Fig. 6 we simulate the ABC fractal-fractional NoV model (36) when both the fractional order " p " and fractal order " q " varies equally. From these figures, we note that changing both " p " and " q " at the same time an epidemic model provides interesting and biologically more feasible results because the infected population is vanishing more significantly as compared to the fractional systems. Thus, from these graphical results, we conclude that utilizing this new idea of FF operator one can observe more accurate results and provide deeper understanding not only for an epidemic model but also to the real world problem arising in science and engineering. In this paper, we propose an epidemic model for NoV transmission in considering environmental noise and fractal-fractional with vaccination effects. Based on the proposed model, there exist a unique time-global solution for any given positive initial value. The thresholds governing the extinction and propagation of the epidemic sickness are determined. The sufficient requirements for extinction of the disease and existence of ergodic stationary distribution of the stochastic system are then obtained from Theorems 2 and 3 utilizing Hasminskii theory and Lyapunov analysis methods. Finally, we applied the concept of FF calculus in the ABC sense to obtain the proposed model. In addition, numerical simulation are given to describe the solution behaviour of our theoretical result. We find that if the noise intensity is large, then the disease will go to extinction. The graphical simulations reveal that the fractal-fractional concept provides better and biologically more reliable results than the classical fractional and ordinary derivatives. In the future, the new approach of modeling known as fractal-fractional operator can be used confidently to study the dynamics of various infectious diseases including the novel COVID-19. Data sharing is not applicable to this article as no new data were created or analyzed in this study. A systematic review and meta-analysis of the global seasonality of norovirus The dynamics of norovirus outbreak epidemics: Recent insights Norovirus seasonality and the potential impact of climate change Norovirus infections and seroprevalence of genotype GII. 4-specific antibodiesin a Spanish population Stochastic optimal control for norovirus transmission dynamics by contaminated food and water Epidemiological study of prevalence of genogroup II human calicivirus (Mexico virus) infections in Japan and Southeast Asia as determined by enzyme-linked immunosorbent assays Duration of immunity to norovirus gastroenteritis Norovirus disease in the United States A norovirus outbreak in a nursing home: norovirus shedding time associated with age Prolonged norovirus shedding in infants 6 months of age with gastroenteritis A new mathematical model of COVID-19 using real data from Pakistan Mathematical epidemiologic and simulation modelling of first wave COVID-19 in Malaysia A novel method for a fractional derivative with non-local and non-singular kernel Prediction of treatment response for combined chemo-and radiation therapy for nonsmall cell lung cancer patients using a bio-mathematical model -functions as exact solutions for Caputo type mass spring damper system under Sumudu transform Analysis and numerical solution of the generalized proportional fractional Cauchy problem Mathematical modelling of the impact of treating latent tuberculosis infection in the elderly in a city with intermediate tuberculosis burden Analysis of series RL and RC circuits with time-invariant source using truncated M, Atangana beta and conformable derivatives Stationary distribution extinction and optimal control for the stochastic hepatitis B epidemic model with partial immunity The complex dynamics of Hepatitis B infected individuals with optimal control Hepatitis B optimal control model with vertical transmission Threshold behaviour of a stochastic SIR model Fractional numerical dynamics for the logistic population growth model under conformable Caputo: a case study with real observations Laplace transform method for economic models with constant proportional Caputo derivative Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system Delayed hepatitis B epidemic model with stochastic analysis Lévy noise impact on a stochastic hepatitis B epidemic model under real statistical data and its fractal-fractional Atangana-Baleanu order model Extinction and stationary distribution of a stochastic SIRS epidemic model with standard incidence rate and partial immunity Alphen aan den Rijn, Netherlands, Sijthoff and Noordhoff (Monographs and Textbooks on Mechanics of Solids and Fluids) New concept in calculus: piecewise differential and integral operators New fractional derivatives with non-local and non-singular kernel: theory and application to heat transfer model Modeling attractors of chaotic dynamical systems with fractal-fractional operators On some new properties of fractional derivatives with Mittag-Leffler kernel Is there a population awareness of microbes: the dynamics of Norovirus transmission In this section, we will illustrate our analytical results by some examples with the help of numerical simulations firstly. We use the first order stochastic iterative techniques of the fourth order Range kutta method to estimate the system solution (1). We adopt stochastic methodologies to simulate the double stochastic integrals because the system is driven by five independent noises dW i (t) for i = 1, 2, 3, 4, 5 . These integrals are approximated utilising the pseudo periodicity of the brownian bridge and their similarity to the deterministic Fourier series approach in our numerical method. For this, we use fourth-order Range kutta stochastic iterative techniques to accomplish the resulting discretization-transformation of the model (1),Here ξ k,i (k = 1, 2, 3, 4) are four free Gaussian general variables with N(0, 1) and �t > 0 time-increment.The values of the parameters listed in Table 2 .Numerical simulations for stochastic stability. Now, we will use numerical simulation to investigate the numerical approximation and biological feasibility of the system (1). As a result, we used the parameters and noise intensities value from Table 2 . The initial value of the Individuals susceptible H(t), vaccinated V(0), asymptomatic U(0), symptomatic A(t), and recovered R(t) are presented in Table 2 for t ∈ [0−150].(48) min T.C.: Conceptualization, data curation, validation, formal analysis, writing-original draft. P.L.: Supervision, Project administration, Funding acquisition, Visualization, review editing. A.D.: Visualization, Writing, Methodology, Software, review editing. All authors reviewed the manuscript. The authors declare no competing interests. Correspondence and requests for materials should be addressed to P.L. or A.D.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.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:// creat iveco mmons. org/ licen ses/ by/4. 0/.