研究小流域位于陕西省延安市湫沟小流域,属于燕沟流域的一条支沟,地理位置为36°29'96″N,109°31'11″E,海拔1 131 m; 淤地坝控制面积为5.23 km²。该淤地坝始建于2003年,由于淤积时间较短,该淤地坝仍然处于有水状态,且并未出现放水现象。因此,所有小流域侵蚀泥沙全部都被拦蓄在该淤地坝内。该地区的气候特征为半干旱大陆性季风气候,多年平均降水量为423 mm,主要集中在7—9月份,且多为暴雨。研究区的土壤侵蚀类型主要是水力侵蚀和重力侵蚀。流域内土壤类型以黄绵土为主,主沟道最底层为裸露的基岩。

1.2.1 样点布设与土壤样品采集 为了研究淤地坝降雨沉积过程,选择延安市的湫沟小流域的湫沟大坝,由于淤积时间较短,该淤地坝仍然处于有水状态。经过前期调查,于2012年选择淤地坝淤积泥沙最厚的点进行淤积泥沙样品的采集,采样使用内径为42 mm的食品级不锈钢空心管进行水上采样。在采样的过程中,将空心管垂直打入,直到原始河床(即基岩位置),随后垂直取出,同时根据采样点所处位置的水面深度,当场对钢管进行焊接和截取,以确保采样钢管内土柱样品的完整性,密封带回实验室。在采样的过程中进行GPS定位,记录采样点的位置信息。

1.2.2 土壤样品处理 将带回的空心管土柱样品取出上部密封物体后,垂直放置,自然风干,待空心管内样品成半干燥状态后沿着管道中心线垂直剖开(该过程要特别小心,尽量减少对空心管内土柱的影响)。将剖开的土柱样品进行0.5 cm/样的切割取样。随后将土样处理后装入密封袋内进行保存,用于土壤颗粒组成的分析,并注明采样时间以及采集的土样所处的剖面深度。

将土柱分割所得的土壤样品平铺在通风口自然风干、除去草石杂物、研磨、过筛(2 mm的筛子)。在100 ml的烧杯中加入2～4 g土壤样品,并加入过氧化氢(H₂O₂)溶液和10%的盐酸(HCL)溶液以消除对土样中的有机质和碳酸盐。最后给每个土壤样品加入六偏磷酸钠(NaPO₃)₆溶液,进行分散处理。对处理的土样使用Mastersizer 2000激光粒度仪进行粒度分析,每个土样重复测量两次,结果取平均值。土壤粒径分类采用1987年的中国制标准,分别为黏粒(<0.002 mm)、细粉粒(0.002～0.01 mm)、粗粉粒(0.01～0.05 mm)、细沙粒(0.05～0.25 mm)、粗沙粒(0.25～1 mm)。本文采用土壤分形维数和土壤质地粗化度综合描述颗粒分布特征。

土壤分形维数可表征土壤粒径的大小组成,还能反映质地的均一程度,分形维数随着土壤黏粒和粉粒含量的增加,其值越来越大; 随着土壤粒径组成中沙粒含量的增大而逐渐减小。计算公式如下:

(V(r<R))/(V_t)=(R/λ_V)^3-D

式中:D表示土壤粒径分形维数; V(r<R)表示粒径小于R的累积体积; V_t表示颗粒总体积; R表示相邻粒级的平均直径; λ_V表示最大粒级土粒平均直径。

土壤质地粗化度指的是坝地剖面沉积旋回层中土壤的沙粒含量(>0.05 mm)与黏粒含量和粉粒含量之和(<0.05 mm)的比值。粗化度越大,则粗颗粒泥沙在土壤粒级组成中的成分越多; 粗化度越小,则土壤粒级组成中的粗颗粒含量越少,土壤受侵蚀程度越大。

本试验原始数据的处理和制图使用Excel 2010软件完成,差异显著性检验使用SPSS 19.0软件完成。

流域坝地侵蚀类型主要以水力侵蚀与重力侵蚀为主,每一场侵蚀性降雨都会形成一次泥沙搬运沉积现象,同时形成一个泥沙沉积旋回层。从图1A坝前沉积泥沙剖面中沙粒含量百分比随深度变化的曲线来看,沙粒在坝前沉积泥沙垂直剖面上呈现高值低值交替出现的周期性变化规律,具体表现为:粗细相间分布,由上到下依次呈现细、粗、细、粗的分布规律,并且粗细分布的厚度差异较大。这主要是由于历次侵蚀性降雨的水动力条件分选沉积作用而造成的,大水大沙下,淤积层厚度大,差异也大,水沙较小时,淤积厚度和差异减小。且由于降雨强度不同,流域内的侵蚀产沙量也不尽相同,进而导致剖面各沉积旋回层厚度的有所差异。

图1 淤地坝坝前沉积泥沙沙粒含量变化特征(A)及沉积旋回厚度变化特征(B)

结合研究小流域降雨特性,根据坝前泥沙粗颗粒(沙粒)的周期变化特征^[22],从沙粒百分含量最小开始到沙粒百分含量逐渐增大到一定程度便开始突然减少为止即为一个完整的沉积旋回。因此,在本研究小流域淤地坝淤积期间,坝内沉积泥沙共存在19个沉积旋回层(图1B)。旋回层的沉积厚度表示单次洪水事件过后,沉积在淤地坝前的泥沙总厚度。在该研究中,不同旋回层的沉积厚度差异显著,平均厚度为12.7 cm,其中,第3旋回层厚度最小(1.5 cm),第17旋回层厚度最大值(37 cm)。不同旋回层的沉积厚度标准差为10.48,变异系数为0.82,属于中等强度变异。从沉积厚度变化趋势分析发现,泥沙沉积较厚的层次均分布在建坝初期,随着沉积时间,除去第1、第4和第5旋回层,其他旋回层厚度均低于10 cm,即坝前沉积泥沙层厚度随着沉积时间显著降低。

2.2.1 不同旋回土壤粒径的变化特征 由表1可知,坝地剖面沉积旋回层中各粒径组成以粗粉粒含量最多,所占比例达到51.73%; 其次为黏粒,所占比例为19.88%; 粗沙粒含量最少,所占比例仅仅为0.25%。研究小流域坝地剖面各粒径含量由大到小分别是粗粉粒(51.73%)、黏粒(19.88%)、细粉粒(17.30%)、细沙粒(10.85%)、粗沙粒(0.25%); 各粒径变异系数由大到小依次为粗沙粒>细沙粒>细粉粒>黏粒>粗粉粒,其中粗砂粒变异系数可达387.04%,而黏粒、细粉粒和粗粉粒仅仅分别为28.93%,36.16%和13.36%,说明粉粒与黏粒在研究小流域的输移、沉积过程中分选性较差。

表1 研究区坝地剖面沉积旋回层各粒径的统计分析

进一步从研究区坝地剖面不同沉积旋回层各粒径的统计分析(图2)可以看出,整个剖面共形成19个沉积旋回层,且每个旋回层中各级沉积泥沙在坝内的空间分布变化较大,变化差异显著(p<0.05)。各沉积层粒径组成均以粗粉粒为主,所占比例在42.85%～62.60%; 不同旋回层中,粗沙粒仅在3个旋回层中出现(第1,9,12旋回层); 细沙粒(3.86%～27.02%)、细粉粒(8.20%～24.22%)和黏粒(12.05%～27.01%)的体积分数接近。而旋回层中粉粒与沙粒极值的空间分布则大致相反,即粉粒含量出现最大值时沙粒含量通常出现最小值,反之亦然。分析表明,退耕还林具有降低土壤细颗粒物质流失的作用,且年际变化显著。

图2 不同沉积旋回层粒径的组成特征

2.2.2 不同旋回土壤粒径分形维数的变化特征 由表2可以看出,该流域坝地剖面各沉积旋回层的土壤分形维数在2.60～2.78,平均值为2.70。其中第8旋回层分形维数最低(2.60); 第14旋回层最高(2.78); 分形维数的线性回归的决定系数在0.91～0.96,相关性极好,说明所有的分形参数在计算上是可信的。但是不同旋回层次的分形维数随着沉积深度有逐渐增大的趋势,即随着沉积时间不同旋回层次沉积泥沙颗粒的分形维数呈现增大趋势,表明随着退耕时间的增加,小流域侵蚀泥沙的细颗粒物质含量增加。但整个剖面上泥沙粒径的分形维数随沉积深度并无显著性差异,这表明该流域土壤质量良好。

为研究分形维数与土壤各粒级含量之间的关系,本文采用皮尔逊相关分析进行分析(表3),分析结果表明,黏粒与分形维数的相关系数最大为0.99; 粗砂粒的相关系数最小为-0.12。分形维数与黏粒、细粉粒、粗粉粒和细砂粒含量在99%检验水平上均呈现出极显著相关关系,而与粗砂粒不存在显著的相关关系; 且黏粒和细粉粒与分形维数呈正相关关系,而粗粉粒和细砂粒与分形维数呈负相关关系。因此,具有较高的黏粒和细粉粒含量(即细物质颗粒含量),较低的砂粒含量的旋回层具有较高的分形维数值。

表2 小流域坝地剖面泥沙粒径的分形维数的统计分析

表3 小流域坝地剖面泥沙粒径的分形维数与各级粒径相关分析

2.2.3 坝地土壤颗粒剖面粗化度的变化 坝地土壤颗粒组成主要以粗粉粒和细沙粒为主,说明径流对小流域侵蚀泥沙颗粒具有一定的分选性。由于土壤侵蚀会引起粒径的分选,因此采用粗化度可以用来说明沉积泥沙颗粒组成的变化。土壤质地粗化度指的是坝地剖面沉积旋回层中土壤的沙粒含量(>0.05 mm)与黏粒含量和粉粒含量之和(<0.05 mm)的比值。粗化度越大,则粗颗粒泥沙在土壤粒级组成中的成分越多; 粗化度越小,则土壤粒级组成中的粗颗粒含量越少,土壤受侵蚀程度越大。

研究结果表明(图3),该研究区小流域坝地沉积泥沙质地粗化度分布范围在0.04～0.37,平均值为0.16。且该坝地接近地表的剖面的质地粗化度均较大,而靠近最底下基岩位置处的质地粗化度都较小,即其泥沙颗粒粒径组成呈现先细后变粗随后又再次细化的趋势。

图3 坝地沉积泥沙质地粗化度的剖面变化