基于作物系数与水分生产函数的向日葵产量预测

云文丽 侯 琼 李建军 苗百岭 冯旭宇

1)(内蒙古气象科学研究所，呼和浩特 010051)2)(内蒙古巴彦淖尔市临河观象台，临河 015000)3)(内蒙古生态与农业气象中心，呼和浩特 010051)



基于作物系数与水分生产函数的向日葵产量预测

云文丽1)*侯 琼1)李建军2)苗百岭1)冯旭宇3)

1)(内蒙古气象科学研究所，呼和浩特 010051)2)(内蒙古巴彦淖尔市临河观象台，临河 015000)3)(内蒙古生态与农业气象中心，呼和浩特 010051)

利用河套灌区向日葵2012年田间水分、分期播种试验数据和两个站点的农业气象历史资料，研究基于向日葵作物系数和水分生产函数的产量预测方法。结果表明：向日葵标准作物系数在生育期内的变化规律是前期小、中期大、后期小,最高值为1.21,出现在开花期。标准作物系数与出苗后日数和大于0℃积温有很好的二次和三次多项式关系，拟合优度在0.93以上。在分析相对叶面积指数和作物系数关系的基础上，提出标准作物系数的相对叶面积指数订正方法，得出河套灌区向日葵作物系数的动态计算式，为水分生产函数中实际蒸散量的计算提供支撑。建立以Jensen模型为基础的向日葵水分生产函数，得到对水分亏缺的敏感顺序从高到低是开花期、花序形成期、成熟期、苗期。综合应用向日葵作物系数方程和水分生产函数模型计算分期播种产量，与实际产量分别相差4.4%和4.1%，初步证明该文提出的方法对产量预测较为理想，在该地区具有很好的适用性。

作物系数； 相对叶面积指数订正； 耗水量； 水分生产函数； 产量预测

引 言

近年来，河套灌区向日葵种植面积占到总播面积的四分之一以上，是经济效益十分显著的作物。河套灌区用于保障粮食生产、农业丰收以及生态用水占到内蒙古自治区引黄用水量的80%以上，而且对水资源的需求量逐年加大，对于光、温资源丰富的该区域，水分成为限制农业发展的主要因子[1]。因此，同时考虑光、温、水3个因子的作物产量气象预报方法[2-4]和考虑参数较多的作物生长模式[5-7]不适合河套灌区，需主要考虑限制因子水分对产量预测的影响，即依据水分生产函数模型的建立预测产量。水分生产函数的主要计算因子是需水量和农田实际蒸散量，而作物系数是计算二者的重要参数。它反映了作物本身的生物学特性、产量水平、土壤耕作条件等对作物耗水量的影响，其最合理的确定方法是采用当地的实测资料[8]。现有研究主要是针对小麦、水稻、玉米等作物给定特定生长阶段的作物系数[9-13]，在向日葵上仅见到文献[14-15]对油用向日葵的阶段作物系数进行了确定，但对食用向日葵尚未见报道，因此，有必要进行食用向日葵全生育期作物系数的动态计算，通过结合参考作物蒸散量，进而确定实际蒸散量，该方法目前最为普遍[16]。实际蒸散量的动态计算为水分生产函数的建立和产量计算奠定了基础。国内外对水分生产函数的研究很多，也建立了一些模型，从这些研究可以看出,不同地区、同一作物的数学模型存在较大差别[17-18]。对于依靠黄河灌溉的河套地区，合理、正确选用满足当地的向日葵水分生产函数模型是非常必要的。本文旨在建立向日葵的作物系数动态计算式，为准确计算需水量和实际蒸散量提供方法。根据产量与不同生育阶段的耗水量，采用乘法和加法模型，选取、建立适合河套灌区向日葵的水分生产函数模型。通过作物系数和水分生产函数二者的应用，试图进一步预测不同水分水平对作物产量的影响，更好地为粮食安全和预警服务。

1 试验与方法

1.1 试验材料与设计

试验设在内蒙古巴彦淖尔市农业气象试验站(40°45′N，107°25′E，海拔为1039.3 m)。该区属于干旱半干旱气候带，多年平均降水量为150.0 mm左右，年平均气温为8.0℃。土壤为壤土，0～1 m土壤容重平均值为1.36 g·cm-3， 田间持水量平均为26.8%，地下水埋深约5.8 m。

试验于2012年6—10月进行。试验水分控制指标设4类，即严重缺水、缺水、水分适宜和偏湿,分别占田间持水量(0,40%]，(40%，55%]，(55%，70%]，(70%，90%]。水分控制时段分3类：二对真叶期-花序形成期、花序形成期-开花期、开花期-成熟期，其他时段水分按55%～70%供给，各处理设2个重复。试验每小区面积为4 m×6 m，小区之间设置有隔离层，试验场设有大型活动式防雨棚。灌水为井灌，灌水量通过水表记录。

供试向日葵品种为LD5009，属中熟种，生育期平均为110 d。6月2日播种，9月18 日收获。大小行种植，种植密度为47250株/hm2左右，田间管理等均以当地大田管理为准。

本文同时利用临河(2005—2011年)、乌拉特前旗(1987—2011年)向日葵农业气象观测的历史资料和2012年分期播种数据，包括土壤水分、发育期和产量等观测资料。

1.2 测定项目与方法

土壤水分测定：采用CPN-503中子土壤水分仪(美国)定点观测，测深0～200 cm，每隔10 cm为1层，测定时间每隔5 d测定1次。表层(0～30 cm)采用土钻取土烘箱烘干法进行测定。利用0～50 cm 的土壤含水量确定小区是否处于控水范围。

叶面积测定：每个小区选择有代表性的植株5株作为叶面积的定株观测样本。二对真叶期开始人工测定叶片最长处的长度与宽度，每5 d测量1次，开花后每7 d观测1次，直至收获。测定向日葵13～14叶期和近开花期的平均叶面积系数为0.64。

产量测定：收获时，每个小区选10株进行考种，分析产量结构。

1.3 研究方法

本文首先基于向日葵田间水分试验数据和农业气象历史资料，采用农田水分平衡方程计算作物实际蒸散量 (见式(1))，当农田充分供水时的实际蒸散量可以近似地认为是作物标准蒸散量，利用作物系数计算方法，建立向日葵全生育期作物系数动态计算式。同时基于阶段实际蒸散量的计算结果和需水量的确定，选取并建立了向日葵水分生产函数模型。最后应用本文建立的向日葵作物系数计算方程和水分生产函数模型计算了分期播种试验中的阶段实际蒸散量和产量，以分期播种中实际观测数据对计算结果进行验证。

1.3.1 实际蒸散量计算

利用农田水分平衡方程计算作物实际蒸散量[19]:

G-R0-Lg。

(1)

式(1)中，ETa为一定时间间隔的实际蒸散量(单位：mm),本文为整个生长季和4个发育时段(出苗期-二对真叶期、二对真叶期-花序形成期、花序期-开花期、开花期-成熟期)，n为取土层次，hi和di分别为第i层土层厚度(单位：cm) 和土壤容重(单位：g·cm-3),W1和W2分别为时段开始和结束时的第i层土壤湿度(计算中分别乘100)，P和I为时段内降水量和灌溉量(单位：mm)；G为地下水补给量[19]：

(2)

计算4—9 月的潜水蒸发量,t为旬数, 4 月第1旬为0；R0和Lg分别为地表径流量(单位：mm) 和渗漏量到计划层以下的水量(单位：mm )，本文将超出计划层田间持水量的水分作为R0和Lg处理。

1.3.2 作物系数计算

作物系数描述了实际作物在农田下垫面反射率、作物冠层特性及空气动力学阻力等方面与参考作物(地表覆盖均匀，土壤水分充足，生长旺盛，株高为12 cm，具有下垫面反射率为0.23和冠层阻力为70 s·m-1)的不同[20]。作物系数计算如下

Kc=ETc/ETo。

(3)

式(3)中，Kc为作物系数，ETc为标准条件下的作物蒸散量(单位:mm)，指在适宜的土壤水分条件下，作物正常发育，获得较高产量的作物蒸腾和棵间蒸发量之和。农田水分平衡方程通常用来计算实际蒸散量ETa,而农田充分供水时的实际蒸散量可以近似地认为是作物标准蒸散量ETc，因此，这种条件下用式(1)计算作物标准蒸散量ETc；ETo为参考作物蒸散量(单位：mm)，采用FAO推荐的Penman-Montieth公式计算，被广泛用于农田参考作物潜在蒸散量的计算[21-22]。

1.3.3 土壤水分修正

(4)

式(4)中，θ为土壤含水量(单位：%)；θp为土壤凋萎湿度(单位：%)；θj为土壤临界含水量(单位：%)，等同于毛管断裂湿度，可取田间持水量的75%[24]。

1.3.4 水分生产函数模型选取

选取乘法模型(Jensen模型、Minhas模型)和加法模型(Blank模型、Singh模型)对河套灌区向日葵水分生产函数进行分析[25]。

(5)

(6)

(7)

(8)

其中，λi，Ai，Bi为作物不同阶段缺水对产量的敏感指数及敏感系数；i为生育期；Ya为处理条件下的实际产量(单位：g·m-2)；Ym为正常灌溉下产量(单位：g·m-2)；ETa为各处理条件下实际蒸散量(单位：mm)；ETm为阶段需水量(单位：mm)，通过水分适宜且获得较高产量的小区作物耗水量(在任何土壤水分条件下的实际蒸散量)确定； 为模型的阶段总数，在模型的计算中n=4。

根据向日葵水分试验实测数据，利用最小二乘法转换为求解线性方程组,得到水分敏感指数(系数)。

2 结果分析

2.1 作物系数确定

2.1.1 水分适宜条件下标准作物系数计算

利用水分适宜小区和降水丰沛年份历史资料，通过式(1)、式(3)和联合国粮食与农业组织(FAO)近期推荐的Penman-Montieth公式计算参考作物蒸散量、实际蒸散量、作物标准蒸散量和作物系数，其中，作物标准蒸散量是农田充分供水时的实际蒸散量，因而也通过式(1)进行计算，计算结果作物系数和作物发育期日序绘图(图1)，二者呈二次多项式关系(表1)，拟合较好，拟合优度(R2)为0.927。向日葵作物系数在生育期内的变化规律是前期小、中期大、后期又小, 其中最高值为1. 21, 出现在播后65 d 左右, 即向日葵开花期。进一步分析大于0℃的积温与作物系数的关系(图2)，呈三次多项式关系(表1)，曲线和实测数据拟合较好，拟合优度为0.949。

图1 标准作物系数与日序的关系Fig.1 Relationship of standard crop coefficient to date of year

图2 标准作物系数与大于0℃的积温的关系Fig.2 Relationship of standard crop coefficient to accumulated temperature greater than 0℃

参数模拟方程式R2样本量标准作物系数Kc=-2E-10T3j+8E-0.7T2j-0.0005Tj+0.26380.92726Kc=-0.0002D2+0.1011D-9.55280.94926叶面积指数LAID=-3E-0.9T3j+2E-0.5T2j-0.0276Tj+14.2360.96524

注：D表示日序，150

对向日葵整个生育期的标准作物系数按照FAO的时段划分进行平均值求算，并与FAO推荐阶段值进行比较[20，26](表2)。由表2可以看出：试验计算的作物系数在初期、前期和后期均较FAO推荐值的阈值范围整体略高，中期和收获期的计算结果较推荐值的阈值范围大，特别是收获期计算的阈值范围为0.24～0.55，FAO推荐值的阈值范围仅为0.35～0.45，整个生育期的计算平均值在推荐值的范围中，即试验计算所得标准作物系数较为合理。

表2 标准作物系数与FAO推荐值的比较Table 2 Compared with FAO recommended values and standard crop coefficient

注：初期:播种-出苗；前期:出苗至13，14片叶；中期：13，14片叶至开花后1周；后期：开花后1周至开花后20 d； 收获期：开花后20 d至成熟。

2.1.2 作物系数叶面积指数订正

当作物受到不同程度的水分亏缺时，叶面积指数成为表征作物生长发育状况的重要指标,其大小与实际蒸散量存在密切关系, 因而可以用作物相对叶面积指数对作物系数进行水分胁迫订正[27-28]。利用水分适宜时生长季逐旬叶面积指数测定值与相应时段大于0℃的积温建立了叶面积指数动态变化方程，用以计算不同时期非水分胁迫条件下的期望叶面积指数(表1)。通过分析不同水分梯度下测定的相对叶面积指数与同期计算的作物系数的关系发现，二者可以用指数回归较好地描述(图3),即

Kci=0.0569e2.73LAI/LAID(R2=0.618，n= 130)。

(9)

式(9)中，Kci为作物系数相对叶面积指数的订正系数,LAID可用表1中的公式计算。

由此，得出河套灌区向日葵的作物系数动态计算式可表示为

Kcx=Kc×Kci。

(10)

式(10)中，Kcx为受水分胁迫影响下的作物系数，Kc为由时间变量或热量变量方程计算的标准作物系数(表1)，Kci为水分胁迫下相对叶面积指数进行订正的作物系数(式(9))。

图3 作物系数与相对叶面积指数关系Fig.3 Relationship of crop coefficient to relative leaf area index

2.2 水分生产函数的建立

2.2.1 总产量与总耗水量关系

利用临河和乌拉特前旗农业气象历史观测资料，结合向日葵水分试验数据，利用式(1)计算整个生育期的实际蒸散量,即作物耗水量,进一步分析向日葵耗水量与产量的关系(样本量为40)，二者呈二次抛物线(图4)，由图4可以看出，在耗水量较低的情况下，随着耗水量的增加，产量上升较快，但当耗水量达到一定数值时，产量达到最大，之后随着耗水量的增加，产量呈现下降趋势，二者之间存在一个适宜耗水界限值[29]。根据拟合方程计算得出向日葵适宜的耗水界限值大约为400～460 mm，产量为496.7～500.6 g·m-2。利用拟合方程对计算的产量结果进行回代，得到图5，可以看出,拟合方程计算的产量和实际产量基本吻合，拟合方程合理。实际产量和回代产量相关系数为0.934，R2为0.9055，平均相对误差为4.37%。

图4 总耗水量和总产量的关系Fig.4 Relationship of total water consumption to total yield

图5 计算的产量和实际产量的比较Fig.5 Comparison of calculated and actual yields

2.2.2 不同生长发育时段需水量确定

选取全生育期耗水量在400～460 mm之间且产量较高的试验小区和历史年份，将向日葵按生育期分成4个阶段,根据式(1)进行生育阶段实际蒸散量即耗水量的计算，由于供水充分,时段的耗水量既为向日葵各生育阶段正常灌溉下的需水量(表3)，对应产量为510 g·m-2。

竹韵有些犯难了，自己跟海力不过是普通主雇关系，上班拿了工资奖金，岂能再收他如此贵重的礼物？从马丽亚眼里的敌视来看，她已经隐隐感觉出了些什么。但是，不收又如何拒绝呢？何不暂时把电脑收下，今后努力工作来回报公司，或者今后用工资付帐也行，其实她也早就想给龙斌买台电脑了。想到这里，她伸手接过电脑，脸上也露出了笑容。

由表3可知, 在当地地理和气候条件下, 整个生育期内向日葵在充分供水时的需水量为450 mm,平均耗水强度为4.09 mm·d-1。向日葵的需水规律表现为苗期需水强度最小, 平均日耗水仅为0.56 mm·d-1，总耗水量为5.0 mm；二对真叶期-花序形成期和开花期-成熟期的日需水强度居中,且日耗水强度相差较小；花序形成期-开花期的需水量和日耗水强度最高，是需水高峰期，这阶段占总需水量的42.22%，平均日耗水达9.5 mm·d-1。

表3 向日葵各生育阶段需水量Table 3 Water requirement of sunflower in different growth stages

2.2.3 水分生产函数建立

根据式(5)～式(8)计算得到4种水分生产函数模型中的敏感指数及检验参数值(表4)。

模型公式表明:λi值越高，缺水后Ya/Ya越低，即由于缺水导致的减产越严重(对缺水越敏感)[30]。Jensen模型和Blank模型λi值从高到低的顺序为开花期、花序形成期、成熟期、苗期，与前期研究所得的向日葵不同生长发育阶段的土壤水分对光合生理特征的影响[31]，以及通过分析向日葵的生长特征得到的整个生育期的适宜土壤含水量和实际的灌溉经验是一致的[1]，即向日葵在花序形成期-开花期阶段耗水量最大，对缺水最为敏感，在苗期向日葵矮小，向日葵对水分要求较低，可见敏感值的顺序合理。而Minhas模型和Singh模型中的敏感指数顺序不符合向日葵的水分生理特性和灌溉实践。进一步分析Jensen模型和Blank模型拟合计算产量和实际产量的相关系数可以看出，均达到0.01显著性水平，但相关性最好的是Jensen模型，相关系数为0.806。同时分析计算产量和实际产量之间的平均相对误差，Jensen模型的误差仅为10.78%，而Blank模型的误差达到32.03%。所以最终应以Jensen模型作为河套灌区向日葵水分生产函数模型为宜。

表4 河套灌区向日葵水分生产函数模型敏感指标及检验参数(n=18)Table 4 Sensitive index and parametric test of sunflower crop water production function in Hetao irrigated area(n=18)

注：*表示达到0.01显著性水平。

根据以Jensen模型为基础建立的向日葵水分生产函数，对临河和乌拉特前旗24年的历史农业气象观测资料进行回代，与实际产量进行对比。由图6可以看出，回代结果较好。分析二者的误差情况，平均相对误差为7.25%，误差在±10%以内的占62.5%，进一步证实了Jensen模型在河套地区的适用性。

图6 Jensen模型计算产量和实际产量比较Fig.6 Comparison of Jensen model calculated and actual yields

2.3 基于作物系数和水分生产函数的产量预测

表5 作物系数方程和水分生产函数模型计算及产量预测Table 5 Calculation of crop coefficient equation and water production function model and yield prediction

3 结论与讨论

本文以内蒙古巴彦淖尔市农业气象试验站2012年向日葵田间水分、分期播种试验数据和两个站的农业气象历史资料为基础，建立向日葵的作物系数动态计算式和水分生产函数模型，应用二者进行产量预测方法适用性研究，得到以下主要结论：

1) 向日葵标准作物系数在生育期内的变化规律是前期小、中期大、后期又小,最高值为1.21,出现在开花期。标准作物系数与出苗后日数和大于0℃积温有很好的二次和三次多项式关系，拟合优度在0.93以上。在分析相对叶面积指数和作物系数关系的基础上，提出标准作物系数的相对叶面积指数订正方法，得出河套灌区向日葵作物系数的动态计算式，为水分生产函数中实际蒸散量的计算提供了支撑。

2) 对比4种水分生产函数模型中的敏感指数顺序，以及计算产量和实际产量之间的平均相对误差，确定以Jensen模型作为河套灌区向日葵水分生产函数模型。向日葵对水分亏缺的敏感顺序从高到低是开花期、花序形成期、成熟期、苗期。

3) 综合应用向日葵作物系数方程和水分生产函数模型计算分期播种产量，与实际产量分别差4.4%和4.1%，初步证明本文提出的方法对产量的预测结果较为理想，在该地区具有很好的适用性。

作物系数的确定受到许多因素的影响。在特定的土壤、气候和管理方式下,当作物不发生水分亏缺时, 作物系数主要受作物生长发育进程和大气蒸发力的影响。文中确定的食用向日葵标准作物系数符合FAO推荐值的阈值范围，但在不同生长阶段均比2011年戴佳信等[14]确定的油用向日葵的双值作物系数偏高，这可能与向日葵品种关系密切。当作物受到不同程度的水分亏缺时，叶面积指数成为表征作物生长发育状况的重要指标,其大小与作物耗水量存在密切关系,因而作物系数可通过作物相对叶面积指数的大小反映[28]。叶面积指数的订正精度与它对实际蒸散量的影响程度有关，同时对整个生长发育期的修正精度会低于对单一生长发育期的精度。

本文通过作物系数对实际蒸散量的计算，结合建立的向日葵水分生产函数的应用，对产量的计算结果和实际产量相差较少，可见该方法不仅可为节水灌溉提供一定的科学依据，同时也为不同阶段缺水对作物产量影响的预评估进行了很好的尝试，但该方法在预评估方面的应用还有待进一步深入。

[1] 云文丽,李建军,侯琼.土壤水分对向日葵生长状况的影响.干旱地区农业研究，2014,32(2)：186-190.

[2] 王叔同,张荣霞,张敏.基于长时段非线性影响分析的冬小麦产量预报.应用气象学报,2000,11(3):377-382.

[3] 魏瑞江,宋迎波,王鑫.基于气候适宜度的玉米产量动态预报方法.应用气象学报,2009,20(5):622-627.

[4] 刘振忠,徐梅.三江平原地区粮食产量预报模型研究.农业工程学报，1999,15(4)：14-18.

[5] 刘春，张春辉，郭萨萨．基于能量模型的水稻生长模型．应用气象学报,2013，24(2)：240-247．

[6] 帅细强,王石立,马玉平,等.基于水稻生长模型的气象影响评价和产量动态预测.应用气象学报,2008,19(1):71-81.

[7] 熊伟.CERES-Wheat模型在我国小麦区的应用效果及误差来源.应用气象学报,2009,20(1):88-94.

[8] 刘钰，Pereira L S.对FAO推荐的作物系数计算方法的验证.农业工程学报，2000,16(5)：26-30.

[9] 王笑影，梁文举，闻大中.北方稻田蒸散需水分析及其作物系数确定.应用生态学报,2005，16(1):69-72.

[10] 刘海军，康跃虎.冬小麦拔节抽穗期作物系数的研究.农业工程学报,2006,22(10):52-56.

[11] 李玉霖,崔建垣,张铜会.奈曼地区灌溉麦田蒸散量及作物系数的确定.应用生态学报,2003,14(6):930-934.

[12] 王宇,周广胜.雨养玉米农田生态系统的蒸散特征及其作物系数.应用生态学报，2010，21(3):647-653.

[13] 陈凤,蔡焕杰,王健，等.杨凌地区冬小麦和夏玉米蒸发蒸腾和作物系数的确定.农业工程学报,2006,22(5):191-193 .

[14] 戴佳信,史海滨,田德龙，等.内蒙古河套灌区主要粮油作物系数的确定.灌溉排水学报,2011,30(3):23-27.

[15] Vu S H，Watanabe H，Takagi K.Application of FAO-56 for evaluating evapotranspiration in simulation of pollutant runoff from paddy rice field in Japan.AgriculturalWaterManagement，2005，76:195-210.

[16] Sentelhas P C，Gillespie T J，Santos E A.Evaluation of FAO Penman-Monteith and alternative methods for estimating reference evapotranspiration with missing data in Southern Ontario，Canada.AgriculturalWaterManagement，2010，97: 635-644.

[17] 刘增进,李宝萍,李远华，等.冬小麦水分利用效率与最优灌溉制度的研究.农业工程学报,2004,20(4):58-63.

[18] 迟道才,王瑄,夏桂敏,等.北方水稻动态水分生产函数研究.农业工程学报,2004,20(3):30-34 .

[19] 华北平原作物水分胁迫与干旱研究课题组.作物水分胁迫与干旱研究.郑州：河南科学技术出版社，1991.

[20] Allen R G，Pereira L S，Raes D，et al.Crop Evapotranspiration：Guidelines for Computing Crop Water Requirements.FAO Irrigation and Drainage Paper No.56.Rome：Food and Agricultural Organization of the United Nations，1998.

[21] Chiew F H S,Kamadalasa N N，Malano H M，et al.Penman-Monteith，FAO-24 reference crop evapotranspiraion and class-a pan data in Australia.AgriculturalWaterManagement，1995，28：9-12.

[22] David M S，Jenifer M J.Utility of Penman-Monteith，Priestley-Taylor，reference evapotranspiraion，and pan evaporation methods to estimate pasture evapotranspiraion.JHydrol，2005,308：81-104.

[23] Jensen E.耗水量与灌溉需水量.马文，译.北京：农业出版社，1982.

[24] 王玉娟，杨胜天，刘昌明，等.植被生态用水结果及绿水资源消耗效用——以黄河三门峡地区为例.地理研究，2009,28(1)：74-84.

[25] 王立坤,单玉芬,马永胜,等.黑龙江省西部半干旱区向日葵水分生产函数试验研究.东北农林大学学报，2009,40(4)：33-36.

[26] 吕厚荃，张玉书，王建林,等.农业干旱等级(审定稿).GB/T—2009，2009.

[27] 侯琼，王英舜，师桂花,等.内蒙古典型草原作物系数的试验研究.草业学报，2011，20(4):34-41.

[28] 刘浩,孙景生,梁媛媛,等.温室滴灌条件下土壤水分亏缺对番茄产量及其形成过程的影响.应用生态学报,2011,22(5):1201-1206.

[29] 侯琼，沈建国.内蒙古主要灌区春小麦、春玉米农田优化灌溉指标研究.干旱地区农业研究，2001,19(1)：33-39.

[30] 云文丽,侯琼,王海梅，等.不同土壤水分对向日葵光合光响应的影响.应用气象学报，2014,25(4):476-482.

[31] 翟胜，梁银丽，王巨媛,等.干旱半干旱地区日光温室黄瓜水分生产函数的研究.农业工程学报，2005,21(4)：136-139.

Yield Prediction of Sunflower Based on Crop Coefficient and Water Production Function

Yun Wenli1)Hou Qiong1)Li Jianjun2)Miao Bailing1)Feng Xuyu3)

1)(InnerMongoliaInstituteofMeteorologicalSciences,Hohhot010051)2)(LinheAstronomicalObservatory,BayannaoerMeteorologicalBureauofInnerMongolia,Linhe015000)3)(InnerMongoliaEcologicalandAgriculturalMeteorologicalCenter,Hohhot010051)

Crop coefficient and water production function are important parameters for water saving irrigation. Through making use of data from the stage sowing test at Bayannaoer (40°45′N，107°25′E，elevation 1039.3 m) of Inner Mongolia in 2012 and historical agro-meteorological data from two monitoring stations, most of the research is carried out in accordance with yield prediction methods based on crop coefficient and water production function. As far as results are concerned, variations in the standard sunflower crop are small in the early stage, large in the medium stage and small again in the late stage. The peak value (1.21) presents itself in the blossom period. There are strong quadratic and cubic polynomial relationships amongst standard crop coefficients days after germination and positive accumulated temperature values (determination coefficient is 0.93). Through comparison with FAO recommended stage values, standard crop coefficient from the test computing is reasonable. In addition, when calibration methods of standard crop coefficient and relative leaf area index are put forward, the actual evapo-transpiration of water production function can be calculated, and dynamic calculation equations of sunflower crop coefficient in the irrigated districts are obtained. There is a quadratic parabola relationship between sunflower water consumption and yield with a suitable water consumption limit. The suitable water consumption threshold is about 400-460 mm, and the yield is 496.7-500.6 g·m-2. Moreover, when water supply is adequate, the water requirement of sunflowers during the entire growth period is 450 mm with an average frequency of 4.09 mm/d. The regular water requirement pattern indicates that minimum water is required during seeding stage, medium water is required during two pairs of true leaves-inflorescence formation stage and blossom-maturity stage, and maximum water is required during inflorescence formation-blossom stage. Additionally, when Jensen model is put forward and established through comparison with 4 sensitive indexes, the order of water deficit from high to low is blossom period, inflorescence formation period, maturity period and seeding period, which is consistent with the regular water requirement pattern. Through integrated utilization of sunflower crop coefficient equation and water production function model, the stage sowing production and production are obtained (504.36g·m-2and 493.83 g·m-2, respectively), which show 4.4% and 4.1% deviations with actual production, respectively. There is preliminary evidence that the prediction method of production proposed is relatively reasonable with a great applicability in this region, and can be further applied to pre-assessment of production affected by water deficit in different stages.

crop coefficient； correction of leaf area index； water consumption； crop water production function； yield prediction

10.11898/1001-7313.20150607

公益性行业(气象)科研专项(GYHY201206021)，内蒙古自治区科技计划项目(20140427)

云文丽,侯琼,李建军,等. 基于作物系数与水分生产函数的向日葵产量预测. 应用气象学报，2015，26(6)：705-713.

2015-01-05收到， 2015-07-22收到再改稿。

* email: 564838918@qq.con