不同施肥条件下土壤水分特征曲线变异及影响因子

　　引言

　　水分是促进植物生长和调节体内外生理生态变化的关键因素,尤其是在黄土高原半干旱区,土壤水分是土壤微生物、植物最主要的水源,是土壤肥力重要的因素之一。土壤水分特征曲线表述了土壤含水率与吸力之间的关系,不仅反映了土壤的持水力,也间接地反映出土壤中孔隙的分布,是模拟土壤水分运动和溶质运移的重要参数,对研究土壤水分的有效性,土壤水分运动溶质运移等有重要的作用。土壤水分特征曲线的影响因素较多、关系复杂。目前尚不能从理论上推求土壤水势与含水率的关系,常采用实验方法测出数据后拟合成经验模型,其中,Van Genuchten 模型因与实测数据拟合程度好而得到广泛应用。因此,通过对土壤水分特征曲线和 Van Genuchten 模型参数空间变异性的研究,可为评估土壤持水能力和释水能力、合理的获取土壤水分运动参数等提供依据。

　　黄土高原生态环境脆弱,降水量低且不均匀,土壤水分严重缺乏,成为影响作物产量的主要限制因素之一。该区域以旱作物为主的雨养农业区,作物供水主要依赖于大气降水,而大气降水对作物的水分供应又需要依赖于土壤调节。因此,探索简单易行且具有良好生态、经济和社会效应的农业管理措施显得尤为迫切和重要。有关黄土高原土壤水分的研究大都集中在土壤水分空间分布、土壤贮水力、农田耗水、水分利用效率等方面,且土壤水分特征曲线方面的研究多集中在林地土壤,而对黄土旱塬区长期施肥条件下土壤水分特征曲线的空间变异性鲜有报道。本文以该地区长期定位实验为对象,深入分析不同施肥条件下土壤水分特征曲线变异及其影响因子,探讨土壤持水、保水和土壤水分的有效性,以期为黄土高原旱区找到可以提高和稳定土地生产力、保证粮食生产的农业管理模式。

　　1、 材料与方法

　　1. 1 实验设计

　　黄土旱原区长期肥料定位实验开始于 1984 年。实验开始前,实验地一直用于种植作物(小麦) 。

　　1984 年,设置空白(CK) 、施氮肥(N,120 kg / hm2) 、氮磷配施(N,120 kg/hm2; P,26 kg/hm2) 、单施有机肥(M,75 t/hm2) 、氮磷有机肥(N,120 kg/hm2;P,26 kg / hm2; M,75 t/hm2) ,按照随机区组,3 次重复的方式布设在田间。种植体系为一年一茬,供试作物为冬小麦(Triticum aestivum L. ,品种 “长武131”系列) 。小区面积 6 m × 4 m,小区间距 0. 3 m,四周保护行 1 m。一般 9 月下旬播种(播量 150 ~190 kg / hm2,行距 20 cm) ,生长期人工及时去除杂草。成熟期人工收割,所有地上部被移出小区,地表残留麦茬高度低于 5 cm。其中,收获整个小区面积的 1/2,用于估算单位面积地上部生物量和籽粒产量。冬小麦收获(次年 6 月) 后土壤休闲,期间圆盘耙机耕(深度 20 cm) 松土蓄墒,播种前再耕翻,准备下一茬冬小麦苗床。作物生长期间及时防治病虫害。氮、磷肥分别为尿素、三料磷肥,在播种前撒施地表后耕翻入土。

　　1. 2 土壤样品采集及土壤水分特征曲线测定

　　2010 年 7 月,在实验地选取小区中用环刀取其原状土,重复 3 次,浸水饱和,测其饱和含水率; 用日本 HITACHI 公司生产的 CR21G 型离心机测定土壤水分特征曲线,通过设置离心机转速以确定相对应吸力,待测试土样达到平衡后,测定该吸力下的质量含水率,离心力增加,依次测定该吸力下的质量含水率,即可得到土壤水分特征曲线。测定的吸力范围为 1 ~1 200 kPa。

　　土壤水分特征曲线是研究土壤持水特性的重要依据之一,揭示了土壤水的能量指标(基质势) 与数量指标(土壤含水率) 之间的内在关系。为了准确地表征各种土壤水分特征曲线的特征,定量地研究土壤的持水能力,土壤物理学界已建立了许多数学模型。在众多模拟水分特征曲线的数学模型中,其中 Van Genuchten 模型因其不仅能够适应多种土壤条件,而且能获得更高程度的模拟效果而被普遍采用,公式为【公式1】

　　本实验选用 RETC 软件,利用 Van Genuchten 模型 (V-G 模型) 模拟不同溶质水溶液的土壤水分特征曲线。

　　1. 3 土壤水分常数及水分有效性划分

　　土壤水分常数和水分有效性可以直接测得,也可由土壤水分特征曲线求得。为了更加直观地显示有机溶质的影响范围和水平,将土壤水势用 pF表示,即土壤水势取 cm 水柱作为单位,并取其绝对值的对数值。

　　依据前人研究结果,一般以 pF为 1. 8 时的土壤含水率为田间持水量,pF为 3. 8 与 4. 2 时的土壤含水率分别为暂时萎蔫系数与永久萎蔫系数。据此计算出田间持水量、暂时与永久萎蔫系数及全有效水(pF= 1. 8 ~ 4. 2) 、速效 水 (pF= 1. 8 ~ 3. 8) 、迟效 水(pF= 3. 8 ~ 4. 2 ) 与无 效水 (pF> 4. 2 ) 的土 壤含水率。

　　1. 4 数据处理

　　统计分析采用 SAS 8. 1 软件和 Excel 2010 进行,当 F 检验显著时,进行各处理的方差分析。

　　2、 结果分析

　　2. 1 不同施肥措施对土壤水分特征曲线的影响

　　不同的施肥措施对土壤持水性能的影响表现出很大的差异性(图 1) 。土壤持水性能表现为有机肥处理(M、NPM) 较高,化肥处理(N、NP) 和不施肥处理(CK) 较低,休闲处理(F) 居中。长期休闲处理(F) 土壤容积密度略有增加(相对 1984 年) ,土壤孔隙度降低,因此,F 处理土壤饱和含水率低于其他施肥处理,其他吸力条件下土壤含水率高于 N、NP 和CK,低于 M、NPM 处理。长期施用化肥(特别是单施氮肥) 破环土壤结构,降低了土壤持水性能。长期单施氮肥处理(N) 在整个吸力范围内,对应的土壤含水率都要小于其他施肥处理。长期施用有机肥(M、NPM) 大大改善了土壤结构,提高了土壤持水性能。

　　因此,在整个吸力范围内,其含水率高于其他处理。

　

　　2. 2 不同施肥措施对土壤水分特征曲线参数的影响

　　为了定量的研究施肥措施对水分特征曲线的影响 ,需要进一步分析V-G模型拟合参数。从表1中可以看出,用 V-G 模型拟合实验结果取得了很好的效果,土壤容积含水率-土壤水势之间的相关关系达到了极显著水平(P < 0. 01) 。在进行拟合过程中,如果剩余含水率(θr) 小于 0. 001 时,RECT 软件程序将以零替代 θr值,拟合结果中各施肥处理均小于 0. 001,故其 θr值均为零。

　　由表 1 可知,长期不同施肥措施对土壤水分特征曲线参数的作用效果也有一定的差异,其中受到影响较为明显的水分特征曲线参数是 θs和 α,对 n的值影响很小,对剩余含水率 θr无明显影响,表明施肥措施对土壤水分特征曲线的位置影响明显,对曲线的斜率无明显影响(P >0. 05) 。

　　长期施肥对土壤理化性质产生很大的影响,这种影响明显地体现在不同施肥对土壤饱和含水率的影响上。休闲处理(F) 最小,有机肥处理(M、NPM)最高,对照(CK) 与化肥处理(N、NP) 居中。与 F 处理相比,小麦的种植不但对土壤结构体不会有明显的破坏作用,且利于结构体内部空气的逸出,使土壤水分的饱和程度有一定的提高。与 CK 相比,化肥处理(N、NP) 无明显改变,有机肥处理(M、NPM) 土壤饱和含水率分别提高了 3%、7%。参数 α 为进气值的倒数,休闲处理(F) 的进气值最大,在小麦连作系统中,氮肥处理(N) 最小,氮磷配施(NP) 最大,其他处理居中。

　　2. 3 不同施肥措施对土壤水分有效性和供水强度的影响

　　土壤释出水能力不仅受土壤基质势支配,还取决于土壤大小孔隙的分配比例。不同施肥处理土壤水分有效性情况见表 2。长期不同施肥处理土壤全有效含水率及速效性有效含水率均发生变化,有机肥处理(M、NPM) 土壤水分的有效性明显高于其他处理,化肥处理(N、NP) 、空白(CK) 、休闲(F) 处理之间无明显差异。在长期定位施肥实验中,CK 处理土壤的速效性有效含水率最低(5. 51) ,与 CK 相比,施肥提高了土壤速效性有效水含量,N、NP、M、NPM分别提高了 7%、5%、10%、19%。有机肥处理全有效含水率为 13. 00 ~ 14. 04,明显高于 CK(11. 75) 。

　　不同施肥处理对土壤的迟效性有效水和缓效性有效水的影响不明显。土壤的速效性有效含水率是突然供水状况优劣的重要标志,长期施有机肥处理土壤的供水能力较强,这与有机肥处理土壤团聚体含量明显增加是一致的。

　　3、 讨论

　　3. 1 土壤持水性与土壤理化特性之间关系分析

　　土壤持水能力高低取决于在一定土壤厚度条件下容积密度土壤容积密度和孔隙的大小]。土壤水分特征曲线是土壤质地、结构和孔隙等土壤物理特性的综合表征。土壤有机质是土壤结构形成和改善的核心物质,土壤在长期施肥过程中,土壤有机质得到不同程度的累积,它们通过改善土壤结构、减低土壤容积密度和增加土壤毛管孔隙度等物理特性对土壤水分和持水性产生作用,其含量的高低与土壤持水功能有重要的关系。本研究中,土壤水分特征曲线参数与土壤理化性状的相关关系分析表明(表 3) ,参数 θs与土壤理化性状有着极显著的相关关系(P < 0. 01) ,参数 α 与土壤理化性状有着显著的相关关系(P < 0. 05) ,参数n仅与有机质有着显著的相关关系。其中,容积密度表现为负相关,有机质含量和孔隙度呈正相关。土壤有机质通过改变土壤其他基本理化特性,进而影响土壤持水特性。因此,有机质对土壤含水率有着显著的控制作用,这与刘效东等研究结果一致。

　　3. 2 不同施肥措施对土壤持水性的影响评价

　　单位基质势的变化引起的含水率变化称为比水容量(C) ,它是土壤释水的量化指标,在评价土壤水分有效性程度方面具有极重要意义。比水容量愈大土壤持水能力愈强,但有效水供给能愈差。

　　从表4中可以看出,不同的施肥处理土壤在5 ~ 20 kPa 段比水容量处于 10- 1级,水分易被植物吸收利用,30 ~ 120 kPa 时,比水容量为 10- 2级,说明此时植物吸收水分难度增加,吸收水量显著减少,水分运动显著减慢,可以根据比水容量 10- 2级出现的吸力值来评价土壤的释水能力和抗旱能力。当土壤比水容量达到 10- 2级时,各施肥处理土壤基质势均出现在 -3. 0 ×104Pa,统计分析结果表明,N 与其他施肥处理间表现出明显的差异性(P < 0. 05) ,其他处理间无明显差异(P > 0. 05) ,但比水容量由强至弱为: M、NPM、NP、CK,其大小代表着土壤供水能力的强弱和土壤水分有效性。韩国军等研究表明,采用草木灰与无机肥两者配施可减少施肥对土壤抗旱性的影响。而本研究有机肥(M、NPM) 处理通过影响小麦地上生物量,增加了土壤水分的消耗,表现出明显的缺水状况。因此,在研究土壤水分亏缺时,应综合考虑其影响因素。

　　4、 结论

　　(1) 长期不同施肥对土壤水分特征曲线表现出明显的差异,但曲线斜率基本一致。化肥处理(N、NP) 破坏了土壤结构,降低了土壤持水性能; 有机肥(M、NPM) 的施入增加了土壤有机质的含量,改善了土壤结构,提高了土壤持水性能。M、NPM 处理饱和含水率 θs最大,F 处理最小,CK 与 N、NP 处理居中,且其间无明显差异。

　　(2) 长期不同施肥对土壤水分的供水强度影响存在一定的差异,有机肥处理(M、NPM) 土壤水分的有效性明显高于其他处理,N、NP、CK 处理之间无明显差异。土壤供水能力由强至弱为 N、CK、NP、NPM、M。

　　参考文献:
　　1. 胡伟,邵明安,王全九. 黄土高原退耕坡地土壤水分空间变异的尺度性研究[J]. 农业工程学报,2005,21(8) : 11 - 16.