崩岗是我国南方丘陵地区普遍存在的严重水土流失类型，它是由水力与重力联合作用形成的特殊地貌类型。崩岗侵蚀产生的大量泥沙会掩埋良田、阻塞河道，对下游人民的生命安全及生产活动造成严重威胁，并造成巨额经济损失^[1]。我国崩岗主要分布于湖北、湖南、安徽、江西、福建、广东、广西等7个省区。2005年水利部调查结果显示，南方7省区崩岗总数共24万个，土壤侵蚀面积达到1 200多km²，年侵蚀量6 000万t，对山下影响范围近200万hm^2[2]。

“崩岗”一词由曾昭璇^①在1960年率先引入到地貌学，它揭示了崩岗的侵蚀方式和侵蚀后的地貌特征，具有发生学与形成学的双重含义^[3-4]。学者们50多年的研究，集中在地质母岩、降雨、植被、坡向等诸多因素对崩岗侵蚀的影响^[5]。

① 曾昭璇.地形学原理(第一册)[M].华南师范学院地理丛书之二, 1960.

许多研究指出崩岗侵蚀的发育具有坡向选择性，但坡向划分方法差异较大，如南北坡崩岗分布研究中，南坡角度有的是8个坡向之一(180°±22.5°)^[6-8]，有的是16个坡向之一(180°±11.25°)^[9]。阴阳坡崩岗分布研究中，南坡与阳坡存在3种关系：南坡等于阳坡^[10]、阳坡包括南坡^{[1, 11-14]}、南坡包括阳坡^[15]。崩岗的坡向确定方法也存在差异，有的^[13-14]用GIS方法统计坡向，大部分用人工方法。这些混乱导致学者们的研究结论虽表述一致，如崩岗分布数量在南坡^②多于北坡^{[1, 6-11, 13]}，阳坡多于阴坡^{[1, 10-13]}；崩岗面积南坡多于北坡^{[6-7, 12]}，阳坡多于阴坡^[12]，但无法准确定量崩岗在不同坡向的分布规律，甚至会质疑崩岗坡向选择性结论的通用性和适用性^[16]。

② 南坡，也称为向南坡或南向坡，阳坡也被称为向阳坡。

不同坡向提取方法、不同坡向划分标准对崩岗分布的坡向选择性有无影响？若已有研究使用相同坡向标准和统计方法，崩岗坡向选择性结论是否成立？不同坡向崩岗分布的定量关系如何？针对上述问题，本文开展试验研究，期望了解坡向的不同提取和划分方法对崩岗坡向研究的适宜性，验证崩岗坡向选择性的结论，明确崩岗坡向选择性的通用性和适用性，定量计算不同坡向的崩岗分布。研究结果将更加深入地揭示崩岗侵蚀发育的规律，可用于指导崩岗侵蚀的防治。

崩岗的坡向角度通过2种方法计算得到。人工计算法是在纸质地形图上根据崩岗所处坡面的主坡向或所属山脉的走势确定崩岗的主方向，再根据水流的方向，确定崩岗的精确方向，手动测量方位角；GIS方法是利用ArcGIS软件通过面积加权确定每个崩岗的坡向角度数值。同时利用GIS方法计算油田镇域的土地坡向角度。

将得到的坡向角度，从0°±22.5°开始，顺时针每45°划分1个坡向，依次为北、东北、东、东南、南、西南、西、西北，得到8个方位(即8分法)，该方法与ArcGIS系统中的坡向划分方法一致。在8分法基础上，一些研究^[13-14]也采用3分法，将东南、南、西南称为阳坡，西北、北、东北称为阴坡，将东、西坡称为半阳坡。在水土保持监测和调查中常用坡向阴阳2分法，根据太阳入射角将北半球的东南坡、南坡、西南坡、西坡称为阳坡，将西北坡、北坡、东北坡、东坡称为阴坡^[17]，坡向的8分法、3分法和2分法的关系详见图 1。

图 1(Figure 1)

图 1 坡向8分法和2分法的对应关系 Figure 1 Corresponding relationship between 8-slope-aspect and 2-slope-aspect

数据来源包含2部分：1)1978年广东省梅州市五华县油田镇1:1万地形图和油田镇边界；2)已发表文献中的各坡向崩岗数据。

五华县位于广东东北部，韩江上游，是全国崩岗分布最密集的地区；油田镇位于五华县中东部，与兴宁市水口镇接壤。E 115°48′48″~115°55′16″，N 23°51′39″~23°58′37″，全镇土地总面积85.6 km²，土壤母质主要是深厚的花岗岩风化物。

首先，借助ArcGIS软件数字化油田镇地形图的等高线、崩岗边界、油田镇边界等要素，根据3种坡向划分方法，计算并统计油田镇各坡向的土地面积、每个崩岗的坡向和面积，并制作崩岗的坡向分布图(图 2)。由于地形图比例尺所限，面积较小的崩岗在图上仅表示为一条线，无法统计面积，因此本文统计的崩岗不包含较小的崩岗。然后，将收集的已有崩岗坡向文献根据前文确定的坡向划分方法，重新统计各坡向崩岗的数量和面积，并计算数量和面积百分比。已有文献共计11篇，研究区分布在广东、福建和江西3省。所有数据使用Excel计算均值和标准差，用SPSS13.0进行方差分析，阴阳坡用t检验，方差分析中3分法用LSD事后检验，8分法用Tukey事后检验。

图 2(Figure 2)

NE: Northeast E: East SE: Southeast S: South SW: Southwest W: West NW: Northwest 图 2 基于2分法(a)、3分法(b)和8分法(c)的油田镇崩岗坡向分布图 Figure 2 Distribution of collapsed gullies by 2-slope-aspect(a), 3-slope-aspect (b), and 8-slope-aspect (c) in Youtian town

根据1:1万地形图利用ArcGIS计算得到，油田镇土地总面积为8 560.49 hm²，崩岗共计293个，崩岗总面积89.79 hm²，崩岗占全镇土地总面积的1.05%，平均每km² 3.4个崩岗，崩岗面积最小305 m²，最大17 065 m²，平均每个崩岗3 064 m²。根据油田镇各坡向面积分布统计(图 2)：2分法中阳坡面积占44.2%，阴坡占55.8%；3分法中阳坡占31.2%，阴坡占45.5%；8分法中北坡面积最大，占19.4%，南坡仅占10.3%。因此油田镇阳坡土地面积大于阴坡，北坡土地面积大于南坡。

比较GIS与人工方法计算的油田镇每个崩岗的坡向角度的误差(表 1)，结果表明，2种方法计算的崩岗坡向中有79个崩岗坡向存在差异，总面积25.11 hm²，占崩岗总数量的27%，占崩岗总面积的28%。坡向角度误差＜30°的崩岗共44个，占崩岗总数的15.0%，占总面积的15.8%，约占总误差数量和面积误差的均值的56%。误差在30°~60°之间的崩岗占崩岗总数的8.2%、占崩岗总面积的9.6%，约占总误差的32%。误差＞60°的崩岗占崩岗总数的3.8%，占崩岗总面积的2.6%，约占总误差的12%。因此，坡向角度误差主要在30°以内。产生误差的崩岗有如下特征：范围大、形状细长且沟道有转折、内部形态破碎、沟头发育越过山脊、处于坡面局部凹地等。GIS方法统计崩岗内部坡向分布，复杂的内部形态可能影响主坡向；而人工目视判读先根据山脉走势确定主坡向，再依水流方向确定精细走向，避免内部变化对主坡向的影响。这是2种方法的计算过程产生的误差。

表 1(Table 1)

表 1 GIS与人工方法计算崩岗坡向的角度差异 Table 1 Angle error of slope aspect of collapsed gully by manual and GIS methods 坡向角度误差范围 Angle error interval of slope aspect θ/(°) 误差崩岗数量 Number of collapsed gullies with error 误差崩岗面积 Area of collapsed gullies with error/m2 占崩岗总个数的比例 Percentage in total number of the collapsed gullies/% 占崩岗总面积的比例 Percentage in total area of the collapsed gullies/% θ≤10 14 38 767 4.8 4.3 10＜θ≤20 15 43 213 5.5 4.8 20＜θ≤30 15 59 719 5.1 6.7 30＜θ≤40 8 26 337 2.7 2.9 40＜θ≤50 7 29 758 2.4 3.3 50＜θ≤60 9 30 105 3.1 3.4 60＜θ≤70 2 4 156 0.7 0.5 70＜θ≤80 2 6 494 0.7 0.7 80＜θ≤90 0 0 0.0 0.0 90＜θ＜180 7 12 607 2.4 1.4 总计Sum 79 251 156 27.0 28.0

表 1 GIS与人工方法计算崩岗坡向的角度差异 Table 1 Angle error of slope aspect of collapsed gully by manual and GIS methods

五华县油田镇各坡向崩岗面积比例的统计结果(表 2)，GIS和人工2种方法计算的各坡向崩岗面积存在一定的误差，2分法中GIS计算的阳坡崩岗面积为65.29 hm²，人工计算是60.55 hm²，相对误差7%，阴坡相对误差17%；在3分法中阳坡相对误差2%，阴坡64%，半阳坡29%。但2种方法计算结果均表明分布在阳坡的崩岗面积远超阴坡，在2分法中阳坡崩岗面积占70.1%，3分法中阳坡占48.1%，阴坡占15.0%。8分法中，2种方法统计的各坡向(北坡、东南坡除外)崩岗面积相对误差在24%~67%之间，多因素方差分析显示，2种方法间、8个坡向间及交互作用下崩岗面积差异均不显著。但2种方法计算的南坡山崩岗面积均明显多于北坡，南坡崩岗面积均值占15.2%，北坡占1.1%。2种方法计算崩岗面积最大的坡向分别是西坡(GIS)和西南坡(人工)，面积最小的均是北坡。根据五华县油田镇各坡向崩岗数量比例统计结果(表 2)，虽然GIS和人工2种方法计算的各坡向崩岗数量存在一定误差，2分法中GIS计算阳坡崩岗数量196个，人工173个，相对误差13%；阴坡相对误差21%；3分法中阳坡相对误差为10%，阴坡为47%，半阳坡为12%。但2种方法的计算结果一致表明阳坡崩岗数量多于阴坡，2分法中阳坡崩岗数量占63.0%，3分法中阳坡占43.0%，阴坡占18.8%。8分法中，2种方法统计的各坡向(北坡除外)崩岗数量相对误差在5%~52%之间。虽然2种方法计算的崩岗数量最多的分别是西坡(67个GIS)、西南和西坡(各56个人工)，但南坡崩岗数量均明显多于北坡，南坡崩岗数量均值占12.5%，北坡占比最小，仅占1.6%。

表 2(Table 2)

表 2 五华油田镇及已有文献各坡向崩岗的面积和数量比例 Table 2 Area and quantity percentage of collapsed gully on different slope aspects in Youtian town of Wuhua county and all former studies % 作者 Author 8分法 8-slope-aspect 3分法 3-slope-aspect 2分法 2-slope-aspect 北 North 东北 NE 东 East 东南 SE 南 South 西南 SW 西 West 西北 NW 阴坡 Shady 阳坡 Sunny 半阳坡 Semi-sunny 阴坡 Shady 阳坡 Sunny A+ 0.0 4.0 17.1 12.7 17.0 17.8 25.2 6.2 10.2 47.5 42.3 27.3 72.7 B+ 2.1 8.0 12.8 12.7 13.4 22.6 18.8 9.6 19.7 48.7 31.6 32.5 67.5 C+ 5.8 0.5 5.8 9.5 34.4 36.5 5.3 2.2 8.5 80.4 11.1 14.3 85.7 F+ — 3.5 — — — 25.6 — — — — — — — L+ 4.2 8.9 6.2 8.7 27.5 20.9 15.8 7.8 20.9 57.1 22.0 27.1 72.9 M+ 12.4 11.3 14.6 11.9 10.1 9.6 15.7 14.4 38.1 31.6 30.3 52.7 47.3 均值+ Mean+ 4.7 b 6.5 b 11.3 ab 11.1 ab 21.3 a 21.5 a 16.2 ab 8.0b 19.5 b 53.1 a 27.5 b 30.8 b 69.2 a 标准差+ SD+ 4.3 4.3 5.1 1.9 9.3 9.8 7.2 4.5 11.8 17.8 11.7 14.0 14.0 A 0.0 5.4 18.8 13.0 14.0 17.0 22.9 8.9 14.3 44.0 41.7 33.1 66.9 B 3.1 9.2 17.8 9.9 10.9 19.1 19.1 10.9 23.2 39.9 36.9 41.0 59.0 D 1.8 0.0 10.0 9.0 56.1 11.8 11.3 0.0 1.8 76.9 21.3 11.8 88.2 E 0.0 0.0 2.5 30.3 1.7 50.4 10.9 4.2 4.2 82.4 13.4 6.7 93.3 F 0.0 0.0 0.0 31.0 15.0 30.0 17.0 7.0 7.0 76.0 17.0 7.0 93.0 G 10.3 8.8 15.6 12.2 20.4 10.6 14.8 7.2 26.3 43.2 30.4 42.0 58.0 H 10.3 8.8 15.6 12.2 20.4 10.6 14.8 7.2 26.3 43.2 30.4 42.0 58.0 J 9.8 8.3 15.3 13.1 22.6 10.4 13.6 6.9 25.0 46.2 28.8 40.3 59.7 K 8.8 6.5 5.7 12.2 24.9 17.6 16.8 7.6 22.9 54.7 22.5 28.6 71.5 L 7.6 7.6 5.3 11.5 23.7 19.1 15.3 9.9 25.1 54.3 20.6 30.4 69.6 M 11.4 12.4 17.2 9.9 10.9 11.1 14.6 12.5 36.3 31.9 31.8 53.5 46.5 N — — — — — — — — — — — 7.6 92.4 均值 Mean 5.7 6.1 11.3 14.9 20.1 18.9 15.5 7.5 19.3 b 53.9 a 26.8 b 28.7 b 71.3 a 标准差 SD 4.7 4.3 6.8 7.9 13.8 12.0 3.4 3.3 10.9 17.0 8.6 16.4 16.4 注：带+表示崩岗的面积比例，不含+表示崩岗的数量比例。A:本文GIS；B:本文人工；C：广东省科考队；D：沈家安等；E：姚清尹等；F：刘瑞华；G：丁光敏；H：阮伏水；J：林敬兰等；K：陈志彪等；L：方芸芸等；M：宋月君；N:吴克刚等。两列小写字母不同代表差异显著(P＜0.05)。Notes: The data with+ refer to the area percentage, those without+ refer to the quantity percentage. A: This article by GIS method. B: This article by artificial method. C: Guangdong Science Exploring Group. D: SHEN Jiaan et al. E: YAO Qingyin et al. F: LIU Ruihua; G: DING Guangmin et al. H: RUAN Fushui. J: LIN Jinglan et al. K: CHEN Zhibiao et al. L: FANG Yunyun et al. M: SONG Yuejun; N: WU Kegang et al. Different lower case letters within the different columns mean significant difference (P＜0.05).

表 2 五华油田镇及已有文献各坡向崩岗的面积和数量比例 Table 2 Area and quantity percentage of collapsed gully on different slope aspects in Youtian town of Wuhua county and all former studies

以本文的坡向划分方法重新统计已有研究中各坡向崩岗面积数据(表 2)，统计分析结果显示不论是2分法还是3分法，阳坡崩岗面积均显著大于阴坡。2分法中阳坡崩岗面积比例均值为69.2%，阴坡为30.8%(t=-4.345，D_f=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分法还是3分法，阳坡崩岗数量均显著多于阴坡。2分法中阳坡崩岗数量比例均值为71.3%，阴坡为28.7%(t=-6.370，D_f=22，P＜0.05)；3分法中阳坡崩岗数量比例均值为53.9%，阴坡为19.3%(F=22.592，P＜0.05)。8分法中，各坡向崩岗数量比例均值是南坡最大，占20.1%；其次是西南、西、东南和东，比例为10%~20%；然后是西北和东北坡，比例＜10%；北坡崩岗数量比例均值最小，仅为5.7%。

1) 虽然GIS和人工2种方法计算的五华县油田镇各个坡向崩岗的数量和面积比例均存在一定误差，但2种方法得到的崩岗分布规律基本一致：阳坡崩岗多于阴坡，南坡崩岗多于北坡，在坡向2分法中阳坡崩岗占66.5%，阴坡占33.5%；在3分法中阳坡占45.0%，阴坡占16.9%；在8分法中南坡占13.8%，北坡占1.3%。油田镇土地面积不论是2分法还是3分法均是阳坡土地面积小于阴坡，8分法中南坡土地面积小于北坡；因此五华油田崩岗分布存在坡向选择性，这一规律与崩岗坡向的提取方法、划分方法以及区域土地坡向分布无关。

2) 以统一的坡向划分方法重新统计已有文献各坡向的崩岗数据，结果表明崩岗分布存在坡向选择性，阳坡崩岗数量和面积均显著多于阴坡，南坡和西南坡崩岗数量和面积均多于北坡，北坡数量和面积均最少。崩岗分布的比例，在2分法中阳坡占70.3%，阴坡占29.7%；3分法中阳坡占53.5%，阴坡占19.4%；8分法中，南坡占20.7%，西南坡占20.2%，北坡占5.2%。这说明坡向是影响崩岗分布的重要因素，本研究实现了崩岗坡向选择性的数值定量。

坡向对崩岗侵蚀的影响是间接的、隐性的，并且非常复杂。南坡接受更多的太阳辐射^[18]和降水^[19]，从而影响局地小气候，造成母岩、土壤和植被等发生变化，引起土壤可蚀性、抗蚀性^[20-21]和剥蚀率^[22]的变化，最终导致崩岗侵蚀更多分布在阳坡和南坡。崩岗侵蚀的发生发育是多个环境因子共同作用的结果。降雨、母岩、植被等因素均在大尺度地带性地影响崩岗侵蚀的分布，而坡向对侵蚀的影响更多是在点的尺度得到体现。

本研究发现面积较大的崩岗形态复杂多变，GIS计算时容易产生较大误差，因此地形复杂的崩岗，需开展专项研究改进坡向提取方法。暂未查到安徽、湖南、湖北、广西等4个省区的崩岗侵蚀分布资料，研究结论有待进一步验证与完善。