林春蕾, 王克勤, 陈若君, 张继辉, 张伯浩, 卢锦钊

(1.西南林业大学 生态与水土保持学院, 云南 昆明650000; 2.重庆大地建设监理有限责任公司, 云南 昆明 650000)

玉溪抚仙湖是中国最大的深水型贫营养淡水湖泊，相应湖容水量2.06×1010m3，占云南省9大高原湖泊总蓄水量的72.8%[1]。其淡水资源是玉溪人民的灌溉、工业用水、生活用水以及工业用水等的主要源泉，并且抚仙湖的优美环境给玉溪人民带来了宝贵的旅游资源，其丰富的渔业资源和旅游资源是当地人民的收入来源之一。但是随着沿湖经济的快速发展，抚仙湖的营养状态明显上升，水质状况呈现下降趋势，湖泊的污染主要来自流域人类活动产生的含氮、磷等的污染物, 尤其是尖山河流域，在暴雨时产生的污染物超标情况比较严重[1-3]。

引起水土流失的因素有人为因素和自然因素，自然因素是水土流失发生的先决条件。自然因素又包括下垫面条件和气候因子，其中降雨是土壤侵蚀的重要影响因子[4]。降雨的各特征因子也不同程度的影响土壤侵蚀，如降雨量、降雨历时、平均降雨强度、最大30 min雨强等。而对于不同土地利用方式下降雨因素与土壤侵蚀的关系，前人的研究主要局限于北方丘陵沟壑区、黄土高原区[5-6]，红壤土区的抚仙湖流域的研究较少。本文拟研究不同雨量级别和不同雨强级别下的侵蚀性降雨特征与产沙关系，以期揭示土壤侵蚀规律，为西南岩溶区的水土流失综合治理提供理论依据。

1 研究区概况与研究方法

1.1 研究区概况

研究区位于云南省玉溪市澄江县尖山河小流域内，位于北纬24°32′00″—24°37′38″，东经102°47′21″—102°52′02″。北接龙街镇广龙村委会，南接禄充管委会，东临抚仙湖，西接晋宁县。最高海拔在流域北部，为2 347.4 m，最低海拔在尖山河入抚仙湖的入口处，为1 722 m，相对高差625.4 m。监测点的监测设施主要布设在流域的下游。流域地处低纬度高原，土壤主要是红紫泥土和红壤，红紫泥土主要分布在尖山大河上游河道顺流左岸方向。多年平均降雨量1 050 mm，雨季为5月下旬至10月下旬，降雨量占全年总降雨量的75%，旱季为11月上旬至次年5月中旬，降雨量占全年降雨量的25%。暴雨基本出现在雨季。

1.2 研究方法

1.2.1 试验区布置 根据水土保持措施类型和项目区自然环境条件，共设坡耕地(2个)、云南松天然次生林、灌草丛、人工幼林和经济林共6个标准径流小区(表1)。

表1 试验小区布置

1.2.2 数据获取和研究方法

(1) 降雨量资料。采用自记雨量计和人工观测结合观测降雨量，在径流小区附近安置自记雨量计(型号：JDZ-1，重庆水问仪器厂)一个，观测降雨量。同时同地安装标准雨量筒一个，用于校对自记雨量计。将雨量计中获取的资料，以6 h为间隔划分次降雨。利用RainRecordSetup 1.06.exe软件分析出次降雨量、降雨历时、降雨强度、I30，I60等降雨特征因子。

(2) 产沙量的获取。泥沙含量采用Model 711悬移质泥沙测定仪测量，在降雨后小区产流时及时进行测定。

2 结果与分析

2.1 次降雨量的年变化特征

由图1和表2可知，尖山河小流域3 a共降雨371场，2014年降雨总量为834.9 mm，侵蚀性降雨量524.8 mm,占全年降雨量的62.86%。2015年的降雨总量是838.4 mm，侵蚀性降雨量为306.8 mm，占全年降雨量的36.59%。2016年总降雨量1035.8 mm,侵蚀性降雨量742.2 mm，占全年降雨量的71.65%。

根据月变化特征图1可以看出，从年内的降雨变化看，3 a中的2月均无降雨事件发生，2014年降雨主要集中于5—9月，3，4，11和12月基本没有降雨，1和2月无降雨事件，大体变化趋势类似正态分布，6月降雨量达最大。2015年的降雨从1月开始，3月无降雨事件，4月的降雨量是全年最小的，全年降雨都很频繁，8月降雨量达到全年的最大值。2016年的降雨主要集中于5—11月，9月的降雨量达到全年最大。从年际变化来看，每年的4月降雨量差异最小，9月的降雨量年际变化幅度较大；就季节变化来看，3 a的降雨量集中于夏秋两季。

图1 不同年限的降雨量年变化

2.2 侵蚀性降雨量等级划分

2.2.1 不同降雨量级别的降雨事件分析 如表2所示，将尖山河小流域3 a的降雨量数据划分为7个降雨量等级，分别为H<5 mm,5 mm≤H<10 mm,10 mm≤H<15 mm,15 mm≤H<25 mm,25 mm≤H<50 mm,50 mm≤H<100 mm和100 mm≤H<150 mm。3 a的降雨中，2#径流小区共侵蚀53次，侵蚀性降雨共1 424.2 mm,3#径流小区共侵蚀58次，侵蚀性降雨共1 529.2 mm，4#径流小区共侵蚀32次，侵蚀性降雨为1 047.2 mm，5#径流小区侵蚀43次,侵蚀性降雨1 356.8 mm，5#径流小区侵蚀28次，侵蚀性降雨620.4 mm。

侵蚀性降雨次数随着降雨量级别的增加，出现跳跃式变化。当降雨量小于50 mm时，降雨量级别与侵蚀性降雨次数和侵蚀性降雨量呈正相关关系，随着降雨量级别的增加，降雨场次减少，但侵蚀性降雨次数增多，造成一次严重土壤侵蚀的降雨为一次短历时的大暴雨。当降雨量大于50 mm时，侵蚀性降雨次数随着降雨量级别增加而减少，这一趋势与降雨特征有关，但是大于50 mm的降雨引起土壤侵蚀的几率几乎是百分之百。对于2#,3#,4#和5#径流小区而言，对土壤侵蚀性贡献最大的降雨量级别为15～25和25～50 mm，4个径流小区的侵蚀性降雨场次分别是26,28,16和23，分别占总的侵蚀性降雨场次的49.06%,44.83%,46.88%和53.49%。而对于6#径流小区而言，当降雨量大于10 mm便开始引起严重的土壤侵蚀。

表2 不同降雨量等级的侵蚀性降雨的分析

2.2.2 不同降雨量级别的产沙量变化 影响坡面产沙的主要因素是下垫面和降雨，而降雨则是产沙的主要来源。如图2所示，从总体上来看，玉米地和烤烟地径流小区的平均产沙模数均显著高于次生林、人工林地和灌草地径流小区。当降雨量小于50 mm时，次生要和灌草地径流小区的产沙量变化幅度不大，而烤烟地径流小区的平均产沙模数随降雨量级别的增加而增大，玉米地径流小区的平均产沙模数大致符合正态分布，随着降雨量级别的变化，人工林地径流小区的平均产沙模数呈现不稳定变化规律，但随降雨量的增大，平均产沙模数呈现增大趋势。综上所述，对于坡耕地而言，大于10 mm的降雨会引起严重土壤侵蚀，而大于50 mm的降雨对灌草地和人工林的土壤侵蚀有显著影响，次生林有较好的水土保持效果。

图2 研究区不同降雨量级别的产沙量变化

从图2可知，在同一降雨，同一坡度条件下，影响土壤侵蚀性的因素主要是地表覆盖类型，不同的地表覆盖对拦截雨水和入渗存在很大的差异性，这样必然导致土壤侵蚀的差异性，顺坡种植烤烟最容易引起土壤侵蚀，其次是玉米坡耕地，顺坡种植不宜拦截雨水，但对泥沙的冲刷作用不是太大。

2.3 侵蚀性降雨强度等级划分

将尖山河小流域3 a的降雨量数据划分为7个降雨强度等级，分别为0～5,5,10～20,20～40,40～60,60～80和80～1 200 mm/h。

由表3可知，随着降雨强度等级的增加，5个径流小区的侵蚀性降雨场次均呈逐渐减少趋势，降雨强度越大，发生土壤侵蚀的频率就越低，但是引起的土壤侵蚀更为严重。雨强级别为0～5 mm/h的侵蚀性降雨次数最多，5个径流小区的侵蚀性次数分别是42,47,24,34和22，分别占总的侵蚀性降雨场次的79.2%,81.03%,75%,79.1%和78.6%。易发生土壤侵蚀的雨强级别是降雨强度大于20 mm/h，20～40的降雨总发生4场，其中1场导致土壤侵蚀，历时仅1 min，侵蚀性降雨量为18.6 mm。不同降雨强度级别的产沙量变化详见图3。

表3 不同降雨强度等级的侵蚀性降雨的分析

图3 研究区不同降雨强度级别的产沙量变化

从图3可以明显的看出，土地利用类型和平均雨强对产沙模数的影响显著，随着侵蚀性降雨强度等级的变化，3#和6#径流小区的平均产沙模数呈现先减少后增加的趋势，符合二次函数的性质，而研究区的降雨类型多为长历时的小型降雨，极少部分高强度的大暴雨。4#径流小区随着平均雨强的增加而增大，2#和5#径流小区的平均产沙模数在雨强的变化下无显著差异。

由此可见，极少的大暴雨引起的土壤侵蚀量高于频繁的小型降雨引起的土壤侵蚀，在制定水土流失防治措施要注意大暴雨对土壤侵蚀的影响。由于坡耕地土壤扰动比较大，因此在不同雨强级别下的产沙量均高于其他土地利用方式。

2.4 产沙与侵蚀性降雨特征指标的相关性分析

利用SPSS分析软件对不同降雨强度级别的产沙和降雨特征指标进行相关性分析，20～40 mm/h的降雨强度等级下，仅发生一次土壤侵蚀，因此剔除这一次数据，不列入进行分析(表4)。

表4 不同降雨强度级别的产沙和降雨特征指标的相关系数

注：**表示极显著相关； *表示显著相关；I30为最大30 min雨强；I60为最大60 min雨强。

从表4看出，4#径流小区与各降雨特征指标无相关关系；当平均降雨强度小于5 mm/h时，2#径流小区的产沙量仅与降雨量和I60呈现显著相关关系，其他降雨特征指标均与产沙无显著相关关系；3#和4#径流小区的产沙与各降雨特征指标都无显著相关关系；5#径流小区与降雨量呈极显著相关关系，6#径流小区与降雨量呈现显著相关关系，这可能是由于6#径流小区长期以顺坡沟垄方式种植烤烟，一旦降雨便会发生水土流失，而3#径流小区种植玉米。随着玉米的生长，植被覆盖度增加，对小型降雨具有一定的拦截作用。当降雨强度为5～10 mm/h时，影响3#径流小区产沙量的因素主要是降雨历时、降雨量、I30和I60，产沙量与降雨历时、降雨量和I60呈极显著正相关关系。

总体上看，在一定降雨强度范围内(0～5 mm/h)，2#径流小区的产沙量与降雨量和最大60 min雨强呈显著正相关关系，6#径流小区的产沙量与降雨量呈显著相关关系。对于3#玉米坡耕地而言，当雨强小于5 mm/h时，影响玉米坡耕地土壤侵蚀的因素是降雨历时，当降雨强度为5～10 mm/h时，产沙量与降雨历时、降雨量和最大60 min雨强呈极显著正相关关系。对于5#径流小区而言，当降雨强度为0～5 mm/h时，产沙量与降雨量呈极显著相关关系，当降雨强度为5～10 mm/h时,产沙量与最大60 min雨强呈显著相关关系，而4#径流小区与各降雨特征指标无相关关系。说明影响4#径流小区的因素可能是下垫面条件或其他因素，由高大乔木组成的人工林具有很强的林冠截留作用，灌草地的拦截泥沙的效应有限，随着降雨量的增大，灌草地的水土保持作用逐渐降低，因此在降雨极为频繁的南方地区而言，乔木林的水土保持效应较灌草地好。

3 讨论与结论

(1) 从3 a的降雨情况来看，大体变化类似正态分布，研究区的降雨主要集中于5—10月，汛期在6—9月；从10月开始，降雨量逐渐减少，而从季节分布来看，降雨集中在夏秋两季，干湿季节明显。这一规律与大陆性季风气候特点有关。

(2) 覆盖度不同的土地利用方式在拦截泥沙和径流等方面存在很大差异，耕作方式和农作物对土壤流失也有不同程度的影响。同一降雨，同一坡度条件下，影响土壤侵蚀性的因素主要是地表覆盖类型，不同的地表覆盖对拦截雨水和入渗存在很大的差异性，这样必然导致土壤侵蚀的差异性，顺坡种植烤烟最容易引起土壤侵蚀，其次是玉米坡耕地，顺坡种植不宜拦截雨水，但对泥沙的冲刷作用不是太大。王占礼等[7]，潘声旺等[8]究结果中，旱地不同土地利用方式下，玉米—油菜轮作的土壤侵蚀量最大。

(3) 极少的大暴雨引起的土壤侵蚀量高于频繁的小型降雨引起的土壤侵蚀，在制定水土流失防治措施要注意大暴雨对土壤侵蚀的影响。由于坡耕地土壤扰动比较大，因此在不同雨强级别下的产沙量均高于其他土地利用方式；而在黄土高原区，大暴雨出现频率较高，造成水土流失的雨型主要是大型暴雨，大多数的降雨一般不产生地表径流[11]。

(4) 许多学者认为降雨量、最大30 min降雨强度和降雨历时对山区坡面产流产沙均存在显著的相关性，并且一些学者认为种植花生的坡耕地与降雨历时无显著相关关系[9-11]，而笔者的结论中，不同的降雨强度级别下，不同土地利用方式的产沙量与降雨特征指标之间的关系也是存在差异性，灌草地的拦截泥沙效应有限，对小型降雨有一定的抗冲能力，而较大暴雨则需要群落层次较多的乔木林，较高的林冠层能够截留一部分雨水，而较低的草木层和枯枝落叶层可降低降雨动能，减小降雨的激溅侵蚀。并且有些植被覆盖度较高的降雨也有一定的抗侵蚀能力，但对于坡度较大的坡耕地而言，仍然建议进行退耕还林措施。

[1] 杨加林,李杰,李经纬等.抚仙湖—星云湖水生生物与水环境研究综述[J].云南地理环境研究,2012,24(2):1001-7852.

[2] 金相灿,屠清瑛.湖泊富营养化调查规范[M].北京:中国环境科学出版社,1990,291-294.

[3] 曹建华,蒋忠诚,杨德生,等.我国西南岩溶区土壤侵蚀强度分级标准研究[J].中国水土保持科学,2008,6(6):1-7.

[4] 马良,左长清,邱国玉.赣北红壤坡地侵蚀性降雨的特征分析[J].水土保持通报,2010,30(1):74-79.

[5] 纪启芳,张兴奇,张科利,等.贵州省喀斯特地区坡面产流产沙特征[J].水土保持研究,2012,19(4):1-5.

[6] 高超,朱继业,朱建国,等.不同土地利用方式下的地表径流磷输出及其季节性分布特征[J].环境科学学报,2005,25(11):1543-1549.

[7] 王占礼,邵明安,常庆瑞.黄土高原降雨因素对土壤侵蚀的影响[J].西北农业大学学报,1998,26(4):101-105.

[8] 潘声旺,雷志华,杨丽娟,等.几种典型边坡植被的产流、产沙特征[J].生态环境学报,2013,22(7):1167-1172.

[9] 林锦阔,李子君,许海超,等.降雨因子对沂蒙山区不同土地利用方式径流小区产流产沙的影响[J].水土保持通报,2016,36(5):7-12.

[10] 郑海金,杨洁,左长清,等.红壤坡地侵蚀性降雨及降雨动能分析[J].水土保持研究,2009,16(3):30-33.

[11] 焦菊英,王万忠,郝小品.黄土高原不同类型暴雨的降水侵蚀特征[J].干旱区资源与环境,1999,13(1):35-41.