周晓晴, 马芊红, 张科利

(北京师范大学 地理科学学部, 北京 100875)

土壤侵蚀是人类面临的限制全球经济社会可持续发展的严峻环境问题，严重的土壤侵蚀会破坏土地生产力，降低生物多样性，威胁区域生态环境安全，加剧山区贫困[1]。受自然因素和人为因素共同作用，陕西省水土流失极其严重，其土壤侵蚀面积和侵蚀强度均高于全国大部分省份[2]。第二次全国遥感调查资料显示，陕西省水土流失面积超过全省土地总面积的60%，且仍有扩大趋势，全省70%的人口和80%的耕地正遭受水土流失危害[3]。陕西省强烈的水土流失是黄河和长江河道泥沙的重要来源，另外水土流失问题不仅是制约陕西省建设成为西部强省的瓶颈，同时也是我国西部大开发战略中亟待解决的问题之一[3]。

地处陕西省北部的黄土高原生态环境脆弱，土壤侵蚀非常强烈，因而成为我国乃至世界土壤侵蚀研究的重点区域[4-5]。目前陕西全省尺度的土壤侵蚀综合研究相对较少，刘康等[6]利用GIS空间叠加技术分析了陕西省水土流失的潜在敏感性和现实敏感性，其中微度到中度敏感性分布面积较大，而改变土地利用方式、破坏植被等人类活动增大了高度现实敏感性区域的面积；李雄飞[7]利用陕西省62个水文站河流泥沙观测资料，推求流域土壤侵蚀模数，结果表明，陕西省土壤侵蚀空间分异特征明显，陕北地区土壤侵蚀十分剧烈，关中和陕南地区侵蚀强度较弱；程琳等[8]通过遥感技术、调查统计、前期积累资料得到陕西省土壤侵蚀各因子值，利用中国土壤流失方程(CSLE)估算土壤侵蚀模数，并结合GIS技术分析土壤侵蚀空间特征，表明陕西省土壤侵蚀整体上呈北高南低的变化趋势。然而，陕西省水土流失规律的已有研究对水土流失影响因子空间分布规律的对比分析尚有不足且均未采用水土流失实地调查资料，而基于地块的野外资料在反映区域水土流失中最直接、最准确、最有效。本研究以第一次水利普查中水土保持专项调查数据为依据，在陕西省抽取9条典型样带(横向7条、纵向2条)，系统分析地形因子和土地利用类型对土壤侵蚀影响的区域差异，探究水力侵蚀空间分异特征，从而为明确陕西省水土流失综合特征，制定合理的水土保持规划提供重要依据。

1 材料与方法

1.1 研究区概况

陕西省地处我国中部腹地(31°42′—39°35′N，105°29′—111°15′E)，南北长达880 km，东西相距160～490 km，面积约20.56万km2。陕西省北部和南部分别为黄土高原和秦巴山地，地势较高，中部为地势低平的关中平原。受地形与季风影响，陕西省降水时空分异明显，陕北、关中、陕南3区年均降雨量分别为279，563，840 mm，自北向南分属中温带半干旱区、暖温带半湿润区和亚热带湿润区[9]，植被类型依次为温带草原、暖温带落叶阔叶林、亚热带常绿阔叶林。陕西省开发历史悠久，土地利用类型以农地、林地、草地为主，退耕还林还草工程使陕北地区林草面积增大，关中平原适宜耕作，陕南地区林地居多[10]。复杂的自然地理条件以及对水土资源的无序利用，使陕西省水土流失面积大、侵蚀强度高、危害严重且难以治理，严重制约着该区生态环境建设与经济社会可持续发展。根据自然环境与水土流失特征的区域差异，本文选取2条纵向样带，7条横向样带，陕北、关中、陕南各3条、1条、3条(图1，表1)。

图1调查样带分布

1.2 数据来源

本研究在陕西省选取的9条样带共包含314个调查单元，1 812个地块，基于2010—2012年开展的水土保持专项调查数据，采用CSLE计算调查区域土壤侵蚀模数[11]。首先将全国划分为4级网格，按照4%的抽样密度，将控制区(5 km×5 km)中心1 km2的基本抽样单元(平原区)或与中心网格相连面积为0.2～3.0 km2的小流域(山地丘陵区)作为实际野外调查单元，土壤侵蚀微弱地区，可降低抽样密度。调查土地利用类型时，先制作野外调查底图和调查表，然后实地到野外调查单元勾绘调查图，填写水蚀野外调查表，依据野外调查单元土地利用现状分类中的二级分类划分土地利用类型。野外调查时，将土地利用类型、水土保持措施和盖度(郁闭度)均相同的连续区域划为同一地块，记录地块基本属性信息，回到室内清绘调查成果草图，根据植被盖度、土地利用、水土保持措施等信息，分别计算水土保持措施因子B，E，T。根据各县1981—2010年30 a逐日降雨资料，计算降雨侵蚀力R；以中国土壤发生分类中的土属信息为依据，确定土壤可蚀性K；利用调查区地形信息提取并计算坡度因子S和坡长因子L[12]。坡度坡长因子计算基于1∶10 000地形图插值生成的DEM，坡度为采用D8算法计算所得最大坡降方向的坡度值，坡长用累积径流算法获取。将各土壤侵蚀因子连乘，即得到各地块和调查单元的多年平均土壤水蚀模数[11]。

表1 调查样带基本信息统计

1.3 数据处理

调查单元和地块基本信息的存储与整理以及各因子值、土壤侵蚀模数的计算在Excel 2010中完成。在SPSS 18.0中对数据进行描述性统计分析，得出各指标平均值、最大值、最小值、标准差等并计算变异系数，对各样带侵蚀模数与坡度及坡长的关系进行相关分析，对不同坡度、坡长、土地利用类型、横向样带上的侵蚀模数均值进行方差分析等。在ArcMap 10.2和Origin 8.5中绘制调查样带分布图、土壤侵蚀模数随坡度坡长变化图、不同土地利用类型土壤侵蚀模数区域差异图等。

2 结果与分析

2.1 坡度与土壤侵蚀

将各样带所有地块坡度平均值作为该样带坡度均值，统计特征表明(表2)，横向样带坡度均值自北向南呈先减小后增大的变化趋势；位于陕南秦巴山地区的H6坡度均值最大，高达48.04°；地处关中平原的H4坡度均值最小，仅为4.03°，其中低于10°的地块面积占该样带总面积的76.22%；纵向样带对比来看，Z2(东侧)坡度均值略高于Z1(西侧)。陕南秦巴山地区坡度变异系数最小，坡度分布集中，其次是陕北黄土高原地区，关中平原为地堑式构造平原，台塬面向河谷平原呈阶梯状倾斜，个别地块坡度远高于均值，提高了该地坡度的变异性。横向样带大部分呈负偏态分布，说明坡度大于均值的地块较多，关中平原样带H4则为正偏态分布，表明大部分地块坡度低于均值4.03°。H4的峰度系数远高于其他样带，坡度分布曲线呈尖峰状。

表2 各样带不同坡度下侵蚀模数均值对比 t/(km2·a)

坡度通过影响土壤质地、累积入渗量、坡面径流水深及流速、坡面流切应力等因素决定土壤侵蚀的强弱，制约土地利用方向及水土保持措施布设。除H4和H5外，其余样带内地块土壤侵蚀模数先随坡度增大而增大，而后随坡度增大后减小，这与众多研究认为土壤侵蚀存在临界坡度的观点一致[13]。可能受主要侵蚀方式影响，不同样带达到最大土壤侵蚀模数的坡度范围有所不同[13]，H4，H5和Z2在坡度15°～20°内土壤侵蚀模数最大，H1，H7土壤侵蚀模数在20°～25°达到最大值，H6与Z1则在25°～30°范围内土壤侵蚀模数最大。H4土壤侵蚀模数随坡度增加整体上呈波动增加趋势，该样带土壤侵蚀模数与坡度呈极显著正相关关系(p<0.01)，相关系数为0.59，高于其他样带；坡度小于30°时，H5土壤侵蚀模数随坡度增大先增后减，在20°～25°范围内侵蚀模数达最大值，坡度大于30°后，土壤侵蚀模数有所增大。整体而言，陕西省土壤侵蚀模数随坡度增加呈先增加后减小的变化特征，在15°～30°土壤侵蚀模数最大(图2)，表明陕西省土壤侵蚀的临界坡度在此范围之内。

图2土壤侵蚀模数随坡度变化特征

从空间分布特征来看，各坡度范围内横向样带土壤侵蚀模数均在黄土高原或秦巴山区南段最高，在关中平原或秦巴山区北段最小。自北向南0°～10°坡度范围内土壤侵蚀模数先增后减，在H2或H3处达最大值；10°～20°坡度范围内土壤侵蚀模数呈增—减—增的变化趋势，H3处达最大值，H5处达最小值；20°～35°坡度范围内整体上土壤侵蚀模数呈先减小后增加的变化特征；>35°后，土壤侵蚀模数总体上由北向南递减。

东西对比来看，0°～20°和>35°坡度范围内，东侧样带Z2土壤侵蚀模数大于西侧Z1；20°～35°范围内，Z1土壤侵蚀模数大于Z2。根据容许土壤流失量和水力侵蚀强度分级[14]可知，黄土高原地区坡度>10°的地块基本上均存在轻度或中度侵蚀，其中样带H1在30°～35°为强烈侵蚀，坡度<20°和>30°时，土壤侵蚀模数均在黄土高原地区达到最大值，充分表明黄土高原地区是我国乃至世界土壤侵蚀极严重的地区之一。关中平原是陕西省的农业基地、经济中心，样带H4在地势起伏地区仍存在轻度侵蚀，需加以关注。陕南地区H7样带山体陡峭，在20°～30°侵蚀强烈，高于其他样带，这与当地人为破坏植被，陡坡开荒有关。

2.2 坡长与土壤侵蚀

坡长是影响坡面径流与侵蚀产沙过程的重要地形因子之一，USLE中规定坡长为地表径流源点到达有沉积出现的地方或明显渠道的水平距离[15]。将各样带所有地块坡长平均值作为该样带坡长均值，统计特征表明(表3)，关中平原地区样带H4坡长均值最小，为18.08 m；陕北黄土高原和陕南秦巴山地的坡长均值相近，约为30 m；东西两条样带坡长均值相似，约28 m。可能是由于人类活动对地表特征的改变较大，关中平原地区坡长变异系数高于黄土高原地区和秦巴山地区。与坡度相比，坡长的变异系数较小，说明陕西省坡长对区域特征、土地利用方式等因素的变化不如坡度敏感；另外坡长与土壤侵蚀模数的相关系数略小于坡度，表明坡长因子对土壤侵蚀模数的影响程度低于坡度。

表3 各样带不同坡长下侵蚀模数均值对比 t/(km2·a)

坡面径流水动力条件、降雨侵蚀力等因素随坡长的变化影响土壤侵蚀过程和强度[15]。各样带土壤侵蚀模数随坡长的变化规律不尽相同，H1，H2，H6土壤侵蚀模数整体上随坡长增大而增大；H3土壤侵蚀模数随坡长增加先增加后减小；H4和H5样带土壤侵蚀模数随坡长呈波动增长，均属微度侵蚀。纵向两条样带，坡长<50 m时，土壤侵蚀模数与坡长呈正相关，坡长>50 m后，土壤侵蚀模数反而有所降低，主要是由于随坡长增加，水体能量多用于搬运坡上侵蚀泥沙，因而侵蚀有所减弱[15]。相关研究表明坡长对土壤侵蚀强度的影响存在3种方式，分别为增长型(土壤侵蚀强度随坡长增加而增大)、增减型(土壤侵蚀强度随坡长增大先增强后减弱)、波动型(坡长增加过程中，土壤侵蚀强度出现两次增加、两次减小)，增长型关系常用于坡面土壤侵蚀模型中[15]。方差分析表明，各坡长范围内土壤侵蚀模数存在显著差异(p<0.01)，40～50 m坡长范围内土壤侵蚀模数最高。总体来看，本研究中坡长与土壤侵蚀模数的关系可归为波动型(图3)，0～30 m和40～50 m土壤侵蚀模数随坡长增加，30～40 m和>50 m土壤侵蚀模数略有下降，在该坡长范围内土壤侵蚀模数降低的原因仍需进一步探究。

图3土壤侵蚀模数随坡长变化特征

各坡长范围段内的土壤侵蚀模数在区域间变化趋势明显，基本上均在黄土高原地区或秦巴山地区达到最大值，而在关中盆地区达最小值。黄土高原地区坡长超过10 m后，土壤侵蚀强度在轻度到强烈间变化，且坡长10～20 m、>30 m后土壤侵蚀模数为所有样带中的最大值，再次表明黄土高原土壤侵蚀问题极其严重。秦巴山地区H6坡长超过30 m后，土壤侵蚀模数高于允许土壤流失量；秦巴山地南部的H7样带所在地区多中深度切割且具岩溶地貌特征的中低山，山高谷深，地势起伏大，基岩疏松破碎，侵蚀较强烈[16]，0～10 m和20～30 m坡长范围内土壤侵蚀模数均高于其他样带。

2.3 土地利用与土壤侵蚀

土地利用方式是人类影响土壤侵蚀的最直接的途径，不合理的人类活动会破坏植被，恶化土壤结构，增加地表径流，从而加剧土壤侵蚀的发生和发展[17]。统计分析调查单元各土地利用类型地块可知(表4)，样带共包括地块1 711个，总面积为23 290 hm2，其中林地面积占比最大，为65.80%，其次为耕地，所占比例为17.34%，草地面积也较大，占调查总面积的9.15%。各土地利用类型侵蚀模数均值从大到小依次为其他土地>耕地>园地>草地>居民点及工矿用地>林地>交通运输用地>水域及水利设施用地。耕地土壤侵蚀模数整体较大，为2 168.88 t/(km2·a)，分别为林地和草地侵蚀模数均值的4.18倍和3.10倍，这主要是由于与林草地相比，耕地受枯落物、地表结皮、根系的保护作用较小，同时人类耕作活动对耕地土壤扰动强烈，土壤结构松散，易被侵蚀。园地中果树冠层在一定程度上削弱了降雨动能，因而侵蚀强度低于耕地；但园地地表常处于裸露状态，缺少植被保护，故土壤侵蚀模数高于林地和草地。

表4 不同土地利用类型土壤侵蚀强度面积占比 %

土地利用类型侵蚀强度分级对比结果表明(表4)，林地和草地几乎全部地区土壤侵蚀强度属轻度侵蚀以下，由此可见，退耕还林还草是保护水土资源，抑制土壤侵蚀，改善区域生态环境的有效生物措施；15.61%的耕地和12.05%的园地受到中度以上土壤侵蚀，5.56%的耕地和4.50%的园地土壤侵蚀强度达强烈以上，加之耕地和园地面积较大，因此耕地、园地两类土地利用类型对陕西省土壤侵蚀的贡献最大。超过98%的居民点及工矿用地和全部的交通运输用地、水域及水利设施用地属微度侵蚀。其他土地利用类型地块偏少，面积较小，土壤侵蚀模数存在较大变异性和不确定性。

耕地、园地、林地、草地4种土地利用类型占地面积最大，且与粮食安全，区域生态环境关系密切。对比区域间土壤侵蚀程度可知(图4)，该4类土地利用方式的土壤侵蚀模数在样带间均存在显著差异(p<0.01)。自北向南耕地土壤侵蚀模数呈增—减—增的变化特征，在H3处达到最大值[8 493.58 t/(km2·a)]，属极强烈侵蚀，其次为秦巴山地区的H7，关中平原H4处达到最小值[356.82 t/(km2·a)]。由北向南园地土壤侵蚀模数也呈增—减—增变化，H2处园地侵蚀模数最大，关中平原地区最低，黄土高原地区高于秦巴山地区；为了便于管理，园地的梯田化率高于耕地，因而整体土壤侵蚀程度低于耕地。林地和草地整体侵蚀程度较低，林草地土壤侵蚀模数均在黄土高原地区的H1处达到最大值，其次为H2，两者对比来看，H1，H2，H6处林地土壤侵蚀模数高于草地，而H4处草地侵蚀模数高于林地，表明黄土高原区与秦巴山地区草地抑制土壤侵蚀的作用更大，关中平原区林地的水土保持效益更佳。纵向两条样带对比来看，西侧Z1耕地侵蚀模数大于东侧Z2，园地、林地、草地侵蚀模数则是Z2大于Z1。

图44种典型土地利用类型侵蚀模数区域差异

2.4 土壤侵蚀区域差异

陕西省黄土高原区、关中平原和秦巴山地区三大自然地理单元的自然环境、人类活动及水土流失特征均存在巨大差异，认识水土流失区域差异可以为因地制宜地布设水土保持措施提供重要依据。整体来看，纵向样带土壤侵蚀模数均值无显著差异，均在1 000 t/(km2·a)左右，这主要是由于两条样带距离较近，自然环境特征、人类活动相似；Z2略高于Z1，主要是由于Z2平均坡度略大于Z1。方差分析结果表明，各样带土壤侵蚀均值存在显著差异(p<0.05)(图5)，随调查样带自北向南变化，土壤侵蚀模数先减小后增加，侵蚀强度从大到小依次为陕北黄土高原区、陕南秦巴山地区、关中平原区。本研究所得土壤侵蚀空间分异特征与李雄飞[7]、程琳等[8]的研究结果一致。

陕北黄土高原土壤侵蚀模数高于其他地区，三条样带侵蚀模数均在2 000 t/(km2·a)左右，接近中度侵蚀。H1到H3土壤侵蚀模数有所降低，这主要是由于样带H1，H2所经地区梁峁起伏，地面支离破碎，H3样带塬面相对平坦[16]，另外由H1到H3植被覆盖度逐渐提高，从而在一定程度上抑制了土壤侵蚀。退耕还林还草工程的实施使黄土高原区土壤侵蚀问题得到有效改善，侵蚀强度有所下降，但土壤侵蚀面积仍占区域面积的60%以上，仍是我国治理水土流失的重中之重[5]。该地应进一步研究退耕后近地表特性变化对土壤侵蚀的潜在影响，从而为退耕还林还草工程的可持续发展和区域水土保持规划提供依据。关中平原地区主要由黄土台塬和河谷平原组成，地形较平缓，适宜农耕，H4是横向样带中侵蚀模数最低的，为187.37 t/(km2·a)。作为陕西省经济、政治、文化中心所在地，关中平原应合理规划利用土地资源，注重防治潜在水土流失。秦巴山地区是陕西省降水量最多的地区，年降雨量达1 000 mm左右，且降雨集中于夏秋两季；另外该地区地势起伏大，山谷相间分布，土层浅薄，多陡坡开荒，人为毁坏植被情况较为严重，因而该区土壤侵蚀强度仅略低于黄土高原地区。该地三条样带侵蚀模数由北向南逐渐增大，样带H7侵蚀模数为横向样带中的最大值。秦巴山地区迫切需要加强坡地治理，对于25°以上的陡坡耕地坚决退耕，将小流域治理与开发承包到户，工矿交通建设中注重保护环境，结合多种途径促进山区农民脱贫致富。

图5横向样带侵蚀模数均值对比

3 结 论

(1) 横向样带坡度均值自北向南呈先减小后增大，在黄土高原地区地块坡度<20°和>30°的土壤侵蚀模数最大，陕南地区H7在坡度20°～30°侵蚀高于其他样带。整体来看，陕西省土壤侵蚀模数随着坡度增大呈先增后减的趋势，土壤侵蚀的临界坡度为15°～30°。

(2) 各样带土壤侵蚀模数随坡长的变化特征复杂，综合来看，土壤侵蚀模数随着坡长的增加呈波动增长趋势。10～20 m和>30 m坡长范围内，黄土高原土壤侵蚀模数最大；<10 m和20～30 m坡长范围内秦巴山地区土壤侵蚀模数最大。

(3) 耕地、园地两种地类对陕西省土壤侵蚀的贡献相对较大，林地和草地整体侵蚀程度较低。耕地与园地土壤侵蚀模数自北向南呈增—减—增变化，最大值分别位于样带H2和H3处。

(4) 纵向样带土壤侵蚀模数均在1 000 t/(km2·a)左右，自北向南横向样带土壤侵蚀模数先减小后增大，侵蚀强度从大到小依次为陕北黄土高原区、陕南秦巴山地区、关中平原区，各地应因地制宜布设水土保持措施，防治水土流失。

参考文献：

[1] 齐清,王天明,寇晓军,等.基于GIS的黄土高原小流域土壤侵蚀定量评价[J].水土保持研究,2009,16(3):1-5.

[2] 李智广,曹炜,刘秉正,等.我国水土流失状况与发展趋势研究[J].中国水土保持科学,2008,6(1):57-62.

[3] 张勇.陕西省水土流失现状、存在问题与防治对策[J].水利发展研究,2004,4(11):47-49.

[4] Li P, Mu X, Holden J, et al. Comparison of soil erosion models used to study the Chinese Loess Plateau[J]. Earth-Science Reviews, 2017,170:17-30.

[5] 孟庆香,刘国彬,杨勤科,等.黄土高原土壤侵蚀时空动态分析[J].水土保持研究,2008,15(3):20-22.

[6] 刘康,康艳,曹明明,等.基于GIS的陕西省水土流失敏感性评价[J].水土保持学报,2004,18(5):168-170.

[7] 李雄飞.陕西省土壤侵蚀模数计算与空间分布特征[J].陕西水利,2014(4):122-123.

[8] 程琳,杨勤科,谢红霞,等.基于GIS和CSLE的陕西省土壤侵蚀定量评价方法研究[J].水土保持学报,2009,23(5):61-66.

[9] 刘闻,曹明明,宋进喜,等.陕西年降水量变化特征及周期分析[J].干旱区地理,2013,36(5):865-874.

[10] 李团胜.陕西省土地利用动态变化分析[J].地理研究,2004,23(2):157-164.

[11] Liu B, Zhang K, Xie Y. An Empirical Soil Loss Equation[C]∥ Yuren J. Proceedings of the 12 th international soil conservation organization conference. Beijing: Tsinghua University Press, 2002.

[12] 国务院第一次全国水利普查领导小组办公室.水土保持情况普查[M].北京：中国水利水电出版社,2010.

[13] 陈晓安,蔡强国,张利超,等.黄土丘陵沟壑区坡面土壤侵蚀的临界坡度[J].山地学报,2010,28(4):415-421.

[14] 陈学兄.基于遥感与GIS的中国水土流失定量评价[D].陕西杨凌:西北农林科技大学,2013.

[15] 孔亚平,张科利,曹龙熹.土壤侵蚀研究中的坡长因子评价问题[J].水土保持研究,2008,15(4):43-47.

[16] 聂树人.陕西自然地理[M].西安：陕西人民出版社,1981.

[17] 赵文武,傅伯杰,吕一河,等.多尺度土地利用与土壤侵蚀[J].地理科学进展,2006,25(1):24-33.