Discrete Dynamics in Nature and Society

Volume 2015 (2015), Article ID 823026, 8 pages

http://dx.doi.org/10.1155/2015/823026

## Dynamical Analysis of a Nitrogen-Phosphorus-Phytoplankton Model

^{1}School of Mathematics and Information Science, Wenzhou University, Wenzhou, Zhejiang 325035, China^{2}Zhejiang Provincial Key Laboratory for Water Environment and Marine Biological Resources Protection, Wenzhou University, Wenzhou, Zhejiang 325035, China^{3}School of Life and Environmental Sciences, Wenzhou University, Wenzhou, Zhejiang 325027, China

Received 7 October 2014; Accepted 23 November 2014

Academic Editor: Garyfalos Papashinopoulos

Copyright © 2015 Yunli Deng 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

This paper presents a nitrogen-phosphorus-phytoplankton model in a water ecosystem. The main aim of this research is to analyze the global system dynamics and to study the existence and stability of equilibria. It is shown that the phytoplankton-eradication equilibrium is globally asymptotically stable if the input nitrogen concentration is less than a certain threshold. However, the coexistence equilibrium is globally asymptotically stable as long as it exists. The system is uniformly persistent within threshold values of certain key parameters. Finally, to verify the results, numerical simulations are provided.

#### 1. Introduction

In marine ecology, phytoplankton play a major role in nutrient cycling, primary production, and global carbon cycling. As is well known, nutrient validity is a necessary factor in phytoplankton population growth. However, the phytoplankton population will breed massively when nutrient input is excessive, eventually leading to eutrophication. Eutrophication can cause water-quality deterioration and fish killing and affect people’s health and recreational activities [1]. It is characterized by frequently recurring algal blooms and reducing species diversity in water bodies at all trophic levels [2]. Therefore, reducing nutrients and eutrophication of water ecosystems is a crucial environmental problem throughout the world.

In 1932, Bertalanffy first proposed the use of mathematical models to study biological systems [3]. Then, some scientific researchers attempted to investigate the biological population dynamics use of mathematical models [4–6]. With the increasing prevalence of eutrophication and algae blooms, the history of mathematical modeling of plankton dynamics and biological eutrophication removal processes is already quite long and has been initiated by the biological sciences. Some approaches have been refined to provide more realistic descriptions of the development of biological natural populations. For instance, ecological models, including impulsive [7–9], diffusion [10, 11], and time delay [12–14], have been taken into account, which can explain certain phenomena in realistic world. In recent years, most efforts have focused on how to control eutrophication, how to predict algae outbreaks, and how to simulate algae spreading tendencies. In this context, many researchers have discussed the dynamic behavior of phytoplankton blooms [15–22]. Huppert et al. [23] presented a simple nutrient-phytoplankton model and explored the dynamics of phytoplankton blooms. Pei et al. [24] investigated a two-zooplankton and one-phytoplankton model with harvesting, which considered the impact of harvesting on the coexistence and competitive exclusion of competitive predators. Mukhopadhyay and Bhattacharyya [25] dealt with a nutrient-plankton model in an aquatic environment in the context of phytoplankton blooms. Zhang and Wang [26] considered a nutrient-phytoplankton-zooplankton model in an aquatic environment and analyzed its global dynamics. Fan et al. [27] proposed a new dynamic nutrient-plankton model and used it to study the relationship between nutritional enrichment and water-quality oscillations. However, researchers have paid little attention to modeling nitrogen-phosphorus-phytoplankton systems.

The Sanyang wetland of Wenzhou is located in a subtropical area and draws on the Wenzhou Economic Development zone and the Longwan zone to the east, linking up with the Chasan Higher Education zone to the south and connecting with the Wenzhou city center area to the northwest; its total area is 13 square kilometers. A certain number of rivers are distributed in a crosswise pattern, forming more than 160 different sizes and shapes of small islands; the proportion of land to water is 1.1 : 1. 47% of the land area is used to grow Mandarin oranges and 15.2% for housing, with the remainder used as agricultural land and fallow land. For the development of Wenzhou, the Sanyang wetland has played an important role in providing water resources, climate regulation, water conservation, flood and drought control, degradation of pollutants, and protection of biodiversity. However, the Sanyang wetland is facing the threat of industrial pollution and sewage, and its formerly crystal-clear water has become a large contaminated area, with local water areas colored red or black and emitting foul odors. Judging from the results of the overall analysis, the quality of the water environment in the Sanyang wetland has been severely damaged, with indicators of nitrogen, phosphorus, and heavy metals seriously exceeding limits. This situation has resulted in frequent nuisance algal blooms, which cause clogging and blocking of filtration systems. An even more frightening threat is that, with the economic development of Wenzhou city, the land and water bodies in the Sanyang wetland are in danger of being heavily invaded by industrial and building land uses as well as agricultural reclamation projects. Therefore, researching on how to enhance the protection of Sanyang wetland natural ecosystems and how to achieve a significant ecological effect on this environment is particularly important and urgent.

Nitrogen and phosphorus are necessary nutrient for plants to live. When small quantities of nutrients flow into a wetland, phytoplankton start to grow. If the process is allowed to proceed, blooms will break out. Laukkanen and Huhtala [1] stated that nitrogen and phosphorus are the primary factors limiting algae blooms, and therefore these two nutrients are considered in the present model. This paper presents a nitrogen-phosphorus-phytoplankton model and uses it to study the interaction between nutrient runoff and phytoplankton growth.

This paper considers a nitrogen-phosphorus-phytoplankton model with Holling type II functional response. The basic model is described by the following ordinary differential equations: where and are, respectively, the density of total nitrogen (mg/L) and total phosphorus (mg/L) at time , is the biomass of the phytoplankton population (mg/L) at time , , are, respectively, the runoff of nitrogen and phosphorus into the wetland, , are, respectively, the natural removal of nitrogen and phosphorus from the water, , are the maximum uptake rates, , are biomass conversion constants, is the natural death rate of the phytoplankton population, , are half-saturation constants, and the terms and represent the response function for nutrient uptake by phytoplankton. The system satisfies the following initial conditions: , , .

This paper aims to obtain a theoretical result for which the values of the bifurcation parameter, the phytoplankton-eradication equilibrium point, and the coexistence equilibrium point are asymptotically stable. To verify the results, numerical simulation has been used to study the controlling relations on the bifurcation parameter .

This paper is organized as follows. In the next section, the theorem governing the positivity and boundedness of solutions is analyzed. In Section 3, the conditions for existence and stability of equilibria are obtained. In Section 4, the uniform persistence of system (1) is examined. Finally, numerical simulations are described in Section 5, and discussion and conclusions are given in Section 6.

#### 2. Positivity and Boundedness of the Solution

In this section, assuming the biologically meaningful initial conditions , , the theorem for the positivity and boundedness of system (1) is stated and proved.

Theorem 1. *Under the given initial conditions, all solutions of system (1) are positive and uniformly bounded.*

*Proof. *From the first equation of system (1),
Hence, .

Similarly, from the second equation of system (1),
From the third equation of system (1), it can be obtained that
Now let us define a Lyapunov function: .

Let .

Then,
The right-hand side of the inequality is bounded , and therefore
Moreover, as . Hence, by the definition of , there are three positive constants , and and such that , for . This completes the proof.

#### 3. Existence and Stability of Equilibria

In this section, the existence of all possible nonnegative equilibria is first discussed.

Obviously, the phytoplankton-eradication equilibrium exists in system (1). Next, in order to research the coexistence of three populations, let us consider the existence of the positive equilibrium.

Theorem 2. *There exists a coexistence equilibrium of system (1) if
*

*Proof. *If the coexistence equilibrium exists, it must satisfy the following three equations:
From (8), we must have and if the coexistence equilibrium exists.

From the first and the second equations of system (8), the following second degree equation in respect of is obtained:
It is easy to calculate that
Hence, there exist two roots in (9) and , as follows:
where . Obviously, and in terms of expression (11). Thus, we just need to take into consideration because of . Let . Then,

Due to , obviously, . Therefore, is a strictly increasing and continuous function of . Moreover, . From (9), we have when ; thus, is continuous at the point . In addition, ; hence, satisfy (9) for .

From the third equation of system (8), it can be obtained that

Because , is a strictly decreasing and continuous function of .

Furthermore,
The vertical asymptote is
In addition,

Based on the analysis of (14), (15), and (16), we know that and may intersect only if .

Let , where . Because , when the condition holds, then . Whereas the condition is inconsistent with the coexistence equilibrium , the condition holds.

Obviously, if (7) holds. Thus, there must be unique , and such that the expression (8) holds, because is a strictly increasing continuous function and is a strictly decreasing continuous function, where , , and is defined as the unique solution of .

Therefore, when holds, there exists a unique coexistence equilibrium in system (1). The proof is complete.

*Next, the stability of the phytoplankton-eradication equilibrium will be discussed, where and .*

*Theorem 3. The phytoplankton-eradication equilibrium is locally asymptotically stable if .*

*Proof. *For , the characteristic equation is
The roots of the characteristics equation are
It is easy to verify that if .

Therefore, the phytoplankton-eradication equilibrium is locally asymptotically stable if . The proof is complete.

*Theorem 4. The phytoplankton-eradication equilibrium is globally asymptotically stable if .*

*Proof. *To investigate the global stability of system (1) at , let us construct the following Lyapunov function:
where are constants to be determined in subsequent steps and , .

Along the trajectories of system (1),
Now, choosing and ,

Note that and are increasing functions and that the solution of (1) is positive. Meanwhile, it can be determined that is equal to . Therefore, each item on the right-hand side of (21) is nonpositive; that is, and if and only if , , and. Based on the Lyapunov-LaSalle theorem, if , then the phytoplankton-eradication equilibrium is globally asymptotically stable. Hence, the proof is complete.

*Finally, the stability of the coexistence equilibrium will be discussed.*

*Theorem 5. The coexistence equilibrium is locally asymptotically stable if .*

*Proof. *For , the characteristic equation is
That is,
where
Note that and are both positive. Based on the Routh-Hurwitz criterion, is locally asymptotically stable if and only if . In fact, after some computations,
Hence, the coexistence equilibrium is locally asymptotically stable if exists. The proof is complete.

*Theorem 6. The coexistence equilibrium is globally asymptotically stable if .*

*Proof. *To investigate the global stability of system (1), at , let us construct the following Lyapunov function:
where are constants to be determined in subsequent steps and , .

Along the trajectories of system (1),
Because , + and choosing and , then

The proof of this is similar to the proof of the global asymptotical stability of . It can be established that . Based on the Lyapunov-LaSalle theorem, the coexistence equilibrium is globally asymptotically stable if exists. Hence, the proof is complete.

*4. Uniform Persistence*

*This section will discuss the uniform persistence of system (1). As discussed below, threshold expressions of some key parameters can be obtained under the condition of all-species persistence. The first step is to state the following lemma.*

*Lemma 7. If , , , and , then + for .*

*Theorem 8. If , then system (1) is uniformly persistent.*

*Proof. *From Theorem 1, there are positive constants , , and and such that , , for .

On the other hand, from the first equation of system (1), it can be determined that
Noting that Lemma 7 is satisfied, then
Therefore, there are positive constants , , and such that
From the second equation of system (1), it follows that . Therefore, there are positive constants , , and such that for .

From the third equation of system (1), it follows that . Therefore, it is evident that if , there are two positive constants and such that for .

Let ; then, it can be obtained that , , and for . Hence, the proof is complete.

*5. Simulation Analysis and Results*

*The stability of the phytoplankton-eradication equilibrium point and the coexistence equilibrium point of system (1) has been demonstrated in the previous section. In this section, numerical simulations will be used to analyze the dynamic behavior of system (1). is chosen as the bifurcation parameter. Table 1 provides the values of the other fixed parameters and their units, which have been obtained from previous studies [27–29].*