Review on the Estimation Model of Crop Evapotranspiration in Solar Greenhouse

ZHAO Ming-yu, YAO Ming-ze, LI Bo, WANG Tie-liang

PDF(678 KB)
Water Saving Irrigation ›› 2023 ›› (4) : 73-81. DOI: 10.12396/jsgg.2022323

Review on the Estimation Model of Crop Evapotranspiration in Solar Greenhouse

Author information +
History +

Abstract

Evapotranspiration model is an effective tool for estimating crop water demand, which provides important parameters for determining irrigation schedule and designing irrigation system, and helps agricultural industry to achieve the goal of water saving, high yield and high quality. There are many existing evapotranspiration models with different application conditions, which are not mature enough to be applied in the greenhouse environment. It is necessary to sort out the models and study the improvement methods of models in the greenhouse environment, so as to improve the accuracy of model estimation. In this paper, the main evapotranspiration estimation models and their applicable conditions were introduced. From the perspective of being applicable to the solar greenhouse environment, the methods of improving the Penman Monteith (PM) model were summarized. On this basis, the research trend of improving the accuracy of evapotranspiration model was pointed out, which enriched the research achievements in the field of crop evapotranspiration estimation in solar greenhouse. The analysis showsthat the PM model has the clearest mechanism and the most extensive application. The key to improve the PM model is to improve the parameters of surface resistance ( rs) and aerodynamic resistance ( ra). In addition, the improvement of crop coefficients ( Kc, Kcb, Ke) can also effectively improve the accuracy of model estimation.

Key words

solar greenhouse / evapotranspiration / simulation / PM model

Cite this article

Download Citations
ZHAO Ming-yu , YAO Ming-ze , LI Bo , WANG Tie-liang. Review on the Estimation Model of Crop Evapotranspiration in Solar Greenhouse. Water Saving Irrigation. 2023, 0(4): 73-81 https://doi.org/10.12396/jsgg.2022323

0 引 言

2022年中央一号文件指出,要加快发展设施农业,因地制宜发展塑料大棚、日光温室、连栋温室等设施[1]。与其他设施相比,日光温室造价低廉,节能环保,是世界上作物低碳生产的典范设施,其面积占据了设施蔬菜总面积的三分之一[2, 3]。我国经济社会用水需求最大的产业是农业,例如,2021年农业用水占总用水量比重为61.56%[4],水资源短缺已经成为制约农业高质量发展的重要因素[5]。习近平总书记提出的新时代治水思路“节水优先、空间均衡、系统治理、两手发力”,明确强调将节水放在首位[6]。对于日光温室而言,节水灌溉既可以节约水资源、提高农业用水利用率,也可以改善温室内高湿环境、减少病虫害、提高作物品质。
作物蒸腾( T)与棵间土壤蒸发( E)合并称为蒸腾蒸发或腾发,消耗的水量称为蒸腾蒸发量( ET)或腾发量( ET),通常又把腾发量称为作物需水量,是农田水分消耗的主要途径。准确估算作物腾发量,从而制定合理的灌溉制度,可以提高用水效率,有助于实现农业高效节水,促进农业产业健康可持续发展。
腾发量估算模型种类多、适用性各有不同,在大田环境应用比较成熟。日光温室与大田环境的光照、温度、湿度、风速等条件差异较大,不能将大田环境下使用的估算模型直接运用于温室环境。有必要归纳常见的作物腾发量估算模型,分析对比模型的适用性和优缺点,筛选出适用于日光温室环境的估算模型,梳理提高模型估算精度的改进方法,并提出下一阶段研究趋势,研究对日光温室作物腾发量估算具有重要意义。

1 腾发量估算模型

按照计算程序划分,腾发量估算模型包含直接估算模型和间接估算模型两种。直接估算模型又分为经验型和机理型,其中经验型通过大量数据拟合建立作物蒸腾蒸发和温室内外气象条件、作物生理生态指标之间的关系,所需参数少、计算简便,不足之处是精度偏低、应用受限于作物类型和地域条件[7];机理型立足于水量平衡、水汽扩散等理论,典型代表是稠密植被条件下的PM单层结构模型和稀疏植被条件下的shuttleworth--Wallace(SW)双层结构模型,以及无平流假设条件下提出的Priestley-Taylor(P-T)模型。间接模型首先要确定参考作物腾发量( ET0),再根据实际作物类型、土壤反射率等条件计算出实际作物腾发量,最常见的是联合国粮食及农业组织(FAO)推荐的单、双作物系数法。

1.1 PM模型

1948年Penman利用空气动力学和能量平衡原理建立了基于气象数据计算水面蒸发的理论公式,1965年Monteith在此基础上,考虑冠层的阻力得到了PM模型[8]
λET=ΔRn-G+ρacpVPD/raΔ+γ 1+rs/ra
式中: λET为潜热通量,W/m2 Δ为饱和水汽压随温度变化曲线的斜率,kPa/℃; Rn为净辐射,W/m2 G为土壤热通量,W/m2 ρa为空气密度,kg/m3 cp为空气比热,kPa/℃; VPD为饱和水汽压差,kPa; γ为湿度计常数,kPa/℃; rs为表面阻力,s/m; ra为空气动力学阻力,s/m。

1.2 SW模型

PM模型将植被和土壤假设为大叶,无法区分土壤蒸发和作物蒸腾。1985年,Shuttleworth和Wallace假设作物冠层为均匀覆盖,在PM模型基础上,考虑将太阳辐射按照Beer定律分配到冠层及土壤表面,发展出能够在全生育期不同地表覆盖条件下正确划分土壤蒸发和作物蒸腾的SW模型[9]
λET=λE+λT=CsPMs+CcPMc
PMs=Δ A+ρacpVPD-Δ rasA-As/raa+rasΔ+γ 1+rss/raa+ras
PMc=Δ A+ρacpVPD-Δ racA/raa+racΔ+γ1+rsc/raa+rac
Cs=1+RsRa/RcRs+Ra-1
Cc=1+RcRa/RsRc+Ra-1
Ra=Δ+γraa
Rs=Δ+γras+γrss
Rc=Δ+γrac+γrsc
A=Rn-G
As=Rns-G
Rns=Rnexp-C LAI
式中: λE为土壤蒸发潜热通量,W/m2 λT为冠层蒸腾潜热通量,W/m2 raa ras分别为冠层到参考面和参考面到土壤表面的空气动力学阻力,s/m; rsc rss分别为冠层阻力和土壤表面阻力,s/m; rac为作物冠层的边界层阻力,s/m; A As分别为总能量和到达下垫面的可利用能量; Rns为到达土壤表面的净辐射,W/m2,采用Beer定律计算; C为消光系数;LAI为叶面积指数。

1.3 P-T模型

1972年,Priestley和Taylor假设无平流影响,将PM模型中包含的空气动力学项整合为经验系数 α,形成P-T模型,根据大区域饱和陆面和海洋的气象观测数据,推荐无平流影响的情况下 α的取值为1.26[10, 11]
ET=α Δλ (Δ+γ) Rn-G
式中: ET表示为估算腾发量,mm/d; λ为水汽化潜热,MJ/kg; α为经验系数; Δ γRnG意义同上。

1.4 作物系数法

为了应用方便,FAO-56对PM模型进行了标准化,设置了参考作物腾发量计算方程式(14),将难以获取的表面阻力和空气动力学阻力融入其中,形成了仅需作物特性和气象条件的单作物系数法[式(15)]和双作物系数法[式(16)]。单作物系数法将作物蒸腾和土壤蒸发融合计算,双作物系数法可以分别计算出二者数值。
ET0=0.408 Δ Rn-G+γ 900 u2T+273 VPDΔ+γ 1+0.34 u2
式中: ET0为参考作物腾发量,mm/d; Rn为作物表面净辐射,MJ/(m2d); G为土壤热通量,MJ/(m2d); T为2 m高处的平均气温,℃; u2为2 m高处的风速,m/s; Δ γ VPD意义同上。
ET=Kc ET0
ET=(Ks Kcb+Ke) ET0
式中: Kc为综合作物系数,与作物种类、品种、生育期和作物群体叶面积指数有关,是作物自身生物学特性的反映; Kcb为基础作物系数,被定义为土壤表层土壤干燥而根区平均含水量不构成土壤水分胁迫条件下 ET ET0的比值,还包括由干燥地表以下土壤水分和稠密植被下的土壤水分补给的残余扩散蒸发部分。 Ke为表层土壤蒸发系数,它代表了作物地表覆盖较小的幼苗期和前期生长阶段中,除 Kcb中包含的残余土壤蒸发效果外,在降雨或灌溉发生后由大气蒸发力引起的表层湿润土壤的蒸发损失比。 Ks为土壤水分胁迫系数,主要和田间土壤有效水分有关,本文介绍土壤供水充足时的状况,即 Ks=1
计算 ET0的方法还有很多,比如FAO 24提供的penman方程、辐射模型、蒸发皿蒸发量计算模型,FAO-56 提供的利用蒸发皿数据计算 ET0,hargreaves方程,PrHo模型,Makkink模型[12],Schendel和Hargreaves-M4经验模型[13]等,但FAO-56唯一推荐的方法是式(14) [9]。该模型可适应不同的气候区域,常作为标准公式验证其他公式的适用性[14]

1.5 小 结

P-T模型涉及参数少,计算方便,但是其无平流的条件假设在现实环境尤其是温室环境下难以达到,下垫面的粗糙度不同或者土壤湿度不同均会导致平流的出现[15],需要根据气象数据重新确定经验系数 α。SW模型估算精度较高[16],可以分别计算土壤蒸发和作物蒸腾,但是其所需参数多,计算过程复杂,需要借助于计算机完成。作物系数法可以拆分腾发组分,为研究降低土壤蒸发提高水分利用效率提供数据参考,其模拟精度关键在于参考作物腾发量的计算和作物系数的选取。
PM模型机理明确,稳定性高,是上述3种计算方法的基准,在作物腾发量估算上得到了广泛应用[17-20]。但是,PM模型最早是在大田环境中提出的,而温室不通风时风速极低、太阳辐射较小、高温高湿等条件与大田差异很大,模型在温室独特的小气候环境下的适用性及其改进研究是近阶段专家学者的关注重点。

2 PM模型改进研究进展

PM模型主要应用于两个方面,一种是标准化后的PM模型用来计算 ET0,另一种是将表面阻力( rs)和空气动力学阻力( ra)代入PM模型后直接计算得到 ET,模型在应用过程中得以不断深化研究和改进。

2.1 在计算 ET0方面的改进

2.1.1 风速为零条件下修正 ra计算方法

王建、陈新明[21, 22]等通过测量,认为温室内风速为零,将 u2=0代入式(14)后,式(14)的空气动力学项为0,仅剩辐射项。但是,温室内蒸发和热量输送一直存在,空气动力学项不应该为0,所以,大田环境下适用的式(14)不能直接应用于温室当中。
造成式(14)不适用于温室环境的原因在于,其空气动力学阻力 ra计算等式受中性稳定条件约束,与温室内空气边界层实际情况不符。当认为温室内风速为零时,按照式(17)计算的 ra将出现无穷大的情况。
ra=lnx-d/hc-d lnx-d/z0k2 u
式中: k为卡曼常数,取为0.4; x为参照高度,m;即风速与温湿度的测量高度; d为零平面位移,m; u为参考高度处的水平风速,s/m; z0为动量传输粗糙度长度,m; hc为平均冠层高度,m。
为此,王建、陈新明等将Thom和Oliver计算空气动力学阻力公式(18)引入到温室环境中,得到式(19)
ra=4.72 lnZ-dZ02/1+0.54 u
ra=4.72 lnZ-dZ02
式中: Z为风速测量高度; Z0为地面粗糙度; d为零平面位移长度, ra u意义同上。
对于参考作物,取作物高度 hc=0.12 m, Z0=0.13 hc d=0.64 hc代入式(19)得出 ra=109.40 s/m,进而得出适用于温室 ET0计算的修正公式:
ET0=0.408 Δ Rn-G+γ 1 713T+273 VPDΔ+1.64 γ
闫浩芳[23]、刘浩[24]、毋海梅[25]、赵宝山[15],王林林[26]等均曾采用式(20)开展温室作物腾发量相关研究,并取得了较好的效果。王林林在式(20)基础上,引入作物高度参数,推导出适于温室茶树的 ET0计算方法。目前,在没有强制通风的温室环境中,式(20)被广泛用于计算参考作物腾发量。

2.1.2 风速为零条件下固定 ra取值

为了评估不同 ET0计算方法在温室环境下的适宜性,Fernández[27]在塑料温室内种植多年生牧草,定期修剪使其高度保持在0.10~0.15 m,使其基本满足参考作物的假想条件。通过实测作物的蒸散量,评估了5种计算 ET0的主要方法,结果表明式(14)的估算结果比其他4种准确度稍高一些。
但是,式(14) ra采用的是风速函数[式(17)],在温室内风速极低的条件下, ra数值介于2 080~20 800之间波动,这与实际获得的 ra情况不符。文中指出,对于保持在0.10~0.15 m高度的温室多年生草作物,整个季节的 ra都是恒定的,通过假设季节内 ra值分别为100、150、200和250 s/m,并再次使用PM模型估算温室 ET0,发现这些 ra明显改善了 ET0估值的拟合度,最佳值为295 s/m。为此,得出 ra取为295 s/m时温室 ET0计算的修正公式[15, 28]
ET0=0.408 ΔRn-G+γ 628T+273 VPDΔ+1.24 γ
Fernández的研究对象为地中海地区的塑料大棚,其研究结果不能直接推广应用于其他地区。但是其研究思路非常值得借鉴:通过实测数据验证不同 ET0计算方法的精度,针对精度最高的计算方法进一步优化,形成改进算法或简便算法。

2.2 在计算 ET方面的改进

2.2.1 改进阻抗( ra rs)算法

目前尚无标准的 ET0计算方法,直接计算 ET具有一定意义。由公式(1)可以知道,通过代入 ra rs可以直接利用PM模型计算出任何实际作物的腾发量。选择合理的 ra rs的计算方法是修正PM模型,提高模拟精度的关键。
(1)改进 ra。目前,计算 ra的主要方法有3种。一种是FAO-56中推荐的式(17),一种是美国国家灌溉排水工程手册中采用的式(18),还有一种是通过热传输系数计算空气动力学阻力的方法[29][式(22)]。
ra=ρa cp2 h LAI
式中: h为热传输系数,W/(m2·k);不同对流类型采用的计算方法不同; LAI为叶面积指数,m2/m2 ρa cp意义同上。
Qiu等研究认为,中国西北地区日光温室内空气对流类型主要为混合对流[30];闫浩芳等研究表明,中国南方地区Venlo型温室内空气对流类型主要为混合对流[31]。Qiu、闫浩芳、龚雪文等采用热传输系数计算公式(23)计算得到混合对流类型下的 ra,将其代入PM模型并应用于温室作物 ET的估算,取得了较好的模拟结果。
h=0.37 kcdc Gr+0.692 Re21/4
式中: kc为空气的导热系数,W/(m2·k); dc为叶片的特征长度,m; Gr为高尔夫数; Re为雷诺数。
3种空气动力学阻力的算法中,式(17)起源并广泛应用于大田环境,在温室低风速下估算精度偏低;式(18)式(22)充分考虑了温室内风速接近于零的情况,可以用于构建改进的PM模型,使模型更适用于温室环境。
在实际应用中,温室内的风速和冠层温度等不易获取,而 ra在一定范围内取定值,不会显著影响PM模型的模拟精度。因此,一些的学者通过设置固定的 ra取值[9, 15, 32, 33],并根据其应用于PM模型时的表现,指出在精度要求不高的情况下, ra取某些定值可以简化计算,提高PM模型的实用性。
(2)改进 rs。PM方程最初用于根据气象数据计算水面蒸发量,通过引进阻力参数后才适用于作物表面的蒸发,阻抗参数的准确性直接影响着PM模型估算精度。表面阻力( rs)不是一个可以通过测量获得的真实物理量,包含冠层阻力( rsc)和土壤表面阻力( rss)两部分[34]。以往的研究往往忽略土壤表面阻力,以冠层阻力代替表面阻力。但是,在作物生长初期等植被稀疏情况下,这样的替代会造成 ET估算准确度降低。因此,建立综合土壤表面阻力和冠层阻力的表面阻力模型,并代入于PM方程是提高 ET估算精度的有效手段。
FAO-56给出了冠层阻力的近似公式:
rsc=r1LAIactive
式中: r1为光照好的叶面总气孔阻力,s/m,是单叶气孔阻力的平均值; LAIactive为有效叶面指数[m2(叶面积)/m2(土壤面积)]。
龚雪文根据并联电路原理,引入土壤表面阻力模型,构建了一个表面阻力模型 rs[式(25)],应用结果表明在灌水期间及叶面积指数较低的情况下,表面阻力模型优于冠层阻力模型[9]
rs=1a LAI/rsc+b/rss
式中: a b是经验系数; rsc rss LAI意义同上。
单叶气孔阻力不易测得,其大小决定于叶片气孔的开闭程度。在充分供水时,叶片气孔的开闭程度与气象因素直接相关,为此,赵宝山、闫浩芳等通过分析气孔阻力与温室内气象因子的响应关系来确定其值,研究结果表明,气孔阻力与太阳辐射呈指数相关[15,31]
除了式(24)以外,还可以通过冠层阻力模型计算 rsc。比较经典的冠层阻力模型有Jarvis(JA)[35]、Katerji-Perrier(KP)[36]、Todorovic(TD)[37]、stand(ST)[38]等,被广泛应用于大田环境。一些学者将其中的一个或几个模型引入到PM方程中进行 ET估算,验证了模型的适用性,一定程度上提高了PM方程的估算精度[15, 18, 34, 39-42]
对冠层阻力模型进行改进,同样能够提高PM模型估算精度。Li等引入土壤表面阻力对JA模型进行改进,形成了CO模型,通过与传统JA、KP、ST、TD等模型模拟结果比较发现,CO模型显著提高了模拟精度,并增强了PM模型的可靠性[43]。氮素能调节气孔导度,进而改变作物蒸腾量,基于此,ZHOU等在温室盆栽番茄中开展实验,建立了考虑氮改良的N-JA模型,结果表明,N-JA模型能更准确的预测番茄植株蒸腾[44]
(3)小结。改进阻抗算法( ra rs)是提高PM模型模拟精度的关键,学者们就此开展了大量研究,取得了一定的研究成果。在 ra的选用上基本达成一致,即大田环境采用式(17),温室环境下计算 ET0 ET分别选用式(18)式(22)。在植被稠密状态下,式(24)具有良好的应用效果,但是植被稀疏情况下怎样合理计算土壤蒸发尚未达成共识,具有较大的研究空间。经典冠层阻力计算模型在国内外大田环境下适用性较强,但是在温室环境下表现如何尚待进一步验证。为此,基于经典冠层阻力模型,考虑土壤蒸发阻力的影响,改进PM模型进行日光温室内作物腾发量的模拟将是未来的研究重要和热点。

2.2.2 改进作物系数法

作物系数法因其所需参数较少,模拟结果比较可靠被FAO-56推荐用于估算作物腾发量,在实际使用中也取得了较好的应用效果[45-47]。尤其是双作物系数法,可以有效区分土壤蒸发( E)和作物蒸腾( T),在稀疏冠层和生育初期条件下,模型更贴近实际,模拟效果更好[48, 49]。作物系数受种植区域、气象条件、作物类型的影响,FAO-56的推荐值不能完全适用,需要根据实际情况进行修正[50]。除了前文2.1中阐述的对 ET0的改进以外,还可以通过改进作物系数( Kc Kcb Ke)来提高作物系数法的模拟精度。
(1)改进 Kc Kc综合考虑土壤蒸发和作物蒸腾,而土壤蒸发受灌溉或降雨影响波动较大,一般采用多日平均数值,因此单作物系数法适合时间步长较大的运算[51] Kc将实际作物与参照作物在作物高度、冠层阻力、土壤蒸发、表面反射等方面的不同融合在一起进行体现,机理复杂,难以直接从某一方面建立改进算法。
为了简化计算,一些学者通过探索单一因素与 Kc之间的响应规律来建立作物需水量模型。刘浩建立了滴灌条件下温室番茄 Kc LAI、积温( Ta)的关系曲线分别为式(26)式(27) [24]
Kc=0.153 3 LAI+0.549 7
Kc=-1.025 1×10-6 Ta2+1.788 5×10-3 Ta+0.264 9
据此对温室番茄需水量估算模型进行大幅度简化,形成只包含辐射、相对湿度和气温等常规气象数据的表达式(28),验证结果表明简化模型具有较高的模拟精度,可以比较真实的反映温室番茄的需水量。
ET=Kc ET0=(-1.025 1×10-6 Ta2+1.788 5×10-3 Ta+0.264 9)0.408 Δ Rn-G+γ 1 694T+273 VPDΔ+1.64 γ
式(28) ET0为刘浩采用陈新明[22]等人研究思路,对PM方程中与风速有关的空气动力学项进行修正,推导得出的表达式。
李毅等构建了冬小麦不同生长阶段单作物系数随砂石覆盖量的变化关系式(29),并得到适用于砂石覆盖条件下作物需水量模型式(30),可以实现在已知日气象数据和砂石覆盖量的情况下初估冬小麦不同生长阶段的 ET [52]
Kcs,i=ai MsGA+bi
ET=(ai MsGA+bi) ET0
式中: Kcs,i为第 i个生长阶段的单作物系数; ai bi为第 i个生长阶段的拟合参数; MsGA为砂石覆盖量。
(2)改进 Kcb。FAO-56给出了根据推荐的基础作物系数( Kcb推荐)以及气象条件和作物因素修正 Kcb的方法:
Kcb=Kcb推荐+0.04 u2-2-0.004 (RHmin-45) h30.3
式中: RHmin是作物生长中期或后期的日最小相对湿度,%, 20%RHmin80% h是作物生长中期或后期的平均株高,m; u2意义同上。
FAO-56推荐的基础作物系数为静态值,难以体现作物生长过程的动态变化,因此丁日升[53]、冯禹[54]、龚雪文[55]、zhao[56]、毋海梅[25]、黄松[57]等均引入冠层覆盖度系数 Kcc计算动态 Kcb,取得了良好的应用效果。
Kcb=Kc,min+Kcc(Kcb,full-Kc,min)
Kcb,full=min(1.0+0.1 h, Kmax)+0.04 (u2-2)-0.004 (RHmin-45) h30.3
Kcc=1-exp(-k LAI)
式中: Kc,min为裸土最小作物系数,取值为0.1; Kcb,full为作物完全覆盖地表时的最大基础作物系数; Kmax为作物系数最大值,取1.2; k为太阳辐射的冠层衰减系数,取0.7; h u2 RHmin LAI意义同上。
龚雪文在此基础上,同时考虑番茄叶片的衰老影响,增设叶片衰老因子 fs修正番茄中、后期的 Kcb
Kcb=(1-fs) (Kcb,min+Kcc (Kcb,full-Kcb,min))
式中:叶片衰老因子 fs在幼苗期、花果期和盛果期取值为0,采摘期取值为0.2。
(3)改进 Ke。FAO-56给出了土壤蒸发系数推荐算法如式(36)所示:
Ke=Kr Kcmax-Kcbfew Kcmax
式中: Kcmax为降雨或灌溉后 Kc的最大值; few为最大的土壤蒸发表面所占的百分比; Kr为取决于表层土壤蒸发(或水分消耗)累积深度的蒸发减小系数,无量纲。
few取裸露土壤( 1-fc)和湿润土壤( fw)表面所占比值的较小者, fc为植被覆盖的土壤面积比。为了准确地估算作物冠层下土面蒸发的动态变化,康绍忠[23, 25, 26, 54, 56, 58, 59]等学者应用 LAI计算 fc[式(37)],实现了对土壤蒸发系数的优化。
fc=1.005 1-exp(-0.6 LAI)1.2
康绍忠[56]对蒸发减小系数做了修正,用实际土壤体积含水率代替了FAO-56推荐公式中的田间持水率,修正方法被闫浩芳[23]等采用,取得了较好的模拟结果。
在覆膜种植情况下,李丰琇[60]和马嘉敏[59]以覆膜面积作为考虑因素,修正 ke的表达式来体现覆膜对土壤蒸发的影响。
(4)小结。目前,对作物系数的改进主要考虑的是不同生长阶段 LAI的变化,通过修正和引入某些系数体现 LAI动态变化对蒸腾蒸发的影响。另外,还有一些学者利用机理模型和实测数据,计算出作物系数实际值,建立作物系数与环境因子的相关关系式,作为适应当地条件的经验模型使用[61]

3 结论与展望

作物腾发量估算是确定作物需水量,制定灌溉制度的科学依据,其研究模型从大田环境提出,经过改进后基本适用于温室条件。PM模型及其衍生模型受到了广泛应用,为了提高作物腾发量估算精度,众多学者围绕阻抗、作物系数进行改进研究,取得了较好的效果,表明阻抗及作物系数是影响模型估算精度的关键因素。基于前人的研究成果,未来可以从提高精度和方便使用两个角度继续优化模型。

3.1 从模型机理出发提高估算精度

准确估算 ET0是提高作物腾发量估算的基础,目前普遍采用的模型[式(20)]认为温室内风速为零与实际情况不符,不通风条件下温室内风速较低但并不为零,而自然或强制通风条件下温室对流情况更是复杂多变,这种环境下作物腾发量的估算需要深入探索[62]。随着监测设备精度的提高,温室内空气流动规律及风速逐渐被掌握,据此研究修改空气动力学阻力 ra的算法,可使得PM模型更贴合温室实际环境。Fernández通过种植多年生牧草,获得参照作物腾发量实测值,并据此验证模型精度的做法非常合理,但是所取得的结果[式(21)]不是放之四海而皆准,应考虑在特定区域和环境下如北方温度较低地区日光温室内按照这一思路进行 ET0模型的重新构建。
rs是影响模型估算精度的重要参数之一,在覆膜种植条件下往往忽略土壤蒸发,并以冠层阻力代替表面阻力。但是,覆膜种植因其会对环境造成污染而逐渐被常规种植模式所替代,此时,必须考虑土壤表面阻力,但是相关研究成果并不成熟,需要探索更合理的 rs算法。

3.2 从响应规律出发简化模型计算

前文所述一些学者通过建立气孔阻力和太阳辐射的关系,砂石覆盖度与单作物系数的关系,简化了模型的计算。提示后人在实测数据缺失或者精度要求不高的区域,可以通过研究某些因素之间响应关系,寻找规律并建立关系式,最终形成方便使用的经验模型。龚雪文通过室内外气象关联法和能量平衡法建立了基于室外气象资料估算室内番茄蒸腾量的经验模型,验证了这一思路的可行性[9]

3.3 借助信息技术提高使用便捷度

在信息化手段极速发展的背景下,以PM模型为理论基础研发设计作物腾发量估算软件,是将理论知识应用于实践的重要手段。CROPWAT被FAO推荐用于确定灌溉制度;SIMDualKc被用于模拟土壤蒸发和作物蒸腾[63, 64],邱让建、杜太生用SIMDualKc来估算温室番茄耗水量也取得较好的效果[65];卢晓鹏通过计算机语言,将单、双作物系数法计算玉米需水量进行编程,极大的简化了实验和人力计算成本[66]
此外,随着测定技术的提升和研究的深入,探索水分胁迫、盐分胁迫、病虫害、杂草等情况下作物蒸腾蒸发规律,分析复合因素作用下温室环境、作物生理生态指标的响应机制也是未来的研究趋势之一。

References

1
中共中央国务院关于做好2022年全面推进乡村振兴重点工作的意见[Z].2022-02.
2
闫冬梅,胡林,周长吉.椭圆管单管拱架日光温室结构性能分析[J].农业工程学报,2022(5):217-224.
YAN D M, HU L, ZHOU C J, et al. Analysis of the performance of the oval tube arch of single-tube solar greenhouses[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022,38(5): 217-224. doi:10.11975/j.issn.1002-6819. 2022.05.026 .
3
孙周平,富宏丹,余朝阁.日光温室番茄优质轻简高效生产技术规程[J].中国蔬菜,2022(5):111-119.
4
中华人民共和国水利部.中国水资源公报[M]. 北京: 中国水利水电出版社, 2021.
5
房建恩,马立伟.农业节水技术推广困境及对策研究[J].农业经济, 2022(5):9-11.
FANG J E, MA L W. Study on the difficulties and countermeasures of agricultural water-saving technology popularization[J]. Agricultural Economy, 2022(5):9-11.
6
任亮,董小涛,王崴.供给侧改革视阈下推广合同节水管理的思考[J].人民黄河,2022(44):33-35.
REN L, DONG X T, WANG W. Thoughts on popularizing contractual water saving management from the Perspective of Supply side Reform[J]. Yellow River, 2022(44):33-35.
7
邱让建,杜太生,刘春伟.温室作物蒸散发估算模型研究进展[J].灌溉排水学报,2015(z2):134-139.
QIU R J, DU T S, LIU C W. Research progress of greenhouse crop evapotranspiration estimation model[J]. Joural of Irrigation and Drainage,2015(z2):134-139.
8
ALLEN R G, PEREIRA L S, RAES D, et al. Crop evapotranspiration-Guidelines for computing crop water requirements[Z]. FAO Irrigation and drainage paper 56. Fao, Rome 300, 1998, D05109.
9
龚雪文.温室滴灌条件下土壤_作物系统水热传输与模拟[D].北京:中国农业科学院,2017.
GONG X W. Transmission and simulation of water and heat in soil-plant system under drip irrigation in solar greenhouse[D]. Beijing: Chinese Academy of Agricultural Sciences Dissertation,2017.
10
姚欢欢,景元书,韩丽娟.基于 Priestley-Taylor 模型的低丘红壤区水热通量变化及影响因素[J].中国农村水利水电,2020(9):1-7.
YAO H H, JING Y S, HAN L J. Research on the change of water and heat flux and influencing factors in Low Hilly Red Soil Region based on Priestley-Taylor model[J]. Rural water conservancy and hydropower in China, 2020(9):1-7.
11
赵若男,肖薇,石立新.内陆小型水体水汽传输系数和 Priestley-Taylor 系数的量化[J/OL].地球科学.
ZHAO R N, XIAO W, SHI L X. Quantification of the water vapor transport coefficient and Priestley-Taylor coefficient over small inland water bodies [J/OL].Earth science.
12
Linvolak Pich, 蔡焕杰, Qaisar Saddique, 等.关中地区夏玉米和冬小麦不同蒸发蒸腾量估算方法的研究[J].节水灌溉,2019(8):39-47.
LINVOLAK P, CAI H J, QAISAR S,et al. Study on different methods for estimating evapotranspiration of summer maize and winter wheat in Guanzhong Area[J].Water Saving Irrigation, 2019(8):39-47.
13
肖璐,崔宁博,赵璐,等. 西北地区夏玉米不同生育期蒸发蒸腾量模拟模型适用性评价[J]. 灌溉排水学报201938(S2):20-29.
XIAO L, CUI N B, ZHAO L, et al. Applicability evaluation of evapotranspiration estimate models for different growth stages of summer maize in Northwest China[J]. Journal of Irrigation and Drainage201938(S2):20-29.
14
陈宇翔,郑荣伟,王庆明,等.关中地区参考作物蒸发蒸腾量时空变化特征及主要因子分析[J].灌溉排水学报,2022,41(4):76-83.
CHEN Y X, ZHENG R W, WANG Q M,et al. Spatiotemporal variation in reference evapotranspiration over Guanzhong Region and its determinants[J].Journal of Irrigation and Drainage,2022,41(4):76-83.
15
赵宝山.大田及温室条件下作物蒸发蒸腾模型及参数的研究[D].江苏镇江:江苏大学,2019.
ZHAO B S. Parametrization of resistances for modeling crop evapotranspiration in the open field and greenhouse environment[D]. Jiangsu Zhenjiang: Jiangsu University,2019.
16
郭燕,刘春伟,雒福佐,等.祁连山高寒草甸蒸发蒸腾量估算及其影响因素分析[J].节水灌溉,2022(5):39-46.
GUO Y, LIU C W, LUO F Z,et al. Estimation of evapotranspiration of alpine meadow in Qilian Mountain and its influencing factors analyzation[J].Water Saving Irrigation,2022(5):39-46.
17
韩沙,马德新 .Penman-Monteith 方程在茶园节水灌溉中的应用研究[J].湖北农业科学201958(16):121-124,138.
HAN S, MA D. Application of Penman-Monteith equation in water saving irrigation in tea garden[J]. Hubei Agricultural Sciences201958(16):121-124,138.
18
刘春伟,吴心语,邱让建.采用作物系数法和PM模型估算南京地区玉米田蒸发蒸腾量[J].节水灌溉,2016(9):12-17.
LIU C W, WU X Y, QIU R J. Evaluation of evapotranspiration of maize with crop coefficient and Penman-Monteith in Nanjing[J].Water Saving Irrigation, 2016(9):12-17.
19
周始威,胡笑涛,王文娥, 等.春玉米生长中期灌水量差异对其蒸发蒸腾量估算精度的影响[J].排灌机械工程学报,2018,36( 9) : 790-795.
ZHOU S W, HU X T, WANG W E, et al. Effect of irrigation quota on estimation accuracy of evapotranspiration in growth middle stage of spring maize[J].Journal of Drainage and Irrigation Machinery Engineering,2018,36(9): 790-795.
20
董斌,孙宁宁,罗金耀.基于棚内气象数据的冬季大棚番茄蒸腾计算[J].武汉大学学报(工学版),2009(5): 601-604.
DONG B, SUN N N, LUO J Y. Calculation of greenhouse tomato transpiration based on indoor meteorological data in winter[J]. Engineering Journal of Wuhan University,2009(5):601-604.
21
王健,蔡焕杰,李红星,等.日光温室作物蒸发蒸腾量的计算方法研究及其评价[J].灌溉排水学报,2006(6):11- 14.
WANG J, CAI H J, LI H X,et al. Study and evaluation of the calculation methods of reference crop evapotranspi ration in solar-heated greenhouse[J].Journal of Irrigation and Drainage,2006(6):11-14.
22
陈新明,蔡焕杰,李红星,等.温室大棚内作物蒸腾蒸发量计算[J].应用生态学报,2007(2):317-322
CHEN X M, CAI H J, LI H X,et al. Calculation of crop evapotranspiration in greenhouse[J]. Journal of Applied Ecology,2007(2):317-322.
23
闫浩芳,毋海梅,张川,等.基于修正双作物系数模型估算温室黄瓜不同季节腾发量[J].农业工程学报,2018,34(15):117-125.
YAN H F, WU H M, ZHANG C,et al. Estimation of greenhouse cucumber evapotranspiration in different
seasons based on modified dual crop coefficient model[J].Transactions of the Chinese Society of Agricultural Engineering,2018,34(15):117-125.
24
刘浩.温室番茄需水规律与优质高效灌溉指标研究[D]. 北京:中国农业科学院,2010.
LIU H. Water requirement and optimal irrigation index for effective water use and high-quality of tomato in greenhouse[D]. Beijing: Chinese Academy of Agricultural Sciences,2010.
25
毋海梅.温室黄瓜耗水模型及高效灌水指标的研究[D].江苏镇江:江苏大学,2019.
WU H M. Study on cucumber water consumption model and efficient irrigation index of greenhouse conditions[D]. Jiangsu Zhenjiang: Jiangsu University,2019.
26
王林林,马文杰,马德新.基于Penman-Monteith方程的温室茶树蒸腾蒸发模型研究[J].节水灌溉,2017(8):30-34.
WANG L L, MA W J, MA D X. Research on evapotranspiration model of tea plants in greenhouse based on Penman-Monteith equation[J]. Water Saving Irrigation,2017(8):30-34.
27
FERNÁNDEZ. Measurement and estimation of plastic greenhouse reference evapotranspiration in a Mediterranean climate [J]. Irrigation Science, 2010, 28:497-509.
28
GONG X W, QIU R J, SUN J S, et al. Evapotranspiration and crop coefficient of tomato grown in a solar greenhouse under full and deficit irrigation[J]. Agricultural Water Management, 2020, 235: 106154.
29
BAILEY B J, MONTERO J I, BIEL C, et al. Transpiration of Ficus benjamina: comparison of measurements with predictions of the Penman-Monteith model and a simplified version[J]. Agricultural & Forest Meteorology, 1993, 65(3/4): 229-243.
30
QIU R, KANG S, DU T, et al. Effect of convection on the Penman–Monteith model estimates of transpiration of hot pepper grown in solar greenhouse[J]. Scientia Horticulturae, 2013, 160(3): 163-171.
31
闫浩芳,赵宝山,张川,等.Penman-Monteith模型模拟Venlo 型温室黄瓜植株蒸腾[J]. 农业工程学报,2019,35(8):149-157.
YAN H F, ZHAI B S, ZHANG C, et al. Estimating cucumber plants transpiration by Penman-Monteith model in Venlo-type greenhouse[J]. Journal of Agricultural Engineering,2019,35(8):149-157.
32
徐立鸿,肖康俊,蔚瑞华. 基于温室环境和作物生长的番茄基质栽培灌溉模型[J]. 农业工程学报202036(10):189-196.
XU L H, XIAO K J, WEI R H. Irrigation models for the tomatoes cultivated in organic substrate based on greenhouse environment and crop growth[J].Transactions of the Chinese Society of Agricultural Engineering202036(10):189-196.
33
VILLARREAL-GUERRERO F, KACIRA M, FITZ-RODRĺGUE E, et al. Comparison of three evapotranspiration models for a greenhouse cooling strategy with natural ventilation and variable high pressure fogging[J].Scientia Horticulturae,2012,134(2):210-221.
34
林馨贝, 周岗, 郑泽涛, 等. 不同冠层阻力模型在夏玉米蒸散发计算中的优化应用[J]. 灌溉排水学报, 2021, 40(6): 28-35.
LIN X B, ZHOU G, ZHENG Z T, et al. Optimizing the canopy resistance models to calculate evapotranspiration from summer maize fields[J]. Journal of Irrigation and Drainage, 2021, 40(6): 28-35.
35
JARVIS P. The interpretation of the variations in leaf water potential and stomatal conductance found in canopies in the field[J]. Philosophical Transactions of the Royal Society of London, B, Biological Sciences, 1976, 273(927): 593-610.
36
KATERJI N, PERRIER A, RENARD D, et al. Modélisation de l'évapotranspiration réelle ETR d'une parcelle de luzerne: rôle d'un coefficient cultural[J]. Agronomie, 1983, 3(6): 513-521.
37
TODOROVIC M. Single-layer evapotranspiration model with variable canopy resistance[J]. Journal of Irrigation and Drainage Engineering, 1999, 125(5): 235-245.
38
STANNARD D I. Comparison of penman-monteith, shuttleworth-wallace, and modified priestley-taylor evapotranspiration models for wildland vegetation in semiarid rangeland[J]. Water Resources Research, 1993, 29(5):1 379-1 392.
39
闫浩芳,周裕栋,张建云,等. 苏南地区典型作物冠层阻力参数及潜热通量的模拟[J]. 农业工程学报202238(9): 101-107.
YAN H F, ZHOU Y D, ZHANG J Y, et al. Parametrization of canopy resistance and simulation of latent heat fluxes for typical crops in southern Jiangsu Province[J]. Transactions of the Chinese Society of Agricultural Engineering, 2022, 38(9): 101-107.
40
郭仲英,吴英楠,李巧珍,等. 六种冠层阻力模型在冬小麦蒸散估算中的应用[J].农业工程学报202238(12):109-117.
GUO Z Y, WU Y N, LI Q Z, et al. Application of six canopy resistance models for estimating winter wheat evapotranspiration[J]. Transactions of the Chinese Society of Agricultural Engineering, 2022, 38(12): 109-117.
41
文建川,景元书,韩丽娟.基于 Penman-Monteith 模型的低丘红壤区稻田蒸散模拟[J].中国农业气象,2020,41(4):201-210.
WEN J C, JING Y S, HAN L J. Simulation of evapotranspiration for paddy rice in Low Hilly Red Soil Region Base[J].Chinese Journal of Agrometeorology,2020,41(4):201-210.
42
李玲.不同冠层阻力模型在半干旱区玉米生态系统蒸散发的应用研究[D]. 兰州:兰州大学,2019.
LI L. Application of different canopy resistance models in evapotranspiration of maize ecosystem in semi-arid areas[D]. Lanzhou: Lanzhou University,2019.
43
LI S, ZHANG L, KANG S, et al. Comparison of several surface resistance models for estimating crop evapotranspiration over the entire growing season in arid regions[J]. Agricultural and Forest Meteorology, 2015, 208:1-15.
44
ZHOU H P, KANG S Z, TONG L, et al. Improved application of the Penman-Monteith model using an enhanced Jarvis model that considers the effects of nitrogen fertilization on canopy resistance. Environmental and Experimental Botany, 2019, 159: 1-12.
45
刘艳伟,朱仲元,朝伦巴根 等.水分胁迫条件下羊草群落生育期双作物系数计算方法的验证[J].中国农村水利水电,2010(4):39-43.
LIU Y W, ZHU Z Y, Chaolunbagen, et al. Validation of estimating dual crop coefficient approach in growing season for the community of Leymus Chinensis under water stress conditions[J].Rural Water Conservancy and Hydropower in China,2010(4):39-43.
46
杨天一, 王军, 张红梅, 等. 基于单作物系数法的华北平原典型农业生态系统蒸散规律研究[J]. 中国生态农业学报 (中英文), 2022, 30(3): 356-366.
YAN T Y, WANG J, ZHANG H M, et al. Evapotranspiration of typical agroecosystems in the North China Plain based on single crop coefficient method[J]. Chinese Journal of Eco-Agriculture, 2022, 30(3): 356-366.
47
杜江涛,张 楠,龚珂宁,等. 基于单作物系数法指导南疆棉花膜下滴灌试验研究[J]. 中国农村水利水电2021(11):141-148.
DU J T, ZHANG N, GONG K N, et al. Experimental research on the irrigation system of mulched-cotton by drip irrigation in south Xinjiang based on the single crop coefficient[J]. China Rural Water and Hydropower2021(11):141-148.
48
余昭君,胡笑涛,冉辉.基于波文比-能量平衡法的半湿润地区葡萄园蒸发蒸腾量估算[J].干旱地区农业研究 2020(4):175-183.
YU Z J, HU X T, RAN H. Estimation of grape evapotranspiration in semi-humid region based on Bowen ratio energy balance method[J]. Agricultural Research in the Arid Areas, 2020(4):175-183.
49
樊引琴,蔡焕杰.单作物系数法和双作物系数法计算作物需水量的比较研究[J].水利学报,2002(3):50-55.
FAN Y Q, CAI H J. Comparison of crop water requirements computed by single crop coefficient approach and dual crop coefficient approach[J]. Journal of Hydraulic Engineering, 2002(3):50-55.
50
耿耘,刘浩,李云峰,等.麦后移栽棉蒸发蒸腾规律和作物系数[J].灌溉排水学报,2022, 41(7): 24-34.
GENG Y, LIU H, LI Y F, et al. Evapotranspiration and crop coefficient of transplanted cotton after wheat harvest[J]. Journal of Irrigation and Drainage, 2022, 41(7): 24-34.
51
王帅兵,李常斌,杨林山,等 .Budyko方程和单作物系数法在区域蒸散发估算中的耦合及应用[J].中国沙漠2015,35(3):683-689.
WANG S B, LI C B, YANG L S, et al. Application of the coupled Budyko Equation and single crop coefficient approach on regional evapotranspiration estimation[J].Journal of Desert Research,2015,35(3):683-689.
52
李毅,付亚亚,唐德秀,等.砂石覆盖条件下冬小麦蒸散量的单、双作物系数法估算[J].农业机械学报,2018,49(3):261-271.
LI Y, FU Y Y, TANG D X, et al. Estimation of evapotranspiration of winter wheat based on single and dual crop coefficient approaches under sand-gravel mulching conditions[J]. Journal of Agricultural Machinery,2018,49(3):261-271.
53
丁日升.干旱内陆区玉米田水热传输机理与蒸散发模型研究[D].北京:中国农业大学,2012.
DING R S. Study on Mechanism of water and heat transport and evapotranspiration model in maize field in arid inland region[D]. Beijing: China Agricultural University,2012.
54
冯禹,崔宁博,龚道枝,等. 基于叶面积指数改进双作物系数法估算旱作玉米蒸散[J]. 农业工程学报201632(9):90-98.
FENG Y, CUI N B, GONG D Z, et al. Estimating rainfed spring maize evapotranspiration using modified dual crop coefficient approach based on leaf area index[J]. Transactions of the Chinese Society of Agricultural Engineering201632(9):90-98.
55
龚雪文,刘浩,孙景生,等.基于双作物系数法估算不同水分条件下温室番茄蒸发蒸腾量[J].农业生态学报,2017,28(4):1 255-1 264.
GONG X W, LIU H, SUN J S, et al. Modeling evapotranspiration of greenhouse tomato under different water conditions based on the dual crop coefficient method[J]. Journal of Agricultural Ecology,2017,28(4):1 255-1 264.
56
ZHAO P, KANG S, LI S, et al. Seasonal variations in vineyard ET partitioning and dual crop coefficients correlate with canopy development and surface soil moisture[J]. Agricultural Water Management, 2018, 197: 19-33.
57
黄松.苏南地区典型作物蒸腾和土面蒸发模型及参数的研究[D]. 江苏镇江:江苏大学, 2020.
HUANG S. Parametrization of resistances for modeling crop transpiration and soil evaporation in the Southern Jiangsu[D]. Jiangsu Zhenjiang: Jiangsu University,2020.
58
戚迎龙,赵举,史海滨,等. 覆土浅埋滴灌玉米田双作物系数模型参数全局敏感性分析[J]. 农业工程学报202036(7):99-108.
QI Y L, ZHAO J, SHI H B, et al. Global sensitivity analysis of model parameters of dual crop coefficient in subsurface drip irrigated maize field[J]. Journal of Agricultural Engineering202036(7):99-108.
59
马嘉敏.覆膜对温室滴灌黄瓜水热通量、生理形态及产量的影响[D]. 江苏镇江:江苏大学, 2021.
MA J M. The effect of film mulching on vapor and heat fluxes, morpho-physiological parameters and yield of cucumber in greenhouse[D]. Jiangsu Zhenjiang: Jiangsu University, 2021.
60
李丰琇,马英杰.基于双作物系数法的新疆覆膜滴灌夏玉米蒸散量估算[J].农业机械学报,2018(11):268-275.
LI F X, MA Y J. Evapotranspiration estimation of summer maize with plastic mulched drip irrigation based on dual crop coefficient approach in Xinjiang[J].Journal of Agricultural Machinery,2018(11):268-275.
61
于文颖,纪瑞鹏,贾庆宇,等 .基于改进的双作物系数法估算辽河三角洲芦苇湿地蒸散量[J].生态学报202040(1) : 325-335.
YU W Y, JI R P, JIA Q Y, et al. Evapotranspiration estimation of phragmites australis wetland in the Liaohe River Delta based on the improved dual crop coefficient method[J]. Journal of Ecology202040(1) : 325-335.
62
刘艳飞.通风与水分调控对温室湿热环境及番茄耗水影响的研究[D]. 郑州:华北水利水电大学,2020.
LIU Y F. Study on the influence of ventilation and water regulation on the humid and hot environment of greenhouse and water consumption of tomato[D]. Zhengzhou: North China University of Water Resources and Hydropower,2020.
63
赵娜娜,刘 钰,蔡甲冰,等. 双作物系数模型 SIMDual_Kc 的验证及应用[J]. 农业工程学报201127(2):89-95.
ZHAO N N, LIU Y, CAI J B, et al. Validation and application of dual crop coefficient model SIMDual_Kc[J]. Journal of Agricultural Engineering201127(2):89-95.
64
石小虎,蔡焕杰,赵丽丽,等. 基于 SIMDualKc 模型估算非充分灌水条件下温室番茄蒸发蒸腾量[J]. 农业工程学报201531(22):131-138.
SHI X H, CAI H J, ZHAO L L, et al. Estimation of greenhouse tomato evapotranspiration under deficit irrigation based on SIMDualKc model[J]. Journal of Agricultural Engineering201531(22):131-138.
65
邱让建,杜太生,陈任强. 应用双作物系数模型估算温室番茄耗水量[J]. 水利学报,2015(6):678-686.
QIU R J, DU T S, CHEN R Q. Application of the dual crop coefficient model for estimating tomato evapotranspiration in greenhouse[J].SHUI LI XUE BAO,2015(6):678-686.
66
卢晓鹏,段顺琼,马显莹,等.单双作物系数法计算玉米需水量的对比研究[J].节水灌溉, 2012(11):18-22.
LU X P, DUAN S Q, MA X Y, et al. A comparative study between single crop coefficient method and double crop coefficient method in calculation of water requirement of maize[J].Water Saving Irrigation, 2012(11):18-22.
PDF(678 KB)

429

Accesses

0

Citation

Detail

Sections
Recommended

/