揣澤堯，任 改，洪 濱，趙 迪，郭晶晶

(生態(tài)環(huán)境部 華南環(huán)境科學(xué)研究所，廣東 廣州 510000)

土壤侵蝕是一個(gè)受地形、氣候、土壤特性、土壤覆蓋和人類活動(dòng)等多種因素影響的復(fù)雜問題[1]，會(huì)導(dǎo)致土壤生產(chǎn)力降低、水質(zhì)污染、洪水風(fēng)險(xiǎn)增加等。土壤侵蝕在全世界廣泛存在，對(duì)環(huán)境的負(fù)面影響很多，被認(rèn)為是最嚴(yán)重的世界性問題之一[2]。為了保護(hù)環(huán)境，減少水土流失對(duì)農(nóng)業(yè)、基礎(chǔ)設(shè)施、水質(zhì)等的負(fù)面影響，有必要開展水土保持專項(xiàng)行動(dòng)[3]。近年來許多學(xué)者致力于如何識(shí)別易受侵蝕地區(qū)的研究，大致分為以下幾類研究：根據(jù)植被覆蓋、土地利用和坡度三個(gè)標(biāo)準(zhǔn)，按保護(hù)的優(yōu)先順序?qū)α饔騼?nèi)的區(qū)域進(jìn)行排序；利用坡度和植被覆蓋兩個(gè)因素對(duì)區(qū)域侵蝕風(fēng)險(xiǎn)進(jìn)行分析。

當(dāng)前侵蝕率的評(píng)估是在不同的尺度水平(田間、流域)下，采用不同的經(jīng)驗(yàn)方法、概念和物理模型進(jìn)行的。這些方法存在一個(gè)缺點(diǎn)——需要大量的數(shù)據(jù)，而這在發(fā)展中國(guó)家的數(shù)據(jù)庫(kù)中較稀缺，存在一定的局限性?；诖?，本文采用層次分析法(AHP)[4]和理想解相似排序法(TOPSIS)[5]，利用坡度、土地利用情況和土壤類型3個(gè)因素對(duì)各子流域的侵蝕脆弱性進(jìn)行排序，以確定需要優(yōu)先保護(hù)的子流域。

1 研究區(qū)概況

選取廣東省韓江流域?yàn)檠芯繀^(qū)域。韓江流域地理位置為東經(jīng)115°13′～117°09′、北緯23°17′～26°05′，是廣東省的第二大流域，流域面積30 112 km2。韓江上游由梅江和汀江匯合而成，梅江為主流，發(fā)源于廣東省紫金縣，汀江發(fā)源于福建省寧化縣的賴家山，梅、汀兩江匯合后稱韓江，由北向南流經(jīng)廣東省的豐順、潮安等縣。以梅江為源頭，干流總長(zhǎng)470 km。將研究區(qū)分為20個(gè)子流域，分別為泗水、高思水、隆文水、小靖河、樟溪、九峰水、沙田水、大圣溪、九村溪、東山溪、荷泗水、永和水、松源水、黃陵水、平安水、龍?zhí)端?、南口水、合溪水、黃岡河、湯東河(下文中用sw-1～sw-20表示)。20個(gè)子流域由于受到侵蝕作用，水土流失、洪澇災(zāi)害、河道淤積等災(zāi)害時(shí)有發(fā)生，給當(dāng)?shù)厣鷳B(tài)環(huán)境以及人民群眾的生命財(cái)產(chǎn)安全帶來了嚴(yán)重威脅。因此，識(shí)別需要優(yōu)先保護(hù)的子流域?qū)τ谟械姆攀傅販p少侵蝕具有重要意義。研究區(qū)位置如圖1所示。

圖1 研究區(qū)位置

2 研究方法

2.1 土壤侵蝕脆弱性評(píng)價(jià)因子的識(shí)別

通過地理信息系統(tǒng)中1∶25 000、1∶50 000地形圖和1∶5 000正射影像，確定了20個(gè)子流域的平均坡度、土壤類型和土地利用方式，將其作為預(yù)測(cè)土壤侵蝕脆弱性的參數(shù)指標(biāo)。研究區(qū)主要土地利用方式有耕地、草甸、城市化地區(qū)、森林、退化森林和工業(yè)區(qū)6類，新積土、淋溶土、石質(zhì)土、崩積層土、中性形成土、黑鈣土、變性土、沖積土8種土壤類型。根據(jù)土壤質(zhì)地(粉砂、砂粒和黏粒的含量)和表層(0～30 cm土層)的有機(jī)碳含量，使用EPIC模型中給出的公式估算了這些土壤類型對(duì)侵蝕的敏感性，用K因子表示，詳見圖2。20個(gè)子流域用于評(píng)估土壤侵蝕脆弱性的三個(gè)參數(shù)(子流域平均坡度S、土地利用因子LU和土壤類型因素ST)數(shù)值如表1所示。

圖2 不同土壤類型的土壤可蝕性因子

表1 各子流域參數(shù)數(shù)值

2.2 層次分析法(AHP)

AHP方法具有良好的數(shù)學(xué)特性，是一種穩(wěn)定、靈活的決策工具，可用于求解復(fù)雜的多目標(biāo)決策問題，包括確定保護(hù)措施的優(yōu)先次序、滑坡易感性測(cè)繪或土壤侵蝕風(fēng)險(xiǎn)評(píng)估，且數(shù)據(jù)獲取難度較低，故將其應(yīng)用到子流域的脆弱性評(píng)級(jí)。

AHP方法的應(yīng)用一般有以下4個(gè)步驟：①將問題分解成目標(biāo)、準(zhǔn)則、次標(biāo)準(zhǔn)等層次結(jié)構(gòu)；②建立兩兩比較矩陣A=(aij)n×n，其中n為矩陣大小，aij表示第i個(gè)決策因素對(duì)第j個(gè)決策因素的重要性，aij≥0，aij·aji=1；③使用特征值法計(jì)算決策因素的相對(duì)權(quán)重(優(yōu)先級(jí))；④所有矩陣都必須滿足一致性檢驗(yàn)，即用一致性指數(shù)(CI)和隨機(jī)指數(shù)(RI)得到的一致性比(CR)要滿足小于0.10的要求。

2.3 理想解相似排序法(TOPSIS)

3 結(jié)果和討論

AHP方法的第一步是確定決策過程的總體目標(biāo)。根據(jù)坡度、土地利用和土壤類型等影響因素，將目標(biāo)確定為選擇最易受土壤侵蝕影響的子流域，這些因素被分解為次級(jí)指標(biāo)(6個(gè)土地利用等級(jí)和8個(gè)土壤類型)。

將所有的決策因素構(gòu)造成判斷矩陣，使用1～9度量表進(jìn)行配對(duì)比較。通過特征值法計(jì)算決策元素的權(quán)重，一級(jí)指標(biāo)和次級(jí)指標(biāo)的成對(duì)比較結(jié)果見表2和表3。

表2 一級(jí)指標(biāo)成對(duì)比較結(jié)果

表3 次級(jí)指標(biāo)成對(duì)比較結(jié)果

將備選方案通過平均坡度(S)、土地利用因子(LU)和土壤類型因素(ST)進(jìn)行成對(duì)比較。使用脆弱性指數(shù)ER對(duì)20個(gè)子流域進(jìn)行排序，該指數(shù)的計(jì)算方式是將指標(biāo)權(quán)重(wi)與替代權(quán)重(dj)相乘，ER值越高，說明越容易受到侵蝕。具體計(jì)算結(jié)果見表4。

表4 層次分析法計(jì)算結(jié)果

接下來應(yīng)用TOPSIS方法確定子流域的排序，將通過AHP方法獲得的權(quán)重代入TOPSIS方法自帶的程序，并對(duì)決策矩陣進(jìn)行歸一化處理，構(gòu)造加權(quán)歸一化矩陣確定正理想解和負(fù)理想解，`具體見表5。

表5 正理想解(A+)和負(fù)理想解(A-)

子流域的最終排名在很大程度上取決于賦予主要指標(biāo)的權(quán)重。為了檢查結(jié)果的穩(wěn)定性，考慮了以下三種情況：

再計(jì)算與理想解(Cj+)的相對(duì)距離，根據(jù)該值對(duì)子流域進(jìn)行排序，具體見表6。

情景1：土地利用占大多數(shù)權(quán)重(wLU=0.714 3，wST=wS=0.142 9)。

情景2：所有指標(biāo)具有相同的權(quán)重(wLU=wST=wS=0.333 3)。

情景3：土地利用占較大權(quán)重 (wLU=0.633 4，wST=0.106 2，wS=0.260 5)。

在情景1的情況下，AHP方法將sw-17、sw-4、sw-6、sw-18和sw-9判別為最易受土壤侵蝕影響的子流域，而TOPSIS方法將sw-18、sw-17、sw-4、sw-3和sw-16判別為最易受土壤侵蝕影響的子流域。圖3比較了AHP方法和TOPSIS方法在三種場(chǎng)景下的排序。Spearman相關(guān)系數(shù)結(jié)果表明基于AHP方法的排序與基于TOPSIS方法的排序具有很強(qiáng)的相關(guān)性，相關(guān)系數(shù)分別為：0.732 3(情景1)、0.703 8(情景2)和0.729 3(情景3)。

表6 每個(gè)子流域與理想解的相對(duì)接近度

圖3 三種情境下AHP和TOPSIS方法的子流域排序比較

為驗(yàn)證上述排序結(jié)果的準(zhǔn)確性，通過對(duì)各分區(qū)采取現(xiàn)場(chǎng)考察、資料查閱及結(jié)果檢測(cè)等手段，選取耕地面積和坡度兩項(xiàng)參數(shù)驗(yàn)證了排序結(jié)果，具體見表7。

表7 排序結(jié)果驗(yàn)證

通過上述調(diào)查研究，發(fā)現(xiàn)sw-17流域耕地面積占比高達(dá)73%，耕地面積占比最低的sw-4流域也達(dá)到了53%。這些需要進(jìn)行優(yōu)先保護(hù)的子流域(sw-17、sw-18、sw-6、sw-4和sw-9)均具有耕地面積占比較大(耕地面積占比超過50%)和邊坡陡峭(坡度超過14°)的特點(diǎn)，同時(shí)這也印證了上述排序結(jié)果的準(zhǔn)確性。

4 結(jié) 論

本研究采用AHP和TOPSIS方法，根據(jù)現(xiàn)有土地利用、土壤類型以及平均坡度對(duì)子流域的侵蝕脆弱性進(jìn)行了排序，兩種方法結(jié)果一致，將小靖河、九峰水、九村溪、南口水及合溪水流域列為需要優(yōu)先保護(hù)的脆弱性地區(qū)。同時(shí)通過實(shí)地驗(yàn)證、查閱資料等方式，從耕地面積占比和平均坡度兩個(gè)指標(biāo)驗(yàn)證了排序的準(zhǔn)確性。識(shí)別這些地區(qū)和對(duì)其進(jìn)行排序，使研究人員能夠更合理、更可持續(xù)地實(shí)施長(zhǎng)期保護(hù)戰(zhàn)略，為自然資源規(guī)劃管理提供重要依據(jù)。下一步的工作計(jì)劃是針對(duì)土壤侵蝕特點(diǎn)提出適宜的措施。