阅读说明：本技术 一种基于Ecopath模型的生态系统评估方法 (Ecosystem assessment method based on Ecopath model ) 是由 王赛 王旭 王团团 王晓迪 王定莹 于 2021-08-25 设计创作，主要内容包括：本发明涉及生态系统研究技术领域,公开了一种基于Ecopath模型的生态系统评估方法,包括如下步骤：A)在研究地区设置多个采集点,对各个采集点的物种进行样品采集,根据所述物种的种类和数量的变化,进行区域划分；B)将所有所述物种进行功能组划分,并将各所述功能组的输入参数输入至Ecopath模型,获取不同区域的食物网结构；C)根据所述食物网结构,对各营养级间能量传输量进行分析,获取各所述营养级之间的传输效率。该方法能掌握不同区域生态系统结构与功能的基本信息,也能为其管理与修复提供针对性的理论支持。(The invention relates to the technical field of ecological system research, and discloses an ecological system evaluation method based on an Ecopath model, which comprises the following steps: A) setting a plurality of collection points in a research area, collecting samples of species of the collection points, and dividing the area according to the variety and the quantity of the species; B) dividing all the species into functional groups, and inputting input parameters of all the functional groups into an Ecopath model to obtain food net structures in different areas; C) and analyzing the energy transmission quantity among all the nutrition levels according to the structure of the food net to obtain the transmission efficiency among all the nutrition levels. The method can master the basic information of the structure and the function of the ecosystem of different areas and also can provide targeted theoretical support for the management and the restoration of the ecosystem.)

1. An ecosystem assessment method based on an Ecopath model is characterized by comprising the following steps:

A) setting a plurality of collection points in a research area, collecting samples of species of the collection points, and dividing the area according to the variety and the quantity of the species;

B) dividing all the species into functional groups, and inputting input parameters of all the functional groups into an Ecopath model to obtain food net structures in different areas;

C) and analyzing the energy transmission quantity among all the nutrition levels according to the structure of the food net to obtain the transmission efficiency among all the nutrition levels.

2. The Ecopath model-based ecosystem assessment method according to claim 1, wherein in the step C), further comprising:

and acquiring mixed nutrition influence values of different regions according to direct or indirect predation or competitive relations of each functional group.

3. The Ecopath model-based ecosystem assessment method according to claim 1, wherein in the step C), further comprising:

and acquiring the criticality index of each function group, and analyzing the key function groups in different areas.

4. The Ecopath model-based ecosystem assessment method according to claim 1, wherein in the step C), further comprising:

and acquiring the total system flux, the total primary production/total respiratory capacity, the total primary production/biomass and the system circulation index of different regions, and judging the maturity of the ecosystems of the different regions.

5. The eco-system assessment method based on the Ecopath model according to any one of claims 1 to 4, wherein in the step A), when the research area is a river, the collection points are arranged along the trend of the river, and the division into regions is performed according to the variation of the species and the number of fishes.

6. An ecosystem assessment method according to claim 5, wherein in step B) the functional group includes at least one debris group.

7. The Ecopath model-based ecosystem assessment method according to claim 5, wherein in step B), the input parameters are biomass, a production quantity to biomass ratio, a consumption quantity to biomass ratio, an ecological nutrition efficiency and an unassimilated real object ratio.

8. The Ecopath model-based ecosystem assessment method according to claim 7, wherein the ecological nutrition efficiency has a numerical range of 0.01 to 0.99, and the ratio of the production amount to the consumption amount has a numerical range of 0.05 to 0.3.

9. The Ecopath model-based ecosystem assessment method according to claim 8, wherein in step C), further comprising:

and assigning values to each data according to the source mode of each input parameter to obtain the Pedigree indexes of different areas.

10. The Ecopath model-based ecosystem assessment method according to claim 8, wherein in step C), further comprising:

and generating graphs according to the changes of the biomass, the ratio of the production capacity to the biomass and the ratio of the consumption to the biomass, and observing whether the biomass, the ratio of the production capacity to the biomass and the ratio of the consumption to the biomass in different areas are in an exponential function relationship.

Technical Field

The invention relates to the technical field of ecological system research, in particular to an ecological system evaluation method based on an Ecopath model.

Background

The river ecosystem refers to a unified body of interaction between a biological community in a river and a river environment, and is an important water body ecosystem. The river ecosystem consists of two major parts, biological and non-biological environments. The abiotic environments mainly comprise riverbed landforms, substrates, flow, depth, inorganic substances (such as nitrogen and phosphorus) participating in substance circulation and organic substances (such as carbohydrate, humus and the like) connecting organisms and abiotic bodies, and the abiotic environments provide substance conditions on which various organisms in a river ecosystem live. The biological part can be mainly classified into producers, consumers and decomposers. With the rapid development of economy, the development of river resources by human beings also enters a peak period, so that the river ecosystem bears huge pressure.

In recent years, with the rapid development of the economy of the Zhujiang Delta area, the influence of human activities such as sewage discharge, hydraulic engineering construction, artificial fishing and the like on the ecological system of the Zhujiang river is more and more serious. The urbanization process is accelerated, so that the utilization types of soil along the banks of the Zhujiang river are changed, the land used for buildings is greatly increased, the cultivated land is sharply reduced, and the hydrological situation of the local river is changed. As a driving force of a river ecosystem, river hydrology is directly related to the environment of a river biological habitat, and the change of the river hydrology has important influence on the structure and the function of the ecosystem.

The traditional river channel is long, the types and the degrees of interference of artificial activities are different in different river sections, the influence and the impact effect on an ecosystem are obviously different, the structure and the function of the ecosystem are different, the existing river condition research is mainly based on water quality assessment, the effect of organisms which are important components of the ecosystem is ignored, the river ecology is researched by more biological diversity indexes or single indicator species organisms afterwards, and the overall characteristics of the river ecosystem are difficult to reflect comprehensively.

Disclosure of Invention

The technical problem to be solved by the invention is as follows: the ecosystem assessment method based on the Ecopath model for researching species by dividing regions is provided, the energy flow process of the ecosystem structure in different regions is analyzed, and the development condition and the attribute of the ecosystem are assessed in a more targeted manner.

In order to solve the above technical problem, the present invention provides an ecosystem assessment method based on an Ecopath model, which includes the following steps:

A) setting a plurality of collection points in a research area, collecting samples of species of the collection points, and dividing the area according to the variety and the quantity of the species;

B) dividing all the species into functional groups, and inputting input parameters of all the functional groups into an Ecopath model to obtain food net structures in different areas;

C) and analyzing the energy transmission quantity among all the nutrition levels according to the structure of the food net to obtain the transmission efficiency among all the nutrition levels.

Preferably, step C) further comprises:

and acquiring mixed nutrition influence values of different regions according to direct or indirect predation or competitive relations of each functional group.

Preferably, step C) further comprises:

and acquiring the criticality index of each function group, and analyzing the key function groups in different areas.

Preferably, step C) further comprises:

and acquiring the total system flux, the total primary production/total respiratory capacity, the total primary production/biomass and the system circulation index of different regions, and judging the maturity of the ecosystems of the different regions.

Preferably, in step a), when the research area is a river, the collecting points are arranged along the trend of the river, and the area division is performed according to the change of the type and the number of the fishes.

Preferably, in step B), the functional group comprises at least one group of debris.

Preferably, in step B), the input parameters are biomass, production to biomass ratio, consumption to biomass ratio, ecological nutritional efficiency and unassimilated real-object ratio.

Preferably, the ecological nutritional efficiency has a value ranging from 0.01 to 0.99, and the ratio of the production amount to the consumption amount has a value ranging from 0.05 to 0.3.

Preferably, step C) further comprises:

and assigning values to each data according to the source mode of each input parameter to obtain the Pedigree indexes of different areas.

Preferably, step C) further comprises:

and generating graphs according to the changes of the biomass, the ratio of the production capacity to the biomass and the ratio of the consumption to the biomass, and observing whether the biomass, the ratio of the production capacity to the biomass and the ratio of the consumption to the biomass in different areas are in an exponential function relationship.

Compared with the prior art, the ecosystem assessment method based on the Ecopath model has the beneficial effects that:

the invention divides different regions by collecting samples of species at different collection points to indicate the variation dividing limits of species types and quantities. And performing functional group division on all the species, and inputting the sampled or calculated input parameters into the Ecopath model, so that the food net structures in different areas are obtained, and the whole energy flow process of an ecological system can be covered. And the maturity of the ecological systems in different areas is judged by analyzing the transmission efficiency of energy among all the nutrition stages, so that not only can basic information of the structures and functions of the ecological systems in different areas be mastered, but also targeted theoretical support can be provided for the management and the repair of the ecological systems.

Drawings

Fig. 1 is a block flow diagram of an ecosystem assessment method based on an Ecopath model according to a preferred embodiment of the present invention.

Fig. 2 is a schematic diagram of the division of the pearl river region according to the preferred embodiment of the invention.

Fig. 3 is a schematic diagram of criticality indexes of functional groups of the agricultural area (Z1) of the Zhujiang river in the preferred embodiment of the invention.

Fig. 4 is a schematic diagram of criticality index of each functional group of the industrial zone of the pearl river (Z2) according to the preferred embodiment of the invention.

Fig. 5 is a schematic diagram of criticality indexes of functional groups in the urban area (Z3) of the pearl river according to the preferred embodiment of the invention.

Fig. 6 is a schematic diagram of criticality indexes of functional groups of a river estuary region (Z4) of a pearl river according to a preferred embodiment of the invention.

Detailed Description

It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings.

The present application is described in further detail below with reference to specific examples, which should not be construed as limiting the scope of the invention as claimed.

As shown in fig. 1, the present invention provides an ecosystem assessment method based on an Ecopath model, which includes the following steps:

A) setting a plurality of collection points in a research area, collecting samples of species of the collection points, and dividing the area according to the variety and the quantity of the species;

B) dividing all the species into functional groups, and inputting input parameters of all the functional groups into an Ecopath model to obtain food net structures in different areas;

C) and analyzing the energy transmission quantity among all the nutrition levels according to the structure of the food net to obtain the transmission efficiency among all the nutrition levels.

The ecosystem assessment method based on the technology is characterized in that samples of species are collected at different collection points to indicate the variation division limit of the species and the number, and different regions are divided. And performing functional group division on all the species, and inputting the sampled or calculated input parameters into the Ecopath model, so that the food net structures in different areas are obtained, and the whole energy flow process of an ecological system can be covered. And the maturity of the ecological systems in different areas is judged by analyzing the transmission efficiency of energy among all the nutrition stages, so that not only can basic information of the structures and functions of the ecological systems in different areas be mastered, but also targeted theoretical support can be provided for the management and the repair of the ecological systems.

The present embodiment is illustrated by using a pearl river ecosystem as an example, and it is understood that the ecosystem assessment method of the present invention can also be applied to other ecosystems.

The river mouth section of the Zhujiang river passes through the urban area of Guangzhou city with dense population, and the ecological system of the Zhujiang river is greatly influenced by human activities. Along the river trend, the fish samples collected from 14 points from the upstream to the downstream river reach of the river mouth of the pearl river show that the types and the number of the fishes have certain changes, more clastic fishes exist in the middle upstream area, the downstream carnivorous fishes are remarkably increased, and the saline fishes such as the seven-thread anchovies, the sea bass and the like gradually appear in the middle downstream area. The species and the quantity of river fishes are closely related to the located ecological system, the coastal zone environment is reflected to a certain extent, the feeding habits of the same kind of fishes in different habitats are different, the ecological system is researched by dividing the pearl estuary into 4 areas by taking the fishes as indicator organisms, and the ecological system is divided into an agricultural area (Z1), an industrial area (Z2), an urban area (Z3) and a estuary area (Z4) according to the coastal zone environment of each area, as shown in figure 2, the periphery of the agricultural area (Z1) is mostly cultivated land, and the greening area is large; the industrial area (Z2) is mainly a river of stone wells which converges into a Zhujiang section and a section from the downstream to the front section of a rear channel, peripheral plants are more and seriously polluted, and fishes are mainly resistant to pollution seeds; the urban area (Z3) extends from the university city to the east river and merges into the Zhujiang river segment, the river segment directly passes through the Guangzhou urban area, and the number of surrounding residents is very large; the estuary area (Z4) is the place where the Zhujiang enters the sea, the saline water and the freshwater are staggered, the habitat diversity of the water area is high, and the types and the number of fishes are extremely large.

The biological groups of fish in the sequin river network ecosystem evolved, and therefore the fish species were divided into "species" level, according to this division method, the number of functional groups in the Z1 region was 39 (27 fish functional groups, 7 macroinvertebrate functional groups, 3 producer functional groups and 1 each of zooplankton functional groups and debris functional groups), the number of functional groups in the Z2 region was 20 (10 fish functional groups, 5 macroinvertebrate functional groups, 3 producer functional groups and 1 each of zooplankton functional groups and debris functional groups), the number of functional groups in the Z3 region was 36 (24 fish functional groups, 7 macroinvertebrate functional groups, 3 producer functional groups and 1 each of zooplankton functional groups and debris functional groups), the number of functional groups in the Z4 region was 47 (36 fish functional groups, 6 macroinvertebrate functional groups, producer functional group 3 and zooplankton functional group and crumb functional group 1 each). Wherein the primary producers mainly comprise phytoplankton, aquatic plants and epiphytic algae, and the invertebrates mainly comprise macrobenthos including shrimps, crabs, insects, polychaetes, oligochaeta, bichaeta and gastropodae.

The specific values of the input parameters in the Z1-Z4 areas of the Yangtze river mouth are not detailed here. The functional groups are arranged in the order of fish, benthos, producer and debris, and the interior of the functional groups is arranged according to the order of the nutritional level. The table lists biomass (B), production to biomass ratio (P/B), consumption to biomass ratio (Q/B), ecological nutritional efficiency (EE), and non-assimilated physical proportions (UC) for each functional group. Wherein the ecological nutrition efficiency of each functional group is 0.01-0.99.

After the original input parameters are input, the operation result generally needs to be checked, if the result which does not accord with the ecological rule appears, the relevant parameters should be adjusted, the operation is carried out repeatedly, and unreasonable factors are minimized, so that the model can objectively reflect the actual condition of the ecological system while achieving the steady state. Among all the input parameters, the ecological nutritional efficiency (EE) is a less accessible parameter, which reflects the fraction of the functional group that is ingested and exported in the system, and therefore ranges between 0 and 1, and furthermore the production/intake (P/Q), also called the conversion of raw food (GE), is generally between 0.05 and 0.3, i.e. the intake is about 3 to 10 times the production.

After the model is built, the operation software evaluates the quality of the model through the Pedigree index, and assigns values to each data based on the data source mode to finally obtain the Pedigree index of the whole model. The Pedigree indexes of the four areas of Zhujiang Z1-Z4 are 0.59, 0.58 and 0.57 respectively, which shows that the data reliability is high and the overall quality of the model is high. Model Pre-balance diagnosis (Pre-balance) is adopted to observe the mass balance of model balance, the Pre-balance program reduces the uncertainty of input parameters by using a basic ecological theory, and the B, P/B, Q/B change condition is visually compared through a simple graph method. Generally, biomass is gradually reduced as the nutritional grade is increased. Each functional group parameter was considered biologically reasonable if an exponential relationship was observed between the increase in 3 parameters and the trophic level. The experimental result can show that B, P/B, Q/B in the four areas of the river mouth of the Zhujiang river has better index correlation with the nutrition grade, which indicates that the model has better quality balance relation.

The food net structure of the Zhujiang Z1-Z4 area can be obtained after Ecopath model operation, and the nutrition level ranges of the four areas are 1-3.36, 1-3.26, 1-3.47 and 1-3.95 respectively. The first nutritional level is phytoplankton, aquatic plants and epiphytic algae, the second nutritional level is mainly benthos, debris-feeding fish, filter-feeding fish, herbivorous fish and omnivorous fish, and the third nutritional level is carnivorous fish. Wherein, the types of carnivorous fishes in the four areas are respectively 8, 1, 6 and 14, and the fishes with the highest nutrition level are respectively erythroculter hainanensis, catfish, erythroculter hainanensis and odontobutis obscurus. The functional groups with the highest nutritional level in the upstream agricultural area (Z1) and the midstream urban area (Z3) are Erythroculter hainanensis, fish and shrimp are mainly used as food, the food sources in the two areas are relatively rich, and large carnivorous fishes do not prey on the fish and shrimp, so that the nutritional level is highest. Clarias fuscus is a typical contaminant resistant species, and is also the only carnivorous fish in the industrial zone (Z2), with the highest nutritional level in that zone being lower than in other zones due to the lack of mutual predation by other carnivorous fish. The estuary region (Z4) is most complex in predation, and even carnivorous fish will be used as feed for each other, so its maximum nutritional level is close to 4.

The energy flow paths in the four regions are essentially identical and all consist of a food chain starting from phytoplankton, aquatic plants, epiphytic algae and debris: phytoplankton → zooplankton, filter feeding fish → carnivorous fish; aquatic plants → grass carp → carnivorous fish; epiphytic algae → shrimps and crabs → omnivorous fishes → carnivorous fishes; clastics → zooplankton → benthos → omnivorous fish → carnivorous fish. The industrial area (Z2) has the simplest food net structure and the least number of food chains; the agricultural (Z1) and urban (Z3) food nets are relatively complex; the estuary region (Z4) had the highest food web complexity and the highest number of food chains. The four regions are not very different in the number of producer functional groups, benthos and plankton functional groups, and the number of fish species is the main reason for determining the number of food chains. The industrial zone (Z2) mainly causes the species of fish to decline due to relatively serious pollution, destroying the structure of the ecosystem; the pollution conditions of an agricultural area (Z1) and an urban area (Z3) are relatively good, and some top-grade carnivorous fishes which have higher requirements on the environment still inhabit the pollution conditions; the estuary area (Z4) is positioned at the junction of the salt and the fresh water, so that the habitat is richer, the inclusion is stronger, and both the fresh water fish and the salt water fish can use the area as the habitat.

On the basis of a food net structure, the energy transmission quantity among all nutrition levels is analyzed, the energy flow among all nutrition levels in the Zhujiang Z1-Z4 area is decomposed, and the energy flow mainly occurs in the first five nutrition levels, so that the change of the transmission efficiency among the first five nutrition levels is mainly analyzed in the embodiment. The values of the efficiency of the interstage transport of nutrients between the food chain of the primary producer, the food chain of the scraps and the whole food network system (which integrates the primary producer and the scraps) of the upstream agricultural sector (Z1) are lower, respectively 5.09%, 5.73% and 5.49%, while those of the industrial sector (Z2) decrease continuously by 3.48%, 3.59% and 3.56%, respectively, with a return increase in the urban sector (Z3), respectively 5.53%, 5.94% and 5.80%, and a maximum in the region of Z4, respectively 9.54%, 9.88% and 9.78%. Comparing the transmission efficiency values between nutrient levels in each region, the transmission efficiency values between nutrient levels in the second nutrient level are higher in the upstream and downstream regions, respectively 9.01% and 11.3%, and relatively lower in the middle of Zhujiang, respectively 5.21% and 4.69%. While the third to fifth nutritional grades all showed a tendency to gradually decrease from the agricultural zone (Z1) to the industrial zone (Z2) and then to the estuary zone (Z4), and all reached a maximum in the estuary zone. In general, the energy flow of the river ecosystem is dominated by the debris food chain, which is especially apparent in the estuary section of the pearl river. The energy influx from the debris in the 4 zones was 60%, 71%, 65% and 69%, respectively, indicating that the debris food chain of the Zhujiang is the primary energy transfer channel. The clastic fish species and biomass are higher than other predatory fishes throughout the whole pearl estuary area, the reason is supposed to be mainly related to the pollution condition of the pearl estuary, and the clastic fish species are generally stronger in adaptability (such as squaliobarbus curriculus and loach) and even resistant to pollution (such as clarias fuscus). Generally, the average nutrient interstage transmission efficiency values of other regions except the downstream estuary region are lower than normal values and lower than other fresh water ecosystems.

Comparing the energy conversion efficiency of the industrial zone (Z2) and the estuary zone (Z4) it can be concluded that the energy conversion efficiency of the industrial zone is lower than that of the estuary zone, mainly because of the lower conversion of primary producers and debris by the second nutrient grade of the industrial zone. A large number of primary producers in the industrial area are not utilized, the main reason is that the biomass of phytoplankton in the industrial area is huge, a large number of unused phytoplankton flow into the debris, and the utilization rate of the debris is low, so that the upward transfer of energy is blocked, and the structure of the food net tends to be simplified.

The mixed nutrient impact can be used to assess the magnitude of the impact of changes in biomass of one functional group on other functional groups. Direct or indirect predation or competition relations exist among the functional groups, the organisms mutually restrict and influence each other to enable the ecological system to be relatively stable, and the nutrition influence matrix output by the model quantifies the direct or indirect action relations among the functional groups. In the ecological system of the pearl estuary, fishing has inhibition effect on most fishes and promotion effect on some benthos and plankton. Due to the inter-species competition relationship, each functional group has a certain inhibition effect on the functional group. Debris, plankton and epiphytic algae contribute to most functional groups, primarily because they are the major food sources for most organisms. The aquatic plants have an inhibitory effect on most functional groups, mainly due to the competitive relationship between the omnivores of few part of the aquatic plants and other organisms.

In the region Z1, the inhibition of the other 37 functional groups of the system by the class Oligochaeta reached a value of-1.28, and the inhibition of the 32 functional groups by the class Gastropoda reached a value of-0.57. The most potent inhibitory effect on the other functional groups was zooplankton, which was-1.70, and the overall feedback effect (sum of absolute values of inhibitory and promoting effects) was also the most potent, 2.51, on the other functional groups. Both shrimps and polychaetes in the region Z2 produced inhibition of 16 other functional groups, the most abundant, with values of-0.43 and-0.57, respectively. The most potent inhibitory effect on other functional groups in the region Z2 was also zooplankton, reaching-2.21, and it had the strongest overall feedback effect, reaching 3.09. The gastropoda class Z3 has an inhibitory effect on 34 functional groups, with a value of-0.72. The inhibition effect on other functional groups is the glossogobius with the strongest inhibition effect, which reaches-3.92, and the comprehensive feedback effect on other functional groups is also 5.36 at most. In the region Z4, the gastropoda class has an inhibitory effect on 40 functional groups, with a value of-0.72. In this region, the zooplankton inhibited the other functional groups as much as-2.58. The comprehensive feedback effect of the tongue sole reaches 4.33 most intensely.

The key functional group is positioned at a fulcrum in the ecological system, has important influence on the energy flow of the whole ecological system, and has important significance for maintaining the energy flow steady state of the food net. The key functional group indexes of different regions of the Zhujiang are shown in figures 3 to 6, wherein the horizontal axis represents the nutrition level of the functional group, the vertical axis represents the criticality index, the higher the value of the index, the higher the criticality of the functional group, and the circle size is in direct proportion to the biomass of the group. As can be seen from the figure, the functional groups with the highest criticality index in the agricultural area (Z1) are erythroculter hainanensis, pelteobagrus fulvidraco, order testudinate, bream and elopichia bambusa respectively; in the industrial area (Z2), the functional groups with the highest key degree index in the industrial area (Z2) are carps, crucians, oligochaeta, clarias fuscus and perches respectively, the functional groups with the highest key degree index in the urban area (Z3) are glossogogogogogogogel, shortgey koi, clarias leather, coilia ectenes and bream respectively, and the functional groups with the highest key degree index in the estuary area (Z4) are tongue sole, globefish, goby tail, goby, goby koch and perch respectively.

The key species in the agricultural (Z1), urban (Z3) and estuary (Z4) regions are top-grade carnivorous fish, while the key functional group in the industrial (Z2) region is carps of omnivorous fish. The key degree of carnivorous fishes in the agricultural area is slightly higher than that of omnivorous fishes, the key degree of the carnivorous fishes in the urban area and the estuary area is obviously higher than that of the omnivorous fishes, and particularly, all the 6 functional groups with the highest key degree are carnivorous fishes in the downstream estuary area. From the level of nutritional grade, the key functional groups are mainly those predators with high nutritional grade, and almost show a trend that the criticality gradually increases along with the increase of the nutritional grade in the middle and downstream areas of the Zhujiang river. The criticality of phytoplankton and aquatic plants, etc. at the first trophic level is very low in each area. Generally, the criticality of organisms with low nutritional levels is relatively low, so that the 5 most critical functional groups in 4 regions of the pearl estuary are all fishes. And the key degree index of the oligochaeta in the industrial area (Z2) reaches 3 rd high and is even higher than that of the only carnivorous fish clarias fuscus in the area, mainly because the water quality in the area is poor, and the oligochaeta is mainly limnodrilus which feeds on scraps and is a key functional group for the energy transfer of a scrap food chain.

Through calculation, some ecological parameters about the ecosystem can be obtained, and the scale, maturity, stability and other states of the ecosystem are generally described. The ecological parameters include total consumption, total output, total respiration, total influent debris, total system flux, total production, nutrient level of the fish harvest, net primary production, total primary production/total respiration, total primary production/biomass, total biomass/total throughput, total biomass (including debris), total fishing, connectivity index, system omnivory index and aroma diversity index, the detailed values of the different regions not being described in detail herein. The total system flux of the 4 areas is 2660t/km respectively according to the result²/year、2594t/km²/year、4972t/km²/year、6039t/km²And/year. Total system throughput is respectively845t/km²/year、525t/km²/year、1469t/km²/year、1652t/km²And/year. The nutrient levels of the fish catch are respectively 2.18, 2.20, 2.14 and 2.30. The total fishing amount is respectively 8.97t/km²/year、5.09t/km²/year、10.83t/km²/year、17.21t/km²And/year. The total system flux shows that the system is the largest in estuary areas and relatively smaller in other areas. The value of total primary production/total breath may reflect the degree of stability of the system, with values closer to 1 indicating more stable systems. The maximum total primary production/total respiration in the Z2 region and the minimum total primary production/total respiration in the Z4 region indicate the relatively least stable Z2 region and the relatively most stable Z4 region. The system omnivorous index reflects the complexity of the system, the closer the value to 1 indicates the more complex the system, and the food net structure in the Z2 area is relatively simplest and the system in the Z4 area is more complex. As can be seen from the system connectivity index, region 0.2054 in region Z2 is significantly higher than other regions by about 0.14, reflecting that the food web connectivity is tighter in region Z2. The system cycle index represents the proportion of the cycle energy flow in the total system flow, and is a reflection of the system maturity. The systems maturity in the Z3 and Z4 regions is higher, and the Z1 and Z2 regions are relatively lower. From the fragrance diversity index, it can be seen that the Z4 region is highest at 3.47, followed by the Z3 and Z1 regions, respectively at 3.21 and 3.15, and the Z2 region is only 2.56 in species diversity.

The ecosystem of the estuary region is the largest in scale, followed by agricultural and industrial areas upstream of the midstream urban area, respectively, as is the trend of biomass size in each region. On one hand, the downstream river channel is wider, the habitat of organisms is larger, and the contained biomass is more; on the other hand, the industrial area is very polluted by human beings, and the biomass is obviously lower than that of other areas. The system connection index shows that the connection index reaches up to 0.21 in the industrial area, and the connection index reaches about 1.46 in other areas, mainly because the types of fishes in the industrial area are less, particularly carnivorous fishes, so that the structure of the food net is relatively simple, and the possible paths for energy transfer are less. Although the connection tightness of the food net in the industrial area is high, the resistance to the outside is weak, and the energy flow transmission is obstructed by the change of any species, so that the structure and the function of an ecosystem are greatly influenced.

By comprehensively evaluating system parameters such as total primary production/total respiration, total primary production/biomass, systemic omnivorous index and systemic circulation index, the estuary region (Z4) was considered to have the highest system maturity and the industrial region (Z2) was considered to have the lowest system maturity in each of the Zhujiang regions. The main reasons can be attributed to: (1) the ecological damage of the industrial area to the local river reach is serious, so that a large amount of fishes disappear; (2) the habitat of the estuary area is more complex, and the fish has stronger adaptability to the mixing of the salt water and the fresh water and has stronger capability of resisting external interference. The systematic omnivory index in the estuary region reached 0.20, being the highest of the 4 regions, indicating that the predation relationship in estuary regions is more complex and that the system is more mature due to longer food chain length and higher influx of debris.

To sum up, the embodiment of the invention presents the structure and function of the ecosystem among 4 different areas of the Yangtze river estuary section from the integral perspective by constructing the Ecopath model of 4 different areas upstream to downstream of the Yangtze river estuary section, evaluates the development condition and attribute of the Yangtze river ecosystem, provides scientific basis for the protection and sustainable development of the Yangtze river biological resources, provides reasonable suggestions for the management of the Yangtze river, and can mainly provide the following conclusions: (1) the ecological system food nets of 4 areas of the Zhujiang estuary section are similar in structure and comprise 4 food chains with phytoplankton, aquatic plants, epiphytic algae and debris as starting points, and the energy flows of the 4 areas mainly comprise the debris food chains. The highest nutritional level in the industrial area is the lowest, the food chain is relatively short because the area only contains one type of carnivorous fish, the food net structure in the midstream urban area and the upstream agricultural area is relatively complex, and the highest nutritional level in the downstream estuary area is the highest, which is probably because the estuary area is directly influenced by the sea, and the saline carnivorous fish has higher predation intensity, so that the length of the food chain is increased. (2) The energy transmission efficiency of 4 areas at the bead river mouth section is lower. The industrial area is mainly due to the fact that the biomass of fishes is low, and the utilization rate of functional groups with low nutrition level is low. Although the transmission efficiency of the agricultural region is not greatly different from that of the urban region, the agricultural region is mainly due to the high transmission efficiency of aquatic plants, and the urban region is mainly due to the high energy transmission efficiency of high nutrition level, and the biomass of carnivorous fishes is high. The ecological system structure of the downstream estuary area is relatively perfect, and the maturity is high. (3) The total system flux of 4 areas of the Zhujiang river mouth section is 2660t/km2/year, 2594t/km2/year, 4972t/km2/year and 6039t/km2/year respectively, which shows that the size of the system of the Zhujiang river mouth is gradually increased from the upstream to the downstream, and the reason is mainly that the biomass which can be accommodated is higher due to the wider downstream river channel. The total primary production/total respiration value reflects the degree of stability of the system, with values closer to 1 indicating more stable systems. The total primary production/total respiration values of 4 areas from upstream to downstream are 1.95, 0.90, 2.03 and 1.28 respectively, the total primary production/total respiration value of the industrial area is 0.90 less than 1, which indicates that the area is polluted by industry and the ecosystem is in degeneration, and the agricultural area and the urban area have higher total primary production/total respiration values compared with the estuary area, which indicates that the two ecosystems are still in development, and the estuary area ecosystem is mature and is more stable.

The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and substitutions can be made without departing from the technical principle of the present invention, and these modifications and substitutions should also be regarded as the protection scope of the present invention.