王坚桦, 邱 凡, 谢福倩, 韦 石, 梁志鑫, 李桂芳

(1.广西大学 农学院, 南宁 530004; 2.广西壮族自治区水土保持监测站, 南宁 530023)

果树是我国农业产业结构的重要组成部分，种类众多，分布广泛，种植面积超过14万km2，占世界果树种植面积的18.1%，总产量达到2.8亿t，占全球总产量的17.0%，位居世界第一位[1]。我国的果园普遍种植在山区坡地，坡面水土流失较为严重，严重制约着当地经济发展，不利于区域生态环境保护[2]。植被是抑制果园坡面发生水土流失的因素之一[3]，果园坡面覆盖植被可以有效减轻水土流失[4]。研究表明，果园生草技术可以有效减缓果园坡面水土流失，改善土壤的理化性状[5-8]。但果园生草容易形成与果树的养分、水分，尤其是氮素营养的竞争，一定程度上会抑制果树的生长[9]，因此该措施在实际的果园生产中并未大面积推广。清耕法对杂草进行彻底清理[10]，可以减少表面杂草与果树争水争肥，利于果树从土壤中汲取更多的养分，但长期清耕会导致土壤板结、肥力下降，破坏果园生态环境，加剧果园水土流失[11]。在南方红壤区，清耕处理的坡面径流量是植被覆盖处理的4.4倍，侵蚀量是67.1倍[3]。刈割作为一种果园管理方式，要求每年不定期对果树下植被进行刈割除草处理，再覆盖于果树下，达到提高土壤肥力和防治水土流失的效果[12]。

广西地处南亚热带区域，气候条件优越，是我国种植热带水果的优势区域之一，具有良好的地域基础，是我国水果总产量超过千万吨的省区之一[13]。目前，广西广泛种植不同品种水果，因地制宜地发展多种水果种植业，比如金桔、香蕉、杧果等[14]。以杧果为例，杧果是广西除香蕉外最大宗的热带水果，享有“热带水果之王”的美称[15]。广西是我国最大的杧果种植区，种植面积达6万hm2，总产量达71.8万t，栽培面积和产量均位居我国第一位[16]。广西现有杧果园大多数建立在丘陵山地和坡地上，90%以上的果园依赖自然降雨以维持其生长与发育[17]，而果园的清耕除草等田间管理措施，会导致地表覆盖状态发生改变，遇到较强降雨时，将加剧坡面侵蚀[18]。但大规模坡地果园的水土流失问题目前尚未得到足够的重视。

本研究以广西杧果果园坡面为研究对象，基于野外原位径流小区观测试验，通过测定2018—2019年次降雨下清耕和刈割处理果园坡面产流产沙量，分析清耕和刈割处理下果园坡面土壤侵蚀特征，结合降雨等级和降雨特征参数等指标，探讨降雨特征对果园坡面侵蚀的影响。研究结果可为南方红壤丘陵区果园坡面水土流失的防治提供理论依据。

1 材料与方法

1.1 研究区概况

试验地在广西木棉麓水土保持科技示范园，位于广西壮族自治区南宁市良庆区那马镇(108°19′00″—108°19′24″E，22°37′26″—22°37′55″N)。园区属于亚热带湿润季风性气候，气候炎热，雨量充沛，年均降雨量1 232.2 mm，年均气温21.8℃，降雨集中且较多暴雨，80%以上的降雨集中在5—9月[19]。研究区土壤类型是赤红壤，成土母岩以沉积岩为主，多为砂质岩，土层厚度范围是100～150 cm，土壤pH 4.5～5.5，有机质含量2～3%，土壤容重为1.2～1.3 g/cm3。园区内现有自然植被主要是荒草和灌丛，主要荒草类型为鬼针草(Bidensbipinnata)和狼把草(Bidenstripartita)。

1.2 试验设计

1.2.1 试验小区布设 本试验样地采用广西木棉麓水土保持科技示范园的标准径流小区，小区布设在山坡中部，小区内坡面横向平整，坡度为16.5°，坡长20 m、宽5 m，土壤条件一致。试验采用6年生杧果为供试对象，一共设置2个处理，刈割和清耕处理，刈割处理是2018年1月和2019年5月、7月对果园杂草进行刈割，用镰刀清除坡地表面杂草，留茬3～5 cm；清耕处理是使果园坡面长期保持无杂草状态，坡面杂草高于1 cm时，用锄头铲除杂草，不保留杂草根系。

1.2.2 试验方法 采用野外径流小区水文观测法。降雨观测时间为2018年1月1日至2019年12月31日。径流小区下方设有集流桶和分流桶，用于承接并测定降雨产生的径流小区径流泥沙量。每次降雨结束后，利用标尺测定集流桶和分流桶水深，确定小区水土流失体积后，将样品充分混匀并收集起来，每个小区收集3个重复样品，每个重复样品采集1 000 ml，烘干后测定其中泥沙含量，计算径流量和土壤流失量参数。在研究区附近布设有气象站，用以收集气象资料，降雨资料的收集包括降雨量、降雨过程、降雨历时等降雨参数，根据降雨参数计算出I30(最大30 min降雨强度)。

2 结果与分析

2.1 降雨特征分析

研究区2018—2019年1—12月降雨量见图1。2018年总降雨量为1 211.5 mm，降雨主要集中在5—9月，共计1 017.6 mm，占全年总降雨量的84.0%，其中7月的降雨量最大，占全年总降雨量的26.1%。2019年总降雨量为928.6 mm，降雨主要集中在3—8月，达到793.8 mm，占全年总降雨量的85.5%，其中4月和8月的降雨量最大，分别占全年总降雨量的17.5%和24.0%。整体而言，2018年降雨量比2019年增加了30.5%，且2018年降雨相对较为集中(主要集中在5—9月)，2019年降雨在全年分布较为分散。

图1 研究区2018年、2019年降雨分布特征

本研究把能够引起径流小区坡面产流产沙的降雨均视为侵蚀性降雨。由表1可知，2018年有21场侵蚀性降雨(3月1日至9月16日)，侵蚀性降雨总量达到992.4 mm，占全年总降雨量的81.9%。21场侵蚀性降雨的降雨量2.0～125.8 mm，平均雨强0.9～76.8 mm/h，I301.6～89.8 mm/h，有10场降雨的I30超过30 mm/h，占总降雨场次的47.6%。按照中国气象降雨等级划分[20]，暴雨及大暴雨降雨等级的降雨量占侵蚀性降雨总量的72.8%(722.8 mm)，表明2018年研究区的侵蚀性降雨类型主要为暴雨及大暴雨。

2019年有22场侵蚀性降雨(2月17日至10月7日)，侵蚀性降雨总量为772.6 mm，占全年降雨量的83.6%。22场侵蚀性降雨的降雨量4.0～145.0 mm，平均雨强0.3～78.0 mm/h，I304.4～184.1 mm/h，有9场降雨的I30超过30.0 mm/h(40.9%)，其中4月16日大暴雨的I30达到184.1 mm/h。大雨、暴雨和大暴雨的降雨量分别占侵蚀性降雨总量的27.3%，14.8%和36.6%，可知2019年研究区的侵蚀性降雨类型主要是大雨和大暴雨(63.9%)。由上可知，2018年、2019年研究区80%以上的降雨是侵蚀性降雨，且40%以上的侵蚀性降雨I30超过30 mm/h。

表1 2018年、2019年研究区侵蚀性降雨特征

2.2 赤红壤果园坡面土壤侵蚀特征

由表2可知，2018年、2019年不同处理下赤红壤果园坡面年径流量和侵蚀量均表现为清耕>刈割。对于年径流量，2018年、2019年清耕小区坡面年径流量分别是刈割小区的2.4，1.7倍，且坡面径流系数也存在明显差异。对于年侵蚀量，2018年、2019年清耕小区坡面年侵蚀量分别是刈割小区的52.2，9.2倍。表明清耕处理会导致果园坡面产流产沙量显著增加，其中对坡面侵蚀量的增加效果更为明显，这与Pan等[21]、朱冰冰等[22]和张锐波等[23]的研究结果一致。

与2019年相比，刈割与清耕处理下，2018年坡面径流量分别增加0.5倍和1.2倍；2018年刈割处理下坡面侵蚀量减少0.5倍，清耕处理下坡面侵蚀量增加1.7倍。清耕小区在2018年的坡面径流量显著高于2019年，主要是因为2018年、2019年的降雨分布，尤其是侵蚀性降雨特征存在明显不同(表1)。

表2 不同处理下果园径流小区年径流量与侵蚀量

2.3 次降雨下赤红壤果园坡面产流产沙特征

为进一步探究清耕对赤红壤果园坡面产流产沙特征的影响，对2018—2019年次降雨条件下果园坡面产流产沙特征进行分析(图2)。次降雨条件下，2018年果园清耕与刈割处理下坡面径流量变化范围是1.1～384.7，2.4～71.7 m3/hm2，整体上均表现为清耕小区径流量大于刈割小区，且二者随次降雨的变化趋势基本一致。次降雨下，不同处理下3—4月果园坡面径流量较低(<7.46 m3/hm2)且变化较平稳，同时二者差异较小；而5—9月出现了不同程度的波动，且二者差异有所增大。对于2019年，次降雨条件下清耕与刈割处理下果园坡面径流量介于2.0～174.0 m3/hm2，2.0～46.3 m3/hm2，同时在年内存在明显的波动特征，主要表现在2—3月和7—9月波动较大，且存在明显差异；4—7月中旬变化趋于平缓，二者的差异较小。

图2 2018年、2019年次降雨条件下赤红壤果园坡面产流特征

次降雨下，2018—2019年不同处理下果园坡面侵蚀量的变化趋势与径流量相似(图3)，整体上，清耕小区的侵蚀量均大于刈割小区，且二者变化趋势基本一致。2018年清耕与刈割处理下果园坡面侵蚀量变化范围是0.1～6 375.6，0.2～34.0 m3/hm2，次降雨下坡面侵蚀量在3—4月波动小(0.1～8.5 m3/hm2)，且二者差异较小；5—9月清耕处理的侵蚀量波动较大(2.7～6 375.6 m3/hm2)，刈割处理波动幅度小，且二者的侵蚀量存在明显差异。对于2019年，清耕和刈割处理下果园坡面侵蚀量介于0.3～2 446.7 kg/hm2，0.3～103.3 kg/hm2，次降雨下，坡面侵蚀量除在9—10月较小且变化平稳外，其余时段下果园坡面侵蚀量呈现明显的波动性，且清耕处理下果园坡面侵蚀量波动明显大于刈割处理。同时结果还表明，坡面侵蚀量主要由一场或几场侵蚀性降雨导致的，例如，2019年3月31日降雨产生的侵蚀量(2 446.6 kg/hm2)占清耕处理下年侵蚀量的79.2%，这与李桂芳[18]，Pan[21]等的研究结果一致。

图3 2018年、2019年次降雨条件下赤红壤果园坡面输沙特征

通过对2018年、2019年果园坡面径流泥沙量进行分析可知，坡面的产流产沙特征，除受果园清耕及刈割管理措施影响外，也受降雨分布、降雨量、降雨等级、降雨强度等[24]因素影响。由图1可知，2018年的降雨量大且集中，2019年的降雨量小且分散，对于清耕小区，2018年坡面的径流泥沙量显著大于2019年，说明降雨的年内分布不同会影响清耕与刈割处理下果园坡面侵蚀特征，与王昭艳等[25]的研究结果一致。2018年5月前的侵蚀性降雨主要以小雨和中雨为主，降雨量较小，所以果园坡面3—5月的径流量较小且平稳；5—8月的降雨普遍是大雨及以上等级的降雨，降雨量较大，所以坡面径流量波动上升。

此外，由2019年果园坡面输沙特征可知，与清耕处理相比，8月1日的大暴雨下刈割处理下坡面侵蚀量较大(45.5 kg/hm2>20.7 kg/hm2)，可能的原因是(1) 本文的清耕处理是对坡面不定期的除草，而7月、8月是地表植被生长较为旺盛时期，刈割处理下7月上旬刚进行了刈割除草，此时清耕和刈割处理下二者坡面植被覆盖度差异较小；(2) 7月下旬研究区降雨较少，坡面土壤较为干燥；(3) 8月1日大暴雨虽降雨量大，但降雨历时长，平均降雨强度和瞬时降雨强度较小。坡面土壤侵蚀过程受降雨特征、下垫面条件及土壤状况等因素的综合作用，自然降雨下野外径流小区坡面产流产沙是这些因子耦合作用的结果。

2.4 不同处理下果园坡面土壤侵蚀对降雨等级的响应

由2018年不同降雨等级下刈割与清耕处理果园坡面产流量可知(图4A)，清耕处理下暴雨和大暴雨产生的坡面径流量分别占坡面年径流量的63.6%和24.2%；刈割处理下则分别占42.8%和36.6%。对于2019年(图4A)，清耕处理下坡面径流量主要由暴雨和中雨导致的，分别占年径流量的41.1%和26.6%；刈割处理下则主要由中雨、大雨和暴雨产生的，分别占年径流量的34.5%，22.1%和21.9%。对比不同处理下果园坡面产流对降雨等级的响应可知，中雨和大雨降雨等级下，2018年清耕比刈割处理下坡面径流量高8.1%和82.2%，2019年则分别高28.5%和17.0%。暴雨和大暴雨降雨等级下，2018年清耕比刈割处理下坡面径流量高2.6，0.6倍，2019年高2.1，0.8倍。由2018年、2019年不同降雨等级下坡面径流量占比可知，由于降雨的年内分布不同，2018年清耕和刈割处理下坡面径流主要来自暴雨和大暴雨降雨等级，2019年则主要源于中雨和暴雨降雨等级。同时，与清耕处理相比，由于地表覆盖度不同，刈割处理下可以减少暴雨和大暴雨降雨等级下果园坡面产流量，且在暴雨下的降低作用最显著[26]。

对不同降雨等级下坡面侵蚀量而言(图4B)，2018年清耕处理下坡面侵蚀量主要来源于暴雨(占比89.5%)，而刈割处理下主要是由大雨、暴雨和大暴雨(占比16.2%，49.3%和22.5%)产生；2019年清耕处理下暴雨降雨等级下侵蚀量占年总侵蚀量的82.7%，刈割处理下主要由中雨和暴雨导致，分别占比43.8%和29.1%。小雨、中雨、大雨、暴雨和大暴雨下，2018年清耕比刈割处理下坡面侵蚀量高0.2～93.7倍不等，2019年高0.5～26.8倍不等。与清耕处理相比，刈割处理在中雨、大雨和暴雨降雨等级下均能显著降低坡面侵蚀量，且对暴雨降雨等级下的降低作用最为明显，这与Lenka等[27]的研究结果一致。同时，对比2018年、2019年不同降雨等级下果园坡面径流量和侵蚀量可以看出，清耕和刈割处理下果园坡面径流量和侵蚀量对降雨等级的响应存在明显不同。表明自然降雨条件下，除降雨量和降雨强度外，降雨时间、降雨年内分布等均会影响坡面土壤侵蚀特征。因此，在后续的研究中需要综合考虑观测时间序列及区域降雨特征等因素对坡面侵蚀的综合影响。

图4 不同降雨等级下坡面径流量和侵蚀泥沙量占比

2.5 不同处理下影响赤红壤果园坡面产流产沙的降雨特征参数

为了进一步分析不同处理下果园坡面产流产沙特征的影响因素，对果园坡面产流产沙量与降雨特征参数进行相关分析，见表3。可知，刈割和清耕处理下影响坡面产流产沙的降雨特征不同。对刈割处理下，坡面径流量与降雨历时和I30均表现为显著相关关系，与降雨量呈极显著相关关系；清耕处理下坡面径流量与降雨历时呈显著相关关系，与降雨量表现出极显著相关关系。而坡面侵蚀量仅在清耕处理下与降雨历时呈显著相关关系，与其他降雨参数均无显著相关关系。这与Pan等[21]的研究结果一致。

表3 不同处理下坡面产流产沙与降雨特征参数的相关分析

3 结 论

(1) 2018年研究区的总降雨量是1 211.5 mm，主要集中在5—9月(占年降雨量的84.0%)；2019年是928.6 mm，主要集中在3—8月(占年降雨量的85.5%)。2018—2019年研究区80%以上降雨是侵蚀性降雨。

(2) 2018年、2019年果园坡面年径流量、年侵蚀量及次降雨下坡面产流产沙量均表现为清耕小区>刈割小区，且侵蚀量差异更为明显。次降雨下，2018年不同处理下坡面径流泥沙量在3—4月变化平稳，5—9月出现波动，且清耕处理波动幅度大；2019年不同处理下坡面径流量在4—7月变化平缓，侵蚀量在年内整体波动较大。

(3) 刈割和清耕处理下径流量和侵蚀量的差异主要是由几场侵蚀性降雨产生的。暴雨是果园坡面产流产沙的主要因素之一，与清耕处理相比，刈割处理可以减少暴雨和大暴雨降雨等级下果园坡面径流量，中雨、大雨和暴雨降雨等级下果园坡面侵蚀量，且对暴雨的降低作用尤为明显。刈割和清耕处理下坡面径流量与降雨历时呈显著相关关系，与降雨量呈极显著相关关系。