赵福年 ,刘 江 ,张 强,3** ,王润元 ,王鹤龄 ,张 凯 ,赵 鸿,齐 月,陈 斐
(1. 中国气象局兰州干旱气象研究所/甘肃省干旱气候变化与减灾重点实验室/中国气象局干旱气候变化与减灾重点开放实验室 兰州 730020;2. 甘肃农业大学资源与环境学院 兰州 730070;3. 甘肃省气象局 兰州 730020)
从全球范围来看,与其他自然灾害相比,干旱是影响农业生产最为严重的气象灾害[1-2]。而影响农业生产的因素众多,既有气象、土壤等环境因子,又有作物自身品种及栽培管理方式等因素[3],因此,如何准确定量评估干旱对农业生产的影响是当前干旱研究的热点问题。
在气象、资源等相关领域,一般将与农业生产过程相关的干旱定义为农业干旱[4]。农业干旱泛指作物生长期间因土壤水分供给不足,无法满足作物生长发育,导致产量降低的现象[5]。但因缺乏动态、定量评估农业干旱过程的手段与工具,许多研究者用研究区代表性作物产量的波动来表示干旱的严重程度[5]。该评估方法将作物生长发育过程视为黑箱,忽略了作物生长发育与供水间的动态变化关系,很难准确掌握农业干旱过程特征及其对作物的影响[3]。
在农学和植物生理相关学科中,大量水分胁迫与作物生理特征相互关系的试验表明,作物生理生态指标对作物水分供给存在阈值响应[6],即作物受到水分胁迫后,生理生态指标不会立即减小,而是当水分供给降低到一定阈值后,作物生理生态指标才会随水分胁迫迅速降低[7]。根据阈值响应特征,有研究者从作物不同叶片气体交换指标对土壤水分的响应及其调控机制存在差异为基础,提出了生长发育期间作物受旱的监测方法,并确定了作物干旱阶段划分标准[8],这为定量评估作物干旱动态过程提供了参考。同时,也有部分研究者利用与作物叶片气体交换紧密相关的指标,如冠层温度、冠层与空气温度差、作物水分胁迫指数(以冠层温度与空气温度的差值与空气饱和差的相互关系为基础建立的指数)等为依据构建作物干旱监测工具,这为生长期内作物受旱状况判定以及灌溉时间确定提供了理论依据[9-10]。
然而,作物遭受干旱胁迫后作物产量未必下降,这与干旱发生的生育时期、干旱持续时间、作物自适应性等多种因素相关。许多研究结果表明只有在生长关键期遭受干旱胁迫才会导致作物明显减产[11]。Lv 等[12]研究发现与干旱程度相比,干旱持续时间对作物生长发育影响更大。还有研究表明,作物生长前期,一定程度的干旱胁迫可促进作物根系生长,反而有利于作物产量提高[13]。尽管作物叶片气体交换包括与作物生长发育及产量形成紧密相关的光合和蒸腾两个生理过程,但作物生长期因干旱而引起叶片气体交换变化的特征与作物最终产量之间的相互关系并不明确。而以生育期作物叶片气体交换过程为基础,构建干旱指标是否能准确反映作物最终产量因干旱而降低的状况依然存在不确定性。因此,有必要开展作物叶片气体交换对干旱响应特征及其对作物产量形成影响的研究。
不过由于传统田间试验时间跨度长且投入大,缺乏叶片气体交换过程与作物产量长时序同步动态观测,限制了作物叶片气体交换与最终产量因干旱而变化的对比分析。而作物生长发育机理模型为该类研究提供了一种可实现的途径。如,Sadok 等[14]利用作物生长模型分析发现,对干旱敏感性(以作物叶片气体交换对水分的响应阈值大小表示)不同的作物品种在相同气候条件下产量存在明显差异。Sinclair 等[15]以作物生长模型模拟,发现对干旱敏感(作物叶片气体交换对水分的响应阈值相对较大)的大豆(Glycine max)品种在水分限制条件下更容易高产。鉴于此,为解决上述科学问题,本研究以春小麦(Triticum aestivum)为例,掌握其叶片气体交换对干旱的响应特征,构建和应用春小麦生长发育机理模型,设置不同水分供给情景,明确春小麦叶片气体交换与最终产量对干旱响应阈值的差异,从而为量化农业干旱过程及其影响提供理论依据。
叶片气体交换试验于2014 年、2015 年和2017年以桶栽方式在中国气象局兰州干旱气象研究所定西干旱气象与生态环境试验站(104.37°E,35.35°N,海拔1920 m)开展。该站位于甘肃省定西市安定区,年平均气温6.3 ℃,年降水量为386 mm,且季节分布极不均匀,主要集中在7-9 月。土壤为黄绵土,碱性(pH 为7~8),0~50 cm 土层平均容重1.15 g·cm-3,田间持水量26.8%,萎蔫系数5.5%。供试土壤采集于大田0~30 cm 土层,经风干处理后过5 mm 筛,以剔除石块、根茬及其他杂物。试验所用桶直径29 cm,深度45 cm。每桶装土14 kg,同时施一次底肥,分别为尿素9 g,过磷酸钙17 g。供试作物为春小麦,2014 年和2015 年品种为‘定西新24 号’,2017 年品种为‘定西30 号’。春小麦分别于2014 年3 月21 日、2015 年3 月17 日及2017 年3 月25 日播种。试验均设置两个水分处理: 一个为对照,充足灌水处理,保证春小麦生长期间土壤含水量始终保持在田间持水量的70%以上;另一个为干旱胁迫处理,在拔节期(2014 年和2017 年试验)或开花期(2015 年试验)持续不灌水直至小麦叶片萎蔫。试验所设处理各6 组重复,共12 桶,播种后至不同水分处理前保持12 个桶管理方式相同,土壤水分供给一致。试验观测项目包括: 桶播后65~74 d (2014 年)、68~77 d (2017 年)及85~95 d (2015 年)每日土壤水分(换算为有效含水量,即土壤实测含水量与萎蔫系数的差值与田间持水量与萎蔫系数差值的比值,该值为无量纲值)、春小麦叶片光合生理过程(包括净光合速率、蒸腾速率、气孔导度、胞间二氧化碳浓度等指标)。土壤有效含水量的计算公式及光合生理过程测量方法与步骤参考文献[8]。
春小麦产量数据分为两类,一部分数据来自田间定位试验,包括: 1)长期定位试验(1987-2018 年)。该试验于定西农业气象试验站开展(104.62°E,35.58°N,海拔1896.7 m),试验只设1 个雨养处理,4 次重复。每年3 月中旬至下旬播种春小麦,播前一次施足底肥,包括农家肥和氮肥(播种及施肥信息参见表1 与表2)。春小麦生长期间观测发育期天数、植株密度以及地上部生物量和产量。观测方法按照中国气象局农业气象观测规范进行。2)分期播种试验。试验于2010 年在定西农业气象试验站进行,试验设置4个播期,分别在3 月8 日、3 月18 日、3 月28 日及4 月7 日播种,每个处理3 次重复。供试作物品种为‘定西新24 号’,各播期春小麦均为雨养处理,播种时一次施足底肥,农家肥15 000 kg·hm-2、尿素150 kg·hm-2。试验期间,观测春小麦发育期天数、植株密度以及地上部生物量和产量,观测方法按照中国气象局农业气象观测规范进行。3)播前灌水试验。试验于2017 年在中国气象局兰州干旱气象研究所定西干旱气象与生态环境试验站开展。该试验设置6 个不同播前灌水处理,每个处理3 次重复。播前1 周各处理分别一次灌水0 mm、20 mm、40 mm、60 mm、80 mm以及100 mm。春小麦供试品种为‘定西新24 号’,于3 月24 日播种,播前施足底肥,包括农家肥15 000 kg·hm-2、尿素150 kg·hm-2。春小麦生长期间,按照中国气象局农业气象观测规范观测发育期、植株密度以及地上部生物量和产量。4)春小麦不同生育期补水试验。试验于2016 年在宁夏彭阳旱地农业试验站进行(106.39°E,35.55°N,海拔1782 m)。供试春小麦品种为‘定西40 号’,试验于3 月25 日播种。播种前基施有机肥30 000 kg·hm-2、尿素225 kg·hm-2,保证肥料供给充足。试验设6 个水分处理,每个处理3次重复,包括在春小麦雨养条件基础上补充灌水0 mm、30 mm、60 mm、90 mm、150 mm 以及完全遮雨处理。生长期补水处理分3 次进行,每次补充水量为总量的1/3,即每次分别灌水0 mm、10 mm、20 mm、30 mm 和50 mm,分别在5 月4 日(六叶期)、5 月22 日(拔节期)和6 月15 日(灌浆期)进行。试验进行期间记录春小麦发育期,并在收获后测定地上部生物量和产量。长期定位试验和生育期补水试验具体观测步骤及试验详情参考文献[16],分期播种试验详细介绍参见文献[17]。
表1 定西农业气象试验站春小麦不同年份播种及施肥信息Table 1 Planting and fertilizing details for spring wheat at Dingxi Agrometeorological Experimental Station
表2 1987-2017 年试验区春小麦发育期信息(月-日)Table 2 Earliest and latest dates of growth stages of spring wheat from 1987 to 2017 (month-day)
另一部分数据为公开发表文献收集获得的黄土高原西端半干旱区春小麦相关试验数据。采用中国期刊全文数据库(CNKI)和Web of Science 数据库通过检索“春小麦”“半干旱区”“定西”及“黄土高原西端”等关键词获得,检索时间为1995-2015 年。春小麦播种信息、处理方式、生育期、地上部生物量及最终产量等数据直接由检索获得的文献文字或表格中提取,或由GetData 图像处理软件(http://getdatagraph-digitizer.com)从文献图中提取。文献中试验包括不同水分、肥料及覆盖方式的处理,研究仅选不同水分、中等和高水平肥料处理及传统耕作方式的试验结果进行分析。文献与试验信息如表3 所示。
表3 文献所收集的春小麦生长发育及产量数据集信息Table 3 Information of spring wheat growth and yield collected from references
春小麦生长模型根据Sinclair 等[24-26]提出的方法构建。该模型可以模拟3 种产量水平,即: 光温限制下的作物潜在产量、光温水限制下的作物雨养产量以及光温水肥限制下的可获得产量。该模型由发育期、叶面积变化以及干物质形成和分配4 个模块组成。发育进程由春小麦生育期内>0 ℃有效积温决定,春小麦干物质积累过程由太阳辐射和光能利用率确定,生育期逐日蒸腾由蒸腾效率、空气饱和差与春小麦干物质积累间的关系所决定。在该模型中,干旱对春小麦干物质及产量形成的影响,由干旱过程中春小麦叶片气体交换对土壤有效含水量的阈值响应特征,通过构建分段函数所确定的胁迫系数(Ks)计算获得(图1),春小麦气体交换阈值响应特征由本研究叶片尺度试验观测确定。因此,通过试验观测与模型模拟,该模型可以用来比较分析作物叶片气体交换与最终产量对水分响应阈值的区别。春小麦生长发育及产量形成过程,由Window 7 系统下运行Fortran90 编译的模型完成。模型输入变量为逐日太阳辐射、最高气温、最低气温和降水量等,运行步长为日。模型部分参数如表4 所示。
图1 春小麦生长模型胁迫系数计算Fig.1 Calculation of stress index for spring wheat growth model
表4 春小麦生长模型主要参数及其取值Table 4 Parameters and their values for the simulation of spring wheat using the ACM-Wheat model
本研究中文献收集数据、中国气象局兰州干旱气象研究所定西干旱气象与生态环境试验站播前灌水试验以及宁夏彭阳旱农试验站生育期灌水试验数据用于模型参数化校正。定西农业气象试验站长期观测试验和分期播种试验数据用于模型验证。根据调参数据中春小麦发育期和相应的气象资料,计算获得春小麦不同发育阶段所需的有效积温,作为参数输入模型进行模拟,并通过查阅文献、试验观测以及模型校正小麦生长发育的参数(表4)。模型评价指标包括相对均方根误差(relative root mean square error,RRMSE)、一致性指数(index of agreement,d)、决定系数(R2)以及回归直线斜率等,具体计算步骤参见文献[30]。
选择定西农业气象试验站2014 年逐日气象因子(太阳辐射、最高和最低温度)为模型驱动要素(图2a,2b),小麦播种时间与农业气象试验站2014 年实际播种时间相同(3 月21 日)。为使春小麦在生长关键期(拔节-开花)遭受不同程度的干旱,模拟情景设置5 个不同水分处理,分别为供水充足的T1 (在播后31 d、41 d、46 d、51 d、60 d、65 d 及71 d 补水,共计230 mm)、供水相对充足的T2 (在播后31 d、41 d、51 d、60 d、65 d 及71 d 补水,共计165 mm)、轻度干旱的T3 (在播后31 d、41 d、51 d、60 d 及65 d补水,共计115 mm)、中度干旱的T4 (在播后31 d、41 d 及51 d 补水,共计50 mm)及重度干旱的T5 (在播后32~72 d 无灌溉)处理,补水时间及具体补水量参见图2d。其他时段为保证各处理水分供给充足,人为设置相同的降水条件,即播后3 d 降水5 mm,播后10 d 降水10 mm,播后21 d 降水15 mm,播后76 d降水40 mm,播后80 d 降水40 mm,播后86 d 降水30 mm,播后91 d 降水20 mm,播后96 d 和101 d 降水15 mm,以及播后111 d 降水5 mm (图2c)。
图2 模拟春小麦生长发育过程的气象条件及不同处理补水量Fig.2 Weather conditions and water treatments for simulation of spring wheat growth
小麦叶片气体交换指标与产量对土壤有效含水量的响应关系采用以下公式进行拟合:
式中:S表示应变量拟合值;x为自变量;Smax表示观测获得的平均最大值;a和b为拟合系数,其中系数a反映了所拟合要素随自变量变化上升或下降的速率,系数b反映了所拟合要素数值随自变量快速下降过程由外凸向内凹转变的拐点。曲线和直线拟合、显著性等计算过程均由R 语言完成。
春小麦叶片气体交换对干旱过程存在明显的阈值响应特征(图3)。在土壤有效含水量较大时,随着土壤有效含水量的变化,小麦叶片气孔导度、蒸腾速率以及净光合速率围绕最大值上下波动(表5,参数Smax);而当土壤有效含水量降至0.50 时,气孔导度随土壤有效含水量下降开始快速降低(图3a);蒸腾速率与净光合速率则是在土壤有效含水量降至0.40时,随土壤有效含水量下降而快速降低(图3b,3c)。
图3 春小麦叶片气体交换指标对干旱过程的响应Fig.3 Response of leaf gas exchange for spring wheat to drought
表5 小麦叶片光合生理指标及作物产出对土壤有效含水量响应函数的参数值Table 5 Parameters of response functions for leaf photosynthetic indices and production of spring wheat to available soil moisture
由于叶片净光合速率是作物干物质积累的基础,因此在本研究模型中使用净光合速率对土壤有效含水量的阈值(图3c),作为春小麦干物质受干旱影响的起始点(表4,参数WSSG,即干物质积累对水分响应阈值的取值)。同时,2000 年前大部分春小麦品种为中晚熟品种,生育期所需积温基本相同,2000 年后春小麦品种为中熟品种,发育期所需积温也相近,因此本研究以2000 年为界,将春小麦生长模型参数分为两组(表4),只改变与发育期相关的模型参数,其他参数设置为一致。经参数化后的春小麦生长模型模拟开花期和收获期的数值,发现模拟获得的开花期和收获期分别可解释观测值39.6%与47.2%的变率(图4),且统计检验显著,同时模拟值与观测值间的RRMSE 均小于10%,d>0.75,说明模型能够较好地模拟反映出春小麦发育期变化进程。
图4 春小麦生长模型模拟的开花期(a)与收获期(b)模拟值与观测值(播种后天数)Fig.4 Observed and growth model-simulated dates (days after sowing) of flower stage (a) and harvest stage (b) of spring wheat
图5 是春小麦生长模型模拟的地上部生物量和产量与观测值的比较,发现无论是调参数据还是验证数据,春小麦生长模型对地上部生物量和产量的模拟值均能够解释观测值70%以上的变率,且统计检验极显著,表明模型对环境变量的响应效果较好。而模型模拟值与观测值之间的RRMSE 均小于30%,说明模型整体模拟精度较高。而d≥0.85,模拟值与观测值之间的关系斜率介于1.0~1.5,说明模型模拟值与实测值吻合程度较高。综合来看,本研究所使用的作物生长模型能够较为准确地模拟半干旱区春小麦干物质与产量形成。
图5 春小麦生长模型模拟的地上部生物量(a)和产量(b)与实测值Fig.5 Growth modle simulation and observation values of aboveground biomass (a) and yield (b) of spring wheat
不同水分处理间,模型模拟获得的春小麦根系层土壤有效含水量以及叶面积指数随生育期变化如图6a-b 所示。供水充足T1 处理和供水相对充足T2处理下,春小麦整个生育期土壤有效含水量都大于0.40,而T3、T4 以及T5 处理下春小麦生育期中段土壤有效含水量分别降至0.21、0.07 以及0.02 附近(图6a),远低于叶片气体交换试验所确定的阈值点0.40 (图3c)。由于水分供给不足,T4 与T5 处理叶面积指数明显小于其他3 个处理(图6b)。尽管T3 处理土壤有效含水量降至0.21 附近,但春小麦叶面积指数与其他两个供水较为充足的处理差异并不明显。分析根系层最小土壤有效含水量与作物地上部生物量以及产量的关系发现,只有T4 和T5 处理的生物量和产量明显偏小,其中尤以T5 处理的生物量和产量最小,远低于其他4 个处理(图6c,6d)。从拟合曲线来看,地上部生物量与产量变动的临界阈值并非叶片气体交换试验所确定的0.40 附近,而是在0.18左右(表5)。
图6 不同水分供给下模拟获得的小麦叶面积指数、根层土壤有效含水量、生物量以及产量的变化特征Fig.6 Variations of simulated leaf area index,available soil water content,biomass and yield for spring wheat under different water treatments
干旱发展过程中,春小麦不同叶片生理指标对干旱响应的阈值不完全相同。其中气孔导度阈值最大,这是因作物维持生长,叶片为了获得最大的碳吸收,适宜条件下气孔会尽可能增大开放程度[31]。而在干旱时,春小麦气孔快速响应环境条件变化,为减少作物体内水分损失,干旱达到一定程度时,部分气孔关闭,气孔导度逐渐下降。在气孔导度开始下降时,作物叶肉内的二氧化碳与大气二氧化碳浓度比往往不会快速下降,这保证了作物进行光合作用的原料供给,因而作物依然能够在气孔少量关闭时正常进行光合作用[32]。只有当干旱导致气孔导度下降到一定程度,作物叶肉中的光合过程才会明显受到影响[33],这一现象导致净光合速率与气孔导度对干旱的响应阈值并不完全相同,这也是本研究中小麦叶片气孔导度与净光合速率对干旱响应阈值存在差异的原因。
作物叶片气体交换对土壤水分的响应阈值,由作物本身的生理特性所决定,其在一定程度上反映了作物对干旱胁迫的敏感性特征[14,29]。本研究小麦叶片气体交换对干旱的响应过程是在桶中观测获取,与一些大田试验观测获得的结果相一致[34]。桶栽试验中小麦根系无法向下延伸,而且小麦生长所需水分无法由深层土壤向上补给,这与大田相比,发生干旱时小麦受旱进程会显著加快[35]。尽管桶栽试验与大田试验在作物生长环境方面存在较大差异,但作物叶片气体交换过程对干旱响应的调控机制(一般由气孔因素和非气孔因素所决定)和阈值响应特征不会发生改变[33,36]。
目前,不同作物模型在模拟干旱对作物生长发育及产量形成过程时,所采用的方法各不相同,但可大致归为4 类: 蒸散比法、作物根区可吸收水分与潜在蒸腾的比值法、临界土壤水分法及土壤有效含水量阈值法等[37-38]。与作物叶片气体交换过程相联系的光合生理过程是决定作物干物质及最终产量形成的基础,因此,用作物叶片气体交换指标对土壤有效含水量的响应阈值作为表征作物受干旱胁迫开始的标志具有可行性。本研究利用观测获得的春小麦叶片净光合速率对土壤有效含水量的响应阈值作为春小麦生长模型中反映干旱对小麦干物质积累的影响指标,能够准确模拟不同年份春小麦地上部生物量和产量形成。该方法直接将小麦受旱响应过程与小麦干物质积累相联系,省去了推算实际蒸散与潜在蒸散比值的过程,原理明确,简单易行,适宜在作物模型构建中应用推广。
尽管春小麦叶片气体交换指标在土壤有效含水量达到阈值后,快速线性降低,但是观察生物量与产量对土壤有效含水量的阈值,发现其明显小于叶片气体交换指标的阈值,这说明作物生长期内叶片尺度上所界定的干旱开始阈值或干旱等级与实际产量减少所定义的阈值会存在较大差异。而这也与实际生产相符,即作物在生长过程中并不是一缺水就会引起作物减产。这也是为何在实际生产中,通过气象要素观测到有气象干旱发生,但是在农业生产中却并没有旱情出现的缘故。一方面,农业生产过程与特定的作物生长时期相联系,在作物生长过程中,只有作物对水分较为敏感的时期缺水才会对作物生长发育造成严重的影响[11,39],比如小麦生长发育的拔节期与开花期。而另一方面,作物产量形成过程受多种因素影响,与干旱程度相比,干旱持续时期的长短是决定干旱影响大小的一个极为重要的因素[12]。在本研究中,水分供给相对较少的处理(T3)小麦根系层土壤有效含水量最低降至了0.20 左右,在叶片尺度上,该数值会导致叶片气体交换指标较水分充足时明显减小。但通过分析春小麦生长模型模拟结果发现,由于土壤有效含水量小于0.40 的时间较短,即使短期缺水导致土壤有效含水量明显下降,干旱对小麦干物质积累的过程影响也不会很大,这也是为何该处理小麦叶面积指数、生物量和产量并未较其他两个水分供给相对充足的处理(T1 与T2)明显减少的原因。相反,当春小麦水分供给不足且干旱持续时间较长时(T4 与T5 处理),其叶面积指数、生物量和产量较水分充足处理明显减小。
在实际生产中,气象部门多以气象要素构建指标反映农业干旱问题,这会出现诸多“天旱地不旱”的情形,从而导致对农业干旱形势的误判。而农学与植物生理学研究中所确定的作物干旱指标(如叶片温度、以冠层气温差值构建的作物水分胁迫指数、气体交换指标如气孔导度和净光合速率,以及一些光谱指数等[10,40])在界定农业干旱程度时也可能会存在一些问题。此类工具依据作物生理指标对水分的响应特征来反映农业干旱,虽然在一定情形下水分是其变动的决定性因子,但这些指标还受其他诸多因素的影响,而且这些指标多属于瞬时值[41-42],忽略了干旱持续时间对作物最终产量的影响,因此用这类指标反映农业干旱,并最终量化作物产量降低程度时也会存在问题。而且作物受旱还是一个累积的过程,瞬时值难以反映干旱累计对作物产量的影响。此外,在作物生长过程中,作物可能不仅是在某一个特定时期受旱,往往会在不同时期遭遇不同程度的干旱胁迫[5,11],这种不同时期不同程度的干旱能否通过叶片或冠层尺度观测获得的瞬时值指标反映,依然是一个亟须论证的问题。而当前众多田间控制试验多以某个特定生育期观测获得的瞬时指标与最终产量建立关系,反映作物受旱问题,忽略了实际生产中干旱发生的随机性。而多个不同时期干旱的叠加及其累计对作物产量的影响,如何用生理指标的瞬时观测值反映依然存在不确定性。
通过控制试验和作物生长模型探讨了小麦叶片气体交换指标与最终产量对干旱的响应阈值差异。发现小麦叶片气体交换指标与最终产量对干旱的响应阈值不完全相同。小麦叶片气体交换指标对干旱反应较快,对干旱胁迫更敏感,而小麦最终产量只有在干旱时间持续达到一定程度才会明显下降。而实际农业生产中,干旱可能在作物生长的任意时期反复发生,因此,与作物叶片气体交换相关生理指标的瞬时观测值,会出现难以量化作物受旱过程,并无法准确反映作物最终产量变化的情况。