Lateral water flow and volumetric water content with rainfall for soils in a suburban Quercus acutissima forest in Nanjing

实验区在南京市东善桥林场铜山分场，31°35′~31° 39′ N，118°50′~118°52′E，年均气温为15.1 ℃，年降水量为1 100 mm，年日照2 199 h，无霜期229 d。试验地选择45年生的麻栎Quercus acutissima林，郁闭度为0.8，密度425株·hm^-2，平均树高16.8 m，平均胸径24.3 cm，平均冠幅9.7 m。林地面坡度为18°，坡长62 m，土壤为黄棕壤。

采用ECH₂O土壤含水率监测系统，ECH₂O传感器通过测量土壤的介电常数来计算土壤体积含水率。本研究着眼于土壤表层5，15，30，40，60，100 cm等6个层次土壤含水率在不同降雨强度条件下的变化特征。在监测点坡面上坡不同郁闭度50%，80%，100%各设置雨量计1个，设置集水槽3个·郁闭度^-1（200 cm×20 cm×20 cm），集水槽附近除去周围过高的草本植被，使其低于径流槽，将降雨导入称重式自动排液式穿透雨测定系统，最后根据集水槽的面积以及3个集水槽的收集到的穿透雨量换算出林内雨量，降雨量结合雨量筒数据取平均值。于2012年6月29号到12月9号以半时监测1次的频率对观测地的土壤水分进行监测。实验选用翻斗式雨量计远程信息传输系统，1.0 mm降雨量通过卫星平台远程发送信息1次。

环刀规格直径5 cm，选定麻栎林代表性的测定地点，采用环刀法取5，15，30，40，60，100 cm等层次原状土后，用水浸泡一定时间使其达到水饱和称量（M₁），放入烘干箱烘干至恒量（M₂），即可算出各层次的最大体积含水率。最大体积含水率（%）＝[（M₁-M₂）/（1×V）]×100。

土壤中侧向流动的径流成为侧向流，因为无地表径流产生，所以侧向流近似为单位时间、单位体积下土壤含水量变化量与实际降雨量的差值。I_in＝ΔW-P。其中：I_in为侧向流，mm；ΔW为含水率变化量所需要的降雨量，mm；P为实测降雨量，mm。

以时刻为横坐标，以降雨强度（mm·h^-1）和土壤瞬时体积含水率（%）为纵坐标，得到5，15，30，40，60，100 cm等6个不同土壤层次土壤含水率对小雨（24 h内降雨量不超过10 mm）、中雨（24 h内降雨量为10~25 mm）、大雨（24 h内降雨量为25~60 mm）的响应变化曲线。

以2012年10月26日为例（图 1），在降雨之前，初始含水率最大和最小的土壤层次分别是30 cm（11.53%）和60 cm（6.67%）。累计降雨量6.8 mm时，5 cm，15 cm，30 cm等3个层次的土壤含水率曲线发生了明显的变化，40 cm，60 cm，100 cm等3个层次的土壤含水率无响应，15 cm和30 cm土壤层次含水率曲线同时达到最高峰，而5 cm土壤层次含水率曲线达到最高峰的时间滞后2 h。有响应的土壤层次中，土壤含水率变化量最大和最小的层次分别是30 cm（2.90%）和5 cm（1.48%）；土壤含水率上升速率最大和最小的分别是30 cm（1.45%·h^-1）和5 cm（0.37%·h^-1）。降雨结束，15 cm和30 cm层次土壤含水率曲线有明显的下降过程，而5 cm层次的土壤含水率曲线无响应。降雨开始，0~30 cm层次土壤含水率表现出随深度增加变化梯度增大趋势，响应曲线表现出土壤水在垂直表面非饱和入渗，侧渗壤中流集中于15~30 cm土壤层次，40~100 cm无响应特征。

Figure 1.  Response of VSWC at 6 depths to rainfall intensity on the day 26th, Oct., 2012 for a small rainfall event

以2012年8月29号为例（图 2），在降雨之前，初始含水率最大和最小的土壤层次分别是15 cm（17.57%），60 cm（8.50%）。在累计降雨量在11.8 mm时，5，15，30，40，60 cm等5个土壤层次土壤含水率都有明显的变化特征。5 cm，15 cm，30 cm等3个土壤层次土壤含水率同时达到最高峰，40 cm层次土壤达到最高峰滞后2 h，60 cm层次土壤含水率缓慢上升。有响应的土壤层次中，土壤含水率变化量最大和最小的层次分别是15 cm（5.72%）和60 cm（1.037%）；土壤含水率上升速率最大和最小的分别是15 cm（5.72%·h^-1）和60cm（0.52%·h^-1）。在达到最高峰之后，5 cm和15 cm层次土壤含水率有显著下降趋势，30 cm和40 cm层次土壤含水率则趋于稳定，60 cm和100 cm层次土壤含水率有很缓慢的上升趋势。中雨开始时，0~15 cm层次的土壤含水率表现出高速增长趋势，曲线显示垂直表面非饱和入渗和测渗集中于0~30 cm土壤层次垂直面上。

Figure 2.  Response df VSWC at 6 depths to rainfall intensity on the day 29th, Aug., 2012 for middle rainfall event

以2012年8月21日为例（图 3），在降雨之前，初始含水率最大和最小的分别是15 cm（15.65%），60 cm（7.24%）。第1次降雨，曲线达到最高峰累计降雨量36.8 mm，各层次土壤含水率曲线出现明显上升，降雨3.0 h时5，15，30，100 cm等4个层次土壤含水率达到最大值，降雨3.5 h时60 cm层次土壤含水率达到最大值，降雨4.5 h时40 cm层次土壤含水率达到最大值，峰值最大和最小的是15 cm（23.4%）和60 cm（9.87%），土壤含水率变化量最大和最小的分别是100 cm（9.79%）和60 cm（2.63%）；土壤含水率上升速率最大和最小的分别是100 cm（4.39%·h^-1）和60 cm（0.87%·h^-1）。

Figure 3.  Rsponse of VSWC at 6 depths to rainfall intensity on the day 21st, Aug., 2012 for storm rainfall event

第1次降雨缓慢结束，5 cm，15 cm，100 cm等3个层次土壤含水率曲线有明显的下降过程，30 cm，40 cm，60 cm等3个层次土壤含水率曲线趋于稳定。第2次降雨，曲线达到最高峰累计雨量18.0 mm，5，15，30，40 cm等4个层次土壤含水率曲线同时达到最高峰，第2次降雨结束后，5，15，30，40 cm曲线呈下降趋势，而60和100 cm层次土壤含水率曲线呈上升趋势，60，100 cm层次土壤含水率曲线达到最高峰，相对滞后1.0~2.0 h，峰值最大和最小的分别是15 cm（26.24%）和60 cm（9.87%）；土壤含水率变化量最大和最小的分别是15 cm（10.59%）和60 cm（5.39%）；退水速率最大和最小的分别是5 cm（3.78%·h^-1）和40 cm（0.63%·h^-1）。

表 1为各个层次最大体积含水率。由图 3降雨过程中在各上层土壤远未达到最大体积含水率时，下层土壤含水率有明显的变化特征，表明是土壤水是非饱和下渗，在上层土壤没有达到土壤饱和含水率，存在土壤水下渗过程。

以时刻为横坐标，雨量（mm）为纵坐标，得到5，15，30，40，60，100 cm等6个层次土壤含水率变化理论所需雨量、实际观测点降雨量随时间的变化曲线，理论所需雨量减去实际观测降雨量等于侧向流，按照1.4节公式I_in＝ΔW-P计算得到图 4~6不同级别降雨不同层次土壤侧向流响应过程。

Figure 4.  Response of lateral flow at 6 depths to rainfall intensity on the day 26th, oct., 2012 for small rainfall event

Figure 5.  Response of lateral flow at 6 depths to rainfall intensity on the day 29th, Aug., 2012 for middle rainfall event

Figure 6.  Response of lateral flow at 6 depths to rainfall intensity on the day 21st, Aug., 2012 for storm rainfall event

对土壤入渗有一定影响的因素很多，如降雨强度的大小^[15]，在图 4降雨量极小时，60~100 cm层次的土壤含水率变化微小，看不到明显的曲线特征，而5~40 cm的土壤含水率变化与降雨过程有良好的同步性，此层次显著地出现了侧向流。图 5~6中，中雨和大雨时，0~100 cm层次的土壤含水率变化与降雨过程有一个很明显的响应曲线特征，侧向流对降雨的响应很显著。

由图 4~6显示：麻栎林土壤5，15，30，40，60，100 cm等6个层次土壤含水变化量变化随降雨量增大而变化显著，而在小雨和中雨时侧向流较小，但在大雨情况下侧向流呈高速增长特征。反映了土壤垂直方向上有侧向流的侵入使土壤含水率变化的加剧，小雨过程中产生最大侧向流2.1 mm，中雨过程中产生最大侧向流2.4 mm，大雨过程中产生最大侧向流28.7 mm。

降雨入渗是坡地水文循环的中心环节，入渗过程也是水文循环的最重要的组成部分之一^[16]，水分入渗影响了降雨地表径流、壤中流、地下水的在分配，一直以来是土壤侵蚀、非点源污染、水资源管理等科学研究的重点内容^[17]。研究表明入渗的直接测量是费力费时昂贵的，而且往往涉及大量的空间和时间变化^[17]，所以弄清楚一场降雨的入渗过程是非常不易的。

本试验结果与鲍彪等^[18]在晋西黄土区刺槐Robinia pseudoacacia林地的研究不同，由于麻栎林根系发达和坡面的复杂性在小雨和大雨时南京城郊麻栎林土壤含水率只出现上升期和退水期，没有出现平台期。麻栎林土壤初始含水率普遍出现15 cm土壤层次高，而60 cm土壤层次低的现象。随着雨量的增大，15 cm土壤层次达到最大土壤含水率，上升的速率快。土壤含水率从降雨前到降雨后呈现低-高-缓慢下降的变化趋势，一般其变异程度随着深度的增加而逐渐减小^[14]，但在此次研究中发现，大雨时100 cm土壤层次土壤含水率的变异程度就比40 cm和60 cm土壤层次高，这可能是由于浅层含量较高的土壤水分随着时间变化水分垂直向下再分布所致。在降雨过程中深层次的土壤含水率变化表现出很明显的滞后关系。研究表明：各层次土壤含水率对降雨强度的响应很明显，随着深度的增大，土壤含水率对降雨的响应程度也逐渐降低。

本研究表明：降雨过程中土壤层次不是达到最大含水率而开始下渗，土壤水是非饱和下渗，在上一土壤层次未达到饱和含水率的情况下就开始下渗。每场降雨垂直面上都有不同程度的侧向流流入，造成垂直面上的土壤含水率变化所需降雨量大于实际观测降雨量。

Bodman等^[19]在考察匀质土层下渗过程中土壤水分坡面变化时发现，土壤下渗以表面饱和到水分传递带、湿润带、湿润锋，逐层次下渗。这些是在平面表层水分积水5 cm条件下的土壤渗透过程，和自然降雨过程条件差别巨大。由于林地土壤孔隙特征表现优先流、侧向流、垂直入渗并存的混合流特征，导致林地土壤非饱和入渗在0~100 cm具有普遍性，和以往特定条件下的入渗特征具有显著不同的规律性，表现为小雨—中雨为无饱和界面，大雨以深层次（100~60 cm）逐渐饱和过程，饱和带上方全部为水分传递带。

裸地或者林下裸露地，由于无地表覆盖，容易发生击溅侵蚀，结果导致裸露土壤泥浆溅散、孔隙堵塞，土壤失去渗透储水功能，产生地表径流和土壤侵蚀，土壤入渗为饱和入渗特征。林地发挥水源涵养功能的机制与裸地不同，林地地表覆盖阻断了击溅侵蚀的发生，土壤不会失去渗透储水功能，难以产生地表径流和土壤侵蚀，降雨过程实际小于表层土壤入渗能力，从而表现出极大不同的非饱和入渗特征。

从小雨到大雨，降雨量均小于观测点土壤实际蓄水量，无地表径流产生，入渗均表现为非饱和入渗。土壤含水率与降雨强度呈良好的线性关系。

不同土壤层次土壤含水率对降雨响应表现出随降雨量不同差异显著。降雨强度为0~4.0 mm·h^-1，土壤含水率自表层到30 cm变异幅度增大；降雨强度为12.0~30.0 mm·h^-1，土壤含水率自表层到60 cm层次变异幅度降低特征，而100 cm层次土壤含水率却上升明显。

5 cm和15 cm土壤层次土壤含水率变化趋势与降雨量变化趋势有着良好的同步性，随着土壤层次的加深，土壤含水率变化趋势与降雨量同步性逐渐下降。浅层次土壤含水率对降雨响应明显，随时间深层次土壤含水率过程曲线形状与降雨过程有一定的平移和延长，中层次土壤含水率同时受降雨入渗和地下水位变动影响，深层次土壤含水率主要受地下水位变动控制。