溫美麗，陳 瑜，何小武，楊 龍，周 晴

(1.廣東省地理空間信息技術(shù)與應用公共實驗室,廣州地理研究所,510070,廣州；2.江西農(nóng)業(yè)大學國土資源與環(huán)境學院,330045,南昌)

崩崗是我國南方丘陵地區(qū)普遍存在的嚴重水土流失類型,它是由水力與重力聯(lián)合作用形成的特殊地貌類型。崩崗侵蝕產(chǎn)生的大量泥沙會掩埋良田、阻塞河道,對下游人民的生命安全及生產(chǎn)活動造成嚴重威脅,并造成巨額經(jīng)濟損失[1]。我國崩崗主要分布于湖北、湖南、安徽、江西、福建、廣東、廣西等7個省區(qū)。2005年水利部調(diào)查結(jié)果顯示,南方7省區(qū)崩崗總數(shù)共24萬個,土壤侵蝕面積達到1 200多km2,年侵蝕量6 000萬t,對山下影響范圍近200萬hm2[2]。

“崩崗”一詞由曾昭璇①曾昭璇.地形學原理(第一冊)[M].華南師范學院地理叢書之二,1960.在1960年率先引入到地貌學,它揭示了崩崗的侵蝕方式和侵蝕后的地貌特征,具有發(fā)生學與形成學的雙重含義[3-4]。 學者們50多年的研究,集中在地質(zhì)母巖、降雨、植被、坡向等諸多因素對崩崗侵蝕的影響[5]。

許多研究指出崩崗侵蝕的發(fā)育具有坡向選擇性,但坡向劃分方法差異較大,如南北坡崩崗分布研究中,南坡角度有的是 8個坡向之一(180°±22.5°)[6-8],有 的 是 16個 坡 向 之 一(180°±11.25°)[9]。陰陽坡崩崗分布研究中,南坡與陽坡存在3種關(guān)系:南坡等于陽坡[10]、陽坡包括南坡[1,11-14]、南坡包括陽坡[15]。崩崗的坡向確定方法也存在差異,有的[13-14]用GIS方法統(tǒng)計坡向,大部分用人工方法。這些混亂導致學者們的研究結(jié)論雖表述一致,如崩崗分布數(shù)量在南坡②南坡,也稱為向南坡或南向坡,陽坡也被稱為向陽坡。多于北坡[1,6-11,13],陽坡多于陰坡[1,10-13]；崩崗面積南坡多于北坡[6-7,12],陽坡多于陰坡[12],但無法準確定量崩崗在不同坡向的分布規(guī)律,甚至會質(zhì)疑崩崗坡向選擇性結(jié)論的通用性和適用性[16]。

不同坡向提取方法、不同坡向劃分標準對崩崗分布的坡向選擇性有無影響？若已有研究使用相同坡向標準和統(tǒng)計方法,崩崗坡向選擇性結(jié)論是否成立？不同坡向崩崗分布的定量關(guān)系如何？針對上述問題,本文開展試驗研究,期望了解坡向的不同提取和劃分方法對崩崗坡向研究的適宜性,驗證崩崗坡向選擇性的結(jié)論,明確崩崗坡向選擇性的通用性和適用性,定量計算不同坡向的崩崗分布。研究結(jié)果將更加深入地揭示崩崗侵蝕發(fā)育的規(guī)律,可用于指導崩崗侵蝕的防治。

1 數(shù)據(jù)與方法

1.1 坡向的計算與分類

崩崗的坡向角度通過2種方法計算得到。人工計算法是在紙質(zhì)地形圖上根據(jù)崩崗所處坡面的主坡向或所屬山脈的走勢確定崩崗的主方向,再根據(jù)水流的方向,確定崩崗的精確方向,手動測量方位角；GIS方法是利用ArcGIS軟件通過面積加權(quán)確定每個崩崗的坡向角度數(shù)值。同時利用GIS方法計算油田鎮(zhèn)域的土地坡向角度。

將得到的坡向角度,從0°±22.5°開始,順時針每45°劃分1個坡向,依次為北、東北、東、東南、南、西南、西、西北,得到8個方位(即8分法),該方法與ArcGIS系統(tǒng)中的坡向劃分方法一致。在8分法基礎(chǔ)上,一些研究[13-14]也采用3分法,將東南、南、西南稱為陽坡,西北、北、東北稱為陰坡,將東、西坡稱為半陽坡。在水土保持監(jiān)測和調(diào)查中常用坡向陰陽2分法,根據(jù)太陽入射角將北半球的東南坡、南坡、西南坡、西坡稱為陽坡,將西北坡、北坡、東北坡、東坡稱為陰坡[17],坡向的8分法、3分法和2分法的關(guān)系詳見圖1。

圖1 坡向8分法和2分法的對應關(guān)系Fig.1 Corresponding relationship between 8-slopeaspect and 2-slope-aspect

1.2 數(shù)據(jù)來源及處理方法

數(shù)據(jù)來源包含2部分:1)1978年廣東省梅州市五華縣油田鎮(zhèn)1∶1萬地形圖和油田鎮(zhèn)邊界；2)已發(fā)表文獻中的各坡向崩崗數(shù)據(jù)。

五華縣位于廣東東北部,韓江上游,是全國崩崗分布最密集的地區(qū)；油田鎮(zhèn)位于五華縣中東部,與興寧市水口鎮(zhèn)接壤。 E 115°48′48″～115°55′16″,N 23°51′39″～23°58′37″,全鎮(zhèn)土地總面積 85.6 km2,土壤母質(zhì)主要是深厚的花崗巖風化物。

首先,借助ArcGIS軟件數(shù)字化油田鎮(zhèn)地形圖的等高線、崩崗邊界、油田鎮(zhèn)邊界等要素,根據(jù)3種坡向劃分方法,計算并統(tǒng)計油田鎮(zhèn)各坡向的土地面積、每個崩崗的坡向和面積,并制作崩崗的坡向分布圖(圖2)。由于地形圖比例尺所限,面積較小的崩崗在圖上僅表示為一條線,無法統(tǒng)計面積,因此本文統(tǒng)計的崩崗不包含較小的崩崗。然后,將收集的已有崩崗坡向文獻根據(jù)前文確定的坡向劃分方法,重新統(tǒng)計各坡向崩崗的數(shù)量和面積,并計算數(shù)量和面積百分比。已有文獻共計11篇,研究區(qū)分布在廣東、福建和江西3省。所有數(shù)據(jù)使用Excel計算均值和標準差,用SPSS13.0進行方差分析,陰陽坡用t檢驗,方差分析中3分法用LSD事后檢驗,8分法用Tukey事后檢驗。

2 結(jié)果與分析

2.1 油田鎮(zhèn)的坡向及崩崗分布

根據(jù)1∶1萬地形圖利用ArcGIS計算得到,油田鎮(zhèn)土地總面積為8 560.49 hm2,崩崗共計293個,崩崗總面積89.79 hm2,崩崗占全鎮(zhèn)土地總面積的1.05%,平均每km23.4個崩崗,崩崗面積最小305 m2,最大17 065 m2,平均每個崩崗3 064 m2。根據(jù)油田鎮(zhèn)各坡向面積分布統(tǒng)計(圖2):2分法中陽坡面積占 44.2%,陰坡占 55.8%；3分法中陽坡占31.2%,陰坡占45.5%；8分法中北坡面積最大,占19.4%,南坡僅占10.3%。因此油田鎮(zhèn)陽坡土地面積大于陰坡,北坡土地面積大于南坡。

2.2 GIS與人工方法計算崩崗坡向角度的差異分析

比較GIS與人工方法計算的油田鎮(zhèn)每個崩崗的坡向角度的誤差(表1),結(jié)果表明,2種方法計算的崩崗坡向中有79個崩崗坡向存在差異,總面積25.11 hm2,占崩崗總數(shù)量的27%,占崩崗總面積的28%。坡向角度誤差＜30°的崩崗共44個,占崩崗總數(shù)的15.0%,占總面積的15.8%,約占總誤差數(shù)量和面積誤差的均值的56%。誤差在30°～60°之間的崩崗占崩崗總數(shù)的8.2%、占崩崗總面積的9.6%,約占總誤差的32%。誤差＞60°的崩崗占崩崗總數(shù)的3.8%,占崩崗總面積的2.6%,約占總誤差的12%。因此,坡向角度誤差主要在30°以內(nèi)。產(chǎn)生誤差的崩崗有如下特征:范圍大、形狀細長且溝道有轉(zhuǎn)折、內(nèi)部形態(tài)破碎、溝頭發(fā)育越過山脊、處于坡面局部凹地等。GIS方法統(tǒng)計崩崗內(nèi)部坡向分布,復雜的內(nèi)部形態(tài)可能影響主坡向；而人工目視判讀先根據(jù)山脈走勢確定主坡向,再依水流方向確定精細走向,避免內(nèi)部變化對主坡向的影響。這是2種方法的計算過程產(chǎn)生的誤差。

表1 GIS與人工方法計算崩崗坡向的角度差異____________________Tab.1 Angle error of slope aspect of collapsed gully by manual and GIS methods

2.3 GIS與人工方法下廣東五華油田鎮(zhèn)的崩崗分布比較

五華縣油田鎮(zhèn)各坡向崩崗面積比例的統(tǒng)計結(jié)果(表2),GIS和人工2種方法計算的各坡向崩崗面積存在一定的誤差,2分法中GIS計算的陽坡崩崗面積為65.29 hm2,人工計算是60.55 hm2,相對誤差7%,陰坡相對誤差17%；在3分法中陽坡相對誤差2%,陰坡64%,半陽坡29%。但2種方法計算結(jié)果均表明分布在陽坡的崩崗面積遠超陰坡,在2分法中陽坡崩崗面積占70.1%,3分法中陽坡占48.1%,陰坡占15.0%。8分法中,2種方法統(tǒng)計的各坡向(北坡、東南坡除外)崩崗面積相對誤差在24%～67%之間,多因素方差分析顯示,2種方法間、8個坡向間及交互作用下崩崗面積差異均不顯著。但2種方法計算的南坡山崩崗面積均明顯多于北坡,南坡崩崗面積均值占15.2%,北坡占1.1%。2種方法計算崩崗面積最大的坡向分別是西坡(GIS)和西南坡(人工),面積最小的均是北坡。根據(jù)五華縣油田鎮(zhèn)各坡向崩崗數(shù)量比例統(tǒng)計結(jié)果(表2),雖然GIS和人工2種方法計算的各坡向崩崗數(shù)量存在一定誤差,2分法中GIS計算陽坡崩崗數(shù)量196個,人工173個,相對誤差13%；陰坡相對誤差21%；3分法中陽坡相對誤差為10%,陰坡為47%,半陽坡為12%。但2種方法的計算結(jié)果一致表明陽坡崩崗數(shù)量多于陰坡,2分法中陽坡崩崗數(shù)量占63.0%,3分法中陽坡占43.0%,陰坡占18.8%。8分法中,2種方法統(tǒng)計的各坡向(北坡除外)崩崗數(shù)量相對誤差在5%～52%之間。雖然2種方法計算的崩崗數(shù)量最多的分別是西坡(67個GIS)、西南和西坡(各56個人工),但南坡崩崗數(shù)量均明顯多于北坡,南坡崩崗數(shù)量均值占12.5%,北坡占比最小,僅占1.6%。

2.4 已有文獻各坡向崩崗數(shù)據(jù)的再分析

以本文的坡向劃分方法重新統(tǒng)計已有研究中各坡向崩崗面積數(shù)據(jù)(表2),統(tǒng)計分析結(jié)果顯示不論是2分法還是3分法,陽坡崩崗面積均顯著大于陰坡。2分法中陽坡崩崗面積比例均值為69.2%,陰坡為30.8%(t=-4.345,Df=8,P=0.002)；3分法中陽坡崩崗面積比例均值為53.1%,陰坡為19.5%(F=7.783,P=0.007)。8分法中西南和南坡的崩崗面積比例均值最大,分別占21.5%和21.3%,其次是西、東、東南、西北和東北,北坡面積最小,僅占4.7%。方差分析顯示(F=5.586,P＜0.05),南與北、東北和西北的崩崗面積差異顯著,西南亦如此。

圖2 基于2分法(a)、3分法(b)和8分法(c)的油田鎮(zhèn)崩崗坡向分布圖Fig.2 Distribution of collapsed gullies by 2-slope-aspect(a),3-slope-aspect(b),and 8-slope-aspect(c)in Youtian town

表2 五華油田鎮(zhèn)及已有文獻各坡向崩崗的面積和數(shù)量比例Tab.2 Area and quantity percentage of collapsed gully on different slope aspects in Youtian town of Wuhua county and all former studies%

重新統(tǒng)計已有的研究中各坡向崩崗數(shù)量比例(表2),統(tǒng)計分析結(jié)果表明不論是2分法還是3分法,陽坡崩崗數(shù)量均顯著多于陰坡。2分法中陽坡崩崗數(shù)量比例均值為71.3%,陰坡為28.7%(t=-6.370,Df=22,P＜0.05)；3分法中陽坡崩崗數(shù)量比例均值為53.9%,陰坡為19.3%(F=22.592,P＜0.05)。8分法中,各坡向崩崗數(shù)量比例均值是南坡最大,占20.1%；其次是西南、西、東南和東,比例為10% ～20%；然后是西北和東北坡,比例 ＜10%；北坡崩崗數(shù)量比例均值最小,僅為5.7%。

3 結(jié)論與討論

3.1 結(jié)論

1)雖然GIS和人工2種方法計算的五華縣油田鎮(zhèn)各個坡向崩崗的數(shù)量和面積比例均存在一定誤差,但2種方法得到的崩崗分布規(guī)律基本一致:陽坡崩崗多于陰坡,南坡崩崗多于北坡,在坡向2分法中陽坡崩崗占66.5%,陰坡占33.5%；在3分法中陽坡占45.0%,陰坡占16.9%；在8分法中南坡占13.8%,北坡占1.3%。油田鎮(zhèn)土地面積不論是2分法還是3分法均是陽坡土地面積小于陰坡,8分法中南坡土地面積小于北坡；因此五華油田崩崗分布存在坡向選擇性,這一規(guī)律與崩崗坡向的提取方法、劃分方法以及區(qū)域土地坡向分布無關(guān)。

2)以統(tǒng)一的坡向劃分方法重新統(tǒng)計已有文獻各坡向的崩崗數(shù)據(jù),結(jié)果表明崩崗分布存在坡向選擇性,陽坡崩崗數(shù)量和面積均顯著多于陰坡,南坡和西南坡崩崗數(shù)量和面積均多于北坡,北坡數(shù)量和面積均最少。崩崗分布的比例,在2分法中陽坡占70.3%,陰坡占29.7%；3分法中陽坡占53.5%,陰坡占19.4%；8分法中,南坡占20.7%,西南坡占20.2%,北坡占5.2%。這說明坡向是影響崩崗分布的重要因素,本研究實現(xiàn)了崩崗坡向選擇性的數(shù)值定量。

3.2 討論

坡向?qū)Ρ缻徢治g的影響是間接的、隱性的,并且非常復雜。南坡接受更多的太陽輻射[18]和降水[19],從而影響局地小氣候,造成母巖、土壤和植被等發(fā)生變化,引起土壤可蝕性、抗蝕性[20-21]和剝蝕率[22]的變化,最終導致崩崗侵蝕更多分布在陽坡和南坡。崩崗侵蝕的發(fā)生發(fā)育是多個環(huán)境因子共同作用的結(jié)果。降雨、母巖、植被等因素均在大尺度地帶性地影響崩崗侵蝕的分布,而坡向?qū)η治g的影響更多是在點的尺度得到體現(xiàn)。

本研究發(fā)現(xiàn)面積較大的崩崗形態(tài)復雜多變,GIS計算時容易產(chǎn)生較大誤差,因此地形復雜的崩崗,需開展專項研究改進坡向提取方法。暫未查到安徽、湖南、湖北、廣西等4個省區(qū)的崩崗侵蝕分布資料,研究結(jié)論有待進一步驗證與完善。

審稿專家對本文提出重要修改意見,廣州地理研究所的梁國昭老師幫助修改全文,陳聰在數(shù)字化和制圖中給予大力幫助,在此謹表示感謝！

4 參考文獻

[1] 阮伏水.福建省崩崗侵蝕與治理模式探討[J].山地學報,2003,21(6):675.RUAN Fushui.Study on slump gully erosion and its control in Fujian province[J].Journal of Mountain Science,2003,21(6):675.

[2] 張學儉.南方崩崗的治理開發(fā)實踐與前景[J].中國水利,2010(4):17.ZHANG Xuejian.The practice and prospect of hill collapsing improving and development in southern China[J].China Water Resources,2010(4):17.

[3] 吳志峰,鐘偉青.崩崗災害地貌及其環(huán)境效應[J].生態(tài)科學,1997,16(2):91.WU Zhifeng,ZHONG Weiqing.Gravity gully erosion and its influence[J].Ecologic Science,1997,16(2):91.

[4] 史德明.我國熱帶、亞熱帶地區(qū)崩崗侵蝕的剖析[J].水土保持通報,1984,4(3):32.SHI Deming.Collapsed gully erosion in tropical and subtropical regions of China[J].Bulletin of Soil and Water Conservation,1984,4(3):32.

[5] 陳培波,劉希林.中國崩崗研究進展的文獻計量分析[J].熱帶地理,2015,35(6):895.CHEN Peibbo,LIU Xilin.Bibliometric analysis of the progress in Benggang research in China[J].Tropical Geography,2015,35(6):895.

[6] 方蕓蕓,陳志強,陳志彪.根溪小流域崩崗分布特征研究[J].吉林師范大學學報(自然科學版),2014,35(3):145.FANG Yunyun,CHEN Zhiqiang,CHEN Zhibiao.Study on distributed characteristics of slump gully in Genxi small watershed[J].Journal of Jilin Normal University(Natural Science Edition),2014,35(3):145.

[7] 劉瑞華.華南地區(qū)崩崗侵蝕災害及其防治[J].水文地質(zhì)工程地質(zhì),2004,31(198):54.LIU Ruihua.Slope disintegration and its control in south China[J].Hydrogeology&Engineering Geology,2004,31(198):54.

[8] 姚清尹,陸兆熊.花崗巖裂隙構(gòu)造及其對風化與巖體破壞的影響[M]∥中國、加拿大水土保持協(xié)作組.廣東省水土保持研究.北京:科學出版社,1989:24.YAO Qingyin,LU Zhaoxiong.Fissure tectonics of granite and its influence on weathering and rock mass destruction[M]∥University of Toronto&Guangzhou Institute of Geography.Soil erosion and land management in the granitic regions of Guangdong province,South China.Beijing:Science Press,1989:24.

[9] 陳志彪,朱鶴健,劉強,等.根溪河小流域的崩崗特征及其治理措施[J].自然災害學報,2006,15(5):83.CHEN Zhibiao,ZHU Hejian,LIU Qiang,et al.Slump gully characteristic of small watershed of Genxi river and its control measures[J].Journal of Natural Disasters,2006,15(5):83.

[10] 沈家安,姚清尹,張仲英,等.德慶縣水土流失的基本情況及其治理[M]∥中國、加拿大水土保持協(xié)作組.廣東省水土保持研究.北京:科學出版社,1989:169.SHEN Jiaan,YAO Qingyin,ZHANG Zhongying,et al.Basic situation and treatment of soil and water loss in Deqing County[M]∥University of Toronto&Guangzhou Institute of Geography.Soil erosion and land management in the granitic regions of Guangdong province,South China.Beijing:Science Press,1989:169.

[11] 丁光敏.福建省崩崗侵蝕成因及治理模式研究[J].水土保持通報,2001,21(5):10.DING Guangmin.Causes and control patterns of collapsed gully erosion in Fujian province20[J].Bulletin of Soil and Water Conservation,2001,21(5):10.

[12] 廣東省科學院丘陵山區(qū)綜合科學考察隊.廣東山區(qū)水土流失及其治理[M].廣州:廣東科技出版社,1991,86.Comprehensive Scientific Expedition of Hilly and Mountainous Areas of Guangdong Academy of Sciences.Soil erosion and its control in mountainous areas of Guangdong[M].Guangzhou:Guangdong Science Press,1991:86.

[13] 林敬蘭,黃炎和,林金石,等.福建省崩崗侵蝕的地質(zhì)地貌背景分析[J].亞熱帶水土保持,2014,26(4):1.LIN Jinglan,HUANG Yanhe,LIN Jinshi,et al.Background analysis on the geology and land form of collapse erosion in Fujian province[J].Subtropical Soil and Water Conservation,2014,26(4):1.

[14] 宋月君.江西省贛縣崩崗發(fā)育的地貌特征分析[J].亞熱帶水土保持,2015,27(3):29.SONG Yuejun.Analysis on the geo-morphological features of collapsed hill development in Gan county of Jiangxi province[J].Subtropical Soil and Water Conservation,2015,27(3):29.

[15] 吳克剛,陸兆熊,郭鼎.華南花崗巖風化殼崩崗侵蝕:以德慶縣深涌小流域為例[M]∥中國、加拿大水土保持協(xié)作組.廣東省水土保持研究.北京:科學出版社,1989:7.WU Kegang,LU Zhaoxiong,GUO Ding.Collapsed gully erosion of granite weathering crust in South China:A case study of Deqing County Deep-chung small watershed[M]∥University of Toronto&Guangzhou Institute of Geogrphy.Soil erosion and land management in the granitic regions of Guangdong province south China.Beijing:Science Press,1989:7.

[16] 劉希林,連海清.崩崗侵蝕地貌分布的海拔高程與坡向選擇性[J].水土保持通報,2011,31(4):32.LIU Xilin,LIAN Haiqing.Distribution choices of elevation and slope orientation of collapsed hills[J].Bulletin of Soil and Water Conservation,2011,31(4):32.

[17] 郭索彥.水土保持監(jiān)測理論與方法[M].北京:中國水利水電出版社,2010:183.GUO Suoyan.Theory and method of soil and water conservation monitoring[M].Beijing:China Water Conservancy and Hydropower Press,2010:183.

[18] 傅抱璞.論坡地上的太陽輻射總量[J].南京大學學報(自然科學版),1958(2):47.FU Baopu.The solar radiation falling on slope[J].Journal of Nanjing University(Natural Science),1958(2):47.

[19] 陳明,傅抱璞,于強.山區(qū)地形對暴雨的影響[J].地理學報,1995,50(3):256.CHEN Ming,FU Baopu,YU Qiang.Influence of topography on storm rainfall[J].Acta Geographica Sinica,1995,50(3):256.

[20] 張曉明,丁樹文,蔡崇法.干濕效應下崩崗區(qū)巖土抗剪強度衰減非線性分析[J].農(nóng)業(yè)工程學報,2012,28(5):241.ZHANG Xiaoming,DING Shuwen,CAIChongfa.Effects of drying and wetting on nonlinear decay of soil shear strength in slope disintegration erosion area[J].Transactions of the Chinese Society of Agricultural Engineering,2012,28(5):241.

[21] 周紅藝,李輝霞.華南活動型崩崗崩壁土體的崩解特性及其影響因素[J].水土保持學報,2017,31(1):74.ZHOU Hongyi,LI Huixia.The study on soil disintegration characteristics and its influence factors of collapsing wall in the collapsing hill erosion region of southern China[J].Journal of Soil and Water Conservation,2017,31(1):74.

[22] 王秋霞,丁樹文,鄧羽松,等.花崗巖崩崗區(qū)不同土層的侵蝕水動力學特征[J].土壤學報,2017,54(3):570.WANG Qiuxia,DING Shuwen,DENG Yusong,et al.Hydrodynamic characteristics of erosion in different soil layers in granite collapsed region[J].Acta Pedologica Sinica,2017,54(3):570.