Abstract

A nonlinear system describing the interaction between toxin-producing phytoplankton and zooplankton was investigated analytically and numerically, where the system was represented by a couple of reaction-diffusion equations. We analyzed the effect of self- and cross-diffusion on the system. Some conditions for the local and global stability of the equilibrium were obtained based on the theoretical analysis. Furthermore, we found that the equilibrium lost its stability via Turing instability and patterns formation then occurred. In particular, the analysis indicated that cross-diffusion can play an important role in pattern formation. Subsequently, we performed a series of numerical simulations to further study the dynamics of the system, which demonstrated the rich dynamics induced by diffusion in the system. In addition, the numerical simulations indicated that the direction of cross-diffusion can influence the spatial distribution of the population and the population density. The numerical results agreed with the theoretical analysis. We hope that these results will prove useful in the study of toxic plankton systems.

1. Introduction

In marine ecosystems, most aquatic life relies on plankton, which comprises phytoplankton and zooplankton. Phytoplankton comprises most of the primary energy sources in aquatic food webs, and it accounts for a large proportion of the world’s fixed production. Phytoplankton is consumed by zooplankton, which provides food for fish and other aquatic animals. In fact, the phytoplankton can also render very useful service by producing a huge amount of oxygen for other animals after absorbing carbon dioxide from environments [1]. Thus, plankton forms the basis of all aquatic food chains and it has an essential role in the study of marine ecology [2, 3]. However, the biomass of phytoplankton population may be of rapid increase or almost equally rapid decrease, and this phenomenon of rapid change in phytoplankton population is called “bloom” at some fixed time. Because some phytoplankton can produce toxin and the accumulation of high biomass, some of these blooms are known as “harmful algal blooms” [4] which can have toxic effects on marine ecosystems or even human health, thereby causing great socioeconomic damage. Thus, research into bloom dynamics is widespread with a special emphasis on harmful algal blooms.

In the past two decades, there have been major increases in harmful plankton blooms in aquatic ecosystems [57]. Studies have shown that there are at least eight different modes and mechanisms that allow harmful phytoplankton species to cause mortality, physiological impairment, or other negative in situ effects [8]. It is well known that the toxin-producing phytoplankton has important impacts on the growth of the zooplankton, and thus studies of marine plankton are ubiquitous and significant, while the dynamic behavior of interacting species in the marine is also a major topic [9]. However, researchers have paid less attention towards toxin-producing plankton blooms in recent years [10, 11], although some have recognized the role of toxin-producing phytoplankton in reducing the grazing pressure due to zooplankton [1215]. In [14], a three-species model comprised a toxin-producing phytoplankton and zooplankton and fish population was studied by Upadhyay and Chattopadhyay, who employed Holling type I, II, and III functional responses to describe the liberation of the toxin, where they concluded that the toxin-producing phytoplankton made a stabilizing contribution to aquatic systems. In [15], the dynamical behaviors of toxin-producing phytoplankton and zooplankton were investigated, where the phytoplankton was divided into two groups, that is, susceptible and infected phytoplankton. Of course, given the global increase in harmful plankton blooms in aquatic ecosystems, the study of the effects of toxic phytoplankton has begun to become a growing concern in recent years [1619].

In marine ecosystems, the diffusion phenomenon of population exists widely, which can impact the distribution of population. Many factors can influence the diffusion, such as turbulence and foraging. Actually, it is indicated that the diffusion induced by mixing play an important role in the population dynamics of the phytoplankton [20]. And climate models predict that global warming will increase the stability of vertical stratification [21] and reduce vertical mixing [22], which means that climate change will affect the population dynamics. Hence, in order to understand how climate changes affect the population dynamics, the study of the diffusion mechanism of population becomes much more important. In this aspect, the reaction-diffusion equation may be a useful tool to study the spatiotemporal dynamics induced by diffusion [23, 24].

The dynamics of interacting populations with self- and cross-diffusion have been studied widely [2527]. In particular, Dubey et al. [25] analyzed a mathematical model of a predator-prey interaction with self- and cross-diffusion and obtained the criteria for local stability, instability, and global stability. Another study [26] considered a nutrient-plankton model of an aquatic environment in the context of phytoplankton blooms to analyze the diffusion-driven instability and stability, as well as cross-diffusion of zooplankton, under the influence of phytoplankton in the spatial model. The results [26] indicated the influence of cross-diffusion of zooplankton, and thus we can consider the effects of more complex cross-diffusion in marine ecosystems. In the previous work, some researchers have studied the dynamical properties of the toxic plankton system with diffusion [2830]. Jang et al. studied a mathematic model which described the phytoplankton-zooplankton interactions with toxin-producing phytoplankton; the study revealed that passive diffusion of both populations can simplify the dynamics of the interactions and exhibit plankton patchiness [28]. In [29], authors studied a delay-diffusion model of marine plankton ecosystem, which exhibited cyclic nature of blooms. In [30], Roy considered a nontoxic phytoplankton-toxin-phytoplankton-zooplankton model which is described by a reaction-diffusion equation. The study demonstrated that spatial movements of planktonic systems in the presence of TPP generate and maintain inhomogeneous biomass distribution of competing phytoplankton, as well as grazer zooplankton. These works have studied the influence of diffusion on toxin plankton system and obtained some good results. However, the influence of cross-diffusion on the toxin plankton ecosystem was seldom considered. In order to study the influence of cross-diffusion in the toxic plankton ecosystem, we propose a toxin-phytoplankton-zooplankton system with self- and cross-diffusion.

2. Model and Linear Stability Analysis

2.1. Model Analysis

In this study, we consider a toxin-phytoplankton-zooplankton system with Holling type II response as follows:where is the density of the phytoplankton population and is the density of the zooplankton population at time . is the phytoplankton natural birth rate; is the phytoplankton natural mortality rate, and thus is the intrinsic growth rate of the phytoplankton. denotes phytoplankton intraspecies competition. is Holling type II response, is the half-saturation constant for a Holling type II functional response, is the rate of predation for zooplankton, is the zooplankton’s natural mortality rate, is the zooplankton intraspecies competition, is the rate of biomass consumption by zooplankton during its growth, and denotes the rate of toxin liberation by the toxin-producing phytoplankton population. We make a basic assumption that . The parameters , , , , , , , ,  and are all positive constants.

Considering the relationship between the climate and the diffusion of species and the fact on 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 phytoplankton population by studying the spatial dynamics of the diffusion system. The spatial system can be described aswhere and are the constant diffusion coefficients for the phytoplankton and zooplankton, respectively, and and are the cross-diffusion coefficients for the phytoplankton and zooplankton, respectively. Biologically, cross-diffusion implies countertransport and it means that the prey exercised a self-defense mechanism to protect against attack by a predator [10], different from the self-diffusion, and the values of and may be positive or negative. We take the zero-flux boundary condition of the plankton populations and in ,   as follows:and the initial conditions arewhere is a bounded spatial domain with smooth boundary .

In addition, we takeas the diffusive matrix. The determinant of is . We assume that the self-diffusion is stronger than cross-diffusion; thus diffusive matrix is positive definite.

Let and be the right hand sides of the equations in (1). It is obvious that the equilibria and always exist.

By , we can obtain the vertical isoclines . This is a downward parabola, so the axis of symmetry is , and the fellowship with the -axis is . By , we can obtain the horizontal isoclines . It is easy to show that the line is the asymptote of . We can prove that increases monotonically when , and the fellowship with the -axis is . Based on the analysis above, we defineand the following conclusions can be established.

Lemma 1. (1) If , , then there exists one positive equilibrium point in system (1), where .
(2) If , then there exists one positive equilibrium point in system (1); it is obvious that .
(3) If and there at least exists to guarantee but , then there exist three positive equilibria , , and in system (1), where ,  .

2.2. Linear Stability Analysis

For and , we can find that is always unstable and that is locally asymptotically stable when .

Now, let be an arbitrary positive equilibrium, and thus satisfies the equations and . Then, the Jacobian matrix for is given byNote thatFrom , we can find that if , then , and thus and ; in this case, the equilibrium is always locally asymptotically stable. But if , then .

Next, we analyze the stability and instability of system (2). Let () be the eigenvalues of on with the zero-flux boundary condition. Setand consider , where is the eigenspace that corresponds to .

In the condition where , we obtain the following conclusions.

Lemma 2. If , then the equilibrium is always locally asymptotically stable.

Theorem 3. If , , , and , then the equilibrium is uniformly asymptotically stable when , .

Proof. The linearization of system (2) at the positive equilibrium iswhere ().
For each (), is invariant under the operator , and is an eigenvalue of on if and only if is an eigenvalue of the matrixNote thatWhen , then , , and the two eigenvalues and have negative real parts. For any , we find the following.
If , then , and if , since and , thenwhere is independent of .
Thus, a negative constant exists, which is independent of , such that for any . By referring to [31], we can prove that the positive equilibrium of system (1) is uniformly asymptotically stable when , .

Theorem 4. (1) If , then the equilibrium is globally asymptotically stable when .
(2) If , , , , then the positive equilibrium is globally asymptotically stable when and , where

Proof. (1) Define a Lyapunov functionWe note that is nonnegative and if and only if .
Furthermore,By substituting the expressions for and from system (1), we obtainUsing the condition that (17) can be simplified as, and always occurs when .
Therefore, if holds.

Proof. (2) Then, by referring to [32, 33], we choose the following Lyapunov function:and by differentiating with respect to time along the solutions of system (2), we can obtain Using Green’s first identity in the plane, we obtainwhereBased on the analysis above,In the case of increasing , to sufficiently large values and if , meet the conditionsthen , and thus . Therefore of system (2) is globally asymptotically stable when and .

Under the condition of we have the following conclusions.

Lemma 5. If , then the equilibrium is locally asymptotically stable when and the equilibrium is unstable when .

Proof. Define under the condition of Lemma 1, we can find that .
Thusand combined with , we can find that , andHence, if , then , and in this case ; thus, is stable. If , then and in this case, ; thus is unstable.
By analyzing Lemma 5, we can find that is the value where the plankton system populations bifurcate into periodic oscillation. By , then we can obtainwhich is a bifurcation value of the parameter . Now, the characteristic equation of the reduces towhere , and the equation has purely imaginary roots. We take . If the stability change occurs at , it is obvious that , , and using , , we get . Thus, the transversality condition for a Hopf bifurcation is satisfied.

Theorem 6. If , , and , then the criterion for the Turing instability of system (2) satisfies the following condition . Consider

Proof. The linearized form of system (2) corresponding to the equilibrium is given bywhere ,  , and are small perturbations in about the equilibrium . We expand the solution of system (34) into a Fourier series:where is the wave number of the solution. By combining (34) and (35), we obtainHence, we can obtain the characteristic equation as follows:whereIn order to allow the system to achieve Turing instability, at least one of the following conditions must be satisfied: and . However, it is obvious that . Thus only the condition that exists can give rise to instability. Therefore, a necessary condition for the system to be unstable is thatWe define If we let be the corresponding value of for the minimum value of , then Thus, the corresponding minimum value of is In addition, if we let , the sufficient condition for the system to be unstable reduces toNoting. If there exists three positive equilibria , , and in system (1), where , , then a bistable system exists whenWe take , , , , , , , . Table 1 show that the condition for Lemma 1(3) is achievable.

3. Numerical Results

3.1. Numerical Analysis

In this section, system (2) is analyzed using the numerical technique to show its dynamic complexity. However, the existence of a positive equilibrium and its stability in system (1) are given first before this. Table 2 shows the existence of the equilibrium and its stability for different values of the parameter , while the other parameters remain fixed: ,,,,,, and Table 2 indicates that the positive equilibrium is unstable when , but it is locally asymptotically stable when and global asymptotically stable when .

The parameters are taken as follows: , , , , , , , and .

Figure 2 presents a series of one-dimensional numerical solutions of system (2). Figure 2(a) shows that the solution of system (2) tends to a positive equilibrium, where , which implies that positive equilibrium is stable. This means that the self-diffusion does not lead to the occurrence of instability. In fact, as shown in Figure 1, the self-diffusion does not cause instability in system (2) when . However, instability may occur when the cross-diffusion coefficients satisfy some conditions.

Therefore, different values of and are employed to illustrate the instability induced by cross-diffusion. Figure 2(b) shows the phytoplankton density, where the values of and are negative, whereas the values of and are positive in Figure 2(d). It is obvious that the spatial distributions of phytoplankton in Figures 2(b) and 2(d) are markedly different. In ecology, positive cross-diffusion indicates that one species tends to move in the direction with a lower concentration of another species, whereas negative cross-diffusion denotes that the population tends to move in the direction with a higher concentration of another species [30]. In Figures 2(c) and 2(f), , and , , respectively, which means that only the influence of zooplankton cross-diffusion is considered in the system. The spatial distributions of phytoplankton are not obviously different in Figures 2(c) and 2(f), but the amplitudes of the phytoplankton density differ, which means that the direction of zooplankton cross-diffusion can affect the phytoplankton density. In Figure 2(e), coefficient is positive and is negative, and it is not difficult to see that the spatial distribution of phytoplankton differs significantly from the other panels shown in Figure 2. These results demonstrate that cross-diffusion plays a significant role in the stability of the system, but they also show the effects of the various types of cross-diffusion that a plankton system may encounter.

3.2. Pattern Formation

In order to illustrate the spatial distribution of phytoplankton more clearly, we show the patterns formed by system (2) in a two-dimensional space in Figure 3.

Figure 3 shows two different spatial distributions of phytoplankton. Figure 3(a) shows the initial spatial distributions of phytoplankton. Figures 3(b) and 3(c) illustrate different types of pattern formation for different values of and at . When , , Figure 3(b) shows that a spotted pattern and the spots appear to be cluttered with an irregular shape. However, when , , a circular spotted pattern appears, and the distribution of the spots is more regular, as shown in Figure 3(c). As discussed in Section 3.1, cross-diffusion will lead to instability in the system, where both positive cross-diffusion and negative cross-diffusion play important roles in the stability of the system. The spatial distributions of phytoplankton shown in Figures 2(b) and 2(c) agree well with the results given in the previous section.

4. Conclusion and Discussion

The system is simple because it is only an abstract description of a marine plankton system. However, the system exhibited some rich features of the marine plankton system. The results show that the spatial distribution of phytoplankton is relatively uniform and stable under some conditions if we only consider the self-diffusion in system (2). When the influence of cross-diffusion is considered, however, the distribution of phytoplankton will change significantly, as shown in Figure 2. Our results demonstrate that the cross-diffusion can greatly affect the dynamic behavior of plankton system.

In Section 2, we showed that complex equilibria exist in the system and that a bistable system exists when the equilibrium satisfies certain conditions. The results shown in Table 1 verified these conclusions. We also found the condition for Hopf bifurcation. In Table 2, we take parameters values of , ,  , , , , , and , and we illustrated the changes in stability by varying the half-saturation constant for Holling type II functional response in system (1). Although the system is simple, the behavior of the system was found to be very rich.

In contrast to previous studies of marine plankton systems, diffusion system (2) considered the mutual effects of self-diffusion and cross-diffusion. Furthermore, various effects of the saturation of cross-diffusion were considered in our study and we analyzed the influence of differences in real saturation. In the paper, theoretical proof of the stability of the system, we considered positive or negative or zero cross-diffusion and we found that the conditions for system instability were caused by cross-diffusion. In our numerical analysis, we focused on the influence of cross-diffusion on the system. Figure 2 shows the influence of cross-diffusion on the system, where we considered the complexity of the actual situation, and the results are illustrated for different values of and . In Figures 2(a)2(f), the results of our theoretical derivation are verified, and they also reflect the complex form of the system under the influence of cross-diffusion. Turing first proposed the reaction-diffusion equation and described the concept of Turing instability in 1952 [34], and pattern formation has now become an important component when studying the reaction-diffusion equation. Thus, based on Figure 2, we analyzed the patterns formed to further illustrate the instability of the system as in Figure 3.

In addition, it is known that climate change can affect the diffusion of population in marine. Our studies show that the diffusion of phytoplankton and zooplankton can affect the spatial distribution of population under some conditions. Hence, climate change may have a significant effect on the population in marine. We hope that this finding can be further studied and proved to be useful in the study of toxic plankton systems.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 31370381).