高宇婷，高甲榮，顧嵐，錢斌天，王越，郭維

(北京林業(yè)大學(xué)水土保持學(xué)院，水土保持與荒漠化防治教育部重點(diǎn)實(shí)驗(yàn)室，100083，北京)

河流健康狀況是對(duì)河流生態(tài)系統(tǒng)結(jié)構(gòu)、功能特征的評(píng)判，處于健康狀況的河流結(jié)構(gòu)合理、功能健全，包括正常能量流動(dòng)且物質(zhì)循環(huán)沒有受到破壞，對(duì)自然干擾的長期效應(yīng)具有抵抗力和恢復(fù)力，能夠維持自身的組織結(jié)構(gòu)長期穩(wěn)定，并能發(fā)揮其正常的生態(tài)環(huán)境效益［1-2］。國際上通過水質(zhì)理化參數(shù)、河流形態(tài)結(jié)構(gòu)、河流生物指標(biāo)、河流水文特征、河岸帶狀況5類指標(biāo)表征河流健康狀況;然而河流生態(tài)系統(tǒng)中許多變量之間存在的相關(guān)性增加了問題分析的復(fù)雜性［3-4］。K.M.Rowntree 等［5］運(yùn)用 RHP 方法，通過河流無脊椎動(dòng)物、魚類、河岸植被、生境完整性、水質(zhì)、水文、形態(tài)等7類指標(biāo)評(píng)價(jià)河流的健康狀況;J.A.Gore等［6］認(rèn)為，河流生態(tài)評(píng)價(jià)須考慮水文、土地利用、地貌、水質(zhì)、生物與生態(tài)、娛樂、經(jīng)濟(jì)和文化等方面。澳大利亞河流評(píng)價(jià)計(jì)劃(AUSRIVAS)、河流無脊椎動(dòng)物預(yù)測和分類系統(tǒng)(RIVPACS)、生物完整性指數(shù)(IBI)、溪流狀態(tài)指數(shù)(ISC)［7-8］通過建立指標(biāo)體系對(duì)河流生態(tài)健康狀況進(jìn)行評(píng)價(jià)，不但提高了河流管理的質(zhì)量，而且有助于不同河段之間進(jìn)行比較，同時(shí)，還可對(duì)受損河段生態(tài)恢復(fù)的成效進(jìn)行評(píng)估，提高了河流健康狀況管理決策的能力。這些指標(biāo)體系雖能較全面地評(píng)價(jià)河溪的近自然狀況，但有些因子只能進(jìn)行定性描述，帶有一定的主觀性，降低了評(píng)價(jià)結(jié)果的可信度。此外，這些評(píng)價(jià)體系帶有地域性，并不適合我國河流自然性的評(píng)價(jià)［9］。筆者通過調(diào)查北京市湯河21個(gè)河段的17個(gè)指標(biāo)，經(jīng)過因子分析得出影響北京市湯河生態(tài)健康的關(guān)鍵因子［10］，研究結(jié)果可為北京市河流生態(tài)健康的管理提供參考。

1 研究區(qū)概況

研究區(qū)位于北京市懷柔區(qū)湯河中下游河岸帶，處于北京市最北端，降水集中在夏季，多年平均降水量約為644.2 mm。湯河發(fā)源于河北豐寧縣鄧家柵子，于大南溝門進(jìn)入北京市懷柔區(qū)，介于E115°30'～116°37.2'，N40°19.2'～ 41°39.2'之間，河流由北到南流速漸緩，在海拔264 m處形成湯河口盆地［11］，流量約為2.45 m3/s。湯河在懷柔境內(nèi)南北長約52.4 km，流域面積618.3 km2。該流域?qū)倥瘻貛Т箨懶约撅L(fēng)型半濕潤氣候，四季分明，雨熱同期，夏季暖熱濕潤，冬季寒冷少雪。目前耕作制度多為一種一熟制，是全區(qū)籽種、雜糧、淡季蔬菜的主產(chǎn)區(qū)。研究區(qū)為帽山至湯河口段，全長46.5 km。北京市湯河評(píng)價(jià)河段位置見圖1。

2 研究方法

2.1 指標(biāo)設(shè)定

國外學(xué)者對(duì)河流健康評(píng)價(jià)影響因素做了大量理論研究，主要集中在地貌、生態(tài)、水文3方面，其中，水文、河流狀況、生態(tài)功能被認(rèn)為是影響河流健康的最主要因素。結(jié)合北京市湯河實(shí)際狀況，根據(jù)相關(guān)文獻(xiàn)資料，將河流健康的主要影響因素劃分為水量因素、水質(zhì)因素、河岸帶因素、物理結(jié)構(gòu)因素、景觀娛樂因素5類。

圖1 北京市湯河評(píng)價(jià)河段位置Fig.1 Location of evaluated sections of Tanghe River in Beijing

水量因素的主要影響子因素為河寬X1、河深X2、流速X3;水質(zhì)因素的主要影響子因素為溶解氧X4、總氨氮X5、pH值X6;河岸帶因素的主要影響子因素為坡度X7、縱比降X8、橫比降X9、緩沖帶植被寬度X10;物理結(jié)構(gòu)因素的主要影響子因素為水與河底接觸點(diǎn)個(gè)數(shù)X11、河流彎曲度X12、心灘面積X13、心灘形狀系數(shù)X14;景觀娛樂因素的主要影響子因素為河流生物多樣性 X15、河灘地覆蓋率 X16、閘壩個(gè)數(shù)X17。

2.2 指標(biāo)測定方法

在2011年的4—8月進(jìn)行野外調(diào)查，在生態(tài)特征典型代表處共21個(gè)河段布設(shè)調(diào)查樣點(diǎn)，采用定量測定法觀測記錄其相關(guān)特性。其中:河寬、河深、流速用卷尺和測速儀測量;溶解氧、總氨氮、pH值用野外便攜式試劑盒測定;坡度、縱比降、橫比降、緩沖帶植被寬度、水與河底接觸點(diǎn)個(gè)數(shù)、河流彎曲度、心灘面積、心灘形狀系數(shù)、河流生物多樣性、河灘地覆蓋率、閘壩個(gè)數(shù)通過實(shí)地測量、內(nèi)業(yè)計(jì)算獲得。

2.3 指標(biāo)檢驗(yàn)

檢驗(yàn)因子相關(guān)性的方法主要有KMO(Kaiser-Meyer-Olk in measure of Sampling Adequacy)樣本測度法和 Baetlett球形檢驗(yàn)法(巴特利 Test of sphericity)。KMO是用來檢驗(yàn)變量間偏相關(guān)性的，其值愈接近1表明愈適合進(jìn)行因子分析。一般認(rèn)為，KMO值在0.7以上適合做因子分析。KMO和Baetlett球形檢驗(yàn)結(jié)果見表1?？芍?，指標(biāo)間顯著性概率是0，小于1%，說明數(shù)據(jù)具有相關(guān)性，是適宜進(jìn)行因子分析的。同時(shí)在顯示結(jié)果中，KMO值為0.813，說明以所涉及的21個(gè)樣點(diǎn)的5類因素17個(gè)指標(biāo)為基礎(chǔ)的數(shù)據(jù)可以進(jìn)行因子分析。

表1 KMO和Bartlett檢驗(yàn)結(jié)果Tab.1 KMO and Bartlett's test

2.4 因子分析

用SPSS20.0對(duì)各指標(biāo)數(shù)據(jù)進(jìn)行因子分析，具體步驟如下。

1)建立多維數(shù)據(jù)矩陣:將所測數(shù)據(jù)以所測指標(biāo)為行、以樣點(diǎn)為列建立多維數(shù)據(jù)矩陣。

2)相關(guān)系數(shù)矩陣:計(jì)算原始變量的相關(guān)系數(shù)矩陣，用以判斷所選變量有無相關(guān)關(guān)系［12］。

3)提取公共因子:確定因子求解方法和因子個(gè)數(shù)。取方差大于1的因子，按因子累計(jì)方差貢獻(xiàn)率確定因子個(gè)數(shù)［13］。

4)因子旋轉(zhuǎn):正交陣右乘A使旋轉(zhuǎn)后的因子載荷陣結(jié)構(gòu)簡化，每個(gè)變量僅在1個(gè)公共因子上有較大的載荷，而在其余的因子上載荷比較小，使得旋轉(zhuǎn)后所得到的因子載荷陣的總方差達(dá)到最大值。

5)計(jì)算因子得分:求出各指標(biāo)的因子得分值，為后續(xù)分析提供參考［14］。

最后用Excel 2007軟件將所得各因子得分均值與各樣點(diǎn)因子得分作雷達(dá)圖進(jìn)行分析。

2.5 指標(biāo)提取

對(duì)檢驗(yàn)后數(shù)據(jù)作標(biāo)準(zhǔn)化處理，采用主成分分析法計(jì)算出用方差最大正交旋轉(zhuǎn)后的因子載荷矩陣，依據(jù)各共同因素中因素負(fù)荷量的大小排序，根據(jù)貢獻(xiàn)率大于85%且特征值大于1的標(biāo)準(zhǔn)抽取因子。

由于因子分析法［15］產(chǎn)生的前4個(gè)主成分特征值分別為1.167、1.128、1.059、1.044(均大于1)，提取出前4個(gè)主成分即湯河生態(tài)健康特征因子，且其對(duì)應(yīng)的方差貢獻(xiàn)率分別為27.675%、26.112%、17.127%、14.918%，累積方差貢獻(xiàn)率超過85%，可以反映17個(gè)指標(biāo)的絕大部分信息。相關(guān)因子載荷陣見表2。

2.6 指標(biāo)修正

提取到各組因子后，為了便于進(jìn)行數(shù)據(jù)分析及因子作用表達(dá)需將其命名，由表2可知:

1)主成分1中，坡度、縱比降、橫比降、閘壩個(gè)數(shù)載荷很大，是主要反映河流岸坡結(jié)構(gòu)的因子，命名為河流岸坡因子(因子1)，其值越大表明河道岸坡情況對(duì)河流生態(tài)健康影響情況越大;

表2 因子旋轉(zhuǎn)后河流健康指標(biāo)因子載荷陣Tab.2 Factor loading array of Tanghe River health factors

2)主成分2中，pH值、溶解氧、總氨氮、河流生物多樣性載荷很大，是主要反映河流水質(zhì)情況的因子，命名為河流水質(zhì)因子(因子2)，其值越大表明河流水質(zhì)情況對(duì)河流生態(tài)健康影響情況越大;

3)主成分3中，流速、水與河底接觸點(diǎn)個(gè)數(shù)、河寬、河深載荷相對(duì)稍大，命名為河床作用因子(因子3)，其值越大表明河床作用對(duì)河流生態(tài)健康影響情況越大;

4)主成分4中，植被帶寬度、河流彎曲度、心灘面積、心灘形狀系數(shù)、河灘地覆蓋率載荷相對(duì)稍大，命名為河流植被因子(因子4)，其值越大表明河流植被對(duì)河流生態(tài)健康影響情況越大。

修正后的因素指標(biāo)體系見表3。

3 結(jié)果與分析

3.1 湯河因子分析

將表3生成的修正指標(biāo)應(yīng)用于北京市湯河生態(tài)健康評(píng)價(jià)，結(jié)果見圖2。從圖2(a)可以看出，湯河因子1和2、3和4得分較相近，均對(duì)湯河河流生態(tài)健康相關(guān)性顯著;從圖2(b)～(e)可以看出，因子1在喇叭溝門、北八道河、古石溝門得分較大，因子2在二道河?xùn)|溝、長哨營得分較大，因子3在長哨營和湯河口得分較大，因子4在頭道穴得分較大。說明這些河段對(duì)各指標(biāo)因子影響較大。

表3 修正后的湯河河流生態(tài)健康影響因素指標(biāo)體系Tab.3 Tanghe River ecosystem health influencing factor model after correction

3.2 指標(biāo)影響分析

從因子的角度分析影響作用最大河段、最小河段及平均值，其相關(guān)度見表4?？梢钥闯?4個(gè)因子特征值中，因子4的平均值最大，其中頭道穴段所占份額最多，其次是因子3，對(duì)應(yīng)湯河口段所占比例最多，說明因子3、因子4對(duì)湯河總體生態(tài)健康影響較大，集中體現(xiàn)在頭道穴段和湯河口段;4個(gè)因子特征值中，因子2的平均值最小，其次是因子1，均為老西溝門段，所占比例最多，說明老西溝門段為因子1和因子2影響湯河總體生態(tài)健康的關(guān)鍵河段。

另外，從表2中還可以看出，各個(gè)子因素的方差貢獻(xiàn)率中，因子1中坡度(9.627%)，因子2中pH值(10.675%)、溶解氧(12.568%)，因子3中流速(10.390%)、水與河底接觸點(diǎn)個(gè)數(shù)(16.112%)，因子4中心灘面積(8.693%)、心灘形狀系數(shù)(8.113%)、彎曲度(8.010%)、緩沖帶植被寬度(4.540%)與湯河總體生態(tài)健康的相關(guān)度較大。

表4 影響湯河生態(tài)健康的因子特征Tab.4 Characteristics of factors affecting Tanghe River

3.3 河段影響分析

將圖2(b)～(e)中顯著正相關(guān)(邊緣最外2環(huán)分布河段)、負(fù)相關(guān)(在雷達(dá)圖最內(nèi)環(huán)分布河段)相關(guān)河段在各因子中的分布比例整理為表5。

1)將數(shù)值前2列橫向求和，可看出顯著相關(guān)河段在因子4、因子1方面所占比例均達(dá)到40%，且第2列數(shù)字均大于第1列數(shù)字。說明河段對(duì)北京市湯河生態(tài)健康狀況的影響主要體現(xiàn)在因子1和因子4上，且下游對(duì)其影響要高于中游(上游為河北豐寧縣鄧家柵子，非本研究涉及范圍)。

2)顯著正相關(guān)河段在因子1中比例最高，顯著負(fù)相關(guān)河段在因子4中比例最高。說明顯著相關(guān)河段在因子1上對(duì)河流生態(tài)健康的影響是正相關(guān)作用，在因子4方面的影響是負(fù)相關(guān)作用。

4 結(jié)論

1)河流岸坡因子、河流水質(zhì)因子、河床作用因子，植被作用因子可以用來評(píng)價(jià)北京市湯河生態(tài)健康狀況，其中河床作用因子、植被作用因子對(duì)北京市湯河河流生態(tài)健康評(píng)價(jià)影響最大。

2)各河段對(duì)北京市湯河生態(tài)健康狀況各因子方面的影響可以通過因子分析法直觀看到，符合北京市湯河生態(tài)地域性特點(diǎn)對(duì)北京市湯河生態(tài)健康狀況的影響。

圖2 因子分析示意圖Fig.2 Diagram of evaluation factors

表5 各河段對(duì)河流生態(tài)健康相關(guān)度比例Tab.5 Percentage of river ecosystem health-related degrees in individual reaches

［1］Dobson A P，Bradshaw A D，Baker A M.Hopes for the future:restoration ecology and conservation biology［J］.Science，1997，277:515-522

［2］The FederalInteragency Stream Restoration Working Group.Stream Restoration:principles，processes and practices［EB/OL］.(2008-09-23)［2013-05-10］.http:∥www.usda.gov/stream restoration

［3］Cairns J Jr.The status of the theoretical and applied science of restoration ecology［J］.The environmental Professional，1991，13:186-195

［4］Raven P J，Holmes N T H，Naura M，et al.Using river habitat survey for environmental assessment and catchment planning in the UK［J］.HydroBiologia，2000，422/423:359-367

［5］Rowntree K M，Wadeson R A.A hierarchical framework for categorising the geomorphology of selected South African rivers［J］.Final Report to the Water Research Commission，1999(1):497

［6］Gore J A，Shields F D.Can large rivers be restored?［J］.BioScience，1995(45):142-152

［7］吳阿娜，楊凱，車越，等.河流健康狀況的表征及其評(píng)價(jià)［J］.水科學(xué)進(jìn)展，2005，16(4):602-608

［8］李文君，邱林，陳曉楠，等.基于集對(duì)分析與可變模糊集的河流生態(tài)健康評(píng)價(jià)模型［J］.水利學(xué)報(bào)，2011，42(7):775-782

［9］高陽，高甲榮，陳子珊，等.河溪近自然治理評(píng)價(jià)指標(biāo)體系探討以及應(yīng)用［J］.水土保持研究，2007，14(6):404-408

［10］姚國征，丁國棟.基于判別、因子分析的采煤沉陷風(fēng)沙區(qū)土壤質(zhì)量評(píng)價(jià)［J］.農(nóng)業(yè)工程學(xué)報(bào)，2012，28(7):200-207

［11］高陽，高甲榮，李付杰，等.基于河道-濕地-緩沖帶復(fù)合指標(biāo)體系的京郊河溪生態(tài)評(píng)價(jià)體系［J］.生態(tài)學(xué)報(bào)，2008，28(10):5150-5160

［12］張晶，董哲仁.河流健康全指標(biāo)體系的模糊數(shù)學(xué)評(píng)價(jià)方法［J］.水利水電技術(shù)，2010，41(12):16-21

［13］向東進(jìn)，李宏偉，劉小雅.實(shí)用多元統(tǒng)計(jì)分析［M］.武漢:中國地質(zhì)大學(xué)出版社，2005:174-183

［14］彭放，楊瑞琰，羅文強(qiáng)，等.數(shù)學(xué)建模方法［M］.北京:科學(xué)出版社，2007:53-65

［15］任若恩，王惠文.多元統(tǒng)計(jì)數(shù)據(jù)分析:理論、方法、實(shí)例［M］.國防工業(yè)出版社，1997:89-97