黄土高原区夏季集中降雨的气候特征，使得水力侵蚀一直是引起黄土高原地区土壤资源流失的主要原因。薄层水流搬运侵蚀是细沟间侵蚀的一大分支，前人已对薄层水流搬运侵蚀的侵蚀发生过程和发生机理进行了细致的研究^[1]。降雨溅蚀作为细沟间侵蚀的另一分支是降雨前期的主要侵蚀类型，会导致土壤颗粒发生跃迁，将土面压实，进而造成土壤入渗能力降低，土壤产生结皮，降低土壤生产力等现象^[2-5]。溅蚀是由于雨滴打击而造成的土壤颗粒发生分散并产生跃迁的一种侵蚀方式^[6]。原理为雨滴由于在下降过程中受重力作用而具有一定的动能，当雨滴掉落地面与地表发生碰撞时，一部分能量被土粒所吸收转化为了内能，另外一部分转化为土体结构发生破坏的形变势能，导致被破坏的部分土粒具有了动能并向四周发生跃迁。坡度作为影响溅蚀的重要因素，前人对坡度对溅蚀影响方面的进行了相关研究。Ellison^[7]研究表明在10°的坡上，向下溅蚀量为向上溅蚀量的3倍；王双涛^[8]、王贵平^[9]、李叶鑫^[10]和田娜玲等^[11]研究表明溅蚀总分散量与向下坡溅蚀量随坡度的增大而增大；焦银龙^[12]、张慧荟^[13]以及赵晓光和吴发启^[14]的研究结果为向下坡溅蚀量随坡度的增加先增加后减小；张科利和细山田健三^[15]、Morgan^[16]研究发现由于径流流速的影响导致向上坡溅蚀量和溅蚀总分散量随着坡度的增加变化不明显^[17-18]。前人研究多为野外观测和径流小区模拟降雨实验，从大面积尺度上对坡度与溅蚀量的关系进行研究，由于所选用的土壤类型、坡度范围与降雨强度不同，所获结果存在差异性^[19-20]。本研究选用关中地区主要耕作土壤——塿土为对象，通过模拟降雨对单位面积尺度上不同方向和不同距离范围溅蚀量进行处理和分析从而深入探究溅蚀的潜在机理。为控制降雨范围，减小降雨对收集装置的影响，防止二次溅蚀造成的误差，选用针头式模拟降雨机。并给针头式模拟降雨机安装线性震动马达解决雨滴始终滴落在同一点的问题，确保了雨滴掉落坡面的随机性使模拟降雨更加接近于自然降雨，提高了数据的准确性和可靠性。研究结果将为黄土高原溅蚀发生机制提供科学依据。

研究区位于陕西省杨凌高新技术产业示范区，107°59'~108°08'E，34°14'~34°20'N，平均海拔为468 m，总面积为135 km²。该地区位于黄土高原南缘，关中平原腹地，属于温带半湿润大陆性季风气候，常年平均温度为12.9℃，年平均降水量为635.1~646.9 mm，是该地区水资源的主要来源，且60%~70%降水集中于7—10月份。

该区分布有7个土类，11个亚类，15个土属，24个土种。其中塿土占总面积的71.7%，为该区的主要土壤类型。塿土是我国陕西关中平原、晋南、豫西一带主要的农业生产土壤，该土壤是人类长期耕种褐土和施肥熟化所形成的特殊土壤类型，土壤肥沃，非常适宜耕作，土壤呈棕灰色。土壤具体理化性质见表 1。

表 1 Table 1

表 1 供试土壤理化性质 Table 1 Physical and chemical properties of the tested soil 土壤类型Soil type 颗粒组成 Particle composition/% 有机质 Organic matter/% 铵态氮 Ammonium nitrogen/（mg·kg–1） 硝态氮 Nitrate nitrogen/（mg·kg–1） 速效钾 Available potassium/（mg·kg–1） 砂粒 Sand 粉粒 Silt 黏粒 Clay 塿土Lou soil 21.5 50.8 30.3 7.4 37 4 64

表 1 供试土壤理化性质 Table 1 Physical and chemical properties of the tested soil

室内模拟降雨试验于西北农林科技大学资源环境学院水工实验室内进行。采用针头式模拟降雨机如图 1所示。降雨机由恒压水箱、雨滴发生装置及钢筋框架三部分构成。恒压水箱通过马氏瓶原理使水箱保持恒压，通过调节水箱中玻璃管内水深改变水箱的压强使降雨强度恒定为80 mm·h^–1。雨滴发生装置为下滴式，由长方形凹槽和底部的80个针头组成。通过色斑法对不同规格尺寸的注射针头所产生的雨滴粒径进行测量后，发现九号针头所产生的雨滴粒径在雨强稳定后集中在1至2 mm范围内，且雨滴数密度与天然降雨十分接近，可以用来进行模拟降雨试验。针头式模拟降雨机上的线性震动马达是为解决雨滴始终滴落在同一点的问题而设置，由研究者自行改进，这确保了雨滴掉落坡面的随机性，使模拟降雨更加接近于自然降雨。使用长方体装土槽（25 cm×20 cm×15 cm）进行不同坡度单位面积土坡的模拟。使用两块木工板（1 220 mm×2 440 mm）合并在一起作为溅蚀量收集装置。

图 1 针头式模拟机降雨机示意图 Fig. 1 Schematic diagram of the new model of needle spinkling device

该试验共设置三组重复。选取5°、10°、15°和20°四个坡度处理，0°平面作为对照。降雨强度设定为80 mm·h^–1，降雨历时为10 min。土壤溅蚀主要是降雨过程前期的侵蚀现象。重复试验发现，当降雨超过10 min后土壤表面会形成积水，此时土粒发生的跃迁不仅仅由雨滴击溅造成，还有其他途径。为了保证试验数据的准确性，选用降雨历时前10 min进行观测。装土时每装5 cm进行一次压实和打毛，将土壤容重控制在1.3 g·cm^–3。将装土槽放在木工板上进行降雨试验。试验时各个方向均存在大于30 cm距离范围的溅蚀土粒，未收集到超过50 cm距离范围的溅蚀土粒，因此本试验设计的溅蚀盘可收集全部溅蚀量。降雨结束后将溅蚀出来的土粒在木工板上风干1 h，使用事先称好的铝盒对木工板上的土粒进行收集，收集结束后将铝盒放入烘箱烘干并称重。一些溅蚀到木工板上较为细小的土粒收集十分困难，这是试验误差的主要来源。每次试验结束后均冲洗收集板，以消除第一次试验未收集的溅蚀量对第二次试验的影响。该研究针对单位面积上的溅蚀量，因此溅蚀量较少，为提高数据的准确性，称重时使用千分之一天平，记录数据时保留三位有效数字。

为了研究坡度对不同方向溅蚀量的影响。采用向下坡溅蚀量、向上坡溅蚀量、侧坡溅蚀量、溅蚀总分散量及溅蚀净搬运量5个指标进行综合对比分析。向下坡溅蚀量是由雨滴击溅造成向坡面下方发生跃迁的土粒总量，本研究中将坡面前方两个扇形区域收集的土粒总量视为向下坡溅蚀量。向上坡溅蚀量是由雨滴击溅造成向坡面上方发生跃迁的土粒总量，本研究中将坡面后方两个扇形区域收集的土粒总量视为向上坡溅蚀量。侧坡溅蚀量是由雨滴击溅作用造成向坡面两侧发生跃迁的土粒总量，本研究中将坡面两侧四个扇形区域收集的土粒总量视为侧坡溅蚀量。溅蚀净搬运量和溅蚀总分散量是具有实际意义的溅蚀量概念，可以描述溅蚀特性。溅蚀总分散量是区域内发生溅蚀的总量，也可以理解为向上坡溅蚀量、向下坡溅蚀量和侧坡溅蚀量的总和。溅蚀净搬运量是真实发生土壤转移的量，也可以理解为向下坡溅蚀量和向上坡溅蚀量之差。

采用Origin2017进行坡度和不同类型溅蚀量关系图表的绘制，使用SPSS19.0对不同坡度下溅蚀量数据进行显著性差异分析。

图 2为不同坡度下不同距离内溅蚀分散量占总量的百分比。由图可知：0~10 cm距离范围内溅蚀分散量占总分散量的百分比随着坡度的增加呈单调递减的趋势。10~30 cm距离范围内溅蚀分散量占总分散量的百分比随着坡度的增加多次出现峰值，不存在明显的相关性关系。大于30 cm距离范围的溅蚀分散量占总分散量的百分比随着坡度增加呈单调递增的趋势。可以得出0~20°范围内溅蚀距离随着坡度的增大而增大。随着坡度的增大，雨滴动能在正交分解后平行于坡面方向的分力增大，导致发生溅蚀的土粒平行于坡面方向运动的初始加速度增大。因此0~20°范围内土粒溅蚀距离随着坡度的增大而增大。

图 2 不同坡度下不同距离范围内溅蚀分散量占溅蚀总分散量的百分比 Fig. 2 Percentage of the splash ersosion dispersion within a certain range to the total dispersed splash erosion relative to distance and slope gradient

表 2为不同坡度下不同方向的溅蚀量，由表可知向下坡溅蚀量在0~20°内随着坡度的增加而增加，向上坡溅蚀量随着坡度的增加先减小后增加，侧坡溅蚀量上下波动无显著规律。

表 2 Table 2

表 2 不同坡度下不同方向溅蚀量 Table 2 Amount of splash erosion relative to direction and slopes gradients 方向 Direction 溅蚀量Amount of splash erosion/g 0° 5° 10° 15° 20° 向下坡溅蚀量① 0.479±0.006eC 0.870±0.019dB 1.250±0.155cB 1.752±0.039bB 3.553±0.225aC 侧坡溅蚀量② 2.676±0.018aA 2.394±0.026bcA 2.267±0.225cA 2.560±0.018abA 2.370±0.234bcB 向上坡溅蚀量③ 0.636±0.027aB 0.624±0.023abC 0.629±0.033aC 0.575±0.016bcC 0.556±0.037cA   注：① Downward splash erosion，② Sideward splash erosion，③ Upward splash erosion. 注：*表内大写字母为不同方向溅蚀量数据间的显著性差异分析，表内小写字母表示各方向溅蚀量与坡度间的显著性差异分析。Note：The capital letters in the table mean significant difference between splash erosion different in direction，and the small letters in the table significant difference between splash erosions relative to direction and slope gradient.

表 2 不同坡度下不同方向溅蚀量 Table 2 Amount of splash erosion relative to direction and slopes gradients

图 3为不同方向溅蚀量与坡度的关系。由图可知在0~20°范围内向下坡溅蚀量随着坡度的增加呈单调递增趋势，对向下坡溅蚀量与坡度进行线性方程拟合，得出y=0.1411x+0173，且相关性显著。在0~15°范围内向上坡溅蚀量随着坡度的增大而减小，在15~20°范围内向上坡溅蚀量随着坡度的增加而增加。

图 3 不同方向溅蚀量与坡度的关系 Fig. 3 Relationship between splash erosion and slope gradient relative to direction

上述结果可总结为：在0~20°范围内向下坡溅蚀量与坡度呈 y=0.1411x+0.173的线性关系，决定系数R²为0.8565。当坡度小于15°时，向上坡溅蚀量随着坡度的增加而减小，当坡度大于15°时，向上坡溅蚀量随着坡度的增加而增加。0~20°范围内侧坡溅蚀量与坡度之间无显著相关性。雨滴在击打地面时所施加的力可正交分解为垂直于坡面的压力与平行于坡面的剪切力。土壤发生颗粒迁移的重要原因是剪切力破坏了土粒原本土体结构使土粒具备了发生跃迁的动能。将平行于坡面的剪切力分为向上坡方向、向下坡方向和侧坡方向。当坡度增大时雨滴的剪切力在向下坡方向的分力越来越大，向上坡方向的分力越来越小，导致向下坡发生跃迁的土粒数量增加，向上坡方向发生跃迁的土粒数量减少。由于土粒自身重力平行于坡面的分力，对其向下坡运动有促进作用，对其向上坡运动有阻碍作用。随着坡度的增大，重力平行于坡面的分力也会随之增大，导致重力分力对土粒向下坡运动的促进作用和向上坡的阻碍作用更加显著。所以在0~20°范围内随着坡度的增大向下坡溅蚀量越来越大，向上坡溅蚀量越来越少。

由图 3可知，侧坡溅蚀量与坡度间不存在显著性关系。理论上随着坡度的增加，雨滴在侧坡方向的剪切力越来越小，侧坡溅蚀量应该随之减少。但由于坡度对侧坡溅蚀量影响不明显，雨强较大时雨滴击打对溅蚀的促进作用较坡度对溅蚀的减缓作用更为显著。因此侧坡溅蚀量与坡度关系不显著。

图 4为溅蚀总分散量、溅蚀净搬运量与坡度的关系。由图 4可知在0~20°的范围内溅蚀总分散量随着坡度增加而增加，并且增加幅度越来越大，关系趋近于二次函数关系。

图 4 溅蚀总分散量和净搬运量与坡度的关系 Fig. 4 Relationships of total dispersion and net transportation of splash erosion with slope gradient

上述结果可总结为：在0~20°范围内，溅蚀总分散量与坡度间呈y=0.0156x²-0.21x+4.5815的二元函数关系，决定系数R²为0.9983。该函数反映了溅蚀总分散量随着坡度的增大而增大，坡度越大，溅蚀总分散量变化越明显。当坡度增大时平行于坡面的剪切力增大，破坏土粒本身土体结构的效果就会越明显，垂直于坡面方向的压力越来越小，使雨滴对地面的压实效果变弱，土粒更容易发生跃迁。随着坡度的增大土壤含水量降低，土粒自身的黏结力越来越低，导致土粒本身更容易向上坡方向、向下坡方向和侧坡方向发生溅蚀。由于随着坡度的增加，土粒自身的重力在沿着平行于地面方向的分力越来越大，因此存在越来越多土粒发生跃迁。

由图 4可知0~20°的范围内，当坡度为0°时溅蚀净搬运量接近于0，当坡度不为0°时溅蚀净搬运量始终大于0，因此在坡面内存在向下坡溅蚀量恒大于向上坡溅蚀量的现象。宏观上体现为溅蚀在坡面内表现为向下堆积的过程。在0~20°范围内溅蚀净搬运量随坡度的增大而增大，增大幅度也越来越大，呈二次函数关系：y=0.0086x²-0.0273x-0.0355, 决定系数R²为0.9742。该函数反映了坡度越大溅蚀净搬运量的效果越明显，可以理解为随着坡度的增加，降雨过程中由于溅蚀造成的土壤转移量也会增大。随着坡度的增大，雨滴剪切力沿下坡方向分力越来越大，沿上坡方向分力越来越小。土粒自身重力平行于坡面方向分力越来越大。从而导致坡面上更多的土粒向下跃迁。溅蚀净搬运量在坡度为0°时接近于0，是由于当坡度为0°时，土粒不受任何平行于土壤表面剪切力分力的影响，导致土壤发生溅蚀跃迁出该面积的土粒总量与该面积周围区域发生溅蚀跃迁入该面积的土粒总量大体一致，宏观表现为该面积未发生土壤搬运。溅蚀净搬运量随着坡度的增大而增大，且恒大于0。

采用小面积的装土槽研究了单位面积上坡度对溅蚀量的影响。选用5个具有不同坡度（0°，5°，10°，15°，20°）的装土槽进行单位面积土坡的模拟，使用较为先进的新型针头模拟降雨装置有效控制降雨面积，防止雨滴对收集装置造成土粒的二次跃迁，可以较为准确的数据。研究发现，在0~20°坡度范围内，溅蚀总分散量、溅蚀净搬运量以及向下坡溅蚀量均随着坡度的增大而增大，向上坡溅蚀量随着坡度的增加先减少后增加。造成这一现象的主要原因在于雨滴剪切力向下坡方向分力越来越大，向上坡方向分力越来越小，土粒自身重力产生平行于坡面方向的分力越来越大。坡度的改变会影响土壤含水量与黏结力，导致上述现象的发生。本研究所得结论与焦银龙、张慧荟以及赵晓光等研究的部分结果一致。文献中报道的相关研究大多基于野外径流小区试验，本研究基于单位面积室内试验，但由于坡度是影响溅蚀量的主要因素——土粒重力分力以及雨滴沿坡面方向剪切力分力大小均未发生显著变化，因此得出的结论具有可比性。由于得出的试验结论一致，可以肯定无论是微观地形还是野外实地研究，发生改变的其他因素均未导致塿土在坡度发生变化时，溅蚀量发生改变。因此，该研究可为进一步分析黄土高原区塿土溅蚀的内在机理提供参考。