key: cord-0927614-fmwgt0uz authors: Raw, S. N.; Mishra, P.; Sarangi, B. P.; Tiwari, B. title: Appearance of Temporal and Spatial Chaos in an Ecological System: A Mathematical Modeling Study date: 2021-05-07 journal: Iran J Sci Technol Trans A Sci DOI: 10.1007/s40995-021-01139-8 sha: 66aaa58ee14f4fa62d7ced1395756335b34de6c2 doc_id: 927614 cord_uid: fmwgt0uz The ecological theory of species interactions rests largely on the competition, interference, and predator–prey models. In this paper, we propose and investigate a three-species predator–prey model to inspect the mutual interference between predators. We analyze boundedness and Kolmogorov conditions for the non-spatial model. The dynamical behavior of the system is analyzed by stability and Hopf bifurcation analysis. The Turing instability criteria for the Spatio-temporal system is estimated. In the numerical simulation, phase portrait with time evolution diagrams shows periodic and chaotic oscillations. Bifurcation diagrams show the very rich and complex dynamical behavior of the non-spatial model. We calculate the Lyapunov exponent to justify the dynamics of the non-spatial model. A variety of patterns like interference, spot, and stripe are observed with special emphasis on Beddington–DeAngelis function response. These complex patterns explore the beauty of the spatio-temporal model and it can be easily related to real-world biological systems. Ecological systems can be delineated by the interactions of different biological populations within a vacillating natural environment. When a similar kind of species lives in a particular geographical area, then they are identified as populations. Population is related to each other by the food they eat and distributed in different trophic levels of the food chain. Predator-prey interaction is a crucial component of the food chain system. The study of mathematical models is very informative in understanding the dynamic behavior of predator-prey interaction. Predator-prey models have become a very extensive field of research in mathematical ecology following the pioneering work of Lotka (1924) and Volterra (1926) . Interaction between predator and its prey is very crucial factor in predator-prey models and it is well known as functional response in mathematical ecology. Holling (1965) introduced four different functional responses known as Holling I, II, III & IV type functional response. The addition of mutual interference between species devises the foremost distinction between the Beddington-DeAngelis (Beddington 1975; DeAngelis et al. 1975 ) and Holling type II functional response. Sometimes in search of food and habitat, species undergo self-inhibition along with corporation among themselves. Also, It is regarded that at high predator density, the feeding rate is relatively low due to mutual interference in predators and the prey dependent functional form ineffective to explain this. In that situation, the Beddington-DeAngelis function describes the real aspect of the ecological systems. Also, this functional response regards handling and interference as exclusive ventures by species. The only drawback of this functional response is that it explicates singular behaviour at the low population density. Here, both Holling type II and Beddington-DeAngelis functional responses are incorporated for modeling the interaction among populations in this paper. Naji and Balasim (2007) investigated the chaotic dynamics of a prey-predator model using the Beddington-DeAngelis type functional response. The influence of predator interference on the stability of a food chain model is studied by Upadhyay et al. (2013) . A lot of work has been done to observe the dynamics induced by interference between species in predator-prey systems. Dimitrov and Kojouharov (2005) analyzed a two species prey-predator model with Beddington-DeAngelis functional response and stated that mutual interference of predator has a stabilizing effect in presence of linear prey growth. Hwang (2003) proved the global stability of a Beddington-DeAngelis type preypredator model using the divergence criterion whereas Chen et al. (2008) used Fluctuation Lemma to prove the same. Haiyin and Yasuhiro (2011) inscribed a detail comparative study of the density-dependent predator-prey model with time delay and without delay in presence of Beddington-DeAngelis functional response. Cui and Takeuchi (2006) mainly focused on species extinction criteria and permanence of a prey-predator system. Zhang et al. (2005) studied the dynamical behavior of a Holling I prey-predator model with mutual interference in the case of pest-free periodic solutions. Parshad et al. (2016) stated that mutual interference between species can cause population explosion and mathematically termed this as finite time blow up. Mathematical modeling is adapted to demonstrate reallife predicaments and to exert prospective enhancements accordingly. Researchers practice numerous techniques to determine the numerical and analytical solution to these mathematical models. This physical interpretation of reallife problems contributes considerably to society. Epidemiological models have been used to describe HIV infection, cancer, and recently the spread of corona virus. Ali et al. (2020) studied HIV-1 infection of CD4 þ T-cells using fractional mathematical model. Sánchez et al. (2020) illustrated the mathematical pattern of the nervous stomach model and this model is designed using the differential transformation (DT) method. Gao et al. (2020a, b, c) studied the epidemic outbreaks of the noble corona virus. Yel et al. (2020) and Gao et al. (2020d) investigated the popular shallow water equations. Baskonus (2016 Baskonus ( , 2019 examined acoustic wave behaviors of Davey-Stewartson equation using the sine-Gordon expansion method, further elucidated the Boussinesq dynamical model applying the modified exponential function method. In the nineteenth century, most mathematical models in population ecology arrogated with the uniform population distribution of species in their habitat (Legendre and Fortin 1989) . The consequences of their interaction with inflammable dynamics are the subject analysis with reactiondiffusion models. Here the reaction is enthralled to predator-prey class, and diffusion leads to the redistribution of their populations. Hence species movement and interaction often give rise to the formation of a spatial pattern in their homogeneous environment. The concept of Turing instability was introduced by Turing (1852) to dispense the interaction of species in a homogeneous environment. Levin and Jackson (1972) first employed the concept of Turing pattern formation in the spatiotemporal system and later, extended by Levin and Segel (1976) . In the past quarter-century, many researchers (Alonso et al. 2002; Banerjee and Abbas 2015; Banerjee and Volpert 2017; Banerjee and Banerjee 2012; Sabrina and Araujo 2007 ) studied two species model systems for pattern formulation in one and two-dimensional space. Gauss (1935) studied stabilization scenario and the long-term existence of certain species in the spatial domain. He observed that diffusivity of the population has the power to stabilize as well as destabilize the coexistence states. Extension properties, persistence and effect of spatial distribution on the stability of population was studied by Luckinbill (1973) . Liu et al. (2008) studied a spatiotemporal phytoplankton-zooplankton model with additive noise and periodic forcing. Gazi (2012) studied the role of diffusive instability of a marine plankton system. Nowadays a wide variety of literature is available evidencing spatial patterns like the spiral, spot, target and spatiotemporal chaos (Cantrell and Cosner 2003; Malchow 1993; Pascual 1993) . Upadhyay et al. (2009) observed a wave of chaos (WOC) in the plankton fish model and confirm the generation of two-dimensional time-independent spatial patterns due to the wave of chaos. Recently, some researchers (Parshad et al. 2014; Wonlyul and Inkyung 2013; Zhifu 2012 ) studied Turing instability and the formation of patterns for two and more interacting species. Maionchi et al. (2006) studied the complex dynamics in the spatiotemporal three species food chain model due to the diffusion of species. They observed very complex patterns in the food chain model. Baghel and Dhar (2014) performed a high-order stability analysis in a three-species food chain model. Guin et al. (2017) derived essential condition for Hopf and Turing bifurcations in a predator-prey model. Also, the study complied with self and cross-diffusion patterns. Upadhyay et al. (2013) studied dynamics induced by the predator interference with generalist predator and Parsad et al. (2016) extended this model and pointed out some other properties of predator interference. They also showed the spatial chaos and Turing instability for the model. Enlightened by the above literature, we have proposed and analyzed a food chain model including a specialist top predator. We have attempted to study the effect of mutual interference between two predators over the prey population. A very limited number of articles have published focusing on pattern formulation in three species predator-prey system. Therefore, in this study, the spatial complexity of the model is also investigated with the help of pattern formulation. The paper is ascertained as follows: Sect. 2, carries the description of the proposed predator-prey model system. Conditions for boundedness, Kolmogorov and stability of model (3) is presented in Sect. 3. Section 4 deals with stability and Turing instability conditions of spatial model system (2). Numerical results are intricately exhibited in Sect. 5 with the help of MATLAB. Section 6 comprises the summed results of our model. Here we have considered a food chain model with prey and two predators. The prey is assumed as the primary consumer whereas predators are assumed to be secondary and tertiary consumers, respectively. The secondary consumer feed on prey U according to the Holling type II functional response and tertiary consumer W feeds on V according to the Beddington-DeAngelis type functional response. Prey (U) has a logistic growth delineated by the functional expression a 1 U 1 À U K , including intrinsic growth parameter a 1 and carrying capacity K. Secondary consumers (V) has a decline in population density with a death rate of a 2 and c expresses the death rate of tertiary consumer (W). Usually, it is seen that species diffuse or show movement in the exploration of food, mate, habitat, and sometimes in fear of predators. Thus, to explain the movement of species (U, V, W) in the x À y direction, we have induced diffusion to the model system. Now the assumed three species reaction-diffusion food chain model is represented as: Holling type II functional response. D and D 1 symbolizes the assistance of the environment given to prey species where D reveals the advantage of prey to evade predation and D 1 describes the loss in maximum gain of secondary consumer Y. The authenticity of this kind of situation in the environment is also given by Upadhyay et al. (2013) . Ecological importance and units of adapted positive parameter values are expressed in the Table 1 . We have assumed D u ; D v and D w as the diffusion coefficients as the associated species actively diffuses in the x À y direction. n shows the outward normal to oX. Here r 2 ¼ o 2 ox 2 þ o 2 oy 2 is the Laplacian operator in 2D Cartesian coordinate system (x, y). The non-dimensional model of system (1) is represented as: Corresponding to system (2), the initial populations and no-flux boundary conditions transform to uðx; y; 0Þ [ 0; vðx; y; 0Þ [ 0; wðx; y; 0Þ [ 0; for ðx; yÞ 2 X; In the next section, we will analytically examine the behavior of the non-spatial system (3). We have limited the non-dimensionalized system (3) Theorem 1 Assume that condition w 7 ! w 4 holds and let K be the set defined by All positive solution of temporal system (3) is uniformly bounded and as time varies the critical solution eventually converges to the region K. Theorem 1 confirms the boundedness and dissipativeness of the system (3). The proof of this theorem is given in ''Appendix 1''. Now, to justify the feasibility and existence of our proposed food chain model, system (3) need to satisfy the Kolmogorov sufficient and necessary conditions (May 1973) . Kolmogorov system ensures periodic oscillations or chaotic dynamics in 3D system if two subsystem of any 3D system is K-system. Let us consider the reduced subsystem of temporal model system (3), we have In the two species dynamical models, the Kolmogorov theorem assures the existence of a stable limit cycle or stable equilibria and ensures biologically relevant parametric values for the model. Therefore, for the biologically relevant predator-prey model, we assume that the subsystem (4) should qualify as K-system. The subsystem (4) contemplate a K-system with the conditions: Moreover, the Kolmogorov system (4) contains three positive equilibrium points. The sustainability of these equilibria are summed as: 1. The equilibrium E 00 ¼ ð0; 0Þ is exists with eigenvalues 1; Àw 5 . Thus, the nature of the eigenvalue is characterized as saddle point. 2. The equilibrium E 10 ¼ ð1; 0Þ always exist and is stable and eigenvalues are À1; Àw 5 3. The equilibrium E 20 ¼ ðũ;ṽÞ, where always exist and if 2 w 5 w 7 holds then E 20 is locally asymptotically stable. In this subsection, we obtain existence conditions for equilibrium points of the temporal model (3) and the stability for positive equilibria E 3 =ðu H ; v H ; w H Þ is investigated. The equilibrium points E 0 ð0; 0; 0Þ and E 1 ð1; 0; 0Þ always exist for model (3). The equilibrium point E 2 ðũ;ṽ; 0Þ exist in u À v plane, whereũ andṽ are given in Eq. (6). Further, due to deficiency of secondary consumer (v) there is no equilibria exists in the v À w plane. Now to study the existence of equilibrium point E 3 ¼ ðu H ; v H ; w H Þ, we consider following algebraic nonlinear equations: Thus, solving the system (8) w H ¼ ðw 11 À w 10 Þv H À w 9 w 10 w 9 w 10 : For v H [ 0 and w H [ 0 it follows that 0\u H \1; 0\ w 9 w 10 w 11 À w 10 H : Next, we explore the dynamical behavior of the model (3) via local stability around the feasible equilibrium points. The calculated variational matrix of the system at the point (u, v, w) is rendered as: The linear stability around the critical points E i ; i ¼ 0; 1; 2; 3 are distinguished as: (i) The equilibrium E 0 is saddle point and eigenvalues are given by 1; Àw 5 ; Àw 10 . (ii) The equilibrium point E 1 exist and eigenvalues are given by À1; The root of characteristic equation p 3 ðkÞ ¼ 0 of the above variation matrix (11) about E 2 ¼ ðũ;ṽ; 0Þ satisfy the following condition: The equilibrium point E 2 ¼ ðũ;ṽ; 0Þ is stable or unstable in the positive w-direction contingents on parameter w 10 and w 11 . If w 11ṽ ðṽ þ w 9 Þ \w 10 , then E 2 rests on being stable. where Then the characteristic equation of (12) is where A 1 ¼ Àða 11 þ a 22 þ a 33 Þ; A 2 ¼ ðða 11 þ a 22 Þa 33 Þ À ða 12 a 21 À a 11 a 22 Þ À a 23 a 32 ; A 3 ¼ ða 12 a 21 À a 11 a 22 Þa 33 þ a 11 a 23 a 32 ; A 1 A 2 À A 3 ¼ ða 2 11 À a 32 a 23 ÞðÀa 22 À a 33 Þ þ ða 2 33 À a 12 a 21 ÞðÀa 11 À a 22 Þ þ a 2 22 ðÀa 11 À a 33 Þ À 2a 11 a 22 a 33 : Therefore, from Routh-Hurwitz criterion if a 11 ; ða 11 þ a 22 Þ; ða 12 a 21 À a 11 a 22 Þ; then the coefficients of Eq. (13) is strictly positive with A 1 A 2 À A 3 [ 0. Hence, the coexisting equilibria E 3 ðu H ; v H ; w H Þ is locally asymptotically stable. Theorem 2 The condition for the happening of a Hopf bifurcation for system (3) at w 8 ¼ w 8 are following Proof From the condition A 1 A 2 À A 3 , we have A 1 A 2 À A 3 ¼ ða 2 11 À a 32 a 23 ÞðÀa 22 À a 33 Þ þ ða 2 33 À a 12 a 21 ÞðÀa 11 À a 22 Þ þ a 2 22 ðÀa 11 À a 33 Þ À 2a 11 a 22 a 33 ¼ 0: Since A 2 [ 0 at w 8 ¼ w 8 , there exist an interval containing w 8 2 ð w 8 À ; w 8 þ Þ, where [ 0. Therefore for the interval w 8 2 ð w 8 À ; w 8 þ Þ the characteristic equation (13) cannot have real roots contain negative real part. For w 8 ¼ w 8 , Eq. (13) can be written as which has exactly three roots Now for transversality condition, substituting above given value in Eq. (13), we get where where U 1 M 1 þ U 2 M 2 6 ¼ 0 and k 3 ¼ ÀA 1 ðw 8 Þ. Thus, the transversality condition Re h dk j da i w 8 ¼ w 8 6 ¼ 0 holds implying the rise of a Hopf bifurcation at w 8 ¼ w 8 . In the next theorem, we analyze the sustainability of equilibria E à 3 irrespective of initial conditions in the domain K. Theorem 3 Under the conditions: v H \w 4 ðu H þ w 4 Þ; ð18Þ w 6 w 11 v c ðv H þ w 8 Þ\w 9 w 10 ðv H þ w 8 w H þ w 9 Þ; w 2 8 w 2 c \d 2 w 9 ðv H þ w 9 Þ; ð21Þ equilibrium point E H 3 is globally asymptotically stable in the interior of the positive octant K and the set K is defined in Theorem 1. The proof of Theorem 3 is conceded in ''Appendix 2''. Here, we analyze the effect of diffusion for the spatially homogeneous equilibrium E 3 ðu H ; v H ; w H Þ for the model system (2). The interaction among the prey u(x, y, t), secondary consumer v(x, y, t) and tertiary consumer w(x, y, t) is described in system (2). To study the spatial dynamics of model system (2), we linearize the model about the equilibrium E 3 ðu H ; v H ; w H Þ as follows where we introduce small perturbations Here c i , i ¼ 1; 2; 3 are constants. k x and k y are the components of wave number along Cartesian plane (x, y) with wave length k k . The variational matrix of the linearized model system is given by Eq. (22) is defined as J ¼ where k is the wave number given by k 2 ¼ k 2 x þ k 2 y and characteristic equation of (24) is given by where q 3 ¼ d 1 d 2 d 3 k 6 À ½a 11 d 2 d 3 þ a 22 d 1 d 3 þ a 33 d 1 d 2 k 4 þ ½d 1 ða 22 a 33 À a 23 a 32 Þ þ d 2 a 11 a 33 þ d 3 ða 11 a 22 À a 12 a 21 Þk 2 þ A 3 ; with A 1 ; A 2 and A 3 are defined in Eq. (13). Theorem 4 Uniform steady state E à of system (2) is locally asymptotically stable if following conditions satisfy: The proof of this theorem is simple and determined by the Routh-Hurwitz criterion. Hence, the proof is omitted here. Theorem 5 The globally asymptotically stable steady state of the temporal model is also globally stable for the corresponding spatial model. (3) is globally asymptotically stable around the positive equilibrium E 3 ¼ ðu à ; v à ; w à Þ, then corresponding spatial system (2) is also globally asymptotically stable in presence of diffusion. Theorem 5 verification is remitted in ''Appendix 3''. This subsection deals with the derivation of conditions that ensure the occurrence of Turing instability of the model system (2). Turing instability occurs when the stable equilibrium point of the temporal system becomes unstable by adding diffusion. It may occur for a certain range of wave number with an authentic choice of parameter values. However, The Turing instability conditions need not be determined by the geometry of the model system, it determined by the reaction rates and diffusion coefficients (Baek et al. 2013 ). Turing instability may occur if at least one root of Eq. (25) has the positive real part for some k ! 0. Disregarding the sign of q 1 and q 2 in Eq. (25), the equation has a positive real root if q 3 ¼ qðk 2 Þ 0. Therefore, diffusion driven instability occurs when qðk 2 Þ 0. Hence condition for diffusion instability is given by where Iran J Sci Technol Trans Sci p 2 ¼ d 1 ða 22 a 33 À a 23 a 32 Þ þ d 2 a 11 a 33 þ d 3 ða 11 a 22 À a 12 a 21 Þ: The minimum value of qðk 2 Þ occur at k ¼ k cr given by Straight forward algebraic calculations show at least one root of Eq. (25) has positive real part, and thus diffusive instability occur if either p 2 or p 0 and p 2 1 [ 3p 0 p 2 . However, Turing instability does not not exist for model (2) if condition (14) hold. We have discussed the Turing instability conditions of the model system (2) for a particular set parameter values with diffusion coefficients. For example, the following parameter values w 4 ¼ 0:5; w 5 ¼ 1:15; w 6 ¼ 2:93; w 7 ¼ 0:58; In this subsection, we venture to support our analytical findings numerically via MATLAB programming. Also, we perceive the diverse and complex dynamics of the model system in this section. These observations are We testify that the model (3) manifests the global behavior for this parameter set. The model system (3) shows chaotic behaviour at w 9 ¼ 0:08, simulated in the Fig. 2a with time evolution of species populations in Fig. 2b . The Lyapunov exponent for this attractor is calculated as 0.065479, -0.057648, and -0.38053 accumulated in Fig. 2c . Here the appearance of both positive and negative Lyapunov exponent reinforces the existence of chaos for the model system (3). We have simulated the periodic behavior at the parameter value w 9 ¼ 0:12 for the model system (3) in Fig. 3 . We have ascertained the Lyapunov exponent for this attractor as 0.0096048, -0.14618, and -0.16712, and the convergence of the first Lyapunov exponent toward 0 validates the periodic behavior of the model system (3). As we increase the parameter value of w 9 , the system emergence towards stability, and finally at w 9 ¼ 0:16, we witness the stable focus dynamics for the temporal model as given in Fig. 4a , b. Negative Lyapunov exponents -0.022655, -0.031018, and -0.20335 as simulated in Fig. 4c authenticates this argument. The tertiary consumer free equilibrium point E 2 with condition (6) is calculated as u ¼ 0:375; v ¼ 0:547; w ¼ 0, as presented in Fig. 5 . In this section, bifurcation analysis is performed for the non-spatial model. All the bifurcation diagrams are plotted for the parameter values given in Eq. (26). The bifurcation diagrams show intricate and rich dynamics representing various sequences of period-doubling and period-halving bifurcation leading to limit cycles and chaos. Also, we can Fig. 7 , the successive maxima of u, v, w is given for the variation of parameter value w 10 in the range 0:2 w 10 0:4. In Fig. 8 , the successive variations of prey (u), consumer (v) and consumer (w) are taken concerning the parameter w 11 in the interval 0:5 w 11 0:75. In this section, we have found inclusive results for the spatial model (2). We have assumed the following set of parameters for the numerical solution of the spatial model system (2): w 4 ¼ 0:5; w 5 ¼ 1:3; w 6 ¼ 2:93; w 7 ¼ 0:58; w 8 ¼ 0:21; w 9 ¼ 0:16; w 10 ¼ 0:35; This subsection deals with the simulation of the system (2) with one dimensional Laplacian operator r 2 ¼ o 2 ox 2 . Here, the populations uðx; 0Þ [ 0, vðx; 0Þ [ 0 and wðx; 0Þ [ 0, x 2 ½0; R, and zero flux boundary conditions are given by We have plotted space versus population density to scrutinize the spatial dynamics. Here, we have considered the non-monotonic initial conditions as where ðu H ; v H ; w H Þ is steady state for the coexistence of primary consumer (prey), secondary consumer and tertiary consumer populations and ¼ 10 À8 , x 1 ¼ 1200, x 2 ¼ 2800 is the parameter affecting the model dynamics. We have studied the behavior of model system (2) in 1D space for the set of parameter values given in Eq. (27). Spatial dynamics of model (2) is given in the form of phase portrait diagram for different value of w 8 in Fig. 10 . We obtain the spatial dynamics of the model (27) varying time. System is stable at time t ¼ 10 but as we increase time, there are onset chaotic faces at t ¼ 500 and t ¼ 5000. In Fig. 9 , the population density of prey and predator shows some jagged patterns representing the chaotic behavior. This phenomenon is known as Wave of Chaos (see Petrovskii and Malchow 2014; Upadhyay et al. 2009 Upadhyay et al. , 2016 in spatio-temporal system. This section deals with the simulation of the spatial system in two-dimensional space. We study spatio-temporal complexity with the help of pattern formulation. Fig. 7 Bifurcation diagrams corresponding to parameter w 10 for the temporal system (3) We have adopted the forward difference scheme to solve the reaction part and the standard five points explicit finite difference scheme method for the two-dimensional diffusion term. Further, we have exerted nonzero initial and zero-flux boundary conditions with the system size of 50  50. Also, we can implement Euler methods to determine time evaluations. Here, the values of the time and space for simulation of the spatial model system is considered as Mh ¼ 0:25; Mt ¼ 0:001. The dynamics of the spatiotemporal system depends upon the choice of initial conditions (Medvinsky et al. 2002) . The initial distribution of the species for two-dimensional simulation is considered as uðx; y; 0Þ ¼ u H þ 0:1 cos 2 ð10xÞ cos 2 ð10yÞ; vðx; y; 0Þ ¼ v H þ 0:1 cos 2 ð10xÞ cos 2 ð10yÞ; wðx; y; 0Þ ¼ w H þ 0:1 cos 2 ð10xÞ cos 2 ð10yÞ; about steady state. Now, we present the evolutionary process of Turing pattern formulation and complex spatial patterns for the two different set of parameter values. Snapshots of the spatial distribution are shown for different time intervals t ¼ 500; 5000, 10,000 in Fig. 11 . We plot the contour spatial distributions for population u, v, w at different time intervals with 50  50 points for d 1 ¼ 1, d 2 ¼ 0:001, d 3 ¼ 0:01 and rest parameter values defined in Eq. (27). In Fig. 11a , we observe stripe patterns with emergence of interference patterns in the left lower boundary of domain. A close examination on these patterns suggest that as time increases, the spatial distribution changes from one pattern to another chaotically and we get very complex patterns. At t ¼ 5000, mixture of stripe and interference patterns are observed, although interference pattern dominate in this case (see Fig. 11b ). Finally, at time t ¼ 10,000, we see interference patterns dominating in the whole domain (see Fig. 11c ). From Figs. 12, 13, we observe that pattern of spatial distribution depends upon mutual interference parameter w 8 . Here the blue spot signifies the patches at low population density and the red spot exposes the patches at high population density. Also, these patterns are named as spot patterns in literature (see Malchow 1993; Medvinsky et al. 2002) . Increasing the value of w 8 in Hopf-Turing domain effect in spatial distribution and it can be seen that complex spot (or 'leopard') patterns changed into the mixture of spot and stripe pattern. If food is abundant the predator does not interfere with each other but a striking competition occurs between predators for food while the scarcity of prey species happens. The main purpose of this work is to detect the complex dynamics in the tri-tropic food chain model when mutual interference is adopted by the predators during the predation. We incorporate Bedington-DeAngelis functional as a predation function between the predators. The dynamics of the non-spatial model system have been studied with the help of phase portrait, time evolution, and bifurcation diagram. Several studies have been investigated respecting the predator-prey model with Beddington-DeAngelis functional response. The research exposes the diverse dynamics including species extinction, persistence, stable or unstable equilibria, limit cycles, and strange attractors (commonly known as chaos) (see Baghel and Dhar 2014; Dobromir and Dimitrov 2005; DeAngelis et al. 1975; Hwang 2003; Naji and Balasim 2007; Upadhyay et al. 2013) . Boundedness and Kolmogorov conditions are obtained for non-spatial model. Figure 4 shows the locally asymptotically positive equilibria E 3 for model (3). We have simulated the equilibria E 2 in Fig. 5 and the coordinate values where the u-nullclines and v-nullclines intersect is calculated as u ¼ 0:375, v ¼ 0:547 and w ¼ 0. Also, this confirms the stability of equilibrium point E 20 for the Kolmogorov system (4). We observe the non-spatial model goes for Hopf bifurcation at some critical value of the interference parameter. Bifurcation diagrams show the consequences of mutual interference and other important parameters of the predator-prey model. The study shows that mutual interference plays a central role in the dynamics of the predator-prey system. It is observed that a large value of interference parameter (w 8 ) ensures stability, but leads to a fall in the top predator population. In Fig. 6 , this result is simulated and for w 8 ! 0:17, the system is approaching toward stability. Furthermore, the interference parameter (w 8 ) has also the ability to stabilize the spatial model (2) and Fig. 10 witness this result. It is seen that some individual among a group of prey changes location either in search of resources, food or to avoid the threat of predation. During this movement of prey species, predator senses the insufficiency of food and they also follow them toward the prey rich habitat. The biological conformity to these sequential movements of species produces structure distinguished as the travelling wave, modulated chaotic wave or chaotic behavior. Many different aspects like random dispersal of species, self prey growth, conversion rule, mutual interference and death rate of predators are considered. It is observed that these aspects of the reaction-diffusion predator-prey model impersonate the spatial structure for some period of interference during diffusion instabilities. We observe the Turing patterns especially emphasizing on Beddington-DeAngelis functional response. A mixture of dominating interference pattern and stripe patterns has been observed and as time increases, one can see the emergence of the interference pattern (see Fig. 11 ). From the point of view of physics, if the light from two point sources overlaps, the interference pattern maps out how the phase difference between the two waves varies in space. Here this is determined by the separation of the point sources and the wavelength. An example of an interference pattern is a generation of Newton's ring pattern. We observe that patterns are strongly dependent on interference parameter (w 8 ) and patterns change chaotically from one pattern to another. For a small value of w 8 ¼ 0:1, hot and cold spot patterns can be seen. For instance when w 8 is varied from w 8 ¼ 0:1 to w 8 ¼ 0:25, mixture of spot with stripe patterns are observed in Figs. 12 and 13. Hot spot patterns are formed by the localized annular patches with a high concentration of population density of secondary and tertiary consumers. The appearance of such stripe (or 'zebra') patters is due to auto-catalysis saturation for high prey concentration values. Since prey can not increase further the region of high prey concentration become wider. The model system (2) encounters structural transitions and produces a variety of patterns (see Figs. 10, 11, 12, 13) for interference parameter and diffusion coefficients, indicating the significant importance of these parameters on pattern selection. The idea contained in this paper gives an impeccable approach to examine how functional responses of species and parameters influence on the dynamics of a food chain model. We observe that the diffusion coefficients as well as mutual interference parameter between predators are responsible for the emergence of chaotic dynamics in the ecological model. In this study, we observe the emergence of complex patterns like spots, stripes and interference. This study shows that patterns in ecological interaction are an excellent tool to examine the fundamental mechanism of complexity that arises in a spatiotemporal model. From the study, we can conclude that if predators mutually interfere during predation, chaotic dynamics can automatically emerge in both temporal and spatial model. Consider the following inequality from the system (3) ) dr dt ¼ uð1 À uÞ À w 5 w 6 v À uv w 7 À w 4 ðu þ w 4 Þðu þ w 7 Þ ð29Þ Fig. 13 Spatial distribution of species at the time t ¼ 5000 and w 8 ¼ 0:25 in two dimension ðx À yÞ plane defined in Eq. (26) and diffusion coefficients are given by d 1 ¼ 1:5; d 2 ¼ 0:01; d 3 ¼ 0:01 Iran J Sci Technol Trans Sci Analytical and numerical study of the HIV-1 infection of CD4? T-cells conformable fractional mathematical model that causes acquired immunodeficiency syndrome with the effect of antiviral drug therapy Mutual interference between predators can give rise to Turing spatial patterns Pattern formation, outbreaks, and synchronization in food chains with two and three species Pattern formation in a semi ratiodependent predator-prey system with diffusion Pattern formulation in three species food web model in spatiotemporal domain with Beddington-DeAngelis functional response Existence and non-existence of spatial patterns in a ratio-dependent predator-prey model Turing instability and spatio-temporal chaos in ratio dependence Holling-Tanner model Spatio-temporal pattern formation in Rosenzweig-MacArthur model: effect of non-local interactions New acoustic wave behaviors to the Davey-Stewartson equation with power-law nonlinearity arising in UID dynamics Complex surfaces to the fractional (2 ? 1)-dimensional Boussinesq dynamical model with the local M-derivative Mutual interference between parasites or predators and its effect on searching efficiency Spatial ecology via reaction-diffusion equations Stability of the boundary solution of a nonautonomous predator-prey system with the Beddington-DeAngelis functional response Permanence, extinction and periodic solution of predator prey system with Beddington-DeAngelis functional response A model for trophic interaction Complete mathematical analysis of predator-prey models with linear prey growth and Beddington-DeAngelis functional response New investigation of batshosts-reservoir-people coronavirus model and application to 2019-nCoV system A new study of unreported cases of 2019-nCOV epidemic outbreaks Novel dynamic structures of 2019-nCoV with nonlocal operator via powerful computational technique Complex solitons in the conformable (2?1)-dimensional Ablowitz-Kaup-Newell-Segur equation Dynamics of a marine plankton system: diffusive instability and pattern formation Stationary patterns induced by self-and cross-diffusion in a Beddington-DeAngelis predator-prey model The functional response of predators to prey density and its role in mimicry and population regulation Global analysis of the predator-prey system with Beddington-DeAngelis functional response Spatial pattern and ecological analysis Hypothesis for origin of planktonic patchiness Dynamics of the density dependent predator-prey system with Beddington-DeAngelis functional responses Resonance and frequency-locking phenomena in spatially extended phytoplankton-zooplankton system with additive noise and periodic forces Coexistence in laboratory population of paramecium aurelia and its predator Didinium nasutam Chaos and pattern in a spatial tritropic food chain Spatio-temporal patterns formulation in nonlinear equilibrium plankton dynamics Spatio-temporal complexity of plankton and fish dynamics Dynamical behavior of a three species food chain model with Beddington-DeAngelis functional response Turing patterns and long-time behavior in a three-species foodchain model Predator interference effects on biological control: the 'paradox' of the generalist predator revisited Diffusion induced chaos in a spatial predator-prey system Wave of chaos: new mechanism of pattern formation in spatiotemporal population dynamics Analytical and approximate solutions of a novel nervous stomach mathematical model Dissipative stucture: an explanation and an ecological example Wave of chaos in a diffusive system: generating realistic patterns of patchiness in planktonfish dynamics The role of top predator interference on the dynamics of a food chain model Wave of chaos in a spatial eco-epidemiological system: generating realistic patterns of patchiness in rabbit-lynx dynamics Fluctuations in the abundance of a species considered mathematically A diffusive one-prey and twocompeting-predator system with a ratio-dependent functional response: I, long time behavior and stability of equilibria New dark-bright soliton in the shallow water wave model Dynamics analysis of a Holling I predator-prey system with mutual interference concerning pest control Cross-diffusion induced Turing instability for a three species food chain model Acknowledgements This work is supported by ECR Award (File No. ECR/2017/000141) SERB-DST, New Delhi.Author's contribution SNR: Supervision; conceptualization (lead); writing -original draft (lead); review and editing (equal); PM: methodology (lead); formal analysis (lead); BPS: investigation; writing -review and editing (equal); BT: writing -review and editing (equal). Conflict of interest The authors declare that they have no conflict of interest. À vw w 6 ðv þ w 8 w þ w 9 Þ :Here, the solutions of linearized ODE system (3) are locally Lipschitz continuous functions. Then, if w 7 ! w 4 , the equation (30) transforms to dr dt uð1 À uÞ À w 5 w 6 v;holds for all positive u and v. Further, Eq. (31) is converted to dr dt þ w 5 rðtÞ uð1 À uÞ þ w 5 u:But, max fuð1 À uÞg ¼ 1 4 , hence we get dr dt þ w 5 rðtÞ 1 4 þ w 5 :Now by using comparison lemma,Thus, for all t ! 0, we have rðtÞ 1 4w 5 þ 1 ) uðtÞ þ vðtÞ w 6 1 4w 5 þ 1:If gðtÞ ¼ uðtÞ þ vðtÞ w 6 þ dwðtÞ with d ¼ 1 w 6 w 11 , then for w 7 ! w 4 , we get dg dt u À w 5 w 6 v À w 10 dw;) dg dt þ pg ð1 þ pÞ;where p ¼ minfw 5 ; w 10 g. By using comparison lemma, we obtain for all t !T ! 0,IfT ¼ 0, thenAppendix 2: Proof of Theorem 3is the considered Lyapunov function with constant k 1 and k 2 around the equilibria E 3 ¼ ðu à ; v à ; w Ã Þ for the temporal system (3). Time derivative of Lyapunov function S iswhere 6 v c Þ and,Appendix 3: Proof of Theorem 5Consider the function E(t) aswhere S(u, v, w) is a Lyapunov function defined in Eq. (32). For global stability of system (2), it is sufficient to show that dE dt 0. Hence,where From above analysis, it is clear that I 2 0 and we have only to show I 1 0 to prove dE dt 0. Hence it is proved that the globally stable equilibrium point E 3 ¼ ðu H ; v H ; w H Þ for system (3) will remain globally asymptotically stable for the diffusive system (2).