Next Article in Journal
Photo-Catalytic Remediation of Pesticides in Wastewater Using UV/TiO2
Next Article in Special Issue
Distribution, Formation and Human Health Risk of Fluorine in Groundwater in Songnen Plain, NE China
Previous Article in Journal
Abiotic Aggregation of Organic Matter in Coastal and Estuarine Waters: Cases in the Eastern Long Island Sound, USA
Previous Article in Special Issue
Evaluation of the Impact of Ecological Water Supplement on Groundwater Restoration Based on Numerical Simulation: A Case Study in the Section of Yongding River, Beijing Plain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on the Application of Typical Biological Chain for Algal Control in Lake Ecological Restoration—A Case Study of Lianshi Lake in Yongding River

1
School of Municipal and Environmental Engineering, Hebei University of Architecture, Zhangjiakou 075000, China
2
Department of Water Ecological Environment, China Institute of Water Resources and Hydropower Research, Beijing 100038, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(21), 3079; https://doi.org/10.3390/w13213079
Submission received: 4 September 2021 / Revised: 23 October 2021 / Accepted: 24 October 2021 / Published: 2 November 2021
(This article belongs to the Special Issue River Ecological Restoration and Groundwater Artificial Recharge)

Abstract

:
Maintaining the health of lake ecosystems is an urgent issue. However, eutrophication seriously affects lakes’ ecological functions. Eutrophication is also the main target of lake ecological restoration. It is vital to carry out research on lake eutrophication control and energy flow evaluation in ecosystems scientifically. Based on in situ survey results for the aquatic life data for Lianshi Lake from 2018 to 2019, the Ecopath model was used to establish an evaluation index system for the typical biological chain to screen out the key species in the water ecosystem, and the fuzzy comprehensive evaluation (FCE) method was used to screen all the biological chains controlling algae. A combination of the FCE coupled with the Ecopath screening method for typical biological chains for algal control was applied to the Lianshi Lake area; the results show that the typical biological chain for algal control is phytoplankton (Phyt)–zooplankton (Zoop)–macrocrustaceans (Macc)–other piscivorous (OthP). Upon adjusting the biomass of Zoop and Macc in the typical biological chain for algal control to three times that of the current status, the ecological nutrition efficiency of Phyt was increased from 0.308 to 0.906. The material flow into the second trophic level from primary producers increased from 3043 to 8283 t/km2/year. The amount of detritus flowing into primary producers for sedimentation decreased from 7618 to 2378 t/km2/year. Finally, the total primary production/total respiratory volume (TPP/TR) decreased from 9.224 to 3.403, the Finn’s cycle index (FCI) increased from 13.6% to 17.5%, and the Finn’s average energy flow path length (FCL) increased from 2.854 to 3.410. The results suggest that the problem of eutrophication can be solved by introducing Zoop (an algal predator) and Macc to a large extent, resulting in improved ecosystem maturity. The research results can facilitate decision making for the restoration of urban lake water ecosystems.

1. Introduction

The rapid expansion of the world’s population has exacerbated the degradation of global lake ecosystems [1]. According to statistics, more than 60% of the lakes in the world are in different degrees of eutrophication [2]. The increasing eutrophication of lakes has become a global water environment problem [3] (e.g., for Lake Canyon in the United States [4], Lake Geneva in Switzerland [5], and Sugarloaf Lake in Australia [6]).
The interaction between lake water environments and water ecology is complex. It is important to ensure the integrity of the ecosystem while improving the quality of the water environment [7,8,9]. Therefore, how to deal with lake eutrophication and restore aquatic ecosystems has become an urgent problem to be solved in current limnology, environmental science, freshwater ecology and other disciplines, and has been closely followed by international scholars [10,11].
At present, the recognized theories of lake ecological restoration mainly include the nutrient salt concentration limit theory by Scheffer [12], multi-steady state theory by Lewontin [13] and biological manipulation theory by Hrbacek [14]. Among them, the biological manipulation theory has the most promising application prospects and has led to many successful cases of lake ecological restoration. For example, the biological manipulation process implemented by putting silver carp (Silc) and bighead carp (Bigc) in Donghu Lake in Wuhan is one of the main reasons for the disappearance of cyanobacteria blooms in the lake [15]. Olin et al. [16] reported that the removal of common carp (Comc) from 10 lakes in southern Finland effectively reduced the biomass of cyanobacteria and the degree of algal outbreaks. Shapiro et al. [17] adjusted the ratio of other piscivarious (Othp) and plankton-eating fish from 1:1.65 to 1:2.2 in the Round Lake, and the concentrations of total nitrogen (TN), total phosphorus (TP) and chlorophyl-a (Chl-a) in the lake were reduced to varying degrees. However, biological manipulation cannot enable all ecosystems to achieve the expected ecological functions and may lead to changes in the species diversity of the ecosystem, the decline of the average trophic level in the system, and the destruction of habitats [18,19]. For example, Razlutskij et al. [20], through 72 days of outdoor experiments, found that introducing Carassius auratus increased the biomass of Phyt. Recent ecological studies have revealed significant negative effects of crucian carp (Cruc) on the water quality and ecological states of shallow lakes, e.g., increasing nutrient levels, leading to reduced water clarity [21,22]. The current international guidance on the application of biological manipulation technology to the process of lake ecological restoration remains far from sufficient [23].
In light of this, this study intended to establish a screening method for typical biological chains of ecosystems and propose a biomass control strategy for the typical biological chain of alga control. The purpose of this study was to provide strong theoretical support for the management of urban lake eutrophication and water ecological restoration, and it consisted of the following: (1) Key species and target biological chains were screened, based on the Ecopath model; the key species of the ecosystem were analyzed, and the biological chain with phytoplankton (Phyt) as the primary producer and including the key species was selected. (2) Typical biological chains were screened, and an evaluation index that could reflect the efficiency and stability of the biological chain was constructed. The index weight was used to select typical biological chains based on the fuzzy comprehensive evaluation method. (3) The regulatory impact was analyzed, the biomass of key species in the typical biological chain was regulated, and then the impacts on the regulation target and the ecosystem were judged.

2. Materials and Methods

2.1. Overview of the Study Area and Data Sources

2.1.1. Overview of the Study Area

As a river-type lake in the urban section of the plain of the Yongding River, Lianshi Lake is located at the junction of Mentougou and Fengtai District in the southwest of Shijingshan District (116°6′58″~116°9′34″ E, 39°56′3″~39°53′15″ N), which belongs to the temperate continental monsoon climate, with high temperatures and rain in summer and cold and dry conditions in winter. The total length of Lianshi Lake is 5.8 km, the average width of the lake is 376 m, the widest point is about 500 m, the total water area is 106 hm2, and the average depth of the water body is 1.6 m. The phytoplankton in Lianshi Lake are dominated by Pseudanabaena moniliformis (Cyanophinyta) and Limnothrix. Lianshi Lake belongs to the northern freshwater city lake. The distribution of the study area and monitoring points are shown in Figure 1.

2.1.2. Data Sources

The data sources included long-term field sampling survey data, formula calculations and reference comparisons. The main source of the biomass B i of each function group was a field survey of the biomass at six points in Lianshi Lake from 2018 to 2019. The macroinvertebrates were collected using a Surber sampler with a mesh diameter of 40 mesh (500 µm) and a sampling area of 0.09 m2, and 70% alcohol was added to the test tube bottle to be stored for inspection. The zooplankton were collected using a No. 25 plankton net (200 mesh) and immediately fixed with 4% formaldehyde. For the collection of phytoplankton, a 1 L water sample (0.5 m below the surface of the water body) was collected in a plexiglass water collector, and 15 mL of Lugol’s reagent was added for fixation. Sampling areas of 0.5 m × 0.5 m and 1 m × 1 m were used for vascular plants and submerged plants, respectively. In the field, the plant species were distinguished and weighed (wet weight). Fish resources were collected mainly using trawl nets, and the travel distance for each sampling point was no more than 100 m. The biomass of organic detritus was calculated based on the relationship between the primary productivity and water transparency [24] (Equation (1)).
l o g D = 0.954 l o g P P + 0.863 l o g E 2.41
where D is the detrital biomass (g C·m−2), PP is the primary production volume (g C·m−2·a−1), and E is the average transparency (m).

2.2. Ecopath Model

2.2.1. Principle of the Ecopath Model

The Ecopath is an ecological model that can directly determine the structure of an ecological system and describe its energy flow and mass transfer through the principle of nutrition dynamics. The Ecopath model stipulates that each function group ( i ) energy input and output in the ecosystem are balanced. The model uses a set of simultaneous linear equations to define an ecosystem, and each function group is represented by a linear equation [25] (Equation (2)).
B i × ( P B ) i × E E i = j = 1 n B j × ( Q B ) i × D C i j + E X i
In the formula, B i is the biomass of the i -th function group. ( P B ) i is the ratio of the annual average production of the i -th function group to the annual average biomass, which is the biomass turnover rate. E E i is the ecological nutrient conversion efficiency of the i -th function group and is usually obtained by the calculation of the model [26]. ( Q B ) i is the ratio of the consumption of the i -th function group to the biomass. D C i j is the ratio of the i-th the prey group to the total predation of the predator group j . E X i is the output of the i -th function group.

2.2.2. Division of Function Groups

The purpose of the establishment of function groups is to merge populations with highly overlapping niches to simplify the food web. Based on the survey results for the Lianshi Lake ecosystem, characteristics of the community, and survival habits of each species, this work identified organisms with similar ecological characteristics, which were grouped into function groups [24]. A total of 14 function groups were set up (Table 1).

2.2.3. Parameter Setting

Some of the parameters refer to lakes with water environmental conditions similar to those of Lianshi Lake, such as the production/biomass (P/B) and consumption/biomass (Q/B) parameters of fish resources that were obtained by querying (Available online: http://www.fishbase.org) (accessed on 20 May 2021). The plant P/B coefficient refers to related research on Taihu Lake [27]. The P/B coefficients of zooplankton (Zoop), other benthos (Othb) and macrocrustaceans (Macc) were estimated based on the measured data by referring to the research results for Taihu Lake [27], Qiandao Lake [28], Dianshan Lake [29] and Zhushan Bay [30]. The P/Q coefficient refers to recognized data. The values for Zoop, OthB and Macc were 0.05 [31], 0.02 [32] and 0.075 [33], respectively. The food composition ( D C i ) is shown in Table 2. In addition to ecological investigation and research, we also referred to related research results [27,34], such as those for Zhushan Lake [30] and Qiandao Lake [35,36]. According to the cited literature, the G S values for general OthP and herbivorous fish (HerF) were set to 0.2 and 0.41, and the G S values for Zoop, OthB and Macc were set as 0.65, 0.94 and 0.7, respectively [31,32,33,37].

2.2.4. Model Balance Calculation

The Ecopath model requires the input of six basic parameters: B i ,   ( P B ) i , E E i , ( Q B ) i , D C i j and E X i . Model balancing was executed, which was established on the basis of conforming to objective laws. The model parameters could be slightly modified to meet the requirements of the model operation, but it was necessary to avoid changing reliable data sources [25]. The nutritional conversion efficiency ( E E i ) ranged from 0 to 1 [38], the group E E was close to 1 in the face of considerable predation pressure, and the underutilized function group had a lower E E value. The value range of G E ( P Q ) was generally 0.1–0.3. If the balance test had one or more E E > 1, it was necessary to locate which predators caused the problem for specific prey groups in the predation mortality. The level of model confidence is mainly related to the quality and reliability of the acquired data. The accuracy of the Ecopath model was measured using the Pedigree index. The higher the index, the higher the quality of the model. The input and output parameters of the balanced model are shown in Table 3.

2.3. Fuzzy Comprehensive Evaluation Method

The fuzzy comprehensive evaluation method is based on fuzzy mathematics. There is a certain characteristic of n things to be evaluated; these n things comprise the object set X = { x 1 , x 2 , , x n } , the factor set U = { u 1 , u 2 , , u n } and the evaluation set V = { v 1 , v 2 , v m } . Suppose that the weight distribution of the factors is the fuzzy subset A on V, denoted as A = { a 1 , a 2 , , a n } . In the formula, a i is the weight corresponding to the i -th factor u i , and the general rule is i = 1 n a i = 1 .
The evaluation steps for the fuzzy comprehensive evaluation method are as follows:
(1)
Establish a factor set.
Assuming there are a total of i factors of the judged object, the factor universe of the evaluated object U is U = { u 1 , u 2 , u i } .
(2)
Determine the membership function.
Assuming that the evaluation level is divided into i levels, the set is V = { v 1 , v 2 , v 3 , v i } .
(3)
Establish a fuzzy relationship evaluation matrix.
(4)
Establish a weight vector.
Determine the weight vector of the judgement factor A = ( a 1 , a 2 , , a n ) . A is the subordination relationship of each factor in U to the thing being judged; it depends on the focus of people when making fuzzy comprehensive judgments, and it is equivalent to assigning weights according to the importance of each factor in the judgment.
(5)
Establish a fuzzy comprehensive evaluation matrix.
According to the calculated maximum membership value, the program selection or category evaluation is performed.

3. Results

3.1. Screening of Typical Biological Chains of Fuzzy Comprehensive Evaluation Coupled with Ecopath

The Ecopath model was used to screen the key species in the Lianshi Lake ecosystem, identify the types of algae-controlling biological chains containing the key species, and establish a biological chain evaluation index system that could simultaneously express the delivery efficiency and stability. Finally, the fuzzy comprehensive evaluation method was used to screen out the biological chain with the largest weight value.
See Flowchart 2 for specific screening ideas (Figure 2).
(1) Analysis of key species in the ecosystem
The key function group plays an important role in the function of the Lianshi Lake ecosystem. This study used the method developed by Libralato et al. [39] to calculate the keystoneness (KS) index of each function group. The one with the largest KS value is regarded as the key function group. Compared to the calculation method for the KS value proposed by Valls et al. [40] and Power et al. [41], this calculation is simpler and comprehensively considers the effects of biomass and capacity. The equation is K S i = l o g [ ε i ( 1 p i ) ] ,   ε i = j 1 n m i j 2 ε, p i = B i i n B K , where K S I is the criticality index of function group i , ε i is the total impacts of function group i in the ecosystem (total impacts), and m i j is the value of the mixed nutrition effect of function group i on function group j , indicating the mutual relationship between each other’s strengths, and P i is the ratio of the biomass B i of function group i to the biomass of the entire ecosystem k n B k . Since K S i and p i are negatively correlated, the criticality index will not be too high due to the high biomass of the function groups (Table 4).
As shown in Figure 3, the first key function groups in the Lianshi Lake ecosystem are submerged macrophytes (SubM). However, other macrophytes (OthM) are redundant in the food web of the entire ecosystem. Other benthos (OthB), Phyt, other fish (OthF) and Zoop have a relative total impact second only to SubM, and their criticality index ranks in the top four.
(2) Screening of target biological chains for alga control and salt reduction
With the goal of alga control, we identified the biological chain with Phyt and SubM as primary producers, and then determined the biological chain containing key species such as Zoop, OthB and OthF. The results are shown in Figure 4.
(3) Establishment of evaluation indicators for typical biological chains
In order to express the characteristics of the Lianshi Lake ecosystem, the ecotrophic efficiency ( E E ), production/respiration ( P / R ), relative total impact ( R T I ), omnivory index ( O I ) and trophic level ( T L ) are listed as typical indicators for biological chain evaluation (Table 5). The evaluation indicators have the following characteristics: (1) E E stands for the utilization and conversion efficiency for the energy contained in the previous trophic organism. The value of EE ranges from 0 to 1, and it is close to 1 in a group with considerable predation pressure. (2) P / R represents an important indicator of the maturity of the biological chain, which is close to 1 in a mature ecosystem. (3) R T I represents the relative total impact of a single function group on the entire ecosystem. (4) O I indicates that, when the omnivorous index is zero, the corresponding prey is oriented. (5) T L refers to the level of the organism in the food chain of the ecosystem. The higher the trophic level of the organism, the greater the contribution of the food chain to the stability of the ecosystem.
(4) Fuzzy comprehensive evaluation method for screening typical biological chains
The research did not consider the comment-level domain in the fuzzy comprehensive evaluation, and selected the typical biological chain according to the maximum weight (from Equation (3) to Equation (6)):
① The factor domain U , U = ( E E , P / R , R T I , O I , T L ) of the judged object was determined.
② The membership function was established.
The higher the value, the higher the membership function:
r i = { 0                                                           ( x S 1 ) x S 1 S 2 S 1                                               ( S 1 < x < S 2 ) 1                                                             ( x S 2 )
(I) The determination of the membership functions of E E ,   P / R ,   R T I , and   O I As the variation range of E E ,   P / R ,   R T I ,   and O I was between 0 and 1, the membership function was:
U A ( x ) = x
(II) The determination of the TL membership function
According to the literature, the lowest trophic level of primary producers in lake ecosystems is 1, and the trophic level of the function group that is generally at the top is 4, so the membership function is:
U A ( x ) = x 1 3
③ Single factor evaluation was performed, and a fuzzy relationship matrix was established.
The fuzzy evaluation matrix was established according to the following formula:
R = [ r i j ] = [ r 11 r 15 r 51 r 55 ]
In the formula, in the i -th row R i = ( r i 1 , r i 2 , r i m ) , i = 1 , m , i is the degree of the membership of the i -th evaluation factor to the evaluation standards at all levels; for the j -th column R j = ( r 1 j , r 2 j , , r n j ) , j is the degree of the membership of each evaluation factor to the j -th evaluation standard. Taking the nine biological chains containing the Zoop group as an example, the membership degrees corresponding to different evaluation indicators of different biological chains were calculated and are shown in Table 6 below.
The fuzzy relation matrix with the Zoop group:
R   = [ 0.257 0.338 0.260 0.251 0.330 0.297 0.457 0.457 0.457 0.328 0.388 0.396 0.341 0.272 0.271 0.332 0.188 0.193 0.631 0.790 0.759 0.624 0.584 0.572 0.765 0.521 0.533 0.058 0.046 0.062 0.137 0.122 0.076 0.067 0.090 0.132 0.418 0.323 0.360 0.389 0.380 0.439 0.319 0.247 0.279 ]
④ A weight vector was established.
The weight vector A = ( a 1 , a 2 , , a n ) of the judgement factor was determined. A is the subordination relationship of each factor in U to the thing being judged. Due to the importance of the distribution of the weight, this study used the analytic hierarchy process to determine the weight of each evaluation index; the hierarchical structure is shown in Figure 5.
The judgement matrix T   =   [ 1 4 0.33 0.143 0.25 0.25 1 0.25 0.143 0.2 3 4 1 0.5 3 7 7 2 1 8 4 5 0.33 0.125 1 ] . The calculated maximum characteristic value is λ m a x = 5.42 , the random consensus ratio is 0.0944 < 0.1, and the consistency test shows that the weight distribution is reasonable, Therefore, the weight of each evaluation index of a typical biological chain is A = ( 0.0712 , 0.0388 , 0.2224 , 0.0534 , 0.1336 ) .
⑤ The comprehensive evaluation results were analyzed.
B = A R = ( 0.2582 , 0.2825 , 0.2839 , 0.2950 , 0.2798 , 0.2581 , 0.2940 , 0.2368 , 0.2663 )
It can be seen that the typical biological chain containing the Zoop group as the key species is ④, equal to phytoplankton (Phyt)–zooplankton (Zoop)–macrocrustaceans (Macc)–other piscivorous (OthP).
According to the same calculation principle, the typical algal control biological chain containing OthB and OthF as the key species does not exist.

3.2. Analysis of Typical Algal Control Biological Chain and Nutritional Structure and Function Response

Through modeling, it was found that the production of primary producers in the Lianshi Lake ecosystem was 10,661.7 t/km2/year, of which the amount of food consumed was 3043 t/km2/year, and the remaining 71.5% was not consumed by other predators but flowed into the debris, entered recirculation, and became deposited through mineralization. Therefore, the effective removal of plant organisms from the water body plays an important role in controlling the total amount of nutrients, dissolved oxygen and transparency of the water body.
In the Lianshi Lake ecosystem, there are 16 biological chains with Phyt as the primary producer. Among them, the number of biological chains where Phyt is directly grazed by Zoop is the largest; it has nine chains. The transfer efficiency of Phyt is 0.308. The transfer efficiency of Zoop reached 0.606. It can be seen that the types of biological chains that impose significant biological density constraints on algae are mainly Zoop predation on Phyt. Combined with the screening of typical biological chains, a typical biological chain of algal control in the Lianshi Lake ecosystem is phytoplankton (Phyt)–zooplankton (Zoop)–macrocrustaceans (Macc)–other piscivorous (OthP).
This study improved the efficiency of nutrient delivery for primary producers by simulating an increase in biomass of Zoop and Macc. Under the guidance of biological manipulation theory, the key factor in the regulation of the algal control food chain is Zoop. Through the introduction of Daphnia magna, the food chain is opened up through insects eating algae and fish-eating insects. A symbiotic system of “Daphnia magna–underwater forest–aquatic animals–microbes” was constructed, and the “grass-type clean water state” self-purification system was restored.
The proliferation of ecological capacity continuously increases the biomass of the target species. By observing the changes in other function groups such as the eaten organisms in the model, when the ecological nutrition transfer efficiency of a function group in the model E E > 1, the model will become unbalanced; the biomass value of the released species before the imbalance of the model is the ecological capacity. The biomass of Zoop and Macc have been increased to 1.5, 2, 3 and 4 times, respectively. It can be seen from Table 7 that, within the ecological capacity of Zoop and Macc, with increases in the biomasses of both, the ecological nutrition transfer efficiency of Phyt gradually increased, from 0.308 (the current situation) to 0.458 (1.5 times), 0.607 (2 times) and 0.906 (3 times). The flow of material flowing into the second trophic level from primary producers also increased from 3043 (the current situation) to 4353 (1.5 times), 5663 (2 times) and 8283 t/km2/year (3 times).
The amount of debris flowing into the primary producers for deposition continuously decreased, from 7618 (the current situation) to 6308 (1.5 times), 4998 (two times) and 2378 t/km2/year (three times). The total primary production/total respiratory volume (TPP/TR) continued to decrease from 9.224 (the current situation) to 6.461 (1.5 times), 4.972 (two times) and 3.403 (three times). The Finn’s circulation index (FCI) continued to rise from 13.59% (the current situation) to 14.33% (1.5 times), 15.23% (two times) and 17.51% (three times), whereas the Finn’s average energy flow path length (FCL) rose from 2.854 (the current situation) to 2.993 (1.5 times), 3.132 (two times) and 3.410 (three times).
Based on the analysis of the nutritional structure of the Lianshi Lake ecosystem, the artificial introduction of Zoop (Daphnia magna) and Macc can increase the transfer efficiency and the maturity of the ecosystem. To a certain extent, it can solve the problem of excessive primary production.

4. Discussion

4.1. Development Characteristics of Lianshi Lake Ecosystem

The total primary production/total respiratory volume (TPP/TR) is an important evaluation index, which is close to 1 in mature ecosystems, far greater than 1 in developing ecosystems, and less than 1 in polluted ecosystems. The Finn’s cycling index (FCI) is the ratio of the circulation flow to the total flow in the system, and the Finn’s mean path length (FMPL) is the average length of each circulation through the food chain. The higher the ratio of material recycling, the longer the food chain through which the nutrient flows, and the FCI value of a mature ecosystem is close to 1. The current Lianshi Lake TPP/TR is 9.224, the FCI is 13.59%, and the FMPL is 2.854, indicating that the Lianshi Lake ecosystem is immature.
In the typical biological chain of algal control, with an increase in the biomass of Zoop (Daphnia magna) and Macc within the ecological capacity, the amount of phytoplankton as a primary producer flowing into the next trophic level gradually increases, the transfer efficiency gradually increases, the FCI and FMPL gradually increase, and the TPP/TR gradually approaches 1, which reduces the risk of lake eutrophication and increases ecosystem maturity to a certain extent.
Therefore, there may be two reasons for the low maturity of the Lianshi Lake ecosystem: first, the biomass of the function group that plays a key role in the Lianshi Lake ecosystem is much lower than the ecological capacity, resulting in insufficient driving force for the ecosystem to develop to a mature state. Second, the biodiversity of Lianshi Lake is low, and the flow of energy to higher levels is hindered. It is recommended to introduce Zoop (Daphnia magna) and indigenous herbivorous fishes to build a food chain in order to promote material and energy cycles.

4.2. Prospects of Ecopath Model in Lake Ecological Restoration

As early as 1975, Shapiro et al. [17] proposed the biological manipulation theory. Biomanipulation methods have been developed for nearly 50 years; there have been many reports in Western European countries, but this technology has not been widely promoted. This may be because biological manipulation methods involve complex biological networks, and too many factors are affected. Traditional research methods can only study the behavior of individual organisms in simple habitats or competitive environments, and there is very little research on the biological chain that significantly affects the regulation goals and ecosystem conditions. Therefore, it is generally believed that research on complex ecosystems must rely on the guidance of mathematical models or theories [42].
The Ecopath with Ecosim (EWE) model, also known as the ecological channel model, is an ecological model that can assess the true structure of the ecosystem and describe its energy flow and mass balance. The Ecopath model was originally created in 1984 [43]. After years of development, the Ecopath model has become a key tool for ecosystem research. Fuzzy comprehensive evaluation (FCE) is based on fuzzy mathematics [44], which can express the fuzzy relationship between various factors and solve the problem of ambiguity between multi-factor evaluations that cannot be solved by traditional methods, and it has been widely used in the field of policy evaluation and risk assessment. A typical biological chain of alga control combines the advantages of the Ecopath model and FCE method, and the results show that the biomass of phytoplankton in Lianshi Lake has been effectively controlled and that the ecosystem’s maturity has been improved.
Ecopath is a powerful model but is mostly used to assess the condition of the ecosystem and provide scientific management and control solutions. In the future, the Ecopath model will continue to be developed and coupled with other models, such as pollutant diffusion models and ecotoxicology models, which will have important scientific research significance for exploring the restoration of lake ecosystems.

5. Conclusions

Based on the survey results for aquatic organisms in Lianshi Lake from 2018 to 2019, this study established, for the first time, a typical biological chain screening method with fuzzy comprehensive evaluation coupled with Ecopath. Among the organisms, Phyt is the primary producer, and the typical biological chain of alga control with Zoop as the key species is phytoplankton (Phyt)–zooplankton (Zoop)–macrocrustaceans (Macc)–other piscivorous (OthP).
In a typical biological chain with significant biological density constraints on algae, when the biomass of Zoop was increased from 7.85 (the current situation) to 23.55 t/km2/year (three times) and Macc was increased from 1.58 (the current situation) to 4.74 t/km2/year (three times), the results show that the ecological nutrition efficiency of Phyt increased from 0.308 (the current situation) to 0.906 (three times), the material flow into the second trophic level from primary producers increased from 3043 (the current situation) to 8283 t/km2/year (three times), the amount of debris flowing into primary producers for sedimentation decreased from 7618 (the current situation) to 2378 t/km2/year (three times), the total primary production/total respiratory volume (TPP/TR) decreased from 9.224 (the current situation) to 3.403 (three times), the Finn’s cycle index (FCI) increased from 13.59% (the current situation) to 17.51% (three times), and the Finn’s average energy flow path length (FCL) increased from 2.854 (the current situation) to 3.410 (three times). In the typical biological chain of alga control, the artificial release of Zoop (Daphnia magna) and Macc can improve the transfer efficiency of phytoplankton as primary producers to a certain extent, reduce the harm caused by eutrophication to lake ecosystems, and improve the maturity of the lake ecosystem.

Author Contributions

P.Z. and X.C. designed the research content; P.Z., H.L. and W.P. processed the data and analyzed the results. All the authors participated in the writing of the manuscript. Y.G. and X.C. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Major Science and Technology Program for Water Pollution Control and Treatment (No.2018ZX07101005).

Acknowledgments

We thank Min Zhang for providing hydrobiological data and Su Wei for assistance with the fieldwork. We would like to thank Guanglei Qiu for his guidance on the translation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhai, Y.; Xia, X.; Yang, G.; Lu, H.; Ma, G.; Wang, G.; Teng, Y.; Yuan, W.; Shrestha, S. Trend, seasonality and relationships of aquatic environmental quality indicators and implications: An experience from Songhua River, NE China. Ecol. Eng. 2020, 145, 105706. [Google Scholar] [CrossRef]
  2. Wang, S.; Li, J.; Zhang, B.; Spyrakos, E.; Tyler, A.N.; Shen, Q.; Zhang, F.; Kuster, T.; Lehmann, M.K.; Wu, Y.; et al. Trophic state assessment of global inland waters using a MODIS-derived Forel-Ule index. Remote Sens. Environ. 2018, 217, 444–460. [Google Scholar] [CrossRef] [Green Version]
  3. Zhao, Y.; Deng, X.; Zhan, J. Research progress on prevention and control strategies of lake eutrophication. Environ. Sci. Technol. 2010, 33, 92–98. (In Chinese) [Google Scholar]
  4. Davis, K.; Anderson, M.A.; Yates, M.V. Distribution of indicator bacteria in Canyon Lake, California. Water Res. 2005, 39, 1277–1288. [Google Scholar] [CrossRef]
  5. Razmi, A.M.; Barry, D.A.; Lemmin, U.; Bonvin, F.; Kohn, T.; Bakhtyar, R. Direct effects of dominant winds on residence and travel times in the wide and open lacustrine embayment: Vidy Bay (Lake Geneva, Switzerland). Aquat. Sci. 2014, 76, 59–71. [Google Scholar] [CrossRef] [Green Version]
  6. Hipsey, M.R.; Antenucci, J.P.; Brookes, J. A generic, process-based model of microbial pollution in aquatic systems. Water Resour. Res. 2008, 44, 76–87. [Google Scholar] [CrossRef]
  7. Su, W.; Wu, J.; Zhu, B.; Chen, K.; Peng, W.; Hu, B. Health Evaluation and Risk Factor Identification of Urban Lakes—A Case Study of Lianshi Lake. Water 2020, 12, 1428. [Google Scholar] [CrossRef]
  8. Shan, B.; Wang, C.; Li, X. The formulation method and case study of river governance plan based on water quality objective management. Acta Sci. Circumstantiae 2015, 35, 2314–2323. (In Chinese) [Google Scholar] [CrossRef]
  9. Birk, S.; Bonne, W.; Borja, A.; Brucet, S.; Courrat, A.; Poikane, S.; Solimini, A.; van de Bund, W.; Zampoukas, N.; Hering, D. Three hundred ways to assess Europe’s surface waters: An almost complete overview of biological methods to implement the Water Framework Directive. Ecol. Indic. 2012, 18, 31–41. [Google Scholar] [CrossRef]
  10. Jiang, Y.; Zhang, J.; Zhu, Y.; Du, Q.; Teng, Y.; Zhai, Y. Design and Optimization of a Fully-Penetrating Riverbank Filtration Well Scheme at a Fully-Penetrating River Based on Analytical Methods. Water 2019, 11, 418. [Google Scholar] [CrossRef] [Green Version]
  11. Li, Y.; Zhang, Q.; Yao, J. Investigation of Residence and Travel Times in a Large Floodplain Lake with Complex Lake-River Interactions: Poyang Lake (China). Water 2015, 7, 1991–2012. [Google Scholar] [CrossRef] [Green Version]
  12. Scheffer, M.; Straile, D.; van Nes, E.; Hosper, S.H. Climatic warming causes regime shifts in lake food webs. Limnol. Oceanogr. 2001, 46, 1780–1783. [Google Scholar] [CrossRef]
  13. Lewonin, R.C. The meaning of stability: Diversity and stability in ecological systems. Brookhaven Symp. Biol. 1969, 22, 13–24. [Google Scholar]
  14. Hrbáčke, J.; Dvořakova, M.; Kořínek, V.; Procházkóva, L. Demonstration of the effect of the fish stock on the species composition of zooplankton and the intensity of metabolism of the whole plankton association. Verh. Int. Ver. Limnol. 1961, 14, 192–195. [Google Scholar] [CrossRef]
  15. Liu, J.; Xie, P. Uncover the mystery of the disappearance of cyanobacteria blooms in East Lake, Wuhan. Resour. Environ. Yangtze Basin 1999, 3, 85–92. (In Chinese) [Google Scholar]
  16. Olin, M.; Rask, M.; Ruuhijärvi, J.; Keskitalo, J.; Horppila, J.; Tallberg, P.; Taponen, T.; Lehtovaara, A.; Sammalkorpi, I. Effects of Biomanipulation on Fish and Plankton Communities in Ten Eutrophic Lakes of Southern Finland. Hydrobiologia 2006, 553, 67–88. [Google Scholar] [CrossRef]
  17. With, J.S.; Wright, D.I. Lake restoration by biomanipulation: Round Lake, Minnesota, the first two years. Freshw. Biol. 1984, 14, 371–383. [Google Scholar] [CrossRef]
  18. He, H.; Han, Y.; Li, Q.; Jeppesen, E.; Li, K.; Yu, J.; Liu, Z. Crucian Carp (Carassius carassius) Strongly Affect C/N/P Stoichiometry of Suspended Particulate Matter in Shallow Warm Water Eutrophic Lakes. Water 2019, 11, 524. [Google Scholar] [CrossRef] [Green Version]
  19. Rakshit, N.; Banerjee, A.; Mukherjee, J.; Chakrabarty, M.; Borrett, S.R.; Ray, S. Comparative study of food webs from two different time periods of Hooghly Matla estuarine system, India through network analysis. Ecol. Model. 2017, 356, 25–37. [Google Scholar] [CrossRef]
  20. Vladimir, R.; Xueying, M.; Natallia, M.; Elena, S.; Dzmitry, L.; Andrei, M.; Erik, J.; Xiufeng, Z. Omnivorous Carp (Carassius gibelio) Increase Eutrophication in Part by Preventing Development of Large-Bodied Zooplankton and Submerged Macrophytes. Water 2021, 13, 1497. [Google Scholar]
  21. Gu, J.; He, H.; Jin, H.; Yu, J.; Jeppesen, E.; Nairn, R.W.; Li, K. Synergistic negative effects of small-sized benthivorous fish and nitrogen loading on the growth of submerged macrophytes—Relevance for shallow lake restoration. Sci. Total Environ. 2018, 610–611, 1572–1580. [Google Scholar] [CrossRef]
  22. He, H.; Hu, E.; Yu, J.; Luo, X.; Li, K.; Jeppesen, E.; Liu, Z. Does turbidity induced by Carassius carassius limit phytoplankton growth? A mesocosm study. Environ. Sci. Pollut. Res. 2016, 24, 5012–5018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Li, C.; Xian, Y.; Ye, C.; Wang, Y.; Wei, W.; Xi, H.; Zheng, B. Wetland ecosystem status and restoration using the Ecopath with Ecosim (EWE) model. Sci. Total Environ. 2019, 658, 305–314. [Google Scholar] [CrossRef] [PubMed]
  24. Christensen, V.; Pauly, D. ECOPATH II—A software for balancing steady-state ecosystem models and calculating network characteristics. Ecol. Model. 1992, 61, 169–185. [Google Scholar] [CrossRef]
  25. Christensen, V.; Walters, C.J.; Pauly, D. Ecopath with Ecosim: A User’s Guide; Fisheries Centre, University of British Columbia: Penang, Malaysia, 2005; Available online: https://www.researchgate.net/publication/267193103_Ecopath_with_Ecosim_A_User’s_Guide (accessed on 15 April 2021).
  26. Huang, X.; Bing, X.; Zhang, X. Application of EwE model in evaluating fishery water ecosystem. J. Hydroecology 2011, 32, 125–129. (In Chinese) [Google Scholar] [CrossRef]
  27. Song, B. Research on the Ecosystem Model of Taihu Lake Fishery and Environment. Ph.D. Thesis, East China Normal University, Shanghai, China. Available online: https://kns.cnki.net/kcms/detail/detail.aspx?FileName=2004087212.nh&DbName=CDFD2004 (accessed on 2 June 2021). (In Chinese).
  28. Yu, J.; Liu, J.; Wang, L. Analysis of Structure and Function of Qiandao Lake Ecosystem Based on Ecopath Model. Acta Hydrobiol. Sin. 2021, 45, 308–317. (In Chinese) [Google Scholar]
  29. Feng, D. Research on ECOPATH Model of Dianshan Lake Ecosystem Structure and Energy Flow Characteristics. Master’s Thesis, East China Normal University, Shanghai, China. Available online: https://kns.cnki.net/kcms/detail/detail.aspx?FileName=1011131168.nh&DbName=CMFD2011 (accessed on 18 June 2021). (In Chinese).
  30. Xian, Y.; Ye, C.; Li, C. Construction and Analysis of EWE Model of Lake Buffer Zone Wetland Ecosystem in Zhushan Bay. Chin. J. Appl. Ecol. 2016, 27, 2101–2110. (In Chinese) [Google Scholar] [CrossRef]
  31. Park, R.A. A generalized model for simulating lake ecosystems. Simulation 1974, 23, 33–50. [Google Scholar] [CrossRef]
  32. Yan, Y.; Liang, Y. Energy flow of benthic fauna in Biandan Pond. Acta Ecol. Sin. 2003, 3, 527–538. (In Chinese) [Google Scholar]
  33. Halfon, E.; Schito, N.; Ulanowicz, R.E. Energy flow through the Lake Ontario food web: Conceptual model and an attempt at mass balance. Ecol. Model. 1996, 86, 1–36. [Google Scholar] [CrossRef]
  34. Li, Y.; Liu, E.; Wang, H. Analysis of the Structure and Function of Taihu Lake Ecosystem Based on Ecopath Model. Chin. J. Appl. Ecol. 2014, 25, 2033–2040. (In Chinese) [Google Scholar] [CrossRef]
  35. Song, G.; Hu, M.; Liu, Q. Using Stable Isotope Technology to Study the Feeding Habits and Nutritional Levels of Gillnet Catch in Qiandao Lake in Autumn. J. Shanghai Ocean. Univ. 2014, 23, 117–122. (In Chinese) [Google Scholar]
  36. Liu, Q. Qiandao Lake Water Conservation Fishery and Its Impact on Lake Ecosystem. Ph.D. Thesis, East China Normal University, Shanghai, China. Available online: https://kns.cnki.net/kcms/detail/detail.aspx?FileName=2005086209.nh&DbName=CDFD2005 (accessed on 30 May 2021). (In Chinese).
  37. Scavia, D. Documentation of CLEANX a generalized model for simulating the open-water ecosystems of lakes. Simulation 1974, 2, 51–56. [Google Scholar] [CrossRef]
  38. Christensen, V.; Walters, C.; Ahrens, R. Databasedriven models of the world’s Large Marine Ecosystems. Ecol. Model. 2009, 9, 1984–1996. [Google Scholar] [CrossRef]
  39. Libralato, S.; Christensen, V.; Pauly, D. A method for identifying keystone species in food web models. Ecol. Model. 2006, 195, 153–171. [Google Scholar] [CrossRef]
  40. Valls, A.; Coll, M.; Christensen, V. Keystone species: Toward an operational concept for marine biodiversity conservation. Ecol. Monogr. 2015, 85, 29–47. [Google Scholar] [CrossRef] [Green Version]
  41. Power, M.E.; Tilman, D.; Estes, J.A.; Menge, B.A.; Bond, W.J. Challenges in the Quest for Keystones. BioScience 1996, 46, 609–620. [Google Scholar] [CrossRef]
  42. Jiang, X.H.; Zhang, D.Y. Discussion on the functions of model in ecology. Resour. Ecol. Environ. Net. Res. 1992, 2, 30–34. [Google Scholar]
  43. Polovina, J.J. Model of a coral reef ecosystem. Coral Reefs. 1984, 3, 1–11. [Google Scholar] [CrossRef]
  44. Chen, T.; Shen, D.; Jin, Y.; Li, H.; Yu, Z.; Feng, H.; Long, Y.; Yin, J. Comprehensive evaluation of environ-economic benefits of anaerobic digestion technology in an integrated food waste-based methane plant using a fuzzy mathematical model. Appl. Energy 2017, 208, 666–677. [Google Scholar] [CrossRef]
Figure 1. The location of Lianshi Lake and the distribution of monitoring points.
Figure 1. The location of Lianshi Lake and the distribution of monitoring points.
Water 13 03079 g001
Figure 2. Screening flowchart of a typical biological chain.
Figure 2. Screening flowchart of a typical biological chain.
Water 13 03079 g002
Figure 3. Key species in the ecosystem of Lianshi Lake.
Figure 3. Key species in the ecosystem of Lianshi Lake.
Water 13 03079 g003
Figure 4. The biological chain of algal control in Lianshi Lake containing key species.
Figure 4. The biological chain of algal control in Lianshi Lake containing key species.
Water 13 03079 g004aWater 13 03079 g004b
Figure 5. Hierarchical structure model.
Figure 5. Hierarchical structure model.
Water 13 03079 g005
Table 1. Ecosystem function groups and main types for Lianshi Lake.
Table 1. Ecosystem function groups and main types for Lianshi Lake.
No.Function
Group
Abbreviation for CompositionIncluded Types
1Other piscivorousOthPHorsemouth, yellow catfish
2Common carpComcCommon carp
3Crucian carpCrucCrucian carp
4Bighead carpBigcBighead carp
5Silver carpSilcSilver carp
6Herbivorous fishHerFGrass carp, bream
7Other fishOthFWheat ear fish, tortoisefish
8MacrocrustaceansMaccGreen prawns, prawns, Chinese mitten crabs, etc.
9Other benthosOthBHydrophilia, Ceratobranchus, Longbrachium, Fanchus, Chironomus, Chironomidae, etc.
10ZooplanktonZoopProtozoa, rotifers, cladocerans, copepods, etc.
11PhytoplanktonPhytCyanobacteria, green algae, euglena, dinoflagellate, cryptophytes, golden algae, etc.
12Submerged macrophytesSubMPotamogeton, Myriophyllum, Hydrilla verticillata, Ceratophyllum, etc.
13Other macrophytesOthMReeds, cattails, etc.
14DetritusDetrOrganic detritus
Table 2. Food composition data entered into the model.
Table 2. Food composition data entered into the model.
No.Prey/Predator12345678910
1OthP0.007
2Comc0.15
3Cruc0.27
4Bigc
5Silc
6HerF0.07 0.006
7OthF
8Macc0.38
9OthB0.0730.13 0.230
10Zoop 0.820.240.5010.213 0.6200.3500.0050.009
11Phyt 0.048 0.3610.620 0.3000.0220.801
12SubM 0.0020.09 0.9970.003 0.101
13OthM0.01 0.003
14Detr0.04 0.670.1380.167 0.1410.3500.8720.190
15Sum1111111111
Table 3. Parameters of Lianshi Lake ecosystem construction.
Table 3. Parameters of Lianshi Lake ecosystem construction.
Function GroupBiomass (t/km2)Production/
Biomass
Consumption/BiomassEco-Nutrition EfficiencyProduction/
Consumption
Proportion of Unassimilated Food
OthP0.131.6706.10.0260.2740.200
Comc0.50.96010.70.2480.0900.200
Cruc0.51.13012.30.3790.0920.200
Bigc1.80.9906.90.0010.1430.200
Silc1.21.1008.00.0010.1380.200
HerF0.270.9877.10.7780.1390.410
OthF2.32.15511.00.0010.1960.410
Macc1.583.09041.00.0620.0750.700
OthB16.1414.130206.50.0990.0200.940
Zoop7.8520.680413.60.6060.0500.650
Phyt47.42185.000 0.308
SubM14601.250 0.186
OthM641 0.001
Detr3.230 0.272
Table 4. Criticality data table for each function group of Lianshi Lake.
Table 4. Criticality data table for each function group of Lianshi Lake.
No.Function GroupBiomass P i   ε i K S i Relative Total Impact
1OthP0.130.00008090.899516−0.046030.622
2Comc0.500.0003110.147707−0.830730.122
3Cruc0.500.0003110.197666−0.70420.17
4Bigc1.800.001120.054332−1.265430.0531
5Silc1.200.0007470.017978−1.745590.0175
6HerF0.270.0001680.429709−0.36690.425
7OthF2.300.0014310.845322−0.07360.749
8Macc1.580.0009830.361432−0.44240.328
9OthB16.140.0100451.0777140.0281190.867
10Zoop7.850.0048850.882764−0.056280.701
11Phyt47.420.029510.914433−0.051860.844
12SubM14600.908571.014112−1.032831
13OthM640.0398280.007863−2.122050.00724
Table 5. Basic data of evaluation indicators for each function group of Lianshi Lake.
Table 5. Basic data of evaluation indicators for each function group of Lianshi Lake.
Serial NumberFunction Group E E P / R O I T L R T I  
1OthP0.0260.5200.1713.3030.622
2Comc0.2480.1260.0482.9580.122
3Cruc0.3790.1300.1862.2420.17
4Bigc0.0000.2190.2552.5060.053
5Silc0.0000.2080.1712.2150.018
6HerF0.7780.3080.0002.0000.425
7OthF0.0000.4970.1252.8630.749
8Macc0.0620.3360.2322.3530.328
9OthB0.0990.5000.0052.0050.867
10Zoop0.6060.1670.0092.0090.701
11phyt0.308-0.0001.0000.844
12SubM0.186-0.0001.0001
13OthM0.000-0.0001.0000.007
14Detritus0.272-0.2511.000-
Table 6. The membership degrees corresponding to different evaluation indicators of different biological chains with the Zoop group.
Table 6. The membership degrees corresponding to different evaluation indicators of different biological chains with the Zoop group.
Zoop Group Project Evaluation Index
E E 0.2570.3380.2600.2510.3300.2970.4570.4570.457
P / R 0.3280.3880.3960.3410.2720.2710.3320.1880.193
R T I 0.6310.790.7590.6240.5840.5720.7650.5210.533
O I 0.0580.0460.0620.1370.1220.0760.0670.090.132
T L 0.4180.3230.3600.3890.3800.4390.3190.2470.279
Table 7. Comparison table for overall characteristics of the Lianshi Lake ecosystem.
Table 7. Comparison table for overall characteristics of the Lianshi Lake ecosystem.
IndexCurrent
State
1.5 Times2 Times3 Times4 TimesUnit
Zoop biomass7.8511.77515.723.5531.4t/km2/year
Macc biomass1.582.373.164.746.32t/km2/year
Phyt Eco-nutrition efficiency0.3080.4580.6070.9061.204-
SubM Eco-nutrition efficiency0.1860.1860.1860.1860.186-
The amount of material flowing into the second trophic level from primary producers304343535663828310904t/km2/year
The amount of debris flowing into the primary producer7618630849982378−241.9t/km2/year
Total primary production/total respiratory volume (TPP/TR)9.2246.4614.9723.4032.587-
Finn’s cycle index (FCI)13.59%14.33%15.23%17.51%20.38%-
Finn’s average energy flow path length (FCL)2.8542.9933.1323.4103.688-
Connection coefficient (CI)0.2250.2260.2260.2260.226-
System omnivorous degree (SOI)0.0920.0920.0930.0930.094-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, P.; Cui, X.; Luo, H.; Peng, W.; Gao, Y. Research on the Application of Typical Biological Chain for Algal Control in Lake Ecological Restoration—A Case Study of Lianshi Lake in Yongding River. Water 2021, 13, 3079. https://doi.org/10.3390/w13213079

AMA Style

Zhang P, Cui X, Luo H, Peng W, Gao Y. Research on the Application of Typical Biological Chain for Algal Control in Lake Ecological Restoration—A Case Study of Lianshi Lake in Yongding River. Water. 2021; 13(21):3079. https://doi.org/10.3390/w13213079

Chicago/Turabian Style

Zhang, Pengfei, Xiaoyu Cui, Huihuang Luo, Wenqi Peng, and Yunxia Gao. 2021. "Research on the Application of Typical Biological Chain for Algal Control in Lake Ecological Restoration—A Case Study of Lianshi Lake in Yongding River" Water 13, no. 21: 3079. https://doi.org/10.3390/w13213079

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop