張 浩，丁 森，張 遠，賈曉波，孟 偉，郭 彪

(1：河北農(nóng)業(yè)大學海洋學院，秦皇島 066003)

(2：中國環(huán)境科學研究院流域水生態(tài)保護技術(shù)研究室，北京 100012)

(3：天津市渤海水產(chǎn)研究所，天津 300457)

西遼河流域魚類生物完整性指數(shù)評價及與環(huán)境因子的關(guān)系*

張 浩1,2，丁 森2，張 遠2，賈曉波2，孟 偉2，郭 彪1,3**

(1：河北農(nóng)業(yè)大學海洋學院，秦皇島 066003)

(2：中國環(huán)境科學研究院流域水生態(tài)保護技術(shù)研究室，北京 100012)

(3：天津市渤海水產(chǎn)研究所，天津 300457)

采用魚類生物完整性指數(shù)(F-IBI)評估西遼河流域環(huán)境質(zhì)量，利用相關(guān)分析篩選了影響F-IBI的環(huán)境因子，結(jié)合局部加權(quán)回歸散點修勻法(LOWESS)估算了這些環(huán)境因子的保護限值.結(jié)果表明：20個候選參數(shù)指標中，鰍科魚類物種百分比、底層魚類物種數(shù)、耐受性魚類個體數(shù)、雜食性魚類物種數(shù)和筑巢產(chǎn)卵魚類個體百分比適合作為西遼河流域F-IBI構(gòu)建參數(shù).采用比值法統(tǒng)一參數(shù)量綱，以參考點F-IBI值的95%分位數(shù)作為健康參考值，對全流域44個采樣點進行健康評價，得到健康點位5個(11.4%)，亞健康點位9個(20.5%)，一般點位12個(27.3%)，差點位11個(25.0%)，極差點位7個(15.9%).電導率、氨氮濃度、堿度和泥沙比例與F-IBI呈顯著負相關(guān)，坡降和草地比例與F-IBI呈顯著正相關(guān).經(jīng)LOWESS和獨立樣本t檢驗發(fā)現(xiàn)，電導率、氨氮濃度、堿度、泥沙比例和草地比例等環(huán)境因子的保護限值分別為531μS/cm、0.55mg/L、4.4mmol/L、47.2%和37.0%，且限值兩側(cè)的F-IBI分值差異顯著，該結(jié)果可作為西遼河流域環(huán)境管理中魚類群落完整性保護的有益參考.

魚類生物完整性指數(shù)；局部加權(quán)回歸散點修勻法；西遼河流域；環(huán)境因子

水體的生物學評價對評估水體生態(tài)系統(tǒng)完整性和水生生物保護具有重要意義[1].生物完整性指數(shù)(index of biotic integrity, IBI)作為一種較為成熟的生物學評價方法，已得到了廣泛認可和應(yīng)用[2-5].關(guān)于IBI的研究，眾多學者較多注重于對其研究的內(nèi)容不斷改進和完善，如拓展生物類群[6-8]、完善指標體系[9]、構(gòu)建方法和健康評價標準[10]等，并對不同的水體類型進行健康評價，而關(guān)于環(huán)境因子與IBI之間的相互關(guān)系以及環(huán)境因子保護限值的研究鮮有報道.

我國關(guān)于生物學評價的研究多集中于地表徑流連續(xù)性的水體[11-15].西遼河所處區(qū)域氣候干旱少雨，多數(shù)支流常出現(xiàn)季節(jié)性斷流，此外，對西遼河流域開展的水質(zhì)生物學評價研究尚不多見，僅見高欣等[16]以大型底棲動物為研究對象建立了本區(qū)域的水質(zhì)生物評價體系.本研究嘗試構(gòu)建西遼河流域魚類生物完整性指數(shù)(fish index of biotic integrity, F-IBI)，了解流域內(nèi)各調(diào)查樣點的健康狀況，并從水化學[17]、物理生境和土地利用[18]等方面篩選出對該流域F-IBI有重要影響的環(huán)境因子并計算其保護限值，以期為該流域魚類保護和環(huán)境管理提供科學依據(jù).

1 材料與方法

1.1 研究區(qū)域概況

西遼河(41.1°～45.2°N, 116.7°～124.2°E)位于我國北方典型的農(nóng)業(yè)區(qū)與牧業(yè)區(qū)交匯過渡地帶，流域面積約為20.1×104km2[19]，其干流長為827km，流域內(nèi)大的支流有西拉木倫河、老哈河、教來河等.該流域的主要氣候特征為干旱少雨，其水分年蒸發(fā)量遠大于補給量，水土流失和草場荒漠化嚴重，區(qū)域內(nèi)大部分河流處于全年斷流或季節(jié)性斷流狀態(tài).該區(qū)域的環(huán)境壓力主要來自自身和人類干擾.

1.2 調(diào)查方法

2012年8月在西遼河流域設(shè)置44個調(diào)查斷面(圖1).對于可涉水河流使用電魚法進行樣本調(diào)查，對于不可涉水河流采用刺網(wǎng)采樣(孔徑大小為5cm)和電魚法相結(jié)合的調(diào)查方法，每個斷面的調(diào)查距離約為300m，調(diào)查時間為30min.對采集的魚類進行物種鑒定和生物學測量，并最終放歸河流.

物理生境特征指標包括坡降、蜿蜒度和底質(zhì)粒徑(巨礫、鵝卵石、小卵石、砂礫和泥沙的比例)；其中坡降和蜿蜒度用ArcGIS 9.3軟件進行分析，底質(zhì)粒徑組成用篩網(wǎng)進行現(xiàn)場分析，并將粒徑大于128mm的定義為巨礫，粒徑為64～128mm的定義為鵝卵石，粒徑為16～64mm的定義為小卵石，粒徑為4～16mm的定義為砂礫，粒徑小于4mm的定義為泥沙.

土地利用方式包括未利用地、水域、林地、建設(shè)用地、耕地和草地等.本研究分析了調(diào)查樣點河岸帶兩側(cè)各3km和樣點上、下游各3km范圍內(nèi)的土地利用情況，該數(shù)據(jù)使用ArcGIS 9.3軟件分析獲得.

1.3 數(shù)據(jù)處理

1.3.1 F-IBI評價計算 該過程包括4個基本步驟：(1) 參照點位的篩選.參照渠曉東等的標準化方法[10]進行參考點位的篩選，即根據(jù)《地表水環(huán)境質(zhì)量標準》(GB 3838-2002)和生境評分結(jié)果對采樣點進行評價：樣點水質(zhì)為Ⅱ類及以上，生境評分為140分及以上，同時人類活動強度與河岸土地利用類型兩項指標得分(≥30)均較高的為參照點；水質(zhì)為Ⅴ類及劣Ⅴ類，生境評分為120分以下，同時人類活動強度與河岸土地利用類型兩項指標得分(<25)均較低的采樣點為受損點.(2) IBI指標體系的構(gòu)建.綜合參考目前魚類IBI的相關(guān)研究報道[12,15,21-23]，并結(jié)合本次魚類調(diào)查的實際情況，從5個方面進行評價體系構(gòu)建：物種組成與豐度、耐受性和敏感性、營養(yǎng)結(jié)構(gòu)、魚類數(shù)量與健康狀況、繁殖共位體.從上述5個方面共選取對環(huán)境變化較為敏感的20個指標作為評價體系的候選指標(表1).(3) 指標的篩選.參考當前IBI計算中的指標篩選方法[21-22]，對上述指標依次進行分布范圍篩選、判別能力分析、相關(guān)性計算等.最終選取適合作為該區(qū)域F-IBI的構(gòu)建指標.其中分布范圍篩選是指候選指標在超過5%以上樣點得分不為0，否則放棄該評價指標[21]；判別能力分析是指對參考點和受損點在25%～75%分位數(shù)分布范圍(即箱體IQ值)重疊情況的比較，若箱體沒有重疊，則IQ值為3，若箱體部分重疊，但是彼此的中位值均不在對方箱體內(nèi)，則IQ值為2，若箱體重疊且一者中位值在對方箱體內(nèi)，則IQ值為1，若箱體重疊且彼此中位值均在對方箱體內(nèi)，則IQ值為0，在研究中僅對IQ值≥2的情況做進一步分析[22]；相關(guān)性計算指對通過判別能力分析的各候選指標進行Pearson 相關(guān)性分析，對于|R|>0.9的指標僅選其一[22].(4) 西遼河F-IBI的健康評價.對各構(gòu)建指標按比值法進行量綱統(tǒng)一[13,21]，各樣點F-IBI值用百分制表示.以參照點F-IBI值分布的95%分位數(shù)值作為健康評價標準，對小于95%分位數(shù)值的分布范圍進行4等分，分別劃分出亞健康、一般、差和極差的評價標準.

1.3.2 環(huán)境因子保護限值計算 環(huán)境變量同F(xiàn)-IBI進行Spearman相關(guān)性分析，對具有顯著相關(guān)性的環(huán)境變量，利用局部加權(quán)回歸散點修勻法(locally weighted regression scatterplot smoother, LOWESS)進行擬合，在曲線最為平滑的情況下(f=0.9)，將曲線明顯轉(zhuǎn)折處作為環(huán)境變量對F-IBI的保護限值[23]，并利用獨立樣本t檢驗法，對各環(huán)境因子的保護限值進行檢驗.

環(huán)境和魚類狀況等基礎(chǔ)性分析采用Excel 2007軟件完成，參考點與受損點箱體圖繪制采用Oringin 8.5軟件完成，環(huán)境因子保護限值分析和環(huán)境因子間相關(guān)性分析圖采用R 3.03語言分析軟件完成，相關(guān)性分析和獨立樣本t檢驗采用SPSS 19.0軟件完成.

2 結(jié)果

2.1 魚類組成和環(huán)境特征

本次調(diào)查共采集到魚類樣本9899尾，隸屬于5目8科25屬28種.其中，鯉形目魚類最多，有20種；鱸形目其次，有4種；鲇形目較少，有2種；刺魚目和鳉形目最少，分別只有1種.鲇、黃顙魚、鯉、寬鰭鱲和大鱗副泥鰍等物種在該流域內(nèi)分布較少；北方條鰍、達里湖高原鰍、棒花魚、洛氏鱥、鯽和麥穗魚等物種數(shù)量較多；鯽、棒花魚、麥穗魚、北方條鰍和泥鰍等物種分布較廣.通過參考相關(guān)文獻[9,12,24]和實際調(diào)查經(jīng)驗，對該流域各魚類相關(guān)生態(tài)信息進行總結(jié)(表2).

表1 西遼河流域魚類IBI候選指標列表

對西遼河流域44個調(diào)查斷面的27個水化學和物理生境等環(huán)境因子的整體狀況進行分析統(tǒng)計.在水化學方面，參照國家《地表水環(huán)境質(zhì)量標準》(GB 3838-2002)，表明水溫、DO和pH均在合理范圍內(nèi)，CODMn、NH3-N、TN和TP濃度等水質(zhì)評價所選變量的最大值均已達到劣Ⅴ類；物理生境方面，各環(huán)境因子變化幅度均較大(表3).

2.2 F-IBI評價

2.2.1 參考點位選擇 按照上述點位篩選方法，最終篩選出了W11、W24、W55、W69和W76共計5個點位為參照點；W6、W33、W39、W48、W93、W95、W99、W100和W106共計9個點位為受損點.

2.2.2 指標篩選 各候選指標都通過了分布范圍的篩選.在候選指標判別能力的篩選中，M2、M7、M10、M11和M19，共計5個指標在參照點和受損點間有顯著差異(P<0.05)，其中M2、M11和M19的IQ值為3，M7和M10的IQ值為2(圖2).對上述5個指標進行Pearson相關(guān)性檢驗(表4)，所有|r|值均小于0.9，最終確定5個指標作為構(gòu)建西遼河流域IBI評價的參數(shù)指標.

2.2.3 評分及評價 根據(jù)篩選出的構(gòu)建西遼河流域IBI的各參數(shù)指標(表5)及各自對干擾的反應(yīng)，計算出各參數(shù)指標在樣點的得分及IBI值，并確定西遼河流域魚類生物完整性評價標準(表6).

2.2.4 評價結(jié)果 根據(jù)該流域F-IBI的健康評價標準，對全流域樣點進行健康評價，結(jié)果發(fā)現(xiàn)健康點位5個，亞健康點位9個，一般點位12個，差點位11個，極差點位7個，其各自在全流域所占比例依次為11.4%、20.5%、27.3%、25.0%和15.9%.健康和亞健康點位主要集中在河流源頭區(qū)的中上游，較差和極差的點位主要分布在河流的下游區(qū)域(圖3)；同時，河流健康狀況與其所處區(qū)域開發(fā)狀況密切相關(guān).

表2 西遼河流域魚類組成

-表示信息不詳.

表3 西遼河流域環(huán)境特征

圖2 通過判別能力篩選的5個候選指標在參照點和受損點的箱體圖Fig.2 Box-plots of 5 candidate metrics which accepted by the discriminatory power analysis between reference sites and impaired sites

表4 5個候選指標的Pearson相關(guān)性計算

Tab.4 The pearson correlation coefficient among 5 candidate biological metrics

候選指標M2M7M10M11M19M21M7-0．451M10-0．430．661M11-0．340．550．421M19-0．620．570．500．771

表5 F-IBI構(gòu)建參數(shù)

表6 西遼河流域魚類生物完整性評價標準

圖3 西遼河流域F-IBI評價結(jié)果Fig.3 F-IBI assessment results in the Xiliao River Basin

2.3 重要環(huán)境變量篩選及其保護限值計算

分析F-IBI值與環(huán)境因子的相關(guān)關(guān)系，以反映水化學、物理生境和采樣點周邊土地利用類型對河流健康的影響程度.采用Spearman相關(guān)性分析IBI值與環(huán)境因子之間的相關(guān)性.結(jié)果表明：IBI與水化學變量中的Cond(P=0.002)、Alk(P=0.003)呈極顯著負相關(guān)，與NH3-N(P=0.035)呈顯著負相關(guān)；與物理生境變量中的坡降(P<0.001)呈極顯著正相關(guān)，與泥沙比例(P=0.038)呈顯著負相關(guān)；與土地利用方式中的草地比例(P=0.006)呈極顯著正相關(guān)(表7).

西遼河流域F-IBI與Cond擬合曲線的拐點為531μS/cm，與NH3-N擬合曲線的拐點為0.55mg/L和1.35mg/L，與Alk擬合曲線的拐點為4.4mmol/L，與坡降擬合曲線的拐點為10.4%，與泥沙比例擬合曲線的拐點為47.2%，與草地比例擬合曲線的拐點為37.0%(圖4).經(jīng)獨立樣本t檢驗發(fā)現(xiàn)Cond、Alk、泥沙比例、草地比例的拐點兩側(cè)F-IBI得分差異顯著，其P值依次為0.002、0.043、0.044和0.009.NH3-N的保護限值(0.55mg/L)有意義(P=0.026)；NH3-N濃度(1.35mg/L)無意義(P=0.19>0.05).坡降的獨立樣本t檢驗P值為0.801，大于0.05，差異不顯著.

表7 西遼河流域F-IBI與環(huán)境因子之間的相關(guān)性系數(shù)

*表示P<0.05；**表示P<0.01.

圖4 西遼河流域與IBI顯著相關(guān)環(huán)境因子的LOWESS曲線及保護限值(垂直虛線代表環(huán)境保護限值)Fig.4 The LOWESS curves and protection limits of environmental factors that significantly associated with IBI in the Xiliao River Basin(vertical dashed line means the protection limit values of environmental factors)

將對西遼河流域F-IBI有重要影響的環(huán)境因子做相關(guān)性分析，結(jié)果表明：Cond與Alk具有較高的相關(guān)性(r=0.870)，其它環(huán)境因子間的相關(guān)性較低，相關(guān)系數(shù)均未超過0.5(圖5).

圖5 西遼河流域環(huán)境因子間的相關(guān)性分析Fig.5 Correlation coefficients of environmental factors in the Xiliao River Basin

3 討論

3.1 關(guān)于健康標準的劃分

對于確定F-IBI核心參數(shù)的健康評價標準，不同學者所采取的標準稍有不同[1,3,25-26]，這會造成最終的評價結(jié)果有所差異.一般來講，參照點IBI值的25%分位數(shù)常被作為健康河流樣點評價標準，該標準的前提是參照點選取較為容易且?guī)缀醪淮嬖谌藶楦蓴_因素.在參照點選取較為困難的情況下，有學者提出[3,25]以所有樣點IBI值的95%分位數(shù)作為健康評價標準.此外還有學者依據(jù)研究流域本身的特征，提出以參照點50%分位數(shù)作為健康評價標準，以5%分位數(shù)作為差等級評價標準，置于兩標準之間的作為亞健康評價標準[25-26].因此，應(yīng)以所選擇的評價標準能客觀反映出研究流域整體及各個樣點的健康狀況為宜.本研究以參考點F-IBI值的95%作為健康樣點評價標準，對于小于該標準的進行4等分，該評價結(jié)果所反映出的西遼河流域不同區(qū)域健康狀況及整體健康狀況與高欣等[16]利用大型底棲動物生物完整性指數(shù)評價的結(jié)果類似.本研究認為采用參考點F-IBI值的95%作為健康評價標準的依據(jù)適合該流域.

3.2 西遼河環(huán)境特征分析及對F-IBI的影響

西遼河流域是我國北方典型的農(nóng)業(yè)區(qū)與牧業(yè)區(qū)交匯過渡地帶，其所流經(jīng)區(qū)域為我國內(nèi)蒙古高原向遼河平原過渡地帶.該流域上游以山地、林地和草場為主，中下游以農(nóng)田和城鎮(zhèn)用地為主，并伴有水域和荒蕪用地[27-28].根據(jù)國家《地表水環(huán)境質(zhì)量標準》，西遼河流域水體CODMn、NH3-N、TN和TP等參數(shù)變化幅度均較大.該區(qū)域上游海拔較高，人類活動相對較少，水質(zhì)較為清潔，中下游區(qū)域受農(nóng)田灌溉、放牧、城市生活污水及企業(yè)廢水排放的影響，水質(zhì)較差.西遼河的物理生境特征以坡降變化和底質(zhì)粒徑組成的變化最為明顯.坡降變化劇烈的原因與該流域海拔變化幅度有關(guān)，較大的坡降有助于提高水體的氧含量[29].底質(zhì)粒徑組成變化大說明該流域河流的生境差異較大，底質(zhì)粒徑與區(qū)域自然地理地貌、氣候特征及人類活動有關(guān).

關(guān)于水化學對F-IBI影響的探討：第一，通常認為河流中的Cond高低受自然和人為兩方面因素影響[30]，不同區(qū)域Cond背景值差異較大.Kimmel等[31]以美國賓夕法尼亞州一條受煤礦廢水影響的溪流為對象，提出保護該河流魚類群落結(jié)構(gòu)的Cond保護限值為3000～3500μS/cm.本課題組在太子河流域的研究發(fā)現(xiàn)，保護溪流中大型底棲動物群落結(jié)構(gòu)完整的Cond保護限值為277.1μS/cm(未發(fā)表).本研究發(fā)現(xiàn)西遼河流域基于魚類群落的Cond保護限值則為531μS/cm.造成上述研究結(jié)果差異的原因可能與流域環(huán)境本底值以及研究的生物類群不同有關(guān).Kimmel等[31]的研究中Cond保護限值較高與水體長期受到工礦廢水影響有直接關(guān)系，在這種長期壓力條件下，水生生物對環(huán)境干擾的耐受性提高，生物群落結(jié)構(gòu)也逐漸趨于穩(wěn)定，較高的環(huán)境背景壓力下研究得到的保護限值有可能會高于其他低環(huán)境干擾壓力區(qū)域的.西遼河流域處于高原向平原過渡地帶，該流域大部分區(qū)域受到降水、灌溉等因素影響，河道內(nèi)水量不足，而在放牧造成草地退化的同時，流域坡度差異導致水體對河岸侵蝕加劇，地表風化嚴重，最終造成水體Cond背景值整體偏高，其中Alk對Cond的貢獻較大(圖5).另外，有研究發(fā)現(xiàn)[31-32]大型底棲動物對環(huán)境的干擾比魚類更為敏感，不同的研究對象會使環(huán)境因子保護限值結(jié)果存在差異.第二，Alk會對魚類群落產(chǎn)生影響.黎道豐等[33]發(fā)現(xiàn)高鹽堿度會降低魚類種類數(shù)量并導致魚類生長緩慢.本研究也發(fā)現(xiàn)相似的結(jié)果，堿度增加會降低F-IBI評價得分.第三，本研究發(fā)現(xiàn)NH3-N的有效保護限值為0.55mg/L，處于Ⅱ～Ⅲ類水質(zhì)之間，這與國家《地表水環(huán)境質(zhì)量標準》中漁業(yè)水質(zhì)要求相一致.W102點位的NH3-N濃度達到4.07mg/L，但F-IBI得分為74.3，同時該點位草地比例較高(84%)，河岸帶生境質(zhì)量較好.一方面較好的河岸帶生境條件有利于對外源NH3-N輸入起到緩沖作用[18]，另一方面較好的生境條件有利于維持較好的魚類群落結(jié)構(gòu)[34].

關(guān)于物理因素對F-IBI影響的探討：西遼河區(qū)域處于內(nèi)蒙古高原向遼河平原的過渡區(qū)域，坡降變化較大.同時受干旱和人為干擾等因素的影響，草場退化和荒漠化情況不容樂觀，因此確定上述環(huán)境因子的保護限值具有重要意義.坡降指河流水面單位距離的落差，其可以直接影響水體的流速，對于底質(zhì)泥沙含量比較高的區(qū)域，水流速度較大導致水體的濁度升高，不利于光合作用的正常進行，導致藻類產(chǎn)氧量降低，同時泥沙可堵塞魚類鰓絲，對其呼吸造成影響[35]；有研究表明[36-37]泥沙對水體中的NH3-N和重金屬等有毒物質(zhì)有吸附作用，并隨著水流帶到下游區(qū)域，對區(qū)域水質(zhì)的改善具有一定的作用，但會加重下游區(qū)域的污染情況；草地對于提高河流生境質(zhì)量，防止水土流失和土地沙漠化有著重要作用，同時可以緩解其它環(huán)境因子的負面影響[18].因此，在區(qū)域管理中，恢復和保持草地比例，防止土地荒漠化和水土流失，控制河流底質(zhì)中的泥沙比例，最終提高F-IBI狀況具有重要意義.

[1] Aparicio E, Carmona-Catot G, Moyle PBetal. Development and evaluation of a fish-based index to assess biological integrity of Mediterranean streams.AquaticConservation:MarineandFreshwaterEcosystems, 2011, 21(4): 324-337.

[2] Karr JR. Assessment of biotic integrity using fish communities.Fisheries, 1981, 6(6): 21-27.

[3] Barbour MT, Gerritsen J, Snyder BDetal. Rapid Bio-assessment Protocols for use in wadeable streams and rivers: Periphyton, benthic macroinvertebrates, and Fish: 2nd ed. Washington DC: 841-B-99-002. USA: U. S. Environmental Protection Agency, EPA, 1999.

[4] Kane DD, Gordon SI, Munawar Metal. The Planktonic Index of Biotic Integrity(P-IBI): An approach for assessing lake ecosystem health.EcologicalIndicators, 2009, 9(6): 1234-1247.

[5] EC(European Communities). Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 establishing a framework for Community action in the field of water policy.OfficialJournaloftheEuropeanCommunities, 2000, L 327:1-72.

[6] 黃 藝,舒中亞.基于浮游細菌生物完整性指數(shù)的河流生態(tài)系統(tǒng)健康評價——以滇池流域為例.環(huán)境科學,2013,34(8): 3010-3018.

[7] 殷旭旺，張 遠，渠曉東等.渾河水系著生藻類的群落結(jié)構(gòu)與生物完整性.應(yīng)用生態(tài)學報,2011,22(10): 2732-2740.

[8] Lunde KB, Resh VH. Development and validation of a macroinvertebrate index of biotic integrity(IBI) for assessing urban impacts to Northern California freshwater wetlands.EnvironmentalMonitoringandAssessment, 2012, 184(6): 3653-3674.

[9] 鄭海濤.怒江中上游魚類生物完整性評價[學位論文].武漢:華中農(nóng)業(yè)大學,2006.

[10] 渠曉東，劉志剛，張 遠.標準化方法篩選參照點構(gòu)建大型底棲動物生物完整性指數(shù).生態(tài)學報,2012,32(15): 4661-4672.

[11] 蔡 琨，張 杰，徐兆安等.應(yīng)用底棲動物完整性指數(shù)評價太湖生態(tài)健康.湖泊科學,2014,26(1): 74-82.DOI 10.18307/2014.0109.

[12] 陳宏文,張 萌,劉足根等.贛江流域淡水生態(tài)系統(tǒng)完整性與健康狀態(tài)的魚類F-IBI值評價.長江流域資源與環(huán)境,2011,20(9): 1098-1107.

[13] 陳 橋,徐東烔，張 翔等.太湖流域平原水網(wǎng)區(qū)底棲動物完整性健康評價.環(huán)境科學研究,2013,26(12): 1301-1308.

[14] 黃亮亮,吳志強,蔣 科等.東苕溪魚類生物完整性評價河流健康體系的構(gòu)建與應(yīng)用.中國環(huán)境科學,2013,33(7): 1280-1289.

[15] 朱 迪,陳 鋒,楊 志等.基于魚類生物完整性指數(shù)的水源地評價.水生態(tài)學雜志,2012,33(2): 1-5.

[16] 高 欣,張 遠,丁 森等.西遼河大型底棲動物生物完整性評價及同環(huán)境因子關(guān)系的探討.中國環(huán)境科學學會學術(shù)年會論文集(2014),2014: 3153-3163.

[17] 吳東浩,于海燕,吳海燕等.基于大型底棲無脊椎動物確定河流營養(yǎng)鹽濃度閾值——以西苕溪上游流域為例.應(yīng)用生態(tài)學報,2010,21(2): 438-488.

[18] Mantyka-Pringle CS, Martin TG, Moffatt DB. Understanding and predicting the combined effects of climate change and land-use change on freshwater macroinvertebrates and fish.JournalofAppliedEcology, 2014, 51(3): 572-581.

[19] 孫小舟,封志明,楊艷昭等.西遼河流域近60年來氣候變化趨勢分析.干旱區(qū)資源與環(huán)境,2009,23(9): 62-66.

[20] 國家環(huán)境保護總局《水和廢水監(jiān)測分析方法》編委會．水和廢水監(jiān)測分析方法:第4版.北京:中國環(huán)境科學出版社,2002.

[21] 裴雪嬌,牛翠娟,高 欣等.應(yīng)用魚類完整性評價體系評價遼河流域健康.生態(tài)學報,2010,30(21): 5736-5746.

[22] 宋智剛,王 偉,姜志強等.應(yīng)用F-IBI對太子河流域水生態(tài)健康評價的初步研究.大連海洋大學學報,2010,25(6): 480-487.

[23] Celeste PK, Nadanovsky P. Income and oral health relationship in Brazil: is there a threshold?CommunityDentOralEpidemiol, 2009, 37(4): 285-293.

[24] 解玉浩.東北地區(qū)淡水魚類.沈陽:遼寧科學技術(shù)出版社,2007: 29-439.

[25] 池仕運,竺維佳,施練東等.應(yīng)用底棲動物完整性指數(shù)評價水源地水庫溪流健康狀態(tài).水生態(tài)學雜志,2012,33(2): 16-25.

[26] Astin L. Developing biological indicators from diverse data: the potomac basic-wide index of benthic integrity(B-IBI).EcologicalIndicators, 2007, 7(4): 895-908.

[27] 孟 偉,張 遠,鄭炳輝.遼河流域水生態(tài)分區(qū)研究.環(huán)境科學學報,2007,27(6): 911-918.

[28] 王金龍.遼河流域水生態(tài)功能三級區(qū)水生態(tài)服務(wù)功能評價[學位論文].沈陽:遼寧大學,2013.

[29] 劉 婷,樊貴盛,陽 輝.瀟河流域水質(zhì)污染的空間變化特征.人民黃河,2014,36(3): 58-60.

[30] 李姍遲,鄭雄波,高全洲等.東江流域河水電導率形成機制.中山大學學報,2013,52(1): 142-146.

[31] Kimmel WG, Argent DG. Stream fish community responses to a gradient of specific conductance.Water,Air,andSoilPollution, 2010, 206(1/2/3/4): 49-56.

[32] Kennedy AJ, Cherry DS, Currie RJ. Field and laboratory assessment of a coal processing effluent in the Leading Creek Watershed, Meigs County, Ohio.ArchivesofEnvironmentalContaminationandToxicology, 2003, 44： 324-331.

[33] 黎道豐,蔡慶華.不同鹽堿度水體的魚類區(qū)系結(jié)構(gòu)及主要經(jīng)濟魚類生長的比較.水生生物學報,2000,24(5): 493-501.

[34] Pinto BCT, Araujo FG, Hughes RM. Effects of landscape and riparian condition on a fish index of biotic integrity in a large southeastern Brazil river.Hydrobiologia, 2006, 556(1): 69-83.

[35] 白音包力皋,陳興茹.水庫排沙對下游河流魚類影響研究進展.泥沙研究,2012,(1): 74-80.

[36] 王而立,王嗣淇,王 亮.西遼河流域沙土對氨氮的吸附行為研究.環(huán)境科學與技術(shù),2011,34(6G): 102-106.

[37] 易雨君.長江水沙環(huán)境變化對魚類的影響及棲息地數(shù)值模擬[學位論文].北京:清華大學,2008.

Assessment of the fish index of biotic integrity and its relationship with environmental factors in the Xiliao River Basin

ZHANG Hao1,2, DING Sen2, ZHANG Yuan2, JIA Xiaobo2, MENG Wei2& GUO Biao1,3

(1:OceanCollege,HebeiAgriculturalUniversity,Qinhuangdao066003,P.R.China)

(2:LaboratoryofRiverineEcologicalConservationandTechnology,ChineseResearchAcademyofEnvironmentalSciences,Beijing100012,P.R.China)

(3:BohaiFisheriesResearchInstituteofTianjin,Tianjin300457,P.R.China)

The environmental quality of Xiliao River Basin is assessed based on the fish index of biotic integrity(F-IBI). Protection limit values of the environmental factors correlated with the F-IBI are determined by the method of locally weighted regression scatterplot smoother(LOWESS). The results indicated that five indices, including the percentage of the Cobitidae fish, the number of demersal fish species, the number of omnivorous fish species, the number of tolerated fish individuals and the percentage of the nests spawning fish individuals, were suitable for constructing the F-IBI assessment system. The ratio method is used to unify dimensionless parameters. The 95% quantile of the F-IBI in all reference sites is used as the reference value of health conditions. Forty four sampling sites of Xiliao River Basin were assessed by F-IBI. The results found that five sites(11.4%) were in excellent grade, nine sites(20.5%) were in good grade, twelve sites(27.3%) were in normal grade, eleven sites(25.0%) were in poor grade, and seven sites(15.9%) were in very poor grade. Electrical conductivity, ammonia nitrogen, alkalinity and proportion of sand were negatively correlated with F-IBI scores in the Xiliao River Basin, whereas slope and proportion of grassland were positively correlated with F-IBI scores. The results of LOWESS and independent samplettest indicated that F-IBI scores significantly changed when the electrical conductivity, ammonia nitrogen concentration, alkalinity, proportion of sand and proportion of grassland reached the protection limit value, which was 531μS/cm, 0.55mg/L, 4.4mmol/L, 47.2% and 37.0%, respectively. The results of this study could provide beneficial information for the regional environmental management of fish integrity protection in the Xiliao River Basin.

Fish index of biotic integrity; LOWESS; Xiliao River Basin; environmental factors

*國家自然科學基金項目(41401066)和國家水體污染控制與治理科技重大專項(2012ZX07501-001-04)聯(lián)合資助．2014-10-26收稿；2015-02-05收修改稿．張浩(1987～)，男，碩士研究生；E-mail: shenduhao@126.com.

J.LakeSci.(湖泊科學)， 2015， 27(5)： 829-839

DOI 10.18307/2015.0509

?2015 byJournalofLakeSciences

**通信作者；E-mail: dingsen@craes.org.cn.