貴州省降雨侵蝕力時空變化特征研究

戴海倫，苑 爽，張科利，朱 青

（1．北京師范大學 地表過程與資源生態國家重點實驗室，北京100875；2．北京師范大學 地理學與遙感科學學院，北京100875；3．貴州省農科院 土肥研究所，貴陽550006）

降雨侵蝕力R是土壤侵蝕預報模型中的重要因子［1］，其值大小表征了降雨對于土壤的潛在侵蝕能力。研究降雨侵蝕力的空間分布對指導水土流失防治、水土保持區劃具有重要的價值與意義。有關降雨侵蝕力的研究，基本都以USLE［2］模型中的計算方法為前提，以小區實測資料為基礎，驗證不同地區的指標組合形式、分析空間分布規律［3－12］。我國疆域遼闊，降雨特征的地域分異性明顯，降雨侵蝕力研究既要考慮模型的統一性，又要充分體現地域差異性。因此，利用不同類型區徑流小區實測資料來驗證并計算降雨侵蝕力十分必要。

貴州省地處喀斯特地貌強烈發育地區，谷深坡陡，人多地少，陡坡墾殖十分普遍，水土流失極為嚴重。遙感調查數據顯示，水土流失面積占全省總面積的41．62%［13］。受石灰巖地區特殊的地球化學過程影響，貴州省土層厚度普遍較薄。在相同的侵蝕強度下，水土流失的危害遠大于其他地區。日趨嚴重的石漠化過程已經嚴重威脅到該區域的水土資源狀況與農業安全。因此，開展喀斯特地區土壤侵蝕定量估算及防治措施效益評價研究尤為重要。寧麗丹等［7］引用并修正了 Yu和 Rosewell［14－15］的降雨侵蝕力計算公式，采用EI10作為最佳組合指標計算了我國西南地區的降雨侵蝕力。楊子生［16］利用布設在滇東北的小區觀測資料研究了坡耕地土壤流失方程，建議采用EI60組合來計算降雨侵蝕力指標。許月卿等［17］利用寧麗丹等［7］修正的降雨侵蝕力計算方法，用貴州省19個站點的日降雨量資料計算并探討了貴州省降雨侵蝕力的時空分布規律，指出貴州省降雨侵蝕力主要分布在夏季；近50a來貴州省降雨侵蝕力呈增加趨勢。盡管上述降雨侵蝕力的計算及分布研究對USLE在我國、特別是西南地區的推廣應用做出了重要貢獻，但在進行貴州省降雨侵蝕力的計算時仍存在以下問題：（1）不同層次、不同區域尺度的水土流失估算需要不同精度的降雨侵蝕力分布圖。然而，在目前的成果中，研究區域范圍多集中于全國、西南地區，尚無法滿足貴州地區降雨侵蝕力計算的精度要求；（2）在進行西南喀斯特地區以及貴州省的降雨侵蝕力研究時，不同學者采用的計算指標不統一，且缺乏當地徑流小區實測降雨過程及泥沙資料的支持。另一方面，已有研究中所采用的站點偏少、密度較低，難以滿足局部區域精確計算的要求。針對以上問題，本文利用貴州省羅甸縣徑流小區的次降雨過程資料，進行貴州省降雨侵蝕力指標選定及計算方法確定，嘗試利用常規雨量站觀測資料計算降雨侵蝕力的簡易方法，并繪制貴州省降雨侵蝕力分布圖，旨在為區域水土流失預報與治理提供理論依據。

1 材料與方法

1．1 數據來源

本文利用貴州省羅甸縣興隆鄉境內小區1992—1997年的32場實測降雨過程資料和土壤流失量資料，進行降雨侵蝕力指標的評價篩選，通過回歸分析得出降雨侵蝕力的簡易算法。再利用貴州省內分布的264個站點（圖1）1956—2000年的實測數據，通過插值算法，繪制貴州省的降雨侵蝕力等值線圖。羅甸小區為長21m、寬5m，坡度為20°的裸露小區，土壤為黃紅壤。各站點雨量資料來源于貴州省水利廳。

圖1 貴州雨量站分布圖

1．2 指標選取方法

盡管USLE采用EI30指標來表征降雨侵蝕力的大小，但由于降雨特性和土壤侵蝕特征的地區差異，應用于各地的指標不盡相同。如有學者建議在黃土高原地區宜采用EI10［11］，而在東北黑土區和南方紅壤地區分別采用E60I30，EI60更適合［18－19］。雖然表達降雨侵蝕力的最佳指標存在地區差異，但其基本表達形式都是降雨動能（E）與某一時段最大降雨強度（I）的乘積（公式1）。本文利用羅甸小區的降雨過程資料和土壤流失量資料，分別計算每場降雨的總動能，再將每場降雨的總動能與該次降雨的5min最大雨強、10min最大雨強、20min最大雨強、30min最大雨強和60min最大雨強分別相乘，生成5個表征降雨侵蝕力的指標EI5，EI10，EI20，EI30，EI60，即R5，R10，R20，R30，R60。 再 分 別 以 每 場 降 雨 的R5，R10，R20，R30，R60為自變量，與該場降雨的土壤侵蝕量進行相關分析（表1）。根據相關系數大小，確定貴州地區計算降雨侵蝕力的最佳指標組合。其中，單位降雨動能的計算采用Brown和Foster［20］推薦的公式，見公式（2）。再由公式（3）得到次降雨總動能，最后利用公式（1）計算次降雨侵蝕力。

式中：er——每一時段的單位降雨動能［MJ／（hm2·mm）］；ir——時段雨強（mm／h）；Pr——對應時段降雨量（mm）；E——次降雨總動能（MJ／hm2）；Im——一定時長范圍內的最大雨強（mm／h）；R——次降雨侵蝕力［（MJ·mm）／（hm2·h）］。

1．3 構建降雨侵蝕力簡易算法

計算降雨侵蝕力需要降雨過程資料，而一般雨量站只有降雨量記錄。為了充分利用觀測資料并保證地區降雨侵蝕力分布圖的精度，需建立采用常規降雨記錄值計算降雨侵蝕力的簡易方法。以由公式（1）計算出的次降雨侵蝕力值為因變量，利用SPSS回歸分析其與降雨量之間關系，建立根據降雨量計算降雨侵蝕力的方法。由此便可利用貴州省雨量站的月降雨量資料計算全省各站不同時段的降雨侵蝕力，并通過Kriging方法進行插值、繪制貴州省降雨侵蝕力等值線圖。

表1 不同降雨侵蝕力指標與土壤流失量的相關性

2 結果與分析

2．1 指標篩選

相關分析結果（表1）表明，雖然羅甸小區的土壤侵蝕量與不同計算指標的相關關系都不如黃土高原等地區高，但仍以EI30的相關性最好（圖2），相關系數為0．640 2。因此，建議在貴州地區計算降雨侵蝕力時仍采用EI30。

圖2 貴州省羅甸小區降雨侵蝕力與土壤流失量關系

2．2 降雨侵蝕力簡易算法

參考其他地區的研究成果，分別嘗試用R＝a＋bP，R＝a＋blnP，R＝abP＋c和R＝aPb＋c共4種不同模型擬合降雨量與降雨侵蝕力之間的關系（a，b，c分別為擬合參數），結果表明降雨量與降雨侵蝕力之間存在良好的冪函數關系（圖3），并得到最佳擬合方程，見式（4），決定系數R2達到了0．724。

式中：R——降雨侵蝕力R30［（MJ·mm）／（hm2·h）］；P——降雨量（mm）。

用常規降雨量計算降雨侵蝕力，可大大增加數據獲取的范圍，使降雨侵蝕力估算結果更符合實際。但需要說明的是，由于在方程擬合時使用的資料都是侵蝕性降雨，而不論是日降雨還是月降雨資料都會包含一定比例的非侵蝕性降雨。如果直接使用常規降雨量計算降雨侵蝕力，結果勢必會高估。因此，在使用該公式估算前需要進行必要的修訂。為此，從南到北分別選擇有自記雨量計資料的羅甸、龍里、遵義和畢節四站的降雨資料，按年計算出侵蝕性降雨量占全年總雨量的比例系數（表2），并以平均系數作為采用降雨量計算降雨侵蝕力的修訂系數。結合上述簡易公式，利用月降雨量資料計算月降雨侵蝕力，再通過將各月值累加計算，得到年降雨侵蝕力值。

圖3 貴州省羅甸小區降雨量與降雨侵蝕力關系

表2 不同地區侵蝕性降雨占年降雨量的比例系數

2．3 貴州省降雨侵蝕力分布及變化特征

2．3．1 多年平均降雨侵蝕力的空間分布 采用公式（4）計算得到的貴州省1956—2000年多年平均降雨侵蝕力的空間分布如圖4所示。由圖4可見，貴州省多年平均降雨侵蝕力具有明顯的區域差異性，總體表現為由南向北遞減。全省范圍內西南部和東南部降雨侵蝕力較大，東北部次之，中部和西北部相對較低。同時，貴州省降雨侵蝕力還存在若干個高值中心。貴州省降雨量變化范圍為730．8～1 708．5mm，均值為1 184．98mm。多年平均降雨侵蝕力介于3 051．0～8 801．6 （MJ·mm）／（hm2·h），均 值 為 5 299．8（MJ·mm）／（hm2·h）。西南部的盤縣、興義形成高值中心，降雨侵蝕力在6 300（MJ·mm）／（hm2·h）左右；西部的織金、普定附近為另一高值中心，降雨侵蝕力在6 650（MJ·mm）／（hm2·h）左右；另外，東南部的荔波、榕江附近和位于東北部的松桃苗族自治縣也形成了高值中心，降雨侵蝕力在6 300（MJ·mm）／（hm2·h）左右。

2．3．2 貴州省降雨侵蝕力的年內分布 通過計算多年平均月降雨侵蝕力并插值制圖后發現，貴州省降雨侵蝕力的年內變化十分明顯，表現出顯著的干濕季變化。每年10月至翌年5月間為干季，6—9月為濕季。干季的降雨侵蝕力普遍偏低，8個月之和變化于1 501．1～4 632．7（MJ·mm）／（hm2·h）之間，空間分布上表現為從東部、東南部向西部、西北部逐漸遞減的趨勢（圖5），干季8個月的降雨侵蝕力占年均降雨侵蝕力的42．6%。在東南部的宰便和東北部的松桃苗族自治縣、銅仁形成高值中心，干季降雨侵蝕力都在3 220（MJ·mm）／（hm2·h）以上；在北部的習水和西北部的赫章形成低值中心，干季降雨侵蝕力都在1 540（MJ·mm）／（hm2·h）以下。全省濕季降雨侵蝕 力 普 遍 偏 高，集 中 于 1 744．0～5 921．2（MJ·mm）／（hm2·h），在空間分布上表現出西部明顯大于東部的特征（圖6）。濕季4個月的降雨侵蝕力占年均降雨侵蝕力的57．4%。高值中心位于西南部的興義市、盤縣和西部的安順市，濕季降雨侵蝕力在4 200 MJ／（mm·hm2·h）左右；東部的錦屏、三穗等地區則為低值中心，低于2 100（MJ·mm）／（hm2·h）。

圖4 貴州省多年平均降雨侵蝕力空間分布

圖5 貴州省干季月份多年平均降雨侵蝕力等值線

2．3．3 貴州省降雨侵蝕力分布的年際變化 在計算各站點多年平均降雨侵蝕力的基礎上，又分別計算了貴州省1956—1979年和1980—2000年兩個時段的年均降雨侵蝕力，并進行插值，繪制出降雨侵蝕力等值線圖（圖7—8）。對比圖7、圖8可以看出，貴州省降雨侵蝕力及其分布趨勢在兩個時段內未發生顯著變化，高值中心和低值中心位置基本一致，但數值大小有一定變化。將1956—1979年與1980—2000年兩時段的R值進行差值計算（后一段減前一段），再繪制變差等值線圖，如圖9所示。結果發現，相對于1956—1979年，1980—2000年間降雨侵蝕力值多保持穩定或有所降低；但在個別地區，降雨侵蝕力明顯增高。如盤縣、錦屏、余慶等地降雨侵蝕力都有所增高，在荔波西部地區R值增大幅度在280（MJ·mm）／（hm2·h）以上。

圖6 貴州省濕季月份多年平均降雨侵蝕力等值線

圖7 貴州省1956－1979年均降雨侵蝕力等值線

圖8 貴州省1980－2000年均降雨侵蝕力等值線

圖9 貴州省1956－1979年和1980－2000年兩時段多年平均降雨侵蝕力變化

3 結 論

本文基于貴州省羅甸小區的實測降雨過程資料、泥沙資料以及264個站點的雨量資料，計算了貴州省降雨侵蝕力并對其時空分布規律進行了分析，主要結論有：

（1）通過對貴州省羅甸小區1992—1997年32場降雨過程資料進行統計分析，得出貴州省降雨侵蝕力的最佳指標組合為EI30。

（2）得出貴州省降雨侵蝕力的簡易算法為：R＝2．0354P1．2159＋45．5649，擬合方程決定系數達0．724。

（3）采用264個站點1956—2000年的雨量資料，代入擬合方程并修訂后發現，貴州省降雨侵蝕力的空間變化呈現由南向北遞減的趨勢。其多年平均降雨 侵 蝕 力 變 化 范 圍 為 3 051．0～8 801．6（MJ·mm）／（hm2·h），平 均 值 范 圍 為 5 299．8（MJ·mm）／（hm2·h）。說明該省區的降雨對土壤具有極強的潛在侵蝕作用。此外，降雨侵蝕力存在明顯的年內變化，干濕季差異顯著，濕季4個月的降雨侵蝕力約占全年降雨侵蝕力的57%。

（4）從全省來看，近20a來貴州省的降雨侵蝕力變化不明顯，但個別地區的降雨侵蝕力有較大幅度的增加，應注重水土流失監測工作的落實與充分開展。

［1］ 章文波，謝云，劉寶元．利用日雨量計算降雨侵蝕力的方法研究［J］．地理科學，2002，22（6）：705－711．

［2］ Wischmeier W H．Predicting Rainfall－erosion Losses from Cropland East of the Rocky Mountains［M］．USDA Agricultural Handbook，No．292，1965．

［3］ 吳素業．安徽大別山區降雨侵蝕力指標的研究［J］．中國水土保持，1992（2）：36－37．

［4］ 何錫君，呂振平，楊軒，等．浙江省降雨侵蝕力時空分布規律分析［J］．水土保持研究，2010，17（6）：31－34．

［5］ 章文波，付金生．不同類型雨量資料估算降雨侵蝕力［J］．資源科學，2003，25（1）：35－41．

［6］ 馬良，左長清，孫勐，等．山東省降雨侵蝕力空間分布特征及簡易方程的研究［J］．水土保持研究，2010，17（2）：28－31．

［7］ 寧麗丹，石輝．利用日降雨量資料估算西南地區的降雨侵蝕力［J］．水土保持研究，2003，10（4）：183－186．

［8］ 范建容，嚴冬，郭祥．GIS支持下的長江上游降雨侵蝕力時空分布特征分析［J］．水土保持研究，2010，17（1）：92－96．

［9］ 劉正佳，劉占仁，王絲絲，等．基于日降雨的沂蒙山區降雨侵蝕力時空變化研究［J］．水土保持研究，2012，19（2）：34－38．

［10］ 張黎明，林金石，于東升，等．我國南方地區降雨侵蝕力指標R的建立研究：以江西鷹潭地區為例［J］．水土保持研究，2011，18（5）：1－4．

［11］ 王萬中，焦菊英，郝小品，等．中國降雨侵蝕力R值的計算與分布（Ⅰ）［J］．水土保持學報，1995，9（4）：5－18．

［12］ 章文波，謝云，劉寶元．中國降雨侵蝕力空間變化特征［J］．山地學報，2003，21（1）：33－40．

［13］ 周忠發，安裕倫．貴州省水土流失遙感現狀調查及空間變化分析［J］．水土保持通報，2000，20（6）：23－25，41．

［14］ Yu B．Rainfall erosivity and its estimation for Australia′s tropics［J］．Australian Journal of Soil Research，1998，36（1）：143－165．

［15］ Yu B，Rosewell C J．A robust estimator of the R－factor for the Universal Soil Loss Equation［J］．Transactions of the ASAE，1996，39（2）：559－561．

［16］ 楊子生．滇東北山區坡耕地土壤流失方程研究［J］．水土保持通報，1999，19（1）：4－12．

［17］ 許月卿，周巧富，李雙成．貴州省降雨侵蝕力時空分布規律分析［J］．水土保持通報，2005，25（4）：11－14．

［18］ 張憲奎，許靖華，盧秀琴，等．黑龍江省土壤流失方程的研究［J］．水土保持通報，1992，12（4）：1－9，18．

［19］ 周伏建，陳明華，林福興，等．福建省降雨侵蝕力指標R值［J］．水土保持學報，1995，9（1）：13－18．

［20］ Brown L C，Foster G R．Storm erosivity using idealized intensity distributions［J］．Transactions of the ASAE，1987，30（2）：379－386．