红壤干旱指标表达模式与阈值农业工程学报.docx
- 文档编号:29967541
- 上传时间:2023-08-03
- 格式:DOCX
- 页数:14
- 大小:113.83KB
红壤干旱指标表达模式与阈值农业工程学报.docx
《红壤干旱指标表达模式与阈值农业工程学报.docx》由会员分享,可在线阅读,更多相关《红壤干旱指标表达模式与阈值农业工程学报.docx(14页珍藏版)》请在冰豆网上搜索。
红壤干旱指标表达模式与阈值农业工程学报
红壤干旱过程中剖面水分特征与土层干旱指标*
陈家宙,吕国安,王石,张丽丽
(华中农业大学农业部亚热带农业资源与环境重点开放实验室,武汉,430070)
摘要:
农田水分管理以及产量评估都需要对土壤作物受旱状况定量化和指标化。
干旱强度和干旱程度结合才能完整地描述土壤作物干旱状况,土壤干旱强度I是土壤剖面失水速率的函数,干旱影响逐渐累积并增强就构成干旱程度D,据此提出了包含I和D二个指数的土层干旱指标表达模式。
通过红壤小区种植玉米并在抽穗期开始设置连续干旱12d-36d等6个不同的处理,研究了红壤干旱过程中剖面水分特征和干旱指标。
结果表明,供试红壤干旱过程中剖面40cm以下含水量下降幅度很小,玉米主要利用了0-40cm土层的水分,监测30-40cm土层含水量的变动情况可以指示玉米受到干旱胁迫的程度。
连续干旱25d后40cm以下土层含水量明显降低,玉米产量也显著下降,此时0-60cm土层和30-40cm土层的干旱程度D均为0.55,可用此指标作为灌溉的依据。
关键词:
土壤干旱;土层干旱指标;干旱强度;干旱程度;红壤
中图分类号:
S274文献标识码:
A
0引言
土壤水分不足是农田系统常见现象,即使在南方红壤湿润区也存在严重的季节性干旱。
土壤水分减少既是大气干旱的结果,又是引起植物干旱的主要原因,因此在干旱预测评价与灌溉管理中,土壤水分状况很关键,是重要的干旱指标。
目前有许多评价土壤-作物系统干旱的指标[1-5],如气象指标(降水量、降水平距、无雨日数等)、土壤含水量指标(植物可吸收水、土壤含水量或水势等)、作物形态和生理指标(茎直径、叶水势、气叶温差、气孔导度、细胞汁液浓度等)、综合指标(水分供需比、作物水分胁迫指数)等。
有的指标定性为主,适合于评价大区域的干旱状况[2],有的则定量精细,可用于评价单个地块的干旱程度,但能准确识别各级干旱的指标很少。
土壤含水量指标是目前研究比较成熟,且能较好反映作物旱情状况的指标,它主要根据土壤含水量相对于田间持水量的百分数多少来判断土壤作物干旱状况,其具体数值随作物类别、品种及生长阶段而变化[1,3]。
但是静态的土壤含水量包含的信息有限,以相对含水量或绝对含水量作为干旱指标是有欠缺的。
适度的土壤干旱在有些作物上不会导致明显减产,在特定的作物和栽培条件下,土壤干旱程度存在一个临界值[6-11],超过该临界值作物生长受到明显抑制,生理指标和产量显著下降,本文称这种土壤干旱程度为土壤干旱阈值。
确定干旱阈值对指导灌溉具有很重要意义,描述土壤干旱阈值必须有合适的土壤干旱指标。
在众多指标中,由Idso和Jackson等[12,13]提出的一个建立在冠气温差基础上的作物水分胁迫指数(CWSI)使用最广泛,得到了很多研究者的认同[6,14,15],理论上,CWSI=1-(实际蒸散/潜在蒸散),由于其计算因子取值非常困难,实际常应用经验模式[11],CWSI=(ΔT-ΔTmin)/(ΔTmax-ΔTmin),其中,ΔT、ΔTmin、ΔTmax分别是实际的冠气温差(Tc-Ta)、充分供水的冠气温差、极限缺水的冠气温差。
但CWSI也存在一些不足,
(1)即使采用经验公式,ΔTmax、ΔTmin的测量还是困难[15];
(2)CWSI值在一些条件下超出0~1范围[16,17];(3)受气象变化影响严重,例如,在同一天不同时间测量,CWSI变化于0.10~0.35[18],给指标的确定带来困难;(4)偏重于作物对干旱的反应,作物冠层覆盖差或裸土条件下,无法评价土壤干旱状况;(5)不能确定灌水量。
此外,单一的CWSI指数值包含的干旱历史信息少,不能很好地反映过去一段时间的干旱状况。
本文首先以土壤水分变化为中心,提出了一个新的土层干旱指标,包括干旱强度I和干旱程度D,作为对土壤含水量指标的改进和CWSI的补充,来定量表达土壤干旱状况;然后研究了红壤干旱过程中剖面水分动态和玉米产量,确定了玉米对红壤干旱产生阈值反应时的土壤干旱程度D值即土壤干旱阈值,以期为评估土壤作物干旱状况并为灌溉时机的选择提供依据。
1土层干旱指标模式的建立
土壤干旱是一个逐渐发展的过程,土壤水分散失速度就是干旱强度,其大小与作物目前生长状况没有直接固定关系;干旱逐渐积累达到一定水平就是干旱程度,即目前的干旱状况,反映了作物受干旱影响的程度。
因此土层干旱指标包含干旱强度I和干旱程度D二个指数。
一方面,对干旱强度而言,土壤-植物干旱的实质是植物失水速率超过根系吸水速率,它们分别与土壤失水速率和土壤供水速率呈正相关,故土壤-植物干旱强度可用土壤失水与供水关系表示。
对某一个土壤层次,本文定义土层干旱强度I如下:
(1)
其中,土层失水包括蒸发、根系吸收(蒸腾)、向下分配等全部水分减少量,与大气蒸发力、作物生长状况和土壤持水能力相关;土层供水与土壤性质(如导水率)和持水状况相关;f(w)表示土层失水与供水关系函数。
土层失水(无论是向上蒸散还是向下分配)速度越快,表明土壤干旱强度越大,能够反映土层水分损失速度的函数可以作为f(w),本文用经验公式表示如下:
(2)
式中,a为一回归经验参数,其计算过程是:
无水分胁迫时土层有效水贮量为x0(mm,即用田间持水量减去萎蔫含水量,田间持水量时土壤无水分胁迫),土层初始有效贮水量为x1(mm),干旱开始后,土层每日损失水分wi(mm),当日剩余有效水为xi(mm),每日相对失水量为ri=wi/xi,至该日累积相对失水量为yi=∑ri,干旱n天后,得到两个随时间变化的序列,即剩余有效贮水量序列X=(x1,x2,……,xn)和累积相对失水量序列Y=(y1,y2,……,yn),进行回归得到对数方程
,回归系数即为参数a。
试验数据和模拟数据证明,一般情况下回归方程的决定系数R2>0.99,Y与X呈极显著负相关,a<0,I的计算值始终落在[0,1]区间。
进一步分析发现,回归系数a具有以下特征:
(1)a<0;
(2)土层连续n天无失水时,a=-1,失水过程│a│>1,得水过程│a│<1;(3)失水速率越慢(wi越小),│a│越接近1,相反越远离1。
可见,a很好地反映了整个过程的干旱速率。
此外,从回归方程
可以看出:
(3)
因此a值反映了土层失水与供水关系,本质上就是函数f(w),干旱强度I只是对其进行了指数化。
另一方面,对干旱程度而言,土壤-植物系统干旱是一个逐渐发展的过程,植物受干旱胁迫的伤害也是逐渐累积的,因此有理由认为土壤失水与供水的历史决定了当前的水分供需关系和干旱程度,即干旱程度是干旱强度累加值的函数,据此,构造干旱程度指标D如下:
(4)
式中,x0为土层无水分胁迫时有效水贮量,x1为土层初始(开始监测含水量时)有效贮水量为,∑I为累积干旱强度,计算时可以取合理的任意天数。
D的计算值始终落在[0,1]区间,该值越大,表明作物受到的土壤水分干旱胁迫程度和伤害越严重。
在计算I和D时,根据监测的土壤含水量的资料来确定计算土层和时间序列,如果某土层含水量低于萎蔫含水量D则取值1。
土层一般取根层,也可以分层计算;时间可以1d为间隔,也可以多日为间隔。
2材料与方法
为了验证与应用土层干旱指标I和D,2005年5~10月在湖北咸宁进行了田间玉米干旱胁迫试验。
供试土壤为第四纪红色粘土发育的红壤,质地为粘土,pH6.52,表层有机质、速效氮磷钾含量分别为3.11g•kg-1、22.5mg•kg-1、2.94mg•kg-1、143.77mg•kg-1,有机质含量很低,0-60cm各层容重变化于1.14~1.46g•cm-3,田间持水量变化于0.3175~0.3547cm3•cm-3,萎蔫含水量变化于0.2033~0.2431cm3•cm-3。
试验地块平坦,用不透水铝塑板分隔为面积相同的矩形小区,每个小区2.7m×1.2m=3.24m2,铝塑板插入地下110cm,露出地表10cm,小区之间地表和根系层互不串水。
小区内种植玉米,播种前和出苗后一段时间内,所有小区进行相同的管理,保证各小区土壤水分状况相同并且没有水分胁迫。
在玉米抽穗期(7月15日)开始通过遮雨棚防雨和人工控制灌水进行土壤水分胁迫,设置6种干旱水平,即分别连续干旱12d,21d,25d,28d,33d,36d后恢复正常灌水,各处理分别记为D12,D21,D25,D29,D33,D36,重复3次。
在玉米生长过程中,不定期多次跟踪测量玉米形态与生理指标,收获后考种。
每个小区内各埋入1根PR1型(英国Delta-T公司生产)长度100cm的分层测水仪探管,每日上午测量10、20、30、40、50、60cm各层土壤体积含水量,计算各层总贮水量和有效水贮量xi,前一日土壤贮水量与当日贮水量之差为失水量wi,失水量与当日有效水贮量之比为相对失水量ri。
土壤水分特征曲线用压力膜仪测量,并以-30KPa基质势对应的含水量作为田间持水量,以-1.5MPa基质势对应的含水量作为永久萎蔫点,超过萎蔫点的含水量视为有效水含量;其它项目测量用常规方法;数据统计分析和曲线回归用SPSS软件进行。
3结果与分析
3.1土壤不同干旱程度下玉米生物量和产量
试验结果表明,土壤干旱对玉米地上部生长和产量都有明显的影响,随着干旱程度加强,对玉米生长的抑制也加强,玉米对土壤干旱存在明显的阈值反应。
图1所示数据统计分析显示,对生物量干重,处理D12、D21、D25、D29与处理D33、D36之间存在显著差异;对生物量鲜重,处理D12、D21、D25与处理D29、D33、D36之间存在显著差异。
处理D29的生物量干重和生物量鲜重分别为处理D12的89%和83%。
这说明,在土壤干旱25-29d的时候,玉米生物量出现明显的下降,表现出对干旱的阈值反应,其中,鲜重比干重更明显,而且阈值出现更早,这也表明鲜重对干旱的反应比干重对干旱的反应更敏感。
图2结果表明,玉米穗干重和籽粒产量有与生物量类似的明显的干旱阈值反应。
土壤干旱使玉米产量降低,处理D21的穗干重和籽粒分别为处理D12的93.1%和93.3%,而在干旱达到25d之后,很快降到处理D12的80%以下。
统计分析表明,对玉米穗干重和籽粒产量,处理D12和D21之间无显著差异,而与处理D25、D29、D33、D36之间存在显著差异。
结果表明,在本试验条件下,土壤干旱21-25d的时候玉米产量出现阈值反应。
图1土壤不同干旱处理下单株玉米的生物量
Fig.1Individualplantbiomassofcornindifferentsoildroughttreatments
图2土壤不同干旱处理下单株玉米的产量
Fig.2Individualplantproductsofcornindifferentsoildroughttreatments
3.2干旱过程中土壤剖面水分特征
图3列出了不同处理在干旱开始时和结束时0-60cm剖面含水量的分布情况,总的趋势是随着干旱进行,各土层含水量逐渐降低,但仍有一些有意义的特征。
含水量首先从表层开始降低,0-10cm和10-20cm土层含水量降低速度很快;干旱持续21d之后,20-30cm和30-40cm土层含水量开始明显降低;继续干旱持续到29d,土层60cm处含水量几乎稳定在田间持水量(0.35cm3•m-3左右)不变,直到33d才明显降低,但仍然保持了0.30cm3/m3左右的较高含水量水平。
这一结果表明,作物主要利用了0-40cm土层的水分;当30cm处含水量明显降低的时候(处理D25),作物已经受到了干旱胁迫的影响,产量下降;当60cm处含水量明显降低的时候(处理D33),土壤干旱对作物已经造成严重影响。
总之,供试红壤干旱过程中剖面水分最明显特征是40cm以下水分难以利用,这与红壤的持水和供水特征有关。
各土层含水量不同的变动特征对土壤干旱程度评价有不同的意义。
图3显示,20cm处土壤含水量变动快,变化幅度较大,而且出现不稳定的波动,即不同处理的20cm土层含水量之间没有表现出明显的规律;60cm处土层含水量变化迟缓,幅度小,对干旱反应不敏感;只有土层30cm或40cm处含水量的变动情况可以指示作物受到干旱胁迫的程度,因此监测30-40cm含水量比监测表层含水量更有意义。
图3不同处理在干旱开始和结束时土壤剖面含水量分布
Fig.3Soilprofilewatercontentatthebeginningandendofsoildroughtindifferenttreatments
3.3干旱过程中土层干旱指标的特征
为了说明干旱程度D与土层干旱状况的关系,图4列出了处理D36在干旱过程中各土层干旱程度D的变化情况。
总体看来,随着干旱时间延长,各土层干旱程度D几乎呈直线增加。
上层土壤(0-10cm、10-20cm、20-30cm)的D增加很快,干旱持续时间不到36d的时候,这几个土层的干旱程度已经达到最大值1;30-40cm土层在干旱初期由于不失水或失水很少,D值缓慢增加,干旱持续22d左右之后,D值增加速度变快,说明此时这2个土层开始明显失水,与图3结果一致;40-50cm和50-60cm土层的D值变化趋势很相似,都是随干旱持续而缓慢增加,2条曲线几乎平行。
引人注意的是,50-60cm土层的D值明显小于其它土层,并且在长达36d的干旱过程中增加幅度很少,表明该层土壤水分状况受干旱过程影响不大,这一方面说明该土层始终保持了较高的含水量,另一方面说明该土层水分向上运输少,向上层土壤补充水分少,如果作物根系不能达到该层,则该层土壤水分很难被作物吸收利用。
这一结论与许多研究结果(即干旱季节红壤50cm以下含水量并不低但难以利用)一致,这从侧面说明土层干旱程度D很好地表达了土壤的干旱状况。
图4土壤干旱过程中处理D36各土层干旱程度D的变化
Fig.4ThedevelopmentofsoillayerdroughtdegreesintreatmentD36duringsoildrying
各处理在干旱结束时(抽穗期)不同土层的干旱程度D列于表2,结果表明,随着干旱时间延长,各土层干旱程度D都增加,达到最大1为止(此时土层平均含水量低于凋萎系数),D值的大小与实际土壤干旱程度大小一致。
不同土层有不同的表现特征,对于0-10cm表层,D值全部达到了1,实际上在本试验条件下,该层土壤失水很快,连续干旱7d时0-10cm贮水量已经低于凋萎含水量,作物很难再从该层吸收水分,但不排除接近10cm附近的含水量高于凋萎系数,作物仍可吸收该处水分,D反映的是该层平均状况。
对于10-20cm和20-30cm土层,随干旱延长D值增加也很快,分别在干旱33d和36d之后达到1,从而不能继续表达土壤干旱状况;而30-40cm和40-50cm土层的D值则平稳增加,在干旱持续36d后D值也没有达到1,这二个土层在很长的时间、很大的干旱范围内其干旱程度可以用D值表达。
50-60cm土层的D值变化很缓慢,D值变化幅度很小。
上述各土层D值的变化特征与上文描述的各土层含水量变化特征有相似之处,表明干旱程度D值比较真实地反映了土层的水分变动情况。
表1各处理干旱结束时(玉米抽穗期)不同土层干旱程度D
Table1Soillayerdroughtdegreesofsixtreatmentsattheendofsoildrought(maizeheadingperiod)
处理
土层干旱程度D
0-10cm
10-20cm
20-30cm
30-40cm
40-50cm
50-60cm
0-60cm
D12
1.0000
0.3451
0.3009
0.2283
0.1715
0.0915
0.3013
D21
1.0000
0.5724
0.4873
0.4838
0.3292
0.1082
0.4897
D25
1.0000
0.7881
0.6520
0.5516
0.3536
0.1469
0.5502
D29
1.0000
0.7979
0.6866
0.5910
0.3914
0.1842
0.6312
D33
1.0000
1.0000
0.8315
0.6506
0.4472
0.2664
0.8367
D36
1.0000
1.0000
1.0000
0.6768
0.4625
0.3174
1.0000
对于干旱强度I,由于各处理是在相同的条件下进行连续干旱,外部气象条件和种植的作物一致,不同处理之间比较意义不大。
不同干旱处理之间干旱强度I没有明显的规律;但不同土层之间的干旱强度存在差异,从表层往下土壤干旱强度降低,表明越往下层土壤失水速度越慢,如图5所示。
试验结果显示,0-60cm不同土层干旱强度变化于0.016-0.16之间,平均为0.11,这就是在本试验条件下红壤-玉米的干旱强度。
图5干旱处理D36不同土层在干旱期间平均干旱强度I
Fig.5Theaveragesoillayerdroughtintensities(I)duringsoildryingintreatmentD36
4讨论与结论
据估计,各国应用的干旱指标有50多种,农业干旱指标也有很多,但能准确识别各级干旱的并不多。
本文提出用两个指标结合来表达土壤干旱状况,干旱强度I表示的是土层目前向干旱发展的快慢,其大小并不反映作物受旱状况;干旱程度D则表示土层目前已经形成的干旱现状,定量地表达了土壤的干旱程度。
D很好地反映了干旱历史对当前干旱程度的影响,D是过去的干旱程度与当前干旱强度的综合,而且是定量的、连续的,较好地反映了整个干旱过程,其大小反映了作物受旱状况,可以评价土壤作物干旱状况从而用于指导灌溉。
可见,用两个指标结合来定量表达土壤干旱状况比用单一的指标更加符合干旱的实际过程。
试验结果统计分析表明,从玉米抽穗期开始,土壤连续干旱25d之后玉米产量下降明显,可以初步确定土壤干旱25d为玉米抽穗期至成熟期干旱阈值出现的时间,因此为了保证玉米较高产量,必须在连续干旱25d之前灌溉。
但是,以持续干旱的天数作为土壤-作物干旱临界值指标是不恰当的,因为在不同的气象条件下,土壤失水速度不同,在持续相同的干旱时间后,土壤-作物实际干旱状况不同,干旱临界值出现的时间就会发生变化,而采用土层干旱指标I和D可以在很大程度上避免这种情况。
根据不同处理的玉米产量,可以得出的结论是,以0-60cm土层为计算基础,土层干旱强度I为0.11的情况下,红壤-玉米的干旱阈值为连续干旱25d,此时对应的干旱程度D为0.5502,即红壤-玉米干旱阈值D为0.55左右,在此之前灌溉玉米不会显著受到干旱的影响。
因此,将I和D两个干旱指标结合,用于指导灌溉将比含水量更加可靠;而且D的计算也只需测定土壤含水量,是将土壤含水量动态进一步定量化和综合化,指标值更加简洁和直观,实用性没有降低反而有所增强。
从理论上说,整个根系层的水分状况都对作物生长有影响,D的计算应该测量整个根系层土壤含水量,这样监测工作强度较大。
但实际上,从图4和表2可以看到,整个0-60cm土层的D值变化情况与20-30cm和30-40cm土层D值变化情况比较接近,说明20-30cm或30-40cm土层水分变化状况代表了0-60cm整个土层水分变化状况;从图3各个土层水分变动分析得到,土层30cm或40cm处含水量的变动情况可以指示作物受到干旱胁迫的程度。
因此,监测土壤20-30cm或30-40cm层次水分状况可以很好地反映土壤-作物干旱状况,不需要监测更多土层的含水量,这样D值的获取更加容易。
我们认为30-40cm土层更好一些,因为在很长干旱时间后该层的D值都没有达到1,还可以定量表达更加干旱的情况,而且在干旱25d的时候,30-40cm土层的D值(0.5516)与0-60cm土层的D值(0.5502)很接近。
在土层含水量低于凋萎系数后,采用本文指标计算的D值就达到了1,即D不能描述土壤含水量小于凋萎系数时的情况,只限于作物可以生长的情况,而对土壤表层而言,含水量可能会小于凋萎系数,如果将计算过程中的凋萎系数修改为风干系数或者剩余含水量,可以扩大D值的描述范围,对此需要进一步研究。
此外,本文仅考察了土壤持续干旱的情况,干旱和降水交替的自然状况没有讨论,在此情况下D的计算方法有所不同。
参考文献
[1]王密侠,马成军,蔡焕杰.农业干旱指标研究与进展[J].干旱地区农业研究,1998.16(3):
119-125.
[2]齐述华,张源沛,牛铮,等.水分亏缺指数在全国干旱遥感监测中的应用研究[J].土壤学报,2005.42(3):
367-372.
[3]胡晓棠,李明思,马富裕,等.膜下滴灌棉花的土壤干旱诊断指标与灌水决策.农业工程学报[J],2002,18
(1):
49-52
[4]张寄阳,段爱旺,孟兆江,等.不同水分状况下棉花茎直径变化规律研究.农业工程学报[J],2005,21(5):
7-11
[5]张寄阳,段爱旺,孟兆江,等.基于茎直径微变化的棉花适宜灌溉指标初步研究.农业工程学报[J],2006,22(12):
86-89
[6]蔡焕杰,康绍忠,张振华,等.作物调亏灌溉的适宜时间与调亏程度的研究.农业工程学报[J],2000,16(3):
24-27
[7]张喜英,裴冬,由懋正.几种作物的生理指标对土壤水分变动的阈值反应.植物生态学报,2000,24(3):
280-283
[8]Sadras,VOandMilroy,SP.Soil-waterthresholdsfortheresponsesofleafexpansionandgasexchange-Areview[J].FieldCropsResearch,1996.47:
253-266.
[9]斐冬,孙振山,陈四龙,等.水分调亏对冬小麦生理生态的影响.农业工程学报[J],2006,22(8):
68-72
[10]刘晓英,罗远培.水分胁迫对冬小麦生长后效影响的模拟研究.农业工程学报[J],2003,19(4):
28-32
[11]Kacira,M,Ling,PPandShort,TH.Establishingcropwaterstressindex(CWSI)thresholdvaluesforearly,non-contactdetectionofplantwaterstress[J].TransactionsoftheASAE,2002.45(3):
775-780.
[12]IdsoSB,JacksonRD,ReginatoRJ.Remotesensingofcropyields[J].Science,1977,196:
19-25.
[13]JacksonRD,ReginatoRJ,IdsoSB.Wheatcanopytemperature:
Apracticaltoolforevaluatingwaterrequirements[J].WaterResourceResearch,1977,13:
651-656.
[14]Alderfasi,AAandNielsen,DC.Useofcropwaterstressindexformonitoringwater
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 红壤 干旱 指标 表达 模式 阈值 农业工程 学报