王麗媛　楊麗偉　羅文磊　陳詩越　陳影影

摘要：東魚河作為南四湖最長的入湖河流，其水質(zhì)狀況對于南四湖的水環(huán)境質(zhì)量具有重要影響。為分析東魚河搖蚊群落結(jié)構(gòu)與水體富營養(yǎng)化的關(guān)系，于2013年5月對東魚河8個采樣點(diǎn)表層松散2 cm底泥中的搖蚊幼蟲頭殼進(jìn)行了分析。東魚河底泥中共鑒定出30個搖蚊屬種，主要屬種由Glyptotendipes pallens-type（20.42%）、Cricotopus sylvestris-type（18.23%）、Chironomus plumosus-type（15.24%）、Tanytarsus mendax-type（7.35%）、Dicrotendipes nervosus-type（7.12%）、Chironomus anthracinus-type（6.23%）、Polypedilum nubifer-type（5.38%）組成，共占所有搖蚊含量的79.97%。其中Glyptotendipes pallens-type、Cricotopus sylvestris-type和Chironomus plumosus-type是東魚河的優(yōu)勢種，占全流域的53.89%。基于生物多樣性指數(shù)及主成分分析對東魚河水質(zhì)進(jìn)行分析，結(jié)果表明，東魚河水體處于輕度-中度污染等級;冗余分析得到NH3-N是影響東魚河搖蚊組合分布的主要水環(huán)境因子，且不同搖蚊屬種對NH3-N的響應(yīng)不同。此外，結(jié)合搖蚊組合生態(tài)特性綜合分析得，中游、中下游污染較嚴(yán)重，是東魚河水污染治理的重點(diǎn)區(qū)域。

關(guān)鍵詞：搖蚊;生物學(xué)評價;冗余分析;水環(huán)境質(zhì)量;東魚河

中圖分類號：X82文獻(xiàn)標(biāo)志碼：A

Abstract：As the longest river discharging into the Nansi Lake，the Dongyu River plays an important role on the water quality of the Nansi Lake.To better understand the relationship between the community structure of chironomid and eutrophication of the Dongyu River，the head shells of larvae preserved in the upper 2-cm-thick loose bottom sediments were investigated with eight samples collected in May 2013.A total of 30 Chironomus species were identified in the Dongyu River surface sediments.The main genus were composed of Glyptotendipes pallens-type （20.42%），Cricotopus sylvestris-type （18.23%），Chironomus plumosus-type （15.24%），Tanytarsus mendax-type （7.35%），Dicrotendipes nervosus-type （7.12%），Chironomus anthracinus-type （6.23%），and Polypedilum nubifer-type （5.38%），accounting for 79.97% of the total counts.Based on the biodiversity index and principal component of the Dongyu River，the water body was in a mild to moderate pollution level.Redundancy analysis showed that NH3-N was the main water environmental factor resulting in the distribution of the chironomid in the Dongyu River，and the responses of different species of Chironomid to NH3-N were varied. It can also be inferred that the Dongyu River is a medium-eutrophic river，according to the ecological characteristics of the chironomid combination，with serious pollution in the middle and lower reaches，which are the key areas for water pollution control.

Key words：Chironomidea;biological evaluation;redundancy analysis;water quality;Dongyu River

南四湖是南水北調(diào)東線重要的調(diào)蓄湖泊，東魚河作為南四湖最大的入湖河流，其水質(zhì)狀況不僅關(guān)系到河流沿線居民的生產(chǎn)生活和經(jīng)濟(jì)發(fā)展，同時也關(guān)系到到南水北調(diào)東線輸水工程的水質(zhì)安全[1-2]，對入南四湖河口水質(zhì)有著重要影響。東魚河作為一條人工河，橋、涵、閘建筑物108座，河流受人工控制較強(qiáng)，東魚河水質(zhì)主要受到沿岸居民生產(chǎn)生活影響。武周虎等[3]就2010-2011年期間的南四湖的35個入湖河口水質(zhì)監(jiān)測數(shù)據(jù)分析其等標(biāo)污染物負(fù)荷，東魚河負(fù)荷比為7.0%，是南四湖10條主要排污河流之一，也成為南四湖治污的重點(diǎn)。東魚河作為國家水資源監(jiān)控能力建立淮河流域二期項(xiàng)目之南四湖19條入湖河流水量監(jiān)控站點(diǎn)之一[4]，其流域范圍內(nèi)自然條件惡劣，降水量小于蒸發(fā)量，更增大了污染物擴(kuò)散及轉(zhuǎn)化的阻力，加劇水質(zhì)安全的威脅。近些年在人類活動的干擾下，水質(zhì)狀況越來越差，徑流量逐漸降低[5]，隨著當(dāng)?shù)爻擎?zhèn)化進(jìn)程加快、城鎮(zhèn)人口數(shù)量增加，城市生產(chǎn)生活用水量大幅度增加，加快東魚河水環(huán)境治理具有現(xiàn)實(shí)意義。

2.2.3 主成分分析

運(yùn)用SPSS軟件將環(huán)境因子數(shù)據(jù)進(jìn)行主成分分析（PCA）。PCA分析是它可以將多維的水環(huán)境因子納進(jìn)統(tǒng)一系統(tǒng)中進(jìn)行定量化研究[24]，然后利用主成分分析建立水環(huán)境評價模型，可以保證在原始數(shù)據(jù)信息處于最小損失的情況下，用少量變量代替多維變量來簡化數(shù)據(jù)結(jié)構(gòu)[25]，作為多指標(biāo)分析，從眾多環(huán)境變量中剔除相關(guān)因子，篩選出主要的少數(shù)獨(dú)立綜合環(huán)境因子，而且篩選出來的這些因子能夠?qū)ρ芯拷Y(jié)果做出充分及合理的解釋[26]。根據(jù)所計(jì)算的觀測點(diǎn)水質(zhì)污染綜合得分，將水質(zhì)污染程度做定量化描述，得分越大，表示水環(huán)境污染程度越嚴(yán)重，進(jìn)而對不同觀測點(diǎn)水質(zhì)狀況進(jìn)行評價分級[27]。本文采用的水環(huán)境指標(biāo)有SAL、DO、NH3-N、NO3-N、NO2-N、SOP、STP、TP、TOC、Hg和As。

3 結(jié)果與分析

3.1 東魚河搖蚊群落特征

東魚河8個采樣點(diǎn)采集1 572.5個搖蚊幼蟲頭殼，分別隸屬于3亞科20屬，共30種。其中，DYH-4搖蚊種類最多有19種，DYH-2和DYH-8最少均為11種，DYH-1、DYH-3、DYH-5、DYH-6、DYH-7分別為13、18、16、17、17種，其中豐度>2%的有22個，詳見圖2。東魚河搖蚊幼蟲頭殼主要屬種由Glyptotendipes pallens-type（20.42%）、Cricotopus sylvestris-type（18.23%）、Chironomus plumosus-type（15.24%）、Tanytarsus mendax-type（7.35%）、Dicrotendipes nervosus-type（7.12%）、Chironomus anthracinus-type（6.23%）、Polypedilum nubifer-type（5.38%）組成，含量共占所有搖蚊亞化石屬種的79.97%。其中Chironomus plumosus-type、Glyptotendipes pallens-type和Cricotopus sylvestris-type是東魚河的優(yōu)勢種，占全流域的53.89%，幾乎全流域分布。

東魚河8個采樣點(diǎn)搖蚊的平均密度為655 ind./m2，DYH-1的搖蚊密度最低，為80 ind./m2，DYH-4的搖蚊密度最高，達(dá)到峰值，為1 447 ind./m2。由圖3可以看出，DYH-1、DYH-2、DYH-3和DYH-8的搖蚊密度低，其余采樣點(diǎn)搖蚊密度高，不同采樣點(diǎn)搖蚊密度差異很大，整體上的搖蚊群落密度沿流域增加。并且搖蚊密度流域分布不均勻，反映出各采樣點(diǎn)水質(zhì)狀況差別很大。[HJ2.1mm]

3.2 東魚河水質(zhì)評價

東魚河8個采樣點(diǎn)搖蚊Margalef指數(shù)、Shannon-Wiener指數(shù)、Pielou指數(shù)的波動范圍依次為2.43～3.54、1.35～1.99、0.53～0.75，東魚河搖蚊平均值依次為2.76、1.66和0.61。根據(jù)表2所示，Margalef指數(shù)的環(huán)境質(zhì)量評價結(jié)果為25%的輕污染等級，75%中污染等級;Shannon-Wiener指數(shù)環(huán)境質(zhì)量評價結(jié)果為100%的中污染等級;Pielou指數(shù)的環(huán)境質(zhì)量評價結(jié)果為87.5%的輕污染等級，12.5%的中污染等級。Margalef指數(shù)和Shannon-Wiener指數(shù)反映采樣點(diǎn)DYH-2、DYH-4、DYH-5、DYH-6、DYH-7、DYH-8的水質(zhì)狀況都為中度污染。采樣點(diǎn)DYH-7三個評價指標(biāo)反映的水質(zhì)狀況都為中污染等級，其余采樣點(diǎn)三個生物 指標(biāo)評價結(jié)果不完全相同，都處于輕度-中度污染狀態(tài)，這表明東魚河不同流域的水體受污染程度不同，反映的水質(zhì)狀況存在一定的差異，整體屬于“輕污染-中污染”等級。

3.3 東魚河水環(huán)境質(zhì)量綜合評價

DCA分析顯示，各軸中最大梯度長度為1.33（SD<2），表明搖蚊對水環(huán)境因子的響應(yīng)為線性模型，因此采用RDA分析搖蚊與水環(huán)境因子間的關(guān)系較合適。RDA揭示了搖蚊群落結(jié)構(gòu)與水環(huán)境因子之間的相互關(guān)系，通過非限定性蒙特卡洛置換（p<0.05，n=499）檢驗(yàn)，剔除共線性參數(shù)（VIF≥5），篩選出NH3-N（F=5.012，P=0.02）是影響搖蚊組成與分布的顯著因子，所解釋的搖蚊變率為45.5%，說明該區(qū)域的環(huán)境變量與搖蚊組分之間具有極強(qiáng)的相關(guān)性。

由圖4（a）RDA物種-環(huán)境排序圖可知，NH3-N的含量與第一軸關(guān)系密切，沿第一軸的負(fù)向直線增加，尤其是物種Microchironomus tabarui-type、Cladopelma laccophila-type、Parachironomus varus-type、Dicrotendipes nervosus和優(yōu)勢種Chironomus plumosus-type、Glyptotendipes pallens-type和Cricotopus sylvestris-type等幾乎都隨著NH3-N濃度的增加呈線性增多;而Paratanytarsus penicillatus-type、Cricotopus bicinctus-type、Tanpus Meigen、Cladotanytarsus mancus-type和Cricotopus trifascia-type等的分布受低濃度的NH3-N的影響較大。其他豐度較高的搖蚊屬種與NH3-N的相關(guān)度不高，表明東魚河水環(huán)境中不同水質(zhì)要素對搖蚊屬種的響應(yīng)差異顯著，且不同物種對不同水環(huán)境要素污染的響應(yīng)差異顯著。

圖4（b）RDA樣點(diǎn)-環(huán)境排序圖可知，東魚河中游、中下游采樣帶點(diǎn)（DYH-4、DYH-5、DYH-6、DYH-7）散落在NH3-N濃度高的區(qū)域。中游、中下游采樣點(diǎn)搖蚊密度大，屬種Chironomus anthracinus-type、Glyptotendipes pallens-type、Propsilocerus akamusi-type等在此區(qū)域分布較多，其中Propsilocerus akamusi-type是Ⅳ水體（中污染水體）的指示物種，并且Margalef指數(shù)、Shannon-Wiener指數(shù)的水質(zhì)評價結(jié)果顯示該區(qū)域水體屬于中度污染等級，污染較嚴(yán)重;東魚河上游及入湖口處采樣點(diǎn)（DYH-1、DYH-2、DYH-3、DYH-8）散落在NH3-N濃度低的區(qū)域，該區(qū)域搖蚊密度小，主要搖蚊屬種為Tanytarsus mendax-type、Polypedilum nubifer-type、Procladius Skuse，其中Tanytarsus mendax-type是Ⅲ類水體（輕污染水）的指示物種，而Polypedilum nubifer-type是在Ⅲ類水體及Ⅳ類水體中豐度最大，Procladius Skuse為世界廣適種，具有耐酸性。因此東魚河上游水體及入湖口處水體污染程度較輕，中游及中下游水體污染較嚴(yán)重。

根據(jù)東魚河8個觀測點(diǎn)的水質(zhì)指標(biāo)，采用PCA分析（累計(jì)貢獻(xiàn)率大于85%，共選出3個主成分）所得結(jié)果見表3。東魚河上游觀測點(diǎn)DYH-1、DYH-2、DYH-3、DYH-8的水質(zhì)處于輕污染，觀測點(diǎn)DYH-4、DYH-5、DYH-7的水質(zhì)處于中度污染，觀測點(diǎn)DYH-6的水質(zhì)處于重污染。即東魚河上游水質(zhì)狀況較好，中游及下游水質(zhì)污染較嚴(yán)重。主成分分析結(jié)果與生物學(xué)評價結(jié)果較為一致，東魚河整體處于輕度-中度污染等級。

4 討論

搖蚊屬種的組合分布受多種水環(huán)境要素的影響，如溫度、PH、水深、食物、溶解氧等的變化都會是搖蚊組合發(fā)生變化。在淺水流域中，水質(zhì)變化會在一定程度上影響搖蚊組合變化分布，因此搖蚊群落組合變化能很好的反應(yīng)水環(huán)境污染狀況[28]。研究表明Chironomus plumosus-type和Propsilocerus akamusi-type為水體富營養(yǎng)化的指示生物[29]，曹艷敏等[30]根據(jù)太湖梅梁灣T0950巖芯搖蚊亞化石組合變化，揭示了太湖梅梁灣湖區(qū)水體從1940年至今營養(yǎng)水平的演化。本文通過非限定性蒙特卡洛置換檢驗(yàn)，僅篩選出NH3-N（F=5.012，P=0.02）是東魚河搖蚊組合分布的顯著影響因子。根據(jù)張菊等[31]的研究結(jié)果，東魚河主要水環(huán)境污染因子為NH3-N、TP和Hg，張恩樓等[15]已證明TP是影響搖蚊群落結(jié)構(gòu)組成與分布的重要水環(huán)境因子，說明東魚河水環(huán)境中TP對搖蚊群落組合分布的影響小于NH3-N。通過RDA分析得到，Microchironomus tabarui-type、Parachironomus varus-type、Cladopelma laccophila-type、Parachironomus varus-type、Dicrotendipes nervosus為NH3-N的耐污染屬種，但其他豐度較高的搖蚊屬種與NH3-N的相關(guān)度不高，表明東魚河水環(huán)境中不同水質(zhì)要素對搖蚊屬種的響應(yīng)差異顯著，且不同物種對不同水環(huán)境要素污染的響應(yīng)差異顯著。

東魚河主要搖蚊屬種中Chironomus plumosus-type為Ⅴ類水體（重污染水體）的指示物種，主要分布于中游、中下游流域;Cricotopus sylvestris-type是Ⅳ類水體指示物種，下游分布較多;Polypedilum nubifer-type、Tanytarsus mendax-type為Ⅳ類水體的指示物種，Polypedilum nubifer-type、Tanytarsus mendax-type主要分布在東魚河中游;Glyptotendipes pallens-type為嗜酸性屬種，分布于中游、中下游河段;Dicrotendipes nervosus-type為Ⅳ類水體指示物種，在中游、中下游流域分布較多，Tanytarsus mendax-type是Ⅲ類水體的指示物種，而Polypedilum nubifer-type是在Ⅲ類水體及Ⅳ類水體中豐度最大。結(jié)合三種生物指標(biāo)評價結(jié)果及主成分分析結(jié)果，東魚河水體處于輕度-中度污染等級，中游、中下游水環(huán)境污染較嚴(yán)重，域范圍內(nèi)差距很大。東魚河中游、中下游水質(zhì)較差，污染主要來源于沿岸農(nóng)業(yè)區(qū)帶來的農(nóng)業(yè)面源污染及生活污水的排放;僅在1997年全區(qū)域的化肥使用量就高達(dá)100多萬t，農(nóng)藥使用量達(dá)2 500萬t[32]，大量的農(nóng)業(yè)污染物殘留，難以根治，使得沿線水環(huán)境污染難以在短時間內(nèi)治理。東魚河上游水體污染較輕，上游主要為工業(yè)廢水排放引起的水環(huán)境污染，如菏澤、定陶縣、成武縣等設(shè)有大量工廠，但在人為干擾下，加強(qiáng)企業(yè)對污水的治理，使得水體得到相應(yīng)改善?？偠灾祟惢顒邮窃斐蓶|魚河水環(huán)境污染的主要原因，隨著城市化進(jìn)程的加快，東魚河水環(huán)境勢必面臨著更加嚴(yán)峻的挑戰(zhàn)，加強(qiáng)東魚河水環(huán)境治理不僅關(guān)系到南水北調(diào)東線水質(zhì)安全，更關(guān)系到沿線居民的生產(chǎn)生活。

本文基于搖蚊反映的東魚河水質(zhì)狀況與張菊等[30]基于東魚河水環(huán)境數(shù)據(jù)所反映的水質(zhì)狀況相吻合，其中NH3-N含量在中游上升，在下游含量下降，且中游水體污染較嚴(yán)重，這與搖蚊組合分布特征相呼應(yīng)，中游Ⅳ類水體指示屬種密集，水環(huán)境污染較嚴(yán)重，上游Ⅲ類水體指示屬種密集，證明根據(jù)搖蚊群落結(jié)構(gòu)所反映的水環(huán)境狀況結(jié)果準(zhǔn)確可靠。根據(jù)張偉分析的東魚河上游東明縣黃軍營村及入湖口處的魚臺縣西姚村的兩個水質(zhì)采樣點(diǎn)2000-2008年的水質(zhì)狀況，僅有黃軍營村水體中的NH3-N含量在2005-2007年間超標(biāo)準(zhǔn)值，其研究結(jié)果表明東魚河水體在2006-2008年為輕污染河流[33]。其研究設(shè)立采樣點(diǎn)較少，僅監(jiān)測了東魚河上游和入湖口兩個采樣點(diǎn)的水質(zhì)狀況，不能反映東魚河體整體的水質(zhì)情況，其研究只能說明東魚河上游及入湖口處水質(zhì)狀況較好。

5 結(jié)論

東魚河搖蚊的優(yōu)勢屬種由Chironomus plumosus-type、Glyptotendipes pallens-type和Cricotopus sylvestris-type組成，占全流域的53.89%。RDA分析顯示NH3-N是影響東魚河搖蚊群落組合分布的顯著因子，共解釋了45.5%的搖蚊變率，屬種Microchironomus tabarui-type、Parachironomus varus-type、Cladopelma laccophila-type分布在NH3-N污染較大的區(qū)域，它們可視為水體NH3-N污染的指示生物。水環(huán)境生物學(xué)評價結(jié)果顯示東魚河不同流域受污染程度不同，根據(jù)搖蚊群落結(jié)構(gòu)分布特征，東魚河整體處于輕度-中度污染狀態(tài)，中游、中下游水體污染較嚴(yán)重，為治理的重點(diǎn)區(qū)域。

參考文獻(xiàn)（References）：

[1] 翟曉萌.南四湖大型底棲動物研究及水質(zhì)評價[D].濟(jì)南：山東師范大學(xué)，2013.（ZHAI X M.Large benthic research and water quality evaluation in Nansi Lake[D].Jinan：Shandong Normal University，2013.（in Chinese））

[2] 羅輝，周建仁，郭忠.南水北調(diào)對南四湖水環(huán)境影響分析與評估[J].河海大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，33（1）：63-67.（LUO H，ZHOU J R，GUO Z.Analysis and assessment of the impact of South-to-North Water Diversion Project on water environment of Nansihu Lake[J].Journal of Hohai University （Natural Science Edition），2005，33（1）：63-67.（in Chinese）） DOI：10.3321/j.issn：1000-1980.2005.01.016.

[3] 武周虎，任杰，張曉波，等.南四湖入湖河口水質(zhì)綜合分析與改善效果評估[J].環(huán)境污染與防治，2012，34（07）：80-84.（WU Z H，REN J，ZHANG X B，et al.Comprehensive analysis and improvement of water quality in the estuary of Nansi Lake[J].Environmental Pollution and Prevention Journal，2012，34（07）：80-84.（in Chinese）） DOI：10.3969/j.issn.1001-3865.2012.07.017.

[4] 韓俊麗.東魚河截污導(dǎo)流工程水質(zhì)保障體系研究[J].山東水利，2014，16（10）：30-31.（HAN J L.Water quality assurance system of East Fish River sewage Interception diversion project[J].Shandong Water Resources，2014，16（10）：30-31.（in Chinese）） DOI：10.3969/j.issn.1009-6159.2014.10.015.

[5] 張翠榮.氣候變化和人類活動對東魚河流域徑流影響的探析[J].水資源開發(fā)與管理，2017，13（7）：66-69.（ZHANG C R.Analysis on influence of climate change and human activity on runoff in Dongyu River basin[J].Water resources development and Management，2017，13（7）：66-69.（in Chinese）） DOI：10.16616/j.cnki.10-1326/tv.2017.07.019.

[6] ARMITAGE P D，CRANSTON P S，PINDER L C V.The Chironomidae：the biology and ecology of non-biting midges[J].Journal of the North American Benthological Society，1995，14（3）：611-738.DOI：10.2307/5810.

[7] 王元元，宋進(jìn)喜，蔣衛(wèi)威，等.搖蚊幼蟲與霍甫水絲蚓擾動對河流沉積物中氮、磷、氧釋放的影響[J].環(huán)境科學(xué)學(xué)報(bào)，2015，35（8）：2504-2511.（WANG Y Y，SONG J X，JIANG W W，et al.Effect of bioturbation of Chironomid larvae and Limnodrilus hoffmeiteri on the release of nitrogen，oxygen and phosphate in the sediments from a river[J].Acta Scientiae Circumstantiae，2015，35（8）：2504-2511.（in Chinese）） DOI：10.13671/j.hjkxxb.2014.1036.

[8] 商景閣，張路，張波，等.中國長足搖蚊 （Tanypus chinensis） 幼蟲底棲擾動對沉積物溶解氧特征及反硝化的影響[J].湖泊科學(xué)，2010，22（5）：708-713.（SHANG J G，ZHANG L，ZHANG B，et al.Bioturbation effect of Tanypus chinensis larvae on denitrification rate and process in sediments[J].Journal of Lake Sciences，2010，22（5）：708-713.（in Chinese）） DOI：10.18307/2010.0512.

[9] PETER B，DANIELA D，RADOVAN P，et al.Relict chironomid communities surviving in the coldest High Tatra Mountain lakes confirmed by a palaeolimnological survey[J].Biologia，2017，72（8）：965-969.DOI：10.1515/biolog-2017-0102.

[10] [ZK（#]HIRABAYASHI K，YOSHIZAWA K，YOSHIDA N，et al.Progress of eutrophication and change of chironomid fauna in Lake Yamanakako，Japen[J].Limnology，2004，5（1）：47-53.DOI：10.1007/s10201-003-0113-2.

[11] 楊麗偉，張菊，鄧煥廣，等.東平湖底泥重金屬污染及其對搖蚊分布的影響[J].應(yīng)用與環(huán)境生物學(xué)報(bào)，2016，22（4）：680-688.（YANG L W，ZHANG J，DENG H G，et al.Heavy metal pollution and its influence on the distribution of chironomids in the sediments of Lake Dongping，Shandong，China[J].Chinese Journal of Applied & Environmental Biology，2016，22 （4）：680-688.（in Chinese）） DOI：10.3724/SP.J.1145.2015.11023.

[12] ZHANG E L，BEDORRD A，JONES R，et al.A subfossil chironomid-total phosphorus inference model for lakes in the middle and lower reaches of the Yangtze River[J].Chinese Science Bulletin，2006，51（17）：2125-2132.DOI：10.1007/s11434-006-2062-8.

[13] 霍堂斌，劉曼紅，姜作發(fā)，等.松花江干流大型底棲動物群落結(jié)構(gòu)與水質(zhì)生物評價[J].應(yīng)用生態(tài)學(xué)報(bào)，2012，23（1）：247-254.（HUO T B，LIU M H，JIANG Z F，et al.Macrobenthos community structure of macrobenthos and bioassessment of water quality in main stream of Songhua River[J].Chinese Journal of Applied Ecology，2012，23（1）：247-254.（in Chinese）） DOI：10.13287/j.1001-9332.2012.0040.

[14] 仲嘉，王衛(wèi)民，單保慶，等.大清河搖蚊幼蟲群落結(jié)構(gòu)特征及水質(zhì)初步評價[J].中國農(nóng)學(xué)通報(bào)，2015，31（32）：106-116.（ZHONG J，WANG W M，SHAN B Q，et al.Community Structure of Chironomid Larvae and Bioassessment of water quality in Daqing River[J].Chinese Agricultural Science Bulletin，2015，31（32）：106-116.（in Chinese）） DOI：10.11924/j.issn.1000-6850.casb15060094.

[15] 張恩樓，ALAN B，RICHARD J，等.長江中下游地區(qū)典型湖泊搖蚊亞化石-湖水總磷定量模型研究[J].科學(xué)通報(bào)，2006，51（11）：1318-1325.（ZHANG E L，ALAN B，RICHARD J，et al.Typical lake Yangtze River chironomid subfossil - LAKE total phosphorus quantitative model[J].Chinese Science Bulletin，2006，51（11）：1318-1325.（in Chinese）） DOI：10.3321/j.issn：0023-074X.2006.11.013.

[16] 張永利，常家忠，葛兆生.南水北調(diào)東線山東段東魚河流域治污探索[J].山東水利，2007，9（7）：53，59.（ZHANG Y L，CHANG J Z，GE Z S.Shandong Section of Eastern South East Fish River pollution exploration[J].Shandong Water Resources，2007，9（7）：53，59.（in Chinese）） DOI：10.3969/j.issn.1009-6159.2007.07.025.

[17] BROOK S J.Late-glacial fossil midge stratigraphies （Insecta：Diptera：Chironomidae） from the Swiss Alps[J].Palaeogeography，palaeoclimatology，palaeoecology，2000，159（3）：261-279.DOI：org/10.1016/S0031-0182 （00） 00089-4.

[18] 張丹，丁愛中，林學(xué)鈺，等.河流水質(zhì)監(jiān)測和評價的生物學(xué)方法[J].北京師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2009，45（2）：200-204.（ZHANG D，DING A Z，LIN X Y，et al.Monitoring and assessment of water quality in rivers with biological methods.Journal of Beijing Normal University（Natural Science），2009，45（2）：200-204.（in Chinese）） DOI：10.3321/ j.issn：0476-0301.2009.02.022.

[19] RDEATH R.Margalef′s Index[J].Encyclopedia of Ecology，2008：2209-2210.DOI：10.1016/B978-008045405-4.00117-8.

[20] 王晶，焦燕，任一平，等.Shannon-Wiener多樣性指數(shù)兩種計(jì)算方法的比較研究[J].水產(chǎn)學(xué)報(bào)，2015，39（8）：1257-1263.（WANG J，JIAO Y，REN Y P，et al.Comparative study on two computing methods for estimating Shannon-Wiener diversity index[J].Journal of Fisheries of China，2015，39（8）：1257-1263.（in Chinese）） DOI：10.11964/jfc.20150509886.

[21] PIELOU E C.Ecological Diversity[M].New York：John Wiley，1995.

[22] 盧少俊.南四湖大型底棲動物的群落結(jié)構(gòu)及演替規(guī)律[D].濟(jì)寧：曲阜師范大學(xué)，2012.（LU S J.Community structure and succession of macrobenthos in the Nansi Lake[D].Jining：Qufu Normal University，2012.（in Chinese））

[23] TER BRAAK C J F.Canonical correspondance analysis：A new eigenvector technique for multivariate direct analysis[J].Ecology，1986，67（5）：1167-1179.DOI：10.2307/ 1938672.

[24] HUANG B R，F(xiàn)RY J D.Root anatomical physiological and Morphological responses to drought stress for Tall Fescue cultivars[J].Crop Science，1998 38（4）：1017-1022.DOI：10.2135/ cropsci1998.0011183X003800040022x.

[25] SUN L F，QI J Q，ZHAO X L，et al.The establishment of principal component analysis assessment model for drinking water quality of city resource[J].Applied Mechanics & Materials，2014，3547（ 675-677）：960-963.DOI：10.4028/www.scientific.net/AMM.675-677.960.

[26] 吉祝美，方里，張俊，等.主成分分析法在SPSS軟件中的操作及在河流水質(zhì)評價中的應(yīng)用[J].環(huán)境研究與監(jiān)測，2012，25 （4）：68-73，57.（JI Z M，F(xiàn)ANG L，ZHANG J，et al.Principal component analysis of SPSS software and its application in river water quality assessment[J].Environmental Research and Monitoring，2012，25 （4）：68-73，57.（in Chinese）） DOI：10.3969/j.issn.1674-0254.2012.04.013.

[27] 劉德林，劉賢趙.主成分分析在河流水質(zhì)綜合評價中的應(yīng)用[J].水土保持研究，2006（3）：124-125，128.（LIU D L，LIU X Z.Application of principal component analysis to the comprehensive evaluation of water quality in river[J].Research of Soil and Water Conservation，2006（3）：124-125，128.（in Chinese）） DOI：10.3969/j.issn.1005-3409.2006.03.040.

[28] LANGDON P E，RUIZ Z，BRODERSEN K L.Assessing lake eutrophication using chironomids：understanding the nature of community response in different lake types[J].Freshwater Biology，2006，51（03）：562-577.DOI：10.1111/j.1365-2427.2005.01500.x.

[29] 舒鳳月，孔亮，王賽迪，等.山東南四湖搖蚊幼蟲群落結(jié)構(gòu)及其對富營養(yǎng)化過程的指示[J].應(yīng)用與環(huán)境生物學(xué)報(bào)，2013，19（1）：141-146.（SHU F Y，K L，WANG S D，et al.Chironomid larva community and indication to eutrophication progress in Lake Nansi，Shandong，China[J].Chinese Journal of Applied and Environmental Biology，2013，19（1）：141-146.（in Chinese）） DOI：10.3724/SP.J.1145.2013.00141.

[30] 曹艷敏，張恩樓，沈吉，等.太湖梅梁灣搖蚊亞化石沉積記錄及營養(yǎng)演化研究[J].湖泊科學(xué)，2011，23（4）：549-554.（CAO Y M，ZHANG E L，SHEN J，et al.Subfossil chironomid assemblages in sediments of Meiliang Bay，Lake Taihu and the trophic reconstruction[J].Journal of Lake Sciences，2011，23（4）：549-554.（in Chinese）） DOI：10.18307/2011.0409.

[31] 張菊，魯長娟，郭娜，等.南四湖支流東魚河水質(zhì)現(xiàn)狀及健康風(fēng)險(xiǎn)評價[J].水電能源科學(xué)，2017，35 （8）：56-59.（ZHANG J，LU C J，GUO N，et al.Water quality and environmental health risk assessment of Dongyu River[J].Water Resources and Power，2017，35 （8）：56-59.（in Chinese））

[32] 張勤勛，孫輝.淺談東魚河的水質(zhì)污染及對策[J].菏澤師專學(xué)報(bào)，1999，21（2）：59-60.（ZHANG Q X，SUN H.Talking about the water pollution and countermeasures of Dongyu River[J].Journal of Heze University，1999，21（2）：59-60.（in Chinese）） DOI：10.16393/j.cnki.37-1436/z.1999.02.021.

[33] 張偉.南四湖主要入湖河流水質(zhì)分析與評價研究[D].濟(jì)南：濟(jì)南大學(xué)，2011.（ZHANG W.Analysis and evaluation on water quality for main inflow river of Nansihu Lake[D].Jinan：Jinan University，2011.（in Chinese））