Experimental Study on Crop Coefficient of Spring Maize in Hetao Irrigation District of Inner Mongolia
-
摘要: 作物系数曲线是估算作物生长季耗水量变化的重要参数。基于2013年4—9月内蒙古巴彦淖尔市临河区田间水分试验和1994—2013年气象站观测资料,利用水量平衡法反求春玉米作物系数,分析生长季内的变化规律, 建立动态模拟方程,并与联合国粮农组织 (FAO) 分段直线法结果进行比较, 提出胁迫条件下作物系数的叶面积修正方法。结果表明:玉米作物系数随发育进程可用三项式曲线描述,变化趋势与产量水平无关, 但随产量增高而变幅增大;以出苗后相对积温为时间变量建立模拟方程效果较好,决定系数 (R2) 均在0.92以上;模拟计算出各站点最大 (1.30~1.48) 和平均 (0.831~0.919) 作物系数,与FAO分段直线法计算的典型值和区间值基本一致,生长中期平均相对误差为3.4%~7.2%;提出利用相对叶面积指数修正作物系数的计算方法;通过2014年实例检验,土壤水分模拟值与实测值的平均相对误差为6.3%,相对误差小于15%的占95.8%。Abstract: The crop coefficient curve is an important parameter for estimating the change of water consumption in growing season, and it plays an important role in water management, such as the simulation of evapotranspiration, irrigation forecasting and irrigation decision-making. Previous crop coefficient studies mainly concentrate on the average value of growth period, while rarely focus on daily changes. A crop coefficient simulation method is in need to improve spring maze irrigation forecasting business in Hetao area. Based on the field moisture test data and meteorological station historical observations, the crop coefficient of spring maize is calculated with water balance method, and a dynamic simulation equation is established considering the variation during growing season. It's then evaluated using results of the United Nations Food and Agriculture Organization (FAO) piecewise linear method, and daily rolling estimation of crop evapotranspiration is achieved. At the same time, the leaf area correction method is put forward to estimate the soil moisture content under the water stress, which can provide basis for the development of maize irrigation forecast. Results show that the crop coefficient of spring maize can be described by three curves of the development process, and the change trend of crop coefficient has nothing to do with the output level, but the variation range increases with the increase of output. Considering heat index can reduce the influence of geographical factors on crop coefficient, the simulation equation of maize coefficient are established based on relative accumulative temperature as time variable after emergence, and the decision coefficient are all above 0.92. The maximum (1.30-1.48) and average (0.831-0.919) maize coefficients of each site are calculated by simulation, results are basically the same as the 3 typical values and interval values got by FAO segmentation method, and the range of averaged relative error during growth period is 3.4%-7.2%. Through the analysis, it is concluded that Kc and relative leaf area index are better described by exponential function, and the calculation method of the standard leaf area index is proposed, which can calculate the crop coefficient in any production condition. The simulated soil moisture is consistent with the measured value with average relative error of 6.3%, and less than 15% for 95.8% circumstance, indicating good application prospects. As the soil moisture supply below 1 m is not considered yet, the model should be improved in the future to explore the calculation method of the lower layer water supply.
-
表 1 4个发育阶段出现日期和平均株高
Table 1 The emergence date of 4 developmental stages and the average height of the plant
地区 气象观测 生长初期 快速生长期 生长中期 生长末期 临河试验田 时间段 05-06—05-31 06-01—05-07 07-06—08-20 08-21—09-15 平均株高/cm 18 110 260 250 临河区 时间段 05-10—06-04 06-05—07-10 07-11—08-25 08-26—09-20 平均株高/cm 110 275 270 准格尔旗 时间段 05-15—06-09 06-10—07-15 07-16—08-25 08-26—09-25 平均株高/cm 120 240 225 土左旗 时间段 05-08—06-02 06-03—07-08 07-09—08-20 08-21—09-15 平均株高/cm 100 266 260 表 2 水分适宜条件下玉米生长期作物系数模拟方程及主要参数 (n=15)
Table 2 Simulation equations and main parameters of crop coefficient of maize in growing season under water suitable conditions (n=15)
地区 表达式 R2 Kc
最大值Kc
平均值产量水平/
(kg·hm-2)临河试验田 Kcb=-5.14Tx3+3.84Tx2+1.795Tx+0.1461 0.9332 1.48 0.919 13501.1 临河区 Kcb=-4.95Tx3+3.59Tx2+1.635Tx+0.2135 0.9208 1.43 0.900 12753.8 准格尔旗 Kcb=-6.91Tx3+7.18Tx2-0.123Tx+0.2912 0.9356 1.38 0.868 11105.2 土左旗 Kcb=-5.36Tx3+4.82Tx2+0.866Tx+0.2128 0.9478 1.30 0.831 10755.3 注:Kcb为水分适宜时的作物系数;Tx表示生长过程相对积温, 单位:℃·d,计算方法为出苗后逐旬积温与生长季总积温的比值。 表 3 各发育阶段持续时间、出现时段和修正作物系数与FAO值的比较
Table 3 The duration and emergence date of development stages and the comparison between modified crop coefficient and FAO value
地区 生长期/d Kc 生长
初期快速
生长期生长
中期生长
末期生育期 生长
初期快速
生长期生长
中期生长
末期FAO推荐值[2] 30 40 50 30 150 0.15 0.15~1.2 1.2 1.2~0.6 临河试验田 25 35 45 25 130 0.2 0.2~1.27 1.27 1.27~0.69 临河区 25 35 45 25 130 0.2 0.2~1.25 1.25 1.25~0.65 准格尔旗 25 35 40 30 130 0.2 0.2~1.22 1.22 1.22~0.62 土左旗 25 35 42 25 127 0.2 0.2~1.19 1.19 1.19~0.60 表 4 不同生长时段多项式计算Kc结果的比较
Table 4 The comparison results of Kc calculated by polynomials of different development stages
地区 生长初期 快速生长期 生长中期 生长末期 平均值 区间 平均值 区间 平均值 区间 平均值 区间 FAO推荐值[2] 0.15 1.20 1.05~1.4 0.60 临河试验田 0.29 0.15~0.45 0.87 0.67~1.08 1.35 1.20~1.48 0.81 1.10~0.47 临河区 0.36 0.20~0.53 0.82 0.61~1.03 1.34 1.22~1.43 0.83 1.12~0.53 准格尔旗 0.40 0.30~0.51 0.84 0.65~1.03 1.27 1.15~1.38 0.89 1.18~0.59 土左旗 0.35 0.25~0.48 0.81 0.63~0.99 1.23 1.14~1.30 0.78 1.00~0.56 -
[1] 刘钰, 汪林, 倪广恒, 等.中国主要作物灌溉需水量空间分布特征.农业工程学报, 2009, 25(12):6-12. doi: 10.3969/j.issn.1002-6819.2009.12.002 [2] Allen R G, Pereira L S, Raes D, et al.Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56.Rome:FAO, 1998. [3] 毛飞, 张光智, 徐祥德.参考作物蒸散量的多种计算方法及其结果的比较.应用气象学报, 2000, 11(增刊Ⅰ):128-136. http://www.cnki.com.cn/Article/CJFDTOTAL-YYQX2000S1016.htm [4] 刘钰, Pereira L C.对FAO推荐的作物系数计算方法的验证.农业工程学报, 2000, 16(5):26-30. http://www.cnki.com.cn/Article/CJFDTOTAL-NYGU200005006.htm [5] 杜斌, 屈忠义, 于健, 等.内蒙古河套灌区大田作物膜下滴灌作物系数试验研究.灌溉排水学报, 2014, 33(4/5):16-20. http://www.cnki.com.cn/Article/CJFDTOTAL-GGPS2014Z1006.htm [6] 赵娜娜, 刘钰, 蔡甲冰.夏玉米作物系数计算与耗水量研究.水利学报, 2010, 41(8):953-969. http://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201008013.htm [7] 阳伏林, 张强, 王润元, 等.黄土高原半干作物系数特征.应用生态学报, 2013, 24(5):1209-1214. http://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201305005.htm [8] 左余宝, 田昌玉, 唐继伟, 等.鲁北地区主要作物不同生育期需水量和作物系数的试验研究.中国农业气象, 2009, 30(1):70-73. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGNY200901017.htm [9] 彭世彰, 索丽生.节水灌溉条件下作物系数和土壤水分修正系数试验研究.水利学报, 2004(1):17-21. http://www.cnki.com.cn/Article/CJFDTOTAL-SLXB200401004.htm [10] 连彩云, 马忠明, 吕晓东, 等.河西绿洲灌区主要作物需水量及作物系数试验研究.灌溉排水学报, 2012, 31(5):136-138. http://www.cnki.com.cn/Article/CJFDTOTAL-CYXB201105034.htm [11] 戴佳信, 史海滨, 田德龙, 等.内蒙古河套灌区主要粮油作物系数的确定.灌溉排水学报, 2011, 30(3):23-27. http://www.cnki.com.cn/Article/CJFDTOTAL-GGPS201103005.htm [12] 冯建设, 王建源, 王新堂, 等.相对湿润度指数在农业干旱监测业务中的应用.应用气象学报, 2011, 22(6):766-772. doi: 10.11898/1001-7313.20110616 [13] 吕厚荃, 钱拴, 杨霏云.华北地区玉米田实际蒸散量的计算.应用气象学报, 2003, 14(6):722-728. http://qikan.camscma.cn/jams/ch/reader/view_abstract.aspx?file_no=20030691&flag=1 [14] 马晓群, 吴文玉, 张辉.农业旱涝指标及在江淮地区监测预警中的应用.应用气象学报, 2009, 20(2):186-194. doi: 10.11898/1001-7313.20090208 [15] Steele D D, Sajid A H, Pruuty L D.New corn evapotranspiration crop curves for southeastern North Dakota.Trans of the ASAE, 1996, 39(3):931-936. doi: 10.13031/2013.27578 [16] Sammis T W, Mapel C L, Lugg D G, et a1.Evapotranspiration, crop coefficients predicted using growing-degree-days.Trans of the ASAE, 1985, 28(3):773-780. doi: 10.13031/2013.32336 [17] Stegman E C.Corn crop curve comparison for the Central and Northern Plains of the US.Applied Engineering in Agriculture, 1988, 4(3):226-233. doi: 10.13031/2013.26612 [18] 孙景生, 刘祖贵, 张寄阳, 等.风沙区春小麦作物系数试验研究.农业工程学报, 2002, 18(6):55-58. http://www.cnki.com.cn/Article/CJFDTOTAL-NYGU200206013.htm [19] 张淑杰, 周广胜, 李荣平.基于涡度相关的春玉米逐日作物系数及蒸散模拟.应用气象学报, 2015, 26(6):695-704. doi: 10.11898/1001-7313.20150606 [20] 侯琼, 王英舜, 杨泽龙, 等.内蒙古典型草原作物系数的动态模拟与确定.植物生态学报, 2010, 34(12):1414-1423. http://www.cnki.com.cn/Article/CJFDTOTAL-ZWSB201012009.htm [21] 云文丽, 侯琼, 李建军, 等.基于作物系数与水分生产函数的向日葵产量预测.应用气象学报, 2015, 26(6):705-713. doi: 10.11898/1001-7313.20150607 [22] 山仑, 康绍忠, 吴普特.中国节水农业.北京:中国农业出版社, 2004:150-155. [23] 麻雪艳, 周广胜.春玉米最大叶面积指数的确定方法及其应用.生态学报, 2013, 33(8):2596-2603. http://www.cnki.com.cn/Article/CJFDTOTAL-STXB201308032.htm