郭思雅，王京剛，王穎，陳澤豪，嚴俊

1. 北京化工大學化學工程學院，北京 100029 2. 北京師范大學環(huán)境學院，北京 100875

2017年4月，國務院宣布在河北保定市內(nèi)設立雄安新區(qū)，以期解決京津冀地區(qū)發(fā)展不平衡的問題，并為我國經(jīng)濟帶來更多的發(fā)展機遇[1]。作為黨中央的歷史性戰(zhàn)略決策，將其規(guī)劃打造成一座生態(tài)城市也成為雄安新區(qū)重點任務之一。而作為承載雄安新區(qū)的“水核心”——白洋淀，由于其獨特的地理位置和豐富的資源受到了國家的高度重視。白洋淀從20世紀80年代開始，由于干旱以及人類活動等原因，湖泊出現(xiàn)干涸，淀區(qū)數(shù)量銳減，淀區(qū)濕地面積也減少了33.91%，淀區(qū)水質(zhì)以Ⅳ和Ⅴ為主。1980—2010年間，白洋淀底棲動物、沉水植物種類減少，沉水植物優(yōu)勢物種由馬來眼子菜、菹草、輪葉黑藻、穗狀狐尾藻和大茨藻轉(zhuǎn)變?yōu)橐越痿~藻和輪藻，生物多樣性降低[1-2]；浮游植物群落優(yōu)勢物種由指示中營養(yǎng)化水體的硅藻逐漸演替為指示水體富營養(yǎng)化的綠藻和藍藻類[3]。為此，政府加大了對白洋淀環(huán)境治理工作的支持和投入[4-6]。從“十一五”規(guī)劃到目前的“十三五”白洋淀國家水專項實施，政府部門陸續(xù)開展了生態(tài)濕地凈化系統(tǒng)建設、河道清潔專項行動等，建立了生態(tài)補水系統(tǒng)，使白洋淀濕地面積增加了約14%，底棲動物和沉水植物生物多樣性均有所增加[7]，修復工作初見成效。

目前，針對白洋淀生態(tài)系統(tǒng)的科學研究主要集中于生物物種組成、數(shù)量變化或其他時空變化特征等[8-12]。生態(tài)系統(tǒng)中物種間關(guān)系復雜，直接或者間接影響著生態(tài)系統(tǒng)中其他物種的物種豐富度、物質(zhì)輸送和能量流動等，而“十一五”后白洋淀生態(tài)系統(tǒng)的上述變化特征如何，目前鮮有報道。由Polovina[13]提出的Ecopath模型，以食物網(wǎng)為主線，以建模方式量化生態(tài)系統(tǒng)各特征參數(shù)，能夠較為清晰地反映出生態(tài)系統(tǒng)的能量流動和營養(yǎng)關(guān)系等。目前，Ecopath模型廣泛運用于水生生態(tài)系統(tǒng)的結(jié)構(gòu)和功能分析，在我國長江口、黃河口、太湖、渤海和南海等地均有應用[14-18]。從“十一五”到“十三五”，需要對白洋淀生態(tài)系統(tǒng)結(jié)構(gòu)和功能的變化有一個全新的認知，因此，本文通過構(gòu)建2010年和2018年全年白洋淀生態(tài)系統(tǒng)模型，分析生態(tài)系統(tǒng)結(jié)構(gòu)和能量流動的特征變化，以期為白洋淀生態(tài)修復提供依據(jù)。

1 研究區(qū)域介紹與數(shù)據(jù)收集(Introduction and data collection of the research area)

1.1 研究區(qū)域概況

研究區(qū)域位于河北省中部(北緯38°43′～39°02′，東經(jīng)115°38′～116°07′)，地處北京、天津和石家莊三市的三角形中心位置[19-20]。淀區(qū)總面積為336 km2，其中，600 m2以上的淀區(qū)占了近70%，白洋淀承擔著淀區(qū)周邊近百萬人口的“安居樂業(yè)”，是河北省最重要的淡水魚產(chǎn)區(qū)。除此以外，作為海河平原最大的湖泊，白洋淀還有調(diào)洪蓄水、降解污染等極為重要的生態(tài)功能，因此被譽為“華北之腎”[21]。

1.2 數(shù)據(jù)來源

根據(jù)Ecopath模型數(shù)據(jù)需求，2010年數(shù)據(jù)來源于文獻調(diào)研[22-25]，模型所采集的生物樣本取自2010年3個季節(jié)(春季、夏季和秋季)共16個點位的平均值。且于2018年3月、7月和11月進行了3次生物樣本采樣調(diào)研。2010年和2018年2個年份的采樣點基本吻合，且基本覆蓋白洋淀所有區(qū)域，采樣方法基本一致，因此，認為二者具有可比較性。2018年設置的采樣調(diào)研點如圖1所示。

圖1 2018年白洋淀采樣調(diào)研點示意圖Fig. 1 The schematic diagram of sampling and investigation points in Baiyangdian Lake in 2018

2 研究方法(Research methods)

2.1 Ecopath模型原理

Ecopath模型運用熱力學原理，使生態(tài)系統(tǒng)中功能群的輸入與輸出保持平衡。Ecopath模型涉及到了物質(zhì)平衡與能量平衡2種核心方程，最終得到方程式：

Bi×(P/B)i×EEi-Ei=ΣBj×(Q/B)j×DCij

式中：i為被捕食功能組，其中，Bi為被捕食功能組i的生物量(t·km-2)；(P/B)i為被捕食功能組i的生產(chǎn)量與生物量比值；EEi為被捕食功能組i的生態(tài)效率，Ei為凈遷移(遷出-遷入)。j為捕食功能組，其中，Bj為捕食功能組j的生物量；(Q/B)j為捕食功能組j的消耗量與生物量比值；DCij為被捕食組i所占捕食組j食物組成的比例。對于不同年份生態(tài)系統(tǒng)的每個功能組，均需要確定B、P/B和Q/B等參數(shù)，并通過調(diào)整參數(shù)來調(diào)整EE，直至Ecopath模型達到穩(wěn)態(tài)。

2.2 功能組的設置

根據(jù)Ecopath模型建模需求，以及調(diào)研白洋淀獲得的各生物食性和分類類群等特性，本研究將白洋淀生態(tài)系統(tǒng)劃分為12個功能組，各功能組物種組成如表1所示，進而利用Ecopath模型模擬白洋淀湖泊生態(tài)系統(tǒng)。

2.3 參數(shù)確定

采樣內(nèi)容包括浮游動植物、水生動物和沉水植物等，所收集的生物樣本帶回實驗室內(nèi)鑒定分析，模型所用生物量均取3次采樣均值。浮游植物、底棲藻類的采樣和計數(shù)方法參考《湖泊富營養(yǎng)化調(diào)查規(guī)范》[26]，并根據(jù)其中的藻類植物細胞平均濕重指標將浮游植物、底棲藻類的細胞數(shù)量換算為生物量。底棲動物通過采泥器進行采樣，大型底棲動物依據(jù)稱量法[27]測定生物量、小型底棲動物依據(jù)體積換算法[27]測定生物量。浮游動物利用大網(wǎng)采集，依據(jù)稱量法[28-30]測定得出生物量。沉水植物在60～80 ℃溫度下烘干至恒重稱量。魚類采樣及調(diào)查方法參考《湖泊富營養(yǎng)化調(diào)查規(guī)范》[26]和《湖泊生態(tài)系統(tǒng)觀測方法》[27]，根據(jù)2018年白洋淀的實地調(diào)查數(shù)據(jù)換算得出其生物量，年平均魚類捕撈量等數(shù)值通過安新縣水產(chǎn)畜牧局公開的相關(guān)數(shù)據(jù)計算得到。

P/B和Q/B系數(shù)估算：參考相近生態(tài)系統(tǒng)[31]對浮游植物、沉水植物P/B和Q/B的估算值，底棲藻類P/B和Q/B參數(shù)值參照淀山湖生態(tài)系統(tǒng)模型的估算值[32]。估算底棲動物P/B參數(shù)值時，底棲動物生產(chǎn)量(P)值采用Brey的經(jīng)驗公式[33]進行估計，個體質(zhì)量干濕比參照其他模型[31-32,34]，取6∶1。魚類功能組的P/B和Q/B值采用經(jīng)驗公式[35-37]計算獲得，計算所需數(shù)據(jù)來自漁業(yè)資源數(shù)據(jù)庫(www.fishbase.org)。未同化食物比例的取值可參考經(jīng)驗值[23-25]。食物網(wǎng)資料：功能組食物組成矩陣的數(shù)據(jù)主要來源于漁業(yè)資源數(shù)據(jù)庫fishbase(www.fishbase.org)和已有的歷史食性研究資料[28-30]。有關(guān)浮游動物、螺類和蝦類的食性參照其他水體已有的食性研究結(jié)果[31-33,38-40]。

2.4 模型平衡

Ecopath模型的基本輸出是建立在生態(tài)系統(tǒng)能量輸入和輸出保持平衡的過程之上，維持模型平衡的最基本的條件是EE≤1[41-42]。Ecopath模型中的功能組的輸入和輸出相等，且功能組的總效率和EE都<1，因此，在Ecopath模型中輸入原始數(shù)據(jù)，若模型不平衡，則需要反復調(diào)整不平衡功能組的相關(guān)參數(shù)直至模型達到最佳狀態(tài)。

3 結(jié)果與討論(Results and discussion)

3.1 生態(tài)系統(tǒng)結(jié)構(gòu)和能流分布

白洋淀生態(tài)系統(tǒng)的平衡Ecopath模型在2010年和2018年2個時期的輸入和輸出結(jié)果如表2所示，可以看出，相較于2010年模型，甲殼類和軟體動物生物量明顯增加，根據(jù)安新縣水產(chǎn)畜牧局公開發(fā)布的數(shù)據(jù)，2010年白洋淀養(yǎng)殖面積為2 365 hm2，水產(chǎn)養(yǎng)殖品種以草魚、鯉魚、鰱魚和鳙魚為主，到“十二五”末，養(yǎng)殖面積超過2 700 hm2，淀區(qū)內(nèi)蝦蟹螺養(yǎng)殖量增加，期間大約有5次蝦蟹螺增殖放流活動，使淀區(qū)內(nèi)甲殼類和軟體動物生物量增加，這與本調(diào)研結(jié)果一致。2010年和2018年白洋淀生態(tài)系統(tǒng)模型的營養(yǎng)級范圍分別為1～3.025和1～3.081，肉食性魚類、雜食性魚類、濾食性魚類和軟體動物營養(yǎng)級上升了0.064、0.228、0.05和0.1，其他功能組無明顯變化。與李華等[31]評價的“1980—2006年白洋淀生態(tài)系統(tǒng)”對比發(fā)現(xiàn)，1980—2006年間肉食性魚類的營養(yǎng)級由3.11降至3.01，2010—2018年白洋淀生態(tài)系統(tǒng)營養(yǎng)級略有升高，這與期間白洋淀鰱魚、鯉魚和蝦蟹螺增殖放流，以及期間每年長達7個月的“禁漁”政策有關(guān)，這些措施使得白洋淀生態(tài)系統(tǒng)多樣性增加，魚類結(jié)構(gòu)趨于復雜化狀態(tài)。初級生產(chǎn)者沉水植物、底棲藻類和浮游藻類能量轉(zhuǎn)換效率呈升高趨勢，其能量被傳遞到更高營養(yǎng)級，其他功能組能量轉(zhuǎn)換效率無顯著變化。

表1 白洋淀生態(tài)系統(tǒng)Ecopath模型各功能組物種組成Table 1 Species composition of each functional group of Ecopath model of Baiyangdian ecosystem

系統(tǒng)總流量是系統(tǒng)中輸出、呼吸、流向碎屑的量和消費量的加和，可以用來代表整個系統(tǒng)的流量規(guī)模。2010年和2018年白洋淀生態(tài)系統(tǒng)生物量、生產(chǎn)量和總流量的分布情況如表3所示。由表3可知，2010年和2018年白洋淀生態(tài)系統(tǒng)各營養(yǎng)級之間總流量隨營養(yǎng)級增加而減少；其中，2018年生態(tài)系統(tǒng)生物量、生產(chǎn)量、流向碎屑量以及總流量均高于2010年，可以看出2018年白洋淀生態(tài)系統(tǒng)規(guī)模有所增大，2018年各營養(yǎng)級生物量都高于2010年，這與2018年白洋淀調(diào)研結(jié)果相吻合；從生產(chǎn)量來看，2010年和2018年第Ⅰ營養(yǎng)級生產(chǎn)量均占總生產(chǎn)量的97%以上，白洋淀2018年生態(tài)系統(tǒng)總生產(chǎn)量是2010年總生產(chǎn)量的1.5倍。2018年流向碎屑的量占總流量的34.53%，低于2010年的41.58%，可見2018年白洋淀生態(tài)系統(tǒng)比2010年的能量利用高，但仍有大部分能量未被充分利用，同時2018年白洋淀初級生產(chǎn)者中沉水植物較2010年增加了10.77%，浮游藻類和底棲藻類生物量略有減少，說明白洋淀在此期間部分藻型湖泊順向演替為草型湖泊。

表2 2010年和2018年白洋淀生態(tài)系統(tǒng)模型參數(shù)的輸入和輸出Table 2 Inputs and outputs of model parameters of the Baiyangdian ecosystem in 2010 and 2018

3.2 營養(yǎng)級的能量轉(zhuǎn)換效率和營養(yǎng)關(guān)系

2010年和2018年白洋淀營養(yǎng)層級間能流效率如圖2所示，白洋淀生態(tài)系統(tǒng)在2010年和2018年2個時期能量流動均可以合并成為5個整合營養(yǎng)級，由于能量來源不同，白洋淀生態(tài)系統(tǒng)存在2種不同的食物鏈，一種是以沉水植物和其他植物為起點的牧食鏈，另一種是以碎屑為起點的腐食鏈。其中，2010年和2018年白洋淀生態(tài)系統(tǒng)牧食鏈的主要結(jié)構(gòu)可以表示為：從第Ⅰ營養(yǎng)級(沉水植物、浮游植物和底棲藻類)開始，流經(jīng)第Ⅱ營養(yǎng)級(軟體動物、草食性魚類和浮游動物)以及第Ⅲ營養(yǎng)級(雜食性魚類和肉食性魚類)進行能量流動；生態(tài)系統(tǒng)腐食鏈為再循環(huán)有機物轉(zhuǎn)化為碎屑，流經(jīng)食碎屑者(底棲生物和浮游動物)、肉食性魚類進行能量流動。

表3 2010年和2018年白洋淀生態(tài)系統(tǒng)生物量、生產(chǎn)量、流向碎屑量和總流量在各營養(yǎng)級的分布Table 3 Distribution of biomass, production, debritus flow and total throughput in Baiyangdian ecosystem in 2010 and 2018

2018年白洋淀生態(tài)系統(tǒng)各營養(yǎng)級間傳輸效率如表4所示，從表4和圖2可知，2010年和2018年2個時期能量分布主要集中于第Ⅰ、Ⅱ營養(yǎng)級，超過白洋淀全部能量的90%，隨營養(yǎng)級升高，營養(yǎng)級能量減少；2010年營養(yǎng)級Ⅰ到營養(yǎng)級Ⅱ以及營養(yǎng)級Ⅳ到營養(yǎng)級Ⅴ的能量轉(zhuǎn)換效率均高于2018年，2018年期間營養(yǎng)級Ⅱ到營養(yǎng)級Ⅲ以及營養(yǎng)級Ⅲ到營養(yǎng)級Ⅳ的轉(zhuǎn)換效率高于2010年，2018年白洋淀生態(tài)系統(tǒng)的總轉(zhuǎn)換效率為7.587%，略高于2010年總轉(zhuǎn)換效率7.576%。對比國內(nèi)相似湖泊生態(tài)系統(tǒng)發(fā)現(xiàn)，白洋淀生態(tài)系統(tǒng)總轉(zhuǎn)換效率高于太湖生態(tài)系統(tǒng)[43]的4.1%和千島湖生態(tài)系統(tǒng)[44]的6.8%，低于淀山湖生態(tài)系統(tǒng)[32]的11.7%，而湖泊生態(tài)系統(tǒng)普遍低于海洋生態(tài)系統(tǒng)總轉(zhuǎn)換效率的平均值10%。這可能因為，湖泊生態(tài)系統(tǒng)相對于海洋生態(tài)系統(tǒng)而言，水生生物種類較少，物種間捕食與被捕食的關(guān)系簡單，從而導致湖泊生態(tài)系統(tǒng)總轉(zhuǎn)化效率偏低。

表4 2010年和2018年白洋淀生態(tài)系統(tǒng)各營養(yǎng)級間傳輸效率Table 4 Transmission efficiency among nutrient levels of Baiyangdian ecosystem in 2010 and 2018

圖2 2010年和2018年白洋淀生態(tài)系統(tǒng)營養(yǎng)層級間能流效率示意圖注：P表示初級生產(chǎn)者，D表示碎屑。Fig. 2 The schematic diagram of energy flow efficiency between nutrition levels of Baiyangdian ecosystem in 2010 and 2018Note: P means primary producer; D means detritus.

混合營養(yǎng)關(guān)系(mixed nutritional effects)通常被用來描述Ecopath模型不同功能組間的相互影響關(guān)系，2010年和2018年白洋淀生態(tài)系統(tǒng)混合營養(yǎng)效應如圖3所示?？梢钥闯?，同一功能組生物之間為負影響，捕食者與被捕食功能群之間存在一定的相互關(guān)系，捕撈生產(chǎn)會對白洋淀生態(tài)系統(tǒng)中魚類有負影響。例如，作為捕食者，2010年白洋淀肉食性魚類功能群對雜食性魚類功能群有明顯負影響，2018年肉食性魚類功能群對雜食性魚類功能群和濾食性魚類功能群有明顯負影響，對其他功能組負效應不明顯。2010年和2018年Ecopath模型中肉食性魚類功能群幾乎可以捕食所有魚蝦，其對雜食性魚類功能群以及濾食性魚類功能群的捕食可以減輕其他生物被捕食的壓力，而2018年濾食性魚類功能群生物量減少造成肉食性魚類功能群對濾食性魚類功能群捕食壓力增加。在2010年和2018年白洋淀生態(tài)系統(tǒng)模型中，作為被捕食者，浮游植物、底棲藻類對大部分功能群有正影響。

圖3 2010年和2018年白洋淀生態(tài)系統(tǒng)混合營養(yǎng)效應Fig. 3 The mixed nutritional effects of Baiyangdian ecosystem in 2010 and 2018

3.3 生態(tài)系統(tǒng)總體特征

表5為2010年和2018年白洋淀生態(tài)系統(tǒng)總體特征參數(shù)對照表。Ecopath Pedigree指數(shù)(Ecopath Pedigree index)被用來評價數(shù)據(jù)和模型整體的質(zhì)量，Coll等[45]指出，在多個不同的生態(tài)系統(tǒng)所建立的Ecopath模型中，Ecopath Pedigree指數(shù)處于0.164～0.676范圍內(nèi)。2010年和2018年白洋淀生態(tài)系統(tǒng)的Ecopath Pedigree指數(shù)分別為0.500和0.432，這介于0.164～0.676范圍內(nèi)且處于較高水平，該模型研究主要基于白洋淀區(qū)域調(diào)查數(shù)據(jù)，數(shù)據(jù)可信度較高。

表5 2010年和2018年白洋淀生態(tài)系統(tǒng)Ecopath模型總體特征參數(shù)對照表Table 5 Comparison of statistics for Baiyangdian ecosystem Ecopath model in 2010 and 2018

2018年生態(tài)系統(tǒng)凈生產(chǎn)量和總流量達到4 941.52 t·km-2·a-1和24 958.28 t·km-2·a-1，分別為2010年的76%和1.58倍。通常而言，系統(tǒng)總流量與系統(tǒng)規(guī)模成正比，由該地區(qū)較早時期生態(tài)系統(tǒng)的參數(shù)[31]可知，1980—2005年間白洋淀系統(tǒng)總流量由15 429.40 t·km-2·a-1上升至57 015.52 t·km-2·a-1后下降至39 985.24 t·km-2·a-1，因此，1980—2010年間白洋淀生態(tài)系統(tǒng)規(guī)模先上升后下降，2010—2018年間生態(tài)系統(tǒng)規(guī)模呈現(xiàn)上升趨勢；2018年流向碎屑總量為2010年的1.5倍，流向碎屑所占百分比由2010年的41.58%減少到2018年的34.53%，且1980—2005年間白洋淀系統(tǒng)流向碎屑所占百分比由35.91%增加到44.54%后減少至42.87%，生態(tài)系統(tǒng)中未被利用的能量占比呈現(xiàn)先上升后下降趨勢，因此1980—2010年間生態(tài)系統(tǒng)成熟度下降，2010—2018年間成熟度升高。

總初級生產(chǎn)量/呼吸量的比值(TPP/TR)用來描述系統(tǒng)成熟程度，TPP/TR由2010年的9.56減少到2018年的1.51，系統(tǒng)成熟度提高。其中，總初級生產(chǎn)量降低24%，而總呼吸量大幅度增加，導致TPP/TR減小，導致這一變化的原因可能是多方面的。對于總初級生產(chǎn)量，沉水植物、底棲藻類和碎屑量增加19%，但底棲藻類生物量減少47%，所以導致總初級生產(chǎn)量降低；對于總呼吸量，2018年生態(tài)系統(tǒng)中消費者生物量增加了86.74%，其中，肉食性魚類、甲殼類、軟體動物和浮游動物均有明顯增加，因此呼吸作用明顯增強。與2010年相比，2018年后白洋淀水體污染程度有所減輕。根據(jù)本課題組采樣分析的結(jié)果可知，與2010年相比，2018年浮游藻類明顯減少，而且根據(jù)我們基于底棲動物完整性指數(shù)(B-IBI)健康等級評價結(jié)果可知，2018年白洋淀56%的水體等于或優(yōu)于亞健康狀態(tài)，33.3%為一般健康，其余處于微病態(tài)；而2010年12個采樣點中有33%的樣點水體等于或優(yōu)于亞健康狀態(tài)，50%一般，16%微病態(tài)，可見，白洋淀水體健康狀態(tài)有所好轉(zhuǎn)。另外，劉鑫等[46]的研究也表明，與2010年相比，2018年后白洋淀水體污染程度有所減輕，白洋淀水體治理效果顯著。李華等[31]研究發(fā)現(xiàn)，水體的營養(yǎng)條件變化會顯著改變淺水湖泊中底棲—浮游食物網(wǎng)的結(jié)構(gòu)和功能，當水體富營養(yǎng)化趨于嚴重時，容易形成浮游藻類爆發(fā)，能流路徑從以底棲—沉水植物路徑為主轉(zhuǎn)變?yōu)橐愿∮沃参锫窂綖橹?，生態(tài)系統(tǒng)從成熟狀態(tài)逆向發(fā)展；實施了生態(tài)恢復工程后，總氮(TN)和總磷(TP)濃度下降，生態(tài)系統(tǒng)狀況有向成熟狀態(tài)發(fā)展的趨勢。白洋淀實施了一系列生態(tài)恢復工程，如外源污染管控、內(nèi)源污染治理和生態(tài)補水等，2010—2018年間白洋淀水質(zhì)好轉(zhuǎn)，因此TPP/TR下降，白洋淀生態(tài)系統(tǒng)向順向演替發(fā)展。另外，陳作志和邱永松[47]研究發(fā)現(xiàn)，過度捕撈會導致大型捕食者缺乏，使小型魚類和無脊椎動物等種類占據(jù)優(yōu)勢地位，生態(tài)系統(tǒng)成熟度呈降低趨勢，漁業(yè)資源呈小型化、低值化趨勢。因此，引起白洋淀生態(tài)系統(tǒng)變化的另一方面原因，是白洋淀實施了漁業(yè)管制措施，如長期“禁漁”等，使得白洋淀捕撈強度有所下降，生態(tài)系統(tǒng)多樣性增加，漁業(yè)資源向大型化趨勢發(fā)展。且2018年采樣中發(fā)現(xiàn)魚蝦總生物量高于2010年，高營養(yǎng)級生物肉食性魚類和雜食性魚類比2010年高出424%和21%，也映證了白洋淀捕撈強度有所下降。此外，由于不同年份采樣點位數(shù)量不完全一致，樣品分析存在一定的誤差，可能導致食物網(wǎng)模型輸出結(jié)果的數(shù)值具有不確定性，但是因大部分采樣點和分析方法是一致的，且數(shù)據(jù)的來源均準確可靠，因此所得數(shù)據(jù)的變化趨勢是可信的。

連接指數(shù)(CI)和系統(tǒng)雜食指數(shù)(SOI)是Ecopath模型中描述系統(tǒng)穩(wěn)定性的參數(shù)，反映了系統(tǒng)內(nèi)部聯(lián)系的復雜程度，CI和SOI越高，系統(tǒng)越穩(wěn)定。Finn’s循環(huán)指數(shù)(FCI)是系統(tǒng)中循環(huán)流量與總流量的比值，表明有機物流轉(zhuǎn)的速度。由表5可知，白洋淀生態(tài)系統(tǒng)的CI和SOI由2010年的0.233、0.065分別增加為2018年的0.273、0.097，說明系統(tǒng)成熟度增加，與2010年相比，F(xiàn)CI增加，食物鏈變長，功能組間影響增大。對比同地區(qū)早期食物網(wǎng)模型[31]發(fā)現(xiàn)，1980—2005年間CI由0.29降至0.27，SOI由0.03升至0.04，F(xiàn)CI由6.88先降低后升高至1.22。整體而言，1980—2010年白洋淀生態(tài)系統(tǒng)成熟度和穩(wěn)定性均下降，白洋淀湖泊由草型湖泊逆向演替為藻型湖泊，2010—2018年間成熟度和穩(wěn)定性有一定程度提高。

本研究分別利用2018年實地調(diào)研數(shù)據(jù)和2010年文獻調(diào)研數(shù)據(jù)，建立了白洋淀不同時期生態(tài)系統(tǒng)的Ecopath模型，量化并對比分析了白洋淀生態(tài)系統(tǒng)結(jié)構(gòu)和能量流動，得出如下結(jié)論。

(1)2018年和2010年白洋淀生態(tài)系統(tǒng)各營養(yǎng)級之間能量隨營養(yǎng)級增大顯著減小，整體呈現(xiàn)金字塔狀；2018年白洋淀生態(tài)系統(tǒng)生物量、生產(chǎn)量、流向碎屑量以及總流量均高于2010年，生態(tài)系統(tǒng)規(guī)模增大。2018年白洋淀生態(tài)系統(tǒng)能量利用率雖然比2010年的高，但仍有大部分能量未被充分利用。

(2) 2010年和2018年白洋淀生態(tài)系統(tǒng)均存在牧食鏈和腐食鏈，由生態(tài)系統(tǒng)食物網(wǎng)結(jié)構(gòu)圖可知，其食物網(wǎng)絡相對簡單，一個物種數(shù)量的變動或消失都可能危及到整個系統(tǒng)的食物網(wǎng)結(jié)構(gòu)。2010年和2018年2個時期能量分布主要集中于第Ⅰ、Ⅱ營養(yǎng)級，2018年白洋淀生態(tài)系統(tǒng)的功能組生態(tài)營養(yǎng)轉(zhuǎn)換率略有提高。

(3)對比2010年和2018年白洋淀生態(tài)系統(tǒng)模型發(fā)現(xiàn)，2個時期的白洋淀生態(tài)系統(tǒng)處于不成熟的發(fā)育期，但從2010—2018年期間，白洋淀生態(tài)系統(tǒng)從發(fā)展階段轉(zhuǎn)向成熟階段，生態(tài)系統(tǒng)的規(guī)模增加，且生態(tài)系統(tǒng)內(nèi)部的網(wǎng)狀結(jié)構(gòu)有所增強。采取合理有效的保護措施，科學治理白洋淀污染，才能解決白洋淀污染問題，全面提升白洋淀的生態(tài)環(huán)境質(zhì)量。