ISRN Biomathematics

Volume 2013 (2013), Article ID 424062, 12 pages

http://dx.doi.org/10.1155/2013/424062

Research Article

Finite Time Blowup in a Realistic Food-Chain Model

Rana D. Parshad,1,2 Hamid Ait Abderrahmane,² Ranjit Kumar Upadhyay,³ and Nitu Kumari1,4

¹Department of Mathematics and Computer Science, Clarkson University, Potsdam, NY 13699, USA

²King Abdullah University of Science and Technology, Applied Mathematics and Computational Science, Thuwal 23955-6900, Saudi Arabia

³Department of Applied Mathematics, Indian School of Mines, Dhanbad, Jharkhand 826004, India

⁴School of Basic Sciences, Indian Institute of Technology Mandi, Mandi, Himachal Pradesh 175 001, India

Received 23 April 2013; Accepted 19 May 2013

Academic Editors: L. Pezard and A. A. Polezhaev

Copyright © 2013 Rana D. Parshad et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigate a realistic three-species food-chain model, with generalist top predator. The model based on a modified version of the Leslie-Gower scheme incorporates mutual interference in all the three populations and generalizes several other known models in the ecological literature. We show that the model exhibits finite time blowup in certain parameter range and for large enough initial data. This result implies that finite time blowup is possible in a large class of such three-species food-chain models. We propose a modification to the model and prove that the modified model has globally existing classical solutions, as well as a global attractor. We reconstruct the attractor using nonlinear time series analysis and show that it pssesses rich dynamics, including chaos in certain parameter regime, whilst avoiding blowup in any parameter regime. We also provide estimates on its fractal dimension as well as provide numerical simulations to visualise the spatiotemporal chaos.

1. Introduction

Interaction networks in natural ecosystems can be visualized as consisting of simple units known as food chains or food webs that are made up of a number of species linked by trophic interaction [1]. A food chain model essentially comprises of the predator-prey relationship between interacting species in a given ecosystem [2]. In their seminal work [3], Hastings and Powell for the first time demonstrated that the evolution of species participating in a tritrophic relationship might be chaotic [3]. This led to a great deal of research activity in analyzing the dynamical behavior of food chain models. Upadhyay and Rai [4] provided a new example of a chaotic population system in a simple three-species food chain with Holling type II functional response. This model is different from the Hastings and Powell model, in that it considers a generalist top predator, one that can switch its food source, in the absence of its favorite prey. Letellier and Aziz-Alaoui [5] and Aziz-Alaoui [6] revisited the Upadhyay and Rai model and found that the chaotic dynamics is observed via sequences of period-doubling bifurcation of limit cycles which however suddenly break down and reverse giving rise to a sequence of period-halving bifurcation leading to limit cycles. Upadhyay in [7] next proposed a three-species food-chain model, by incorporating mutual interference, in the original model [4], thus generalizing the models in [3, 4, 8]. Parshad and Upadhyay [9] considered the diffusive form of the model proposed in [7]. Under certain restrictions on the parameter space, they proved global existence of solutions as well as the existence of a finite dimensional global attractor, for the diffusive model. The original three-species model for a generalist top predator [4] and its variants are very rich dynamically and have led to a number of works in the literature [4–9]. They continue to be recently investigated, for rich dynamical behavior, such as turing instability [10]. Despite these rich dynamics, the Parshad and Upadhyay model possesses a drawback, in that there is the possibility of finite time blowup in the model. Recall in general we say that the norm of a solution blows up in finite time if

where is a suitable function space where one looks for the solution , to the PDE or ODE at hand. Blow-up phenomenon in physical and biological problems has a rich history and has been studied in the context of gas dynamics, chemotaxis, cooperative, and competitive lotka-Volterra type systems [11–20]. Understanding blowup in a model can help to determine how robust a system is, as well as providing insight on the characteristics of the corresponding physical or biological phenomena.

As the Parshad and Upadhyay model generalises many of the aforementioned three-species models, our result actually implies finite time blowup in all of those models as well. This, however, has not been reported in the literature. It is well documented in the literature [5] that in the absence of the middle predator the top predator becomes extinct if . This essentially means that if the middle predator is absent the top predator has to switch its food source. However, it is probably not able to predate efficiently enough, and thus the population cannot be sustained, unless the rate of sexual reproduction is above a certain threshold. If this is not so, that is, , then the species goes extinct. Also, it grows unboundedly if the inequality reverses. We remark that the situation can be much worse than unbounded growth and that finite time blowup is actually possible.

In the current paper our primary contributions are the following.

(1)We show that the diffusive model of Parshad and Upadhyay [9] can blow up in finite time, for certain parameter regime and large enough initial data, via Theorem 2. This immediately shows the blowup of the diffusive form of the model of Upadhyay and Rai [7]. The temporal model proposed in [4] is also shown to blow up in finite time, in certain parameter regime. (2)We propose a modification to the diffusive model and show global existence of classical solutions and global attractor, for this model, via Theorems 6 and 14. The modified model is well posed and does not blow up in finite time, for any data or in any parameter regime. (3)We reconstruct the attractor of the modified model, via nonlinear time series analysis, and show that it also exhibits rich dynamics, including chaos in certain parameter regime, just like the old model, whilst avoiding finite time blowup in any parameter regime. Numerical simulations are also provided to show the blowup in the classical models, as well as spatiotemporal chaos, in the modified model.

This paper is organized as follows . In Section 2, we recount the diffusive three-species food-chain model [9], that generalises several other known models in the literature. Finite time blowup results are summarized in Section 3. We propose a modification to the general model in Section 4 and show global existence of classical solutions and global attractor for this model. Various estimates made here concerning the middle predator and prey species are very similar to [21]. We reconstruct the chaotic attractor of the modified model, using nonlinear time series analysis, and estimate its fractal dimension. This, as well as some numerical simulations, is presented in Section 5. Some important conclusions are discussed in Section 6.

2. Food Chain Model System

In this section, we recount the model from [9]. Consider a situation where a prey population is predated by individuals of population . The population in turn serves as a favourite food for individuals of population . This interaction is represented by the following system of a simple prey-specialist predator-generalist predator interaction with the inclusion of spatial spread:

where , , , , , , , , , , , , . The parameters for are mutual interference parameters that model the intraspecific competition among predators when hunting for prey [22–27]. The problem is posed on .   is bounded, and is assumed to be smooth. We consider the Dirichlet boundary conditions

We also assume suitable initial conditions. In this model, prey population of size serves as the only food for the specialist predator population of size . This predator population, in turn, serves as a favorite food for the generalist predator population of size . The equations for rate of change of population size for prey and specialist predator have been written following the Volterra scheme; that is, predator population dies out exponentially in the absence of its lone prey. The interaction between this predator and the generalist predator is modeled by the modified version of the Leslie-Gower scheme where the loss in a predator population is proportional to the reciprocal of per capita availability of its most favorite food. is the intrinsic growth rate of the prey population , is the intrinsic death rate of the predator population in the absence of the only food , measures the rate of self-reproduction of generalist predator , and , , , and are the maximum values which per capita growth rate can attain. measures the strength of intraspecific competition among the individuals of the prey species . and quantify the extent to which environment provides protection to the prey and may be thought of as a refuge or a measure of the effectiveness of the prey in evading a predator’s attack. is the value of at which per capita removal rate of becomes . represents the residual loss in population due to severe scarcity of its favourite food . For the coefficient of the third term on the right-hand side of (2) is obtained by considering the probable effect of the density of the prey’s population on predator’s attack rate. If this coefficient is multiplied by (the prey population at any instant of time), it gives the attack rate on the prey per predator. Denote , when , , which is the maximum that it can reach. The third term on the right-hand side of (3) represents the per capita functional response of the invertebrate predator and was first introduced by Holling [28] in the ecological literature. He carried out specially designed experiments to determine the nature of these functional responses. The ecological role of per capita functional response has been well described by May [29]. The interaction terms appearing in the rate equation restore to some extent the symmetry which characterizes the Lotka-Volterra model. The generalist predator , in (4), is sexually reproducing species. We assume that males and females are equal in number and every individual has got equal opportunity to meet an individual of opposite sex. The first term of (4) represents growth rate of the sexually reproducing species in well-mixed conditions. measures the limitation on the growth of the generalist predator by its dependence on per capita availability of its most favorite prey .

The study of ecological problems render’s one to regard communities of individuals as subpopulations affecting the survival of the individuals of other populations. One aspect of the dynamics of community interactions would be mutual interference among the interacting subpopulations, which in our model system are represented by the parameters . It has been shown in [24–26] that mutual interference is generally a “stabilizing” process. The parameters , , and are diffusion coefficients of the populations.

3. Finite Time Blowup

In this section, we show that (2)–(4) blow up in finite time, for certain parameter regime and large enough initial data. In particular, the equation for as described via (4) blows up. To this end we employ the first eigenvalue method for blow up,  [30]. First recall Jensen’s inequality with a probability measure.

Lemma 1. Let be a probability measure on , and let be a nonnegative function. Then the following holds whenever is a convex function:

Note, this is easily generalised to higher dimensions. We next present our blow-up result.

Theorem 2. Consider the three-species food-chain model as described via (2)–(4). For and initial data large enough, that is

there exists a finite time , such, that for any ,

Here, is a function such that on , on . Furthermore, and .

Proof. Consider the equation for rewritten as

Since is assumed positive, we obtain

Thus,

Now we multiply (11) by , consider , and integrate by parts to obtain

Via assumptions on the eigenfunction , it is seen to be a probability measure. That is, we can define a probability measure on and . Now we consider the convex function (note that any works as the in question is always positive). Thus, the application of Jensen’s inequality yields

Thus, substituting the previous equation into (12) yields

Simple integration in time of the previous equation yields that blows up in finite time; that is

where

and .

After some algebra, we obtain that for blowup the data has to satisfy

This proves the theorem.

Remark 3. Note that the blowup of immediately implies that the norm of blows up as follows:

Remark 4. We also show blow-up in the Upadhyay-Rai model proposed in [4], that is, the temporal case with and. Also, for the model in [4], the blow-up time is given by

Remark 5. We note that finite time blowup is also possible in the same weighted norm, for Neumann boundary conditions. In this case, the eigenfunction will have to be differently defined, that is, satisfying Neumann boundary conditions itself.

4. Modification to the Old Model

The aim of the current section is to introduce an improvement to (2)–(4), which is well posed. The key issue in amending the old model, (2)–(4), is to correct (4). Upadhyay and Rai introduced sexual reproduction in (4), by considering . Under this assumption, the number of males and females in the population is assumed to be equal, [5]. Thus, mating takes place as , and consequently the population grows at rate . We point out that this assumption may fail in various structured populations, where due to specific demographics, the numbers of females and males are rarely equal. We modify this by considering the mating term to be , where we restrict . This modification may be interpreted as modeling a population where the number of males is different from the number of females (even if very slightly so). If the numbers of males and females are in fact equal or close to, one can choose . Next, the mutual interference parameter is incorporated, by raising the reaction terms to the power. The primary role of the previously stated modification is that now the reaction term in (4) is , and because we enforce , standard parabolic theory [31] guarantees global existence. Based on the above, we introduce the following three-species food-chain model, with a generalist top predator,

The problem is posed on . is bounded, and is assumed to be smooth. We consider the Dirichlet boundary conditions

Note, the proofs of global existence to follow will work the with Neumann boundary conditions as well. The forms of (2)-(3) are kept essentially the same; however, we assume that , as commonly done in numerical simulations, [5]. We also assume . This facilitates the estimates to follow. Restricting to the range has distinct advantages over the old model (2)–(4). Firstly, finite time blowup is prevented, and this is independent of the sign of . Secondly, in the absence of , if , the dynamics of (2)–(4) is quite boring, and decays to 0, [5, 6]. However, if , even if and , can still persist. Thus, this formulation plausibly models a true generalist . One that does not go extinct, even if its favorite food source does.

4.1. Preliminaries

The system (20)–(22) is an example of a good semilinear parabolic system, and the standard theory for parabolic systems [31] ensures the existence of a nonnegative local solution. Next, comparison arguments can be used to derive global existence for initial data in . We assume that , , and are nonnegative because they are densities. This assumption ensures that ,  ,  ,  and . The nonlinearities are thus well defined. Standard parabolic theory guarantees the uniform boundedness of . We state the following theorem.

Theorem 6. Consider the three-species food-chain model described by (20)–(22). For any positive initial data and , , there exists a global classical solution to the system.

Remark 7. In this section, we demonstrate a stronger version of global existence. That is, we derive uniform bounds on the solution with the initial data only in .

The opposite signed nonlinearities here cause problems in making direct estimates and proving the existence of bounded absorbing sets, for the variable. To circumvent this difficulty, we resort to grouping estimates via appropriate addition of the equations at hand. This technique is very similar to what we adopted in [21]. This condition plays a key role in the dissipation process. If it is satisfied, there is often a global attractor. We will also show the existence of a global attractor for the reaction diffusion system (20)–(22).

Remark 8. Note, in [9] in order to proceed, we need to assume . Thus, our results in [9] are only true in the parameter regime, such that . Furthermore in [9], strong boundedness requirements on the solutions are assumed. Here, we enforce no such requirement, and the methods of analysis are quite different from [9].

Definition 9. Consider a semigroup, , acting on a reflexive Banach space, . Then, the global attractor, , for this semigroup is an object that has the following properties:

(i) is compact in ,(ii) is invariant; that is, ,(iii)if is bounded in , then , .

Next, various preliminaries are presented, detailing the phase spaces of interest and recalling certain standard theories. Let us define our phase spaces of interest as

To prove the existence of a global attractor, we are required to show [32]

(i)the existence of a bounded absorbing set in the phase space; (ii)the asymptotic compactness property of the semigroup in question.

These are defined next.

Definition 10 (bounded absorbing set). A bounded set, , in a reflexive Banach space, , is called a bounded absorbing set if, for each bounded subset of , there is a time, , such that for all . The number is referred to as the compactification time for . This is essentially the time after which the semigroup compactifies.

Definition 11 (asymptotic compactness). The semigroup associated with a dynamical system is said to be asymptotically compact in if, for any bounded in and a sequence of times ,   possesses a convergent subsequence in .

4.2. Global Existence and Existence of Absorbing Sets in

We begin by multiplying (20) by and integrating by parts over . This yields

We then use Hölder’s inequality followed by Young’s inequality to obtain

Next, using the compact Sobolev embedding, , and the positivity of and , we find that

where depends explicitly on and .

Thus, application of Gronwall’s Lemma gives us the following estimate:

Therefore, there exists time given explicitly by

such that, for all , the following estimate holds uniformly:

We now make a local in time estimate for . By integrating (26) in the time interval , we obtain

Thus, via a mean value theorem for integrals, there exists a time such that the following estimate holds:

We now move on to showing the existence of an absorbing set for in . By multiplying (20) by and (21) by , adding the two together, and setting , we obtain

By multiplying (33) by , integrating by parts over , and using Hölder’s inequality followed by Young’s inequality, keeping in mind the positivity of the solutions, we find that

Thus, by using the compact Sobolev embedding, , we obtain

Next, we multiply the previous equation by and integrate in time from to to find

where . This follows via (30).

We now make the following estimate via time integration of (26), after multiplying through by :

This follows via (30).

We now substitute the above into (36) and use the form of to obtain

Therefore, there exists time , given explicitly by

such that, for all , the following estimate holds uniformly:

Therefore,

where is a constant that is independent of the time and the initial data.

We next derive local in time estimates for . We proceed by multiplying (21) by and integrating by parts over to obtain

We now use Hölder’s inequality, positivity of solutions, the embedding (via assumption on ), and the earlier estimate via (40) to yield

By integrating this inequality in the time interval , for , we obtain

The existence of an absorbing set for in follows easily via the form of (22).

Standard estimates imply the existence of a time given explicitly by

such that, for all , the following estimates hold uniformly:

We can now state the following lemma.

Lemma 12. Consider which are solutions to the diffusive three-species food-chain model described via (20)–(22), with , , and . There exists a time and a constant C independent of time and the initial data and dependent only on the parameters in (20)–(22), such that, for any , the following uniform estimates hold:

We now make the estimate. We take the inner product of (20) with and integrate by parts over . Thus, we obtain

We apply the Cauchy-Schwartz inequality, the Cauchy with epsilon, and Young’s inequalities on the last term, along with the embedding (via assumption on ) to obtain

This yields

We now recall the following lemma.

Lemma 13 (the uniform Gronwall lemma). Let , , and be nonnegative functions in . Assume that is absolutely continuous on and the following differential inequality is satisfied:

If there exists a finite time and some such that

for any , where , and are some positive constants, then

We apply the previous lemma to (50) with

and , , and we use the estimates via (31) and Lemma 12 to obtain

The estimates on the other components, and , are obtained similarly.

Recall that in we have via the Sobolev embedding that ,  for all ,

and similarly, for and , we have

Using the previous equation it is easy to show that the reaction terms , , and are in , for . In general, is required [33], but in the current scenario ; thus the condition reduces to . The estimates trivially follow from Hölder’s inequality, applied on the reaction terms along with (56) and (57). This proves the global existence result.

4.3. Existence of Global Attractor in

We can now use the standard theory [32] to prove that the semigroup for (20)–(22) possesses a global attractor. Thus, we can state the following result.

Theorem 14. Consider the three-species food-chain model described by (20)–(22), , and . There exists a global attractor in , for the solution semigroup generated by this model.

Proof. We have shown that the system is well posed via Theorem 6. Thus, there exists a well-defined semigroup . Lemma 12 establishes the existence of bounded absorbing sets in . Thus, given a sequence that is bounded in , we know that, for ,

Here, is the bounded absorbing set in . Now, , for any large enough . Therefore, for , we have

This implies that we have the following uniform bound, via (57):

By standard functional analysis theory [32], a subsequence, still labelled as , exists such that

which implies, via the compact Sobolev embedding of , that

This yields the asymptotic compactness of the semigroup in . Similar analysis is possible for components and . We have the existence of an absorbing set via Lemma 12 and the asymptotic compactness via (62). These in conjunction yield the existence of the global attractor, via standard methods [32].

5. Numerical Simulation and Attractor Reconstruction

We now carry out numerical simulations of (2)–(4). Our goal is to firstly show spatiotemporal chaos, even if . This is indeed the case in certain parameter regime; see Figure 3. We also want to numerically validate the finite time blow-up results in 1D as well as in 2D. For 1D simulation, The system is simulated in MATLAB (R2010) via the PDEPE solver for partial differential equation. The algorithm essentially uses a central difference in space and an implicit time stepping method. All the simulations performed have been refined several times on spatial grids with 300, 600, 900, and 1200 points on the domain of length . These refinements lead to the same general shape and structure of the figures. Table 1 lists the parameter values used in the simulations.

Table 1: Parameters used in simulations.

Figure 1
Figure 2
Figure 3

0.01
1
0.01

0.000001
1
0.000001

0.0001
1
0.0001

2
2
2

0.8
0.8
0.8

0.15
0.15
0.15

0.07
0.07
0.04

2
2
1.97

0.55
0.55
0.55

2
2
2

0.5
0.5
0.5

1.2
1.2
1.2

10
10
10

13
13
13

10
10
10

20
20
20

In order to explore the dynamics of the model in 2D, we use a finite difference method. A forward difference scheme is used for the reaction terms. For the diffusion terms, a standard five-point explicit finite difference scheme is used. The numerical simulation is carried out at different time levels for two dimensional spatial model system. The system of equations is numerically solved over 200 200 mesh points, on a domain of size , with spatial resolution = = 1 and time step . The initial condition used is a small perturbation about , and the boundary conditions used are the no flux Neumann conditions. Values for and are taken to be 200 each. Note this is fine, as the blow-up results hold for Neumann boundary conditions as well. See Figure 1 to visualise the blowup in 1D, and Figure 2 to visualise the blowup in 2D.

(a)

(b)

Figure 1: Blowup is numerically demonstrated for the system (2)–(4), in these plots, when the mutual interference parameters are and. This is exactly the model proposed in [4]. We see blowup in the ODE model in the left, at , and PDE model in the right, at . Note here we take .

Figure 2: Here, finite time blowup for the model system in two dimensional spatial domain is seen. This occurs at time .

(a)

(b)

Figure 3: The densities of the and species are shown as contour plots in the -plane. The corresponding parameters are listed in Table 1. The long-time dynamics (we run up to time ) is seen to be that of spatiotemporal chaos. This shows that one can have interesting dynamics even in the case that , whilst avoiding finite time blowup.

Our next goal is to quantify the chaotic dynamics present in the system. The dynamics of the three-species model described via (2)–(4) is depicted through the reconstruction of the trajectories in phase space with the aide of delay-time method. This method, proposed by Takens and Mane [34] ensures equivalence between the topological properties of actual and reconstructed attractors. Numerical simulations were conducted for three values of the mutual interference parameter . That is , , and . In all cases, was set equal to 2. The system of PDE (2)–(4) was solved again via the MATLAB routine PDEPE. We then fixed certain points in space, and recorded the trajectory of the solution in time, yielding a time series. Very interesting behavior was observed; this consists of suppression of chaos. In fact for a value of , chaos is found when in Table 1 in Upadhyay et al. [35]; see Figure 4(a).

(a)

(b)

(c)

Figure 4: The figures above show the dependence of the dynamics on . For , one sees chaotic dynamics with the largest Lyapunov exponent equal to 2.23. Reducing to 1.9 brings the largest lyapunov exponent equal to 0.26. Further reducing it to 1.2 brings the dynamics to a limit cycle.

We next calculated the dimension of this attractor, using a combination of singular-value decomposition and the Grassberger-Procaccia algorithm, [36]. We essentially want to embed the time series in an dimensional embedding space of embedded vectors, which take the typical form

Here, is called the “delay-time," is the sampling interval, and is the “window length" which represents the time spanned by each embedding vector. This algorithm entails choosing and the delay, so that the window is sufficiently large. We then perform a singular value decomposition and calculate the correlation dimension, from the correlation integral. We continue to refine the embedding dimension till we have maximized the straight line parts in the log-log plot of the correlation integrals. The dimension of this attractor is approximately equal to 2.3 while the value of the largest Lyapunov exponent, estimated with the method proposed by Lai and Chen [37], is equal to 2.23. This clearly indicates that the dynamics is chaotic. Decreasing slightly the value of to 1.9, the dynamics remained chaotic and the dimension of the attractor remained almost unchanged, that is, equal to while the largest Lyapunov exponent had decreased to 0.26; see Figure 4(b). This tells us that the divergence of close trajectories becomes now very slow. The important point here is that the proposed modification to the Parshad-Upadhyay model and subsequently the classical Upadhyay-Rai models still exhibits chaos in certain parameter regime. That is for . Note since , this parameter set is not very realistic biologically. However, when we decrease further the value of to , the dynamics of the system changes radically and becomes oscillatory with single period; this is depicted by a limit cycle in phase space; see Figure 4(c). From this numerical analysis, it seems plausible that the mutual interference parameter plays a crucial role in the dynamics of the three-species model. Setting a value of the parameter (), very different from the classical one , suppresses the chaos predicted by the model when . However, keeping the parameters close to 2, that is , still retains chaotic dynamics, whilst eliminating the possibility of finite time blowup.

6. Discussion

Various nonlinear effects, like structural instability, dissipative structures, dynamical chaos, and finite time blowup, do exist in food chain models. However, the problem of detecting these phenomena in real ecosystems is far from being well understood. Various studies have suggested that the biology of the top generalist predator has an important role to play as far as the dynamics of food-chain models is concerned [10, 38]. The nature of interactions between the populations and communities may be responsible for chaos evading its capture in the wild. In this work, we have shown that a large class of three-species food-chain models [4, 7, 9] can blowup in finite time, for certain parameter regime and large enough initial data. We reiterate that the situation can be worse than unbounded growth if , as reported earlier [5]. We propose a possible correction based on the analysis of the way mating is modeled in (4). That is, if we restrict to the range , then finite time blowup is prevented, and this is independent of the sign of . However, if and , the top predator can still persist, for . Thus, this formulation plausibly models a true generalist predator , one that does not go extinct, even if its favorite food source does. Finally, by varying in the range , one can adapt the new model to populations where the numbers of males and females vary. If the numbers are very close, we take , . One could also model a population where the males far outnumber the females or vice versa. Here one might take , saying we have roughly 10 females per 100 males in a population or vice versa. Thus, making the case for the model in structured populations, where such demographics may be dominant. Also, simulations performed by varying the mutual interference parameter indicate clearly that there is a bifurcation point , in the parameter space, such that, chaotic dynamics is seen for and oscillatory dynamics for . It would be interesting to obtain the actual value of the mating parameter from real field data. Comparing this “real value" to might then be a means to again answer the long standing question of why chaos has evaded its capture in the wild. Lastly, we remark that the blow-up result obtained via Theorem 2 is sufficient but not necessary. It is not optimal in terms of data. That is, it might still be possible for blowup to occur if the data does not meet the largeness condition

For the classical diffusive Upadhyay-Rai model [4], we show blowup numerically; see Figure 1. Here, we have precisely all the parameters previously stated, including the first eigenvalue and eigenfunction of the Dirichlet laplacian. We see that, for the parameters chosen to demonstrate blowup in Figure 1,

Thus blowup takes place, even for data smaller than that required by the largeness condition. It would thus be very interesting to try and sharpen this largeness requirement on the data, perhaps via employing certain other techniques available in the blow-up literature, than the version of the first eigenvalue method that we have resorted to currently.

References

R. K. Upadhyay, “Observability of chaos and cycles in ecological systems: lessons from predator-prey models,” International Journal of Bifurcation and Chaos, vol. 19, no. 10, article 3169, 65 pages, 2009.

J. D. Murray, Mathematical Biology, Springer, New York, NY, USA, 1993.

A. Hastings and T. Powell, “Chaos in a three-species food chain,” Ecology, vol. 72, no. 3, pp. 896–903, 1991.

R. K. Upadhyay and V. Rai, “Why chaos is rarely observed in natural populations,” Chaos, Solitons and Fractals, vol. 8, no. 12, pp. 1933–1939, 1997.

C. Letellier and M. A. Aziz-Alaoui, “Analysis of the dynamics of a realistic ecological model,” Chaos, Solitons and Fractals, vol. 13, no. 1, pp. 95–107, 2002.

M. A. Aziz-Alaoui, “Study of a Leslie-Gower-type tritrophic population model,” Chaos, Solitons and Fractals, vol. 14, no. 8, pp. 1275–1293, 2002.

R. K. Upadhyay, “Chaotic dynamics in a three species aquatic population model with Holling type II functional response,” Nonlinear Analysis. Modelling and Control, vol. 13, no. 1, pp. 103–115, 2008.

R. K. Upadhyay, S. R. K. Iyengar, and V. Rai, “Chaos: an ecological reality?” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 8, no. 6, pp. 1325–1333, 1998.

R. D. Parshad and R. K. Upadhyay, “Investigation of the long time dynamics of a diffusive three species aquatic model,” Dynamics of Partial Differential Equations, vol. 7, no. 3, pp. 217–244, 2010.

N. Kumari, “Pattern formation in spatially extended tritrophic food chain model systems: generalist versus specialist top predator,” ISRN Biomathematics, vol. 2013, Article ID 198185, 12 pages, 2013.

N. Mizoguchi, H. Ninomiya, and E. Yanagida, “Diffusion-induced blowup in a nonlinear parabolic system,” Journal of Dynamics and Differential Equations, vol. 10, no. 4, pp. 619–638, 1998.

J. Morgan, “On a question of blow-up for semilinear parabolic systems,” Differential and Integral Equations, vol. 3, pp. 973–978, 1990.

H. F. Weinberger, “An example of blowup produced by equal diffusions,” Journal of Differential Equations, vol. 154, no. 1, pp. 225–237, 1999.

H. A. Levine, “Role of critical exponents in blowup theorems,” SIAM Review, vol. 32, no. 2, pp. 262–288, 1990.

P. Souplet, “A note on Diffusion-induced blow up,” Journal of Dynamics and Differential Equations, vol. 19, no. 3, pp. 819–823, 2007.

R. Ferreira, A. De Pablo, F. Quirós, and J. D. Rossi, “The blow-up profile for a fast diffusion equation with a nonlinear boundary condition,” Rocky Mountain Journal of Mathematics, vol. 33, no. 1, pp. 123–146, 2003.

M. Chlebik and M. Fila, “On the blow-up rate for the heat equation with a nonlinear boundary condition,” Mathematical Methods in the Applied Sciences, vol. 23, pp. 1323–1330, 2000.

Z. Lin, “Blowup estimates for a mutualistic model in ecology,” Electronic Journal of Qualitative Theory of Differential Equations, vol. 8, no. 1, pp. 1–14, 2002.

G. Whitham, Linear and Nonlinear Waves, John Wiley & Sons, New York, NY, USA, 1974.

T. Hillen and K. J. Painter, “A user's guide to PDE models for chemotaxis,” Journal of Mathematical Biology, vol. 58, no. 1-2, pp. 183–217, 2009.

R. D. Parshad, A. Kasimov, and H. Abderrahmane, “Long time behavior and the Turing instability in a diffusive three species food chain model,” Under Review

J. R. Beddington, “Mutual interference between parasites on predators and its effects on searching efficiency,” Journal of Animal Ecology, vol. 44, pp. 331–340, 1975.

L. H. Erbe and H. I. Freedman, “Modeling persistence and mutual interference among subpopulations of ecological communities,” Bulletin of Mathematical Biology, vol. 47, no. 2, pp. 295–304, 1985.

L. H. Erbe, H. I. Freedman, and V. S. H. Rao, “Three-species food-chain models with mutual interference and time delays,” Mathematical Biosciences, vol. 80, no. 1, pp. 57–80, 1986.

H. I. Freedman, “Stability analysis of a predator-prey system with mutual interference and density-dependent death rates,” Bulletin of Mathematical Biology, vol. 41, no. 1, pp. 67–78, 1979.

H. I. Freedman and V. Sree Hari Rao, “The trade-off between mutual interference and time lags in predator-prey systems,” Bulletin of Mathematical Biology, vol. 45, no. 6, pp. 991–1004, 1983.

M. P. Hassell, “Mutual interference between searching insect parasites,” Journal of Animal Ecology, vol. 40, pp. 473–486, 1971.

C. S. Holling, “The functional role of invertebrate predators to prey density,” Memoirs of the Entomological Society of Canada, vol. 45, pp. 3–60.

R. M. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, NJ, USA, 2001.

L. Evans, Partial Differential Equations, vol. 19 of Graduate Studies in Mathematics, Springer, New York, NY, USA, 2nd edition, 1998.

C. V. Pao, Nonlinear Parabolic and Elliptic Equations, Plenum, New York, NY, USA, 1993.

R. Temam, Infinite Dimensional Dynamical Systems in Mechanics and Physics, vol. 68 of Applied Mathematical Sciences, Springer, New York, NY, USA, 2nd edition, 1997.

D. Henry, Geometric Theory of Semi-Linear Parabolic Equations, vol. 840 of Lecture Notes in Mathematics, Springer, New-York, NY, USA, 1984.

F. Takens and R. Mane, “Dynamical systems and turbulence,” in Warwick, R. Rand and L. S. Young, Eds., vol. 898 of Lecture Notes in Mathematics, Springer, Berlin, Germany, 1981.

R. K. Upadhyay, S. R. K. Iyengar, and V. Rai, “Stability and complexity in ecological systems,” Chaos, Solitons and Fractals, vol. 11, no. 4, pp. 533–542, 2000.

P. Grassberger and I. Procaccia, “Dimensions and entropies of strange attractors from a fluctuating dynamics approach,” Physica D, vol. 13, no. 1-2, pp. 34–54, 1984.

D. Lai and G. Chen, “Statistical analysis of lyapunov exponents from time series: a jacobian approach,” Mathematical and Computer Modelling, vol. 27, no. 7, pp. 1–9, 1998.

V. Rai and R. K. Upadhyay, “Chaotic population dynamics and biology of the top-predator,” Chaos, Solitons and Fractals, vol. 21, no. 5, pp. 1195–1204, 2004.