# A Reaction Diffusion Model to Describe the Toxin Effect on the Fish-Plankton Population.

1. IntroductionUnderstanding the functioning of an ecosystem is a major issue for resource and environmental management. However, this goal remains difficult to attain due to the complexity of natural systems, especially in the aquatic environment, where many processes of all types interact with living organisms [1-3]. Plankton is the basis of all aquatic food chains and phytoplankton in particular occupies the first trophic level and the fluctuations in its abundance determine the production of marine environment. Rapid increase and almost equally rapid decrease separated by periods are the common features of plankton populations. In a broad sense, planktonic blooms can be divided into two types: "spring blooms" and "red tides." Spring blooms occur seasonally for the changes in temperature or nutrient availability which are connected with seasonal changes in thermocline depth, strength, and consequent mixing. Red tides are localized outbreaks and occur due to high water temperature [4-7]. Toxic substances produced by phytoplankton species reduce the growth of zooplankton by decreasing grazing pressure and this is one of the important common phenomena in plankton ecology [2, 6, 8, 9].

Within the broad perspective drawn above, the present paper explores and compares the coupled dynamics of the phytoplankton and the zooplankton in a number of mathematical models. The system phytoplankton-zooplankton has attracted considerable attention from various fields of research [3, 8,10,11]. It is an important issue in mathematical ecology. The literature abounds in models focusing on various aspects of the problem. Recently, the attention has been focused on the role of the space in explaining heterogeneity and the distribution of the species and the influence of the spatial structure on their abundance [2, 6, 8]. However, the very question of the interactions between phytoplankton, zooplankton, and fish depending on space is far from being fully elucidated. As part of our work, we will highlight a threshold value of the toxin released by phytoplankton, below which the effect of the toxin influences less the dynamics of the zooplankton-fish system. To our knowledge, this has never been addressed in the literature. The proposed model consists of three interactive components, zooplankton, herbivorous fish, and toxin-producing phytoplankton, which reduce the growth of fish and zooplankton population. The model studied here is of the reaction-diffusion type, describing the dynamics of the phytoplankton-zooplankton-fish system in the sense of the work of Courchamp et al. [12].

The paper is organized as follows. As far as Section 3 is concerned, we will establish mathematical results such as the existence of a solution, stability of equilibria, and persistence, relating to the constructed model in the Section 2. Section 4 will be devoted to numerical experiments to illustrate the mathematical results. Finally, Section 5 is devoted to the conclusion and perspectives.

2. Mathematical Model

In this section, we propose a model to describe the dynamics of the phytoplankton-zooplankton-fish system in the presence of toxin. We begin our modeling by a general model describing the dynamics of the prey-competitor-predator system, based on the equations with ordinary derivatives. And then we transform this model into a model of reaction-diffusion type while remaining in the logic of the works of Courchamp et al. [12] and Bendahmane [13]. It should be noted that the aim is to take into account the effect of the toxin on the fish population through that of zooplankton.

2.1. Original Model Formulation. If we designate by P the density of the prey, by Z the density of the competitor, and by F that of the predator, according to [12], the general model is written as follows:

[mathematical expression not reproducible], (1)

where

(i) [[phi].sub.2], [g.sub.5], [g.sub.6], [g.sub.7], [g.sub.8], [g.sub.9], [g.sub.10], and [g.sub.11] are positive [C.sup.[infinity]] functions,

(ii) [[phi].sub.2](P) is the growth function of the prey,

(iii) [g.sub.5](P,Z) is the function that measures the mortality due to the competition between the same individuals of the prey,

(iv) [g.sub.6](P, Z) is the function that measures the mortality of the prey due to the competitor's consumption,

(v) [g.sub.7](P, Z) is the competitor's growth function,

(vi) [g.sub.8](P, Z, F) is the function that measures the mortality of the competitor from the release of the substance produced by the prey,

(vii) [g.sub.9](P, Z, F) is the function that measures the mortality of the competitor from the consumption of the predator,

(viii) [g.sub.10](Z, F) is the growth function of the predator,

(ix) [g.sub.11](Z,F) is the function that measures the external mortality of the predator population.

We continue our modeling by fixing the expressions of the functions intervening in model (1). The dynamics of the system are represented by Figure 1.

According to Figure 1, for any time T >0, the dynamics of the phytoplankton (prey)-zooplankton (competitor)-fish (predator) system is governed by the following ODE system:

[mathematical expression not reproducible], (2)

where

(i) [r.sub.p] and [[mu].sub.[epsilon]] denote, respectively, the phytoplankton growth rate and the phytoplankton external mortality rate,

(ii) [b.sub.0] and [[alpha].sub.0] denote, respectively, the mortality due to competition between the individuals of the phytoplankton population and the specific predation rate of the zooplankton on the phytoplankton population,

(iii) [[gamma].sub.0] is the half-saturation constant of zooplankton predation on the phytoplankton,

(iv) [r.sub.z] denotes the external mortality rate of zooplankton,

(v) [[theta].sub.p] denotes the rate of toxin production per phytoplankton species,

(vi) [[alpha].sub.1] and [[alpha].sub.2] represent, respectively, the maximum value that the consumption rate of the phytoplankton by the zooplankton can reach and the maximum value that the reduction rate by the zooplankton can reach,

(vii) [[gamma].sub.2] is the value of zooplankton for which the elimination rate by the individual zooplankton becomes [[alpha].sub.2]/2,

(viii) [r.sub.f] and [[mu].sub.f] represent, respectively, the growth rate and the external mortality rate of the fish,

(ix) [[alpha].sub.3] is the maximum value that the rate of reduction by the fish can reach,

(x) [[gamma].sub.3] represents the residual loss caused by high rarity of the zooplankton of individuals of the fish population.

2.2. A Spatially Structured Model. Considering the relationship between the climate and the diffusion of species and the fact of the existence of diffusion in population, system (1) is developed into a spatial system with diffusion. We expect to explore the effect of climate change on plankton population by studying the spatial dynamics of the diffusion system. We will introduce the concept of spatial structure in the model; that is to say, the population's densities depend now on the time and space. Diffusion models are a simple and reasonable choice for modeling dispersion of populations on a spatial domain [14, 15]. Indeed, let [[delta].sub.1](x), [[delta].sub.2](x), and [[delta].sub.3](x) be, respectively, the diffusivity coefficients of P, Z, and F. Based on the work established in [12], the structured model associated with model (1) can be modeled for x [member of] [OMEGA] and T > 0 as follows:

[mathematical expression not reproducible], (3)

where [OMEGA] [subset] [R.sup.n] (n [greater than or equal to] 2) is the spatial domain in which species occur and zero-flux boundary condition.

[mathematical expression not reproducible], (4)

for Q = P, Z, F, where v is the unit normal vector to [OMEGA] along [partial derivative][OMEGA] and nonnegative initial conditions

Q (x, 0) = [Q.sub.0] (x) > 0, Q = P, Z, F, x [member of] [OMEGA]. (5)

So, model (2) can be written by the following system:

[mathematical expression not reproducible]. (6)

We make the following assumptions:

([H.sub.a]) All demographic parameters of system (6) are positive constants.

([H.sub.b]) The diffusivity coefficients of system (6) are independent of x.

Let us consider [[delta].sub.1](x) = [[delta].sub.1], [[delta].sub.2](x) = [[delta].sub.2], and [[delta].sub.3](x) = [[delta].sub.3] in the remaining of our work.

Taking into account hypotheses ([H.sub.a]) and ([H.sub.b]), model (6) becomes

[mathematical expression not reproducible]. (7)

Remark 1. The model presented here is a model of predator prey system with a self-diffusion but it is also possible to extend these types of system to systems with a cross-diffusion in the same logic of the works of Andreianov et al. [16, 17].

3. Mathematical Results

The mathematical results are based on the works in [3, 16, 18-20].

3.1. Reduction of Model Parameters. To simplify the writing, we will make the following changes of variables:

[mathematical expression not reproducible]. (8)

Thus, we obtain, respectively, systems (2) and (7) of the following forms:

[mathematical expression not reproducible], (9)

[mathematical expression not reproducible]. (10)

3.2. Existence and Boundedness of Solution. Before stating the boundedness of the solution, we give the following theorem concerning the existence and uniqueness of the solution of system (7).

Theorem 2 (see [16, 21-23]). System (7) has a unique local solution ([U.sub.1] (*, t), [U.sub.2] (*, t), [U.sub.3] (*, t)) under the condition 0 [less than or equal to] t < [T.sub.max], where [T.sub.max] depends on nonnegative initial data [U.sub.01] (x), [U.sub.02](x), and [U.sub.03](x).

The following theorem ensures the boundedness of the solution.

Theorem 3. The set D = [0, 1] x [0, c + c/4b] x [0,1/[alpha]a+1/4[alpha]b+ M/[alpha]b] attracts all global solutions ofsystem (7) initiatingin the positive orthant, where [alpha] = 4b/[(b + p).sup.2](4bc + c + 4br).

Proof. By noticing that [U.sub.1] [less than or equal to] [E.sub.1], [U.sub.2] [less than or equal to] [E.sub.2], and [U.sub.3] [less than or equal to] [E.sub.3] and from [10, 19], we have [lim.sub.t[right arrow]+[infinity]][E.sub.1](t) [less than or equal to] 1 and [lim.sub.t[right arrow]+[infinity]] ([E.sub.1]) + (1/c)[E.sub.2](t)) [less than or equal to] 1 + 1/4b. Finally, we get that

[mathematical expression not reproducible]. (11)

The following theorem gives the existence of global solution.

Theorem 4 (see [21-23]). System (10) admits a global solution ([U.sub.1](*,t), [U.sub.2](*,t), [U.sub.3] (*,t)) for any t > 0 and for any regular initial positive functions [U.sub.01](x) [less than or equal to] 1, [U.sub.02](x), and [U.sub.03](x).

Proof. In fact,

(i) on one hand, we have [U.sub.1](x,t) [greater than or equal to] 0, [U.sub.2](x,t) [greater than or equal to] 0, and [U.sub.3](x, t) [greater than or equal to] 0 because 0 is the subsolution of each equation of system (10);

(ii) on the other hand, [U.sub.1] satisfies the following problem:

[mathematical expression not reproducible]. (12)

According to the principle of comparison, we have [U.sub.1] (x, t) [less than or equal to] U(t) [less than or equal to] 1, where U(t) = [U.sub.01][e.sup.t]/(1 - [U.sub.01] + [e.sup.t]) is the solution of the system [mathematical expression not reproducible]. (13)

In the same order, [U.sub.2] satisfies d[U.sub.2]/[partial derivative]t = c[U.sub.1][U.sub.2]/ ([U.sub.1] + a) - b[U.sub.2] - [U.sub.2][U.sub.3]/([U.sub.2] + d) + [[epsilon].sub.2][DELTA][U.sub.2].

According to the principle of comparison, we have (1/c)([partial derivative][U.sub.2]/[partial derivative]t) [less than or equal to] (1/c)(d[E.sub.2]/dt), where [E.sub.2] is the solution of the second equation of system (2) with initial condition [E.sub.2](0) = [max.sub.[bar.[GAMMA]]] [U.sub.02](x).

According to [10, 19], we have (1/c)([partial derivative][U.sub.2]/[partial derivative]t) [less than or equal to] (1/c)(d[E.sub.2]/dt) + d[E.sub.1]/dt. Let us denote [sigma] = (1/c)[E.sub.2] + [E.sub.1]; then (1/c)([partial derivative][U.sub.2]/[partial derivative]t) [less than or equal to] d[sigma]/dt [less than or equal to] 1/4 + b - b[sigma]. Therefore, we have [U.sub.2](x, t) [right arrow] c + c/4b, 0 [right arrow] t < T.

From the third equation of system (10), we find [alpha]([partial derivative][U.sub.3]/[partial derivative]t) = [alpha]p[U.sub.3] - [alpha]q[U.sub.3.sup.2]/([U.sub.2] + r) + [alpha][[epsilon].sub.3][DELTA][U.sub.3]. From the comparison principle, we can deduce that [alpha]([partial derivative][U.sub.3]/[partial derivative]t) [less than or equal to] [alpha](d[E.sub.3]/dt), where [E.sub.3] is the solution of the third equation of system (2) with initial condition [E.sub.3] (0) = [max.sub.[bar.[OMEGA]]][U.sub.3](x). Therefore,

[mathematical expression not reproducible]. (14)

Then,

[mathematical expression not reproducible]. (15)

And thus,

[mathematical expression not reproducible]. (16)

3.3. Analysis of Stationary Solutions. The positive equilibria are the solutions ([U.sub.1](x), [U.sub.2](x), and [U.sub.3](x)) of system (7), where [U.sub.3](x) [greater than or equal to] 0, [U.sub.2](x) [greater than or equal to] 0, and [U.sub.3](x) [greater than or equal to] 0. Suppose that there exist two positive constants [m.sub.1] < 1 and [m.sub.2] such that [U.sub.1] (x) [greater than or equal to] [m.sub.1] > 0, [U.sub.2](x) [greater than or equal to] [m.sub.2] >0, and [U.sub.3](x) > 0. The equilibrium of system (7) satisfies

[mathematical expression not reproducible]. (17)

We make the following assumptions:

[mathematical expression not reproducible]. (18)

Theorem 5. Let assumptions ([H.sub.a]), ([H.sub.1]), ([H.sub.2]), and ([H.sub.3]) hold; system (17) has at least one positive solution.

Proof. We will construct a pair of subsolution and oversolution of problem (17):

[mathematical expression not reproducible], (19)

which verifies

[mathematical expression not reproducible]. (20)

By fixing [bar.[U.sub.1]] = 1 > [[U.bar].sub.1] = [m.sub.1] and [[U.bar].sub.2] = [m.sub.2], according to Theorem 3 and under hypothesis ([H.sub.1]), we have [[bar.U].sub.2] [less than or equal to] c+c/4b, and c + c/4b [less than or equal to] (([m.sub.1] + a)/[m.sub.1])[m /(1 + a) + [m.sub.1](1 - [m.sub.1])]. This implies that [[bar.U].sub.2] [less than or equal to] (([m.sub.1] + a)/[m.sub.1])[[m.sub.2]/(1+a)+[m.sub.1](1- [m.sub.1])]. Thus, the right-hand term of the first inequality of (20) is satisfied. Let [[bar.U].sub.1] = 1, and to satisfy the left term of the second inequality of system (20), we need [[U.bar].sub.3] to satisfy the following inequality: c/(1 + a) - b - [[U.bar].sub.3]/([[bar.U].sub.2] + d) [less than or equal to] 0. Let [[U.bar].sub.1] = [m.sub.1] and [[U.bar].sub.2] = [m.sub.2]. According to Theorem 3 and under hypothesis ([H.sub.2]), we have 0 [less than or equal to] [[bar.U].sub.3] [less than or equal to] 1/[alpha] + c/4b[alpha] + M/[alpha] and 1/[alpha] + c/4b[alpha] + M/[alpha] [less than or equal to] [cm.sub.1]/([m.sub.1] + a) - b - [[bar.U].sub.3]/([m.sub.2] + d).

Therefore,

[[cm.sub.1]/[m.sub.1] + a] - b - [[[bar.U].sub.3]/[m.sub.2] + d] [greater than or equal to] 0. (21)

Thus, the right-hand side of the second inequality of system (20) is satisfied. To satisfy the right term of the third inequality of system (20), we need [[U.bar].sub.3] [less than or equal to] (p/q)([m.sub.2] + r). From Theorem 3 and hypothesis ([H.sub.3]), this last condition is satisfied.

3.4. Stability of Homogeneous Stationary Solutions. Let us consider a solution U(x,t) of system (7) in the following form:

[mathematical expression not reproducible], (22)

where [U.sub.s] (x) is the positive solution of (17). Let us replace U(x, t) by (22) in system (10) and, by identifying all terms that are linearized in J, we have

[mathematical expression not reproducible], (23)

with

[mathematical expression not reproducible]. (24)

We assume the following assumptions:

[mathematical expression not reproducible]. (25)

Theorem 6. Let ([H.sub.4]), ([H.sub.5]), and ([H.sub.6]) hold; then the equilibrium Us is locally stable.

Proof. Concerning the proof of this theorem, we use the results established in [22, 23]. In fact, we consider the following system:

[mathematical expression not reproducible], (26)

where [[chi].sub.j] is the eigenfunction of -[DELTA] in [OMEGA] and the eigenvalue [[beta].sub.j] satisfying 0 < [[beta].sub.0] < ... < [[beta].sub.j] < ....

So, we have

[mathematical expression not reproducible]. (27)

From (27) and (23), we have [[??].sub.j] = [C.sub.j][J.sub.j], where [C.sub.j] is the matrix defined as

[mathematical expression not reproducible]. (28)

[U.sub.s] is stable if and only if each [J.sub.j](t) decreases to zero. So, [C.sub.j] has three eigenvalues with negative real parts [[xi].sub.k], k = 1, 2, 3, such that [[xi].sub.3] + a[[xi].sub.2] + b[xi] + c = 0, with

[mathematical expression not reproducible]. (29)

For the three eigenvalues to admit a negative real part, we must necessarily have

[mathematical expression not reproducible]. (30)

For j = 0, from (30), we have

[mathematical expression not reproducible]. (31)

According to the Routh-Hurwitz criterion [24], sufficient stability conditions for positive stationary solutions can be obtained.

We will establish the stability of the homogeneous equilibrium ([U.sup.*.sub.1], [U.sup.*.sub.2], [U.sup.*.sub.3]) under certain conditions. In fact, we have the following hypothesis:

[mathematical expression not reproducible], (32)

Thus, we have the following theorem.

Theorem 7. Suppose that the hypotheses ([H.sub.a]), ([H.sub.1]), ([H.sub.2]), and ([H.sub.7]) are satisfied; then ([U.sup.*.sub.1], [U.sup.*.sub.2], [U.sup.*.sub.3]) is globally asymptotically stable.

Proof. Let us define the function

[mathematical expression not reproducible]. (33)

We have

[mathematical expression not reproducible], (34)

where

[mathematical expression not reproducible]. (35)

Under hypothesis ([H.sub.7]), the matrix S is positive, and consequently [T.sub.1] < 0.

On the other hand, we have

[mathematical expression not reproducible]. (36)

4. Numerical Experiments

In this section, we perform extensive numerical simulations of the spatial model system (7) in two-dimensional space using the forward finite difference method. The set of parameter values used for the numerical simulation are given in Table 1 [11, 15-17, 25]. Here, the system is studied on a spatial domain: [OMEGA] = [0, 50] x [0, 50]. It is assumed that the fish, zooplankton, and phytoplankton are spread over the whole domain at the beginning.

4.1. Pattern Formation. Here, we will illustrate the mathematical predictions, by numerical simulations, concerning the behavior of the dynamics under hypotheses (H1)-(H7). The qualitative results of different pattern formations due to the variation of [b.sub.0] and [[alpha].sub.0] are shown. Here, we consider the value of released toxin [[theta].sub.p] = 0.13. These numerical results show that, for every strictly positive initial condition, under assumptions ([H.sub.1])-([H.sub.7]), the nonhomogeneous equilibrium is always globally asymptotically stable. Figures 2-7 show the spatial structures formation for the threes species described in (7). These numerical results confirm the mathematical results about the existence and the stability of positive equilibrium according to the values of [b.sub.0] and [[alpha].sub.0]. In this case, we will speak of a subsistence phenomenon of the fish and zooplankton population.

Remark 8. From a biological point of view, these results (Figures 2-7) show that there is coexistence between the three populations despite the release of the toxin into the aquatic environment. This means that, despite the harmful effects of the toxin released by phytoplankton, fish and zooplankton populations persist.

4.2. Analysis of the Dynamics Behavior with Toxin Effect. In this section, we consider that [b.sub.0] = 0.2562 and [[alpha].sub.0] = 0.71. Figures 8-11 show that, after a transitional phase, an equilibrium can be established with a strong distribution of the three populations, which is not far from the results obtained previously. As a biological interpretation, we can say that if the toxin is released below this value, the impact is not significant on the fish and zooplankton populations (Figures 8-10). A less dense distribution of fish and zooplankton than before is observed. This explains the considerable decrease of these species due to the increase in the number of toxic phytoplankton. On the other hand, we observe in Figure 11 a sudden change concerning the spatial distribution of phytoplankton (strong and bad spatial distribution).

Remark 9. Figure 12 shows the estimate of the error by using the standard norm [[parallel] * [parallel].sub.2] on the iterations [16, 17, 26]. These figures show a decrease in error estimation over the course of the iterations. This confirms the convergence of the numerical scheme used.

5. Conclusion and Future Work

Our interest in this paper is the construction of a reaction-diffusion model for modeling the dynamics of fish population under a diet on plankton (zooplankton and phytoplankton), taking into account the effect of the toxin. The model formulation is derived from an ODE system by considering an isotropic distribution as in [12]. It should be noted that we consider diffusion independently of the spatial variable in the construction of the reaction-diffusion model. The mathematical results allowed us to establish conditions of existence of equilibrium which depend on the demographic parameters. We have also given some results about the stability of the stationary equilibria and we have established the conditions on the inexistence of the equilibrium with strictly positive components. We continued our study through numerical experiments in order to confirm our mathematical results. Also, our numerical results have yielded interesting results on the effect of the toxin on the dynamics. This is why we are led to conclude that the release of the toxin under certain conditions in the aquatic space contributes to the regulation of the system. The phytoplankton bloom was observed during our simulations and is in perfect agreement with the biological observations.

Despite important results obtained in this paper, in order to further our study, we can consider for future work subdividing the phytoplankton population into toxic phytoplankton and nontoxic phytoplankton to extend our results to a system of cross-diffusion type.

https://doi.org/10.1155/2018/2037093

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors' Contributions

All the authors read and approved the final manuscript.

References

[1] M. Kot, Elements of Mathematical Biology, Cambridge University Press, Cambridge, UK, 2001.

[2] L. Matthews and J. Brindley, "Patchiness in plankton populations," Dynamics and Stability of Systems, vol. 12, no. 1, pp. 39-59, 1997.

[3] W. Ouedraogo, B. Sangare, and S. Traore, "A mathematical study of cannibalism in the fish-plankton model by taking into account the catching effect," Advanced Modeling and Optimization, vol. 18, no. 2, pp. 197-216, 2016.

[4] E. Briand, Contribution a la comprehention du determinisme de la mise en place des proliferation du cyanobacteries et leur production de toxine [PhD. Thesis], Universite Paris, 2008.

[5] M. L. Rosenzweig and R. H. MacArthur, "Graphical representation and stability conditions of predator-prey interactions," The American Naturalist, vol. 97, no. 895, pp. 209-223, 1963.

[6] J. Steele and E. Henderson, "The significance of interannual variability," in Towards a Model of Ocean Bio-geochemical Processes, G. T. Evans and M. J. R. Fasham, Eds., pp. 237-260, Springer, 1993.

[7] B. Traore, B. Sangare, and S. Traore, "A mathematical model of Malaria transmission with structured vector population and seasonality," Journal of Applied Mathematics, vol. 2017, Article ID 6754097, 15 pages, 2017.

[8] A. M. Edward and J. Brindley, "Zooplankton mortality and the dynamical behavior of plankton population models," Bulletin of Mathematical Biology, pp. 303-339, 1999.

[9] C. Kohlmeier and W. Ebenhoh, "The stabilizing role of cannibalism in a predator-prey system," Bulletin of Mathematical Biology, vol. 57, no. 3, pp. 401-411, 1995.

[10] M. A. Aziz Alaoui and M. Daher Okaye, "Boundness and global stability for a predator-prey model with modified Leslie Gower and Holling-type 2," Applied Mathematical Letters, vol. 16, no. 7, pp. 1069-1075, 2003.

[11] W. Ouedraogo, B. Sangare, and S. Traore, "Some mathematical problems arising in biological models: a predator-prey model fish-plankton," Journal of Applied Mathematics & Bioinformatics, vol. 5, no. 4, pp. 1-27, 2015.

[12] F. Courchamp, M. Langlais, and G. Sugihara, "Rabbits killing birds: modelling the hyperpredation process," Journal of Animal Ecology, vol. 69, no. 1, pp. 154-164, 2000.

[13] M. Bendahmane, "Analysis of a reaction-diffusion system modeling predator-prey with prey-taxis," Networks and Heterogeneous Media, vol. 3, no. 4, pp. 863-879, 2008.

[14] B. Ainseba, F. Heiser, and M. Langlais, "A mathematical analysis of a predator-prey system in a highly heterogeneous environment," Differential and Integral Equations, vol. 15, no. 4, pp. 385-404, 2002.

[15] S. Gaucel, M. Langlais, and D. Pontier, "Invading introduced species in insular heterogeneous environments," Ecological Modelling, vol. 188, no. 1, pp. 62-75, 2005.

[16] B. Andreianov, M. Bendahmane, and R. Ruiz-Baier, "Analysis of a finite volume method for a cross-diffusion model in population dynamics," Mathematical Models and Methods in Applied Sciences, vol. 21, no. 2, pp. 307-344, 2011.

[17] G. Galiano, M. A. Garzon, and A. Jungel, "Semi-discretization in time and numerical convergence of solutions of a nonlinear cross-diffusion population model," Numerische Mathematik, vol. 93, no. 4, pp. 655-673, 2003.

[18] B. I. Camara, "Complexite de dynamiques de modeles proiepredateur avec difusion et applications. Mathematics," Mathematics Universite du Havre, 2009 (French).

[19] M. Daher Okiye and M. A. Aziz-Alaoui, "On the dynamics of a predator-prey model with the Holling-Tanner functional response," in Proceedings of the ESMTB Conference, MIRIAM Editions, pp. 270-278, 2002.

[20] G. Kreisselmeier and R. Engel, "Nonlinear observers for autonomous Lipschitz continuous systems," Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 48, no. 3, pp. 451-464, 2003.

[21] H. Amann, "Dynamic theory of quasilinear parabolic equations. I. Abstract evolution equations," Nonlinear Analysis. Theory, Methods & Applications, vol. 12, no. 9, pp. 895-919, 1988.

[22] D. G. Aronson and H. F. Weinberger, "Nonlinear diffusion in population genetics, combustion, and nerve pulse propagation," Lecture Notes in Math, vol. 446, pp. 5-49, 1975.

[23] J. Smoller, Shock Waves and Reaction-Diffusion Equations, Springer, New York, NY, USA, 1983.

[24] M. M. Fahmy and J. O'Reilly, "A note on the Routh-Hurwitz test," Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 27, no. 2, pp. 483-485, 1982.

[25] J. Dockery, V. Hutson, K. Mischaikow, and M. Pernarowski, "The evolution of slow dispersal rates: a reaction diffusion model," Journal of Mathematical Biology, vol. 37, no. 1, pp. 61-83, 1998.

[26] J. Pousin and J. Rappaz, "Consistency, stability, a priori and a posteriori errors for Petrov-Galerkin methods applied to nonlinear problems," Numerische Mathematik, vol. 69, no. 2, pp. 213-231, 1994.

Wendkouni Ouedraogo, Hamidou Ouedraogo, and Boureima Sangare (iD)

Department of Mathematics and Informatics, UFR/ST, Universite Nazi Boni, 01 BP 1091, Bobo-Dioulasso 01, Burkina Faso

Correspondence should be addressed to Boureima Sangare; mazou1979@yahoo.fr

Received 10 December 2017; Accepted 5 March 2018; Published 10 May 2018

Academic Editor: Kenneth Karlsen

Caption: Figure 1: Typical ecological situation presented by food-chain model.

Caption: Figure 2: Spatial distribution of the system with [b.sub.0] = 0.25 and [[alpha].sub.0] = 0.6.

Caption: Figure 3: Spatial distribution of the system with [b.sub.0] = 0.38 and [[alpha].sub.0] = 0.54.

Caption: Figure 4: Spatial distribution of the system with [b.sub.0] = 0.4032 and [[alpha].sub.0] = 0.342.

Caption: Figure 5: Spatial distribution of the system with [b.sub.0] = 0.2562 and [[alpha].sub.0] = 0.71.

Caption: Figure 6: Spatial distribution of the system with [b.sub.0] = 0.61 and [[alpha].sub.0] = 0.69.

Caption: Figure 7: Spatial distribution of the system with [b.sub.0] = 0.69 and [[alpha].sub.0] = 0.595.

Caption: Figure 8: Dynamics behavior for [[theta].sub.p] = 0.18.

Caption: Figure 9: Dynamics behavior for [[theta].sub.p] = 0.29.

Caption: Figure 10: Dynamics behavior for [[theta].sub.p] = 0.46.

Caption: Figure 11: Dynamics behavior for [[theta].sub.p] = 0.935.

Caption: Figure 12: Corresponding convergence for error estimation in function of the iterations, where the simulated time is t = 50: (a) for P (b) for Z, and (c) for F.

Table 1: Numerical values for the parameters of system. Parameter Description Value [[gamma].sub.0] Half-saturation constant of the Z 0.7 predation on P [[gamma].sub.2] The value of Z (elimination rate by the 0.4 individual of Z) [[gamma].sub.3] The residual loss caused by high rarity 0.9 of Z of individuals of F [r.sub.p] Growth rate of the phytoplankton 1.6 population [r.sub.z] External mortality rate of the 0.4 zooplankton population [[gamma].sub.1] Diffusivity coefficient of the 0.03 phytoplankton population [[gamma].sub.2] Diffusivity coefficient of the 0.03 zooplankton population [[mu].sub.[epsilon]] External mortality rate of 0.9 phytoplankton population [[delta].sub.3] Diffusion coefficient of the fish 0.05 population [[alpha].sub.1] Consumption rate of the phytoplankton 0.75 population by Z [[alpha].sub.2] Consumption rate of the zooplankton 0.4 population by F [[alpha].sub.3] Mortality rate of the fish population 0.122 due to the competition [[mu].sub.f] External mortality rate of the fish 0.02 population [r.sub.f] Growth rate of the fish population 0.4

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Ouedraogo, Wendkouni; Ouedraogo, Hamidou; Sangare, Boureima |

Publication: | Journal of Mathematics |

Date: | Jan 1, 2018 |

Words: | 5023 |

Previous Article: | On the Characterization of Some Classes of Four-Dimensional Matrices and Almost B-Summable Double Sequences. |

Next Article: | Around Chaotic Disturbance and Irregularity for Higher Order Traveling Waves. |

Topics: |