覆膜条件下对AquaCrop模型冬小麦生长动态和土壤水分模拟效果的评价分析

刘匣，丁奠元，张浩杰，褚晓升，余坤，冯浩,4

（1西北农林科技大学中国旱区节水农业研究院，陕西杨凌 712100；2扬州大学水利与能源动力工程学院，江苏扬州 225009；3西北农林科技大学水利与建筑工程学院，陕西杨凌 712100；4中国科学院水利部水土保持研究所，陕西杨凌 712100）

覆膜条件下对AquaCrop模型冬小麦生长动态和土壤水分模拟效果的评价分析

刘匣1,3，丁奠元2,3，张浩杰1,3，褚晓升1,3，余坤1,3，冯浩1,3,4

（1西北农林科技大学中国旱区节水农业研究院，陕西杨凌 712100；2扬州大学水利与能源动力工程学院，江苏扬州 225009；3西北农林科技大学水利与建筑工程学院，陕西杨凌 712100；4中国科学院水利部水土保持研究所，陕西杨凌 712100）

【目的】通过评价AquaCrop模型对覆膜条件下冬小麦的生长发育、土壤水分、产量以及水分利用效率的模拟效果，为AquaCrop模型在覆膜条件下的校准和应用提供科学的方法和理论依据。【方法】试验设置不覆盖（CK）和白色地膜覆盖（PM）两个处理，于2013年10月至2016年6月年在陕西杨凌进行田间试验，利用2014 —2015年度试验数据对AquaCrop模型进行参数校准，利用2013—2014年度和2015—2016年度的冬小麦观测数据对AquaCrop模型进行验证。【结果】AquaCrop模型较好地模拟了冬小麦冠层覆盖度，冠层覆盖度模拟值和实测值之间的决定系数（R2）为0.86—0.99，均方根误差（RMSE）为2.1%—8.1%。AquaCrop模型也较好地模拟了冬小麦生物量和土壤贮水量，其中地上部生物量的模拟值和实测值之间的R2均大于0.95，RMSE为0.814—1.933 t·hm-2；CK土壤贮水量模拟值和实测值间的相关系数均大于 0.85，PM土壤贮水量模拟值和实测值间的相关系数均大于0.75，CK和PM土壤贮水量模拟值和实测值的均方根误差表现为9.2 mm＜RMSE＜17.6 mm，标准均方根误差（NRMSE）小于5.5%。冬小麦产量实测值和模拟值相对误差（RE）为-4.4%—9.0%，PM产量实测值和模拟值的平均值较CK分别提高40.5%和40.3%，表现出较好的一致性，处理间成显著性差异。水分利用效率实测值和模拟值RE 为-10.4%—-1.5%，PM水分利用效率实测值和模拟值的平均值较CK分别提高54.1%和47.5%，同样表现出较好的一致性，处理间成显著性差异。在冠层覆盖度、地上部生物量、产量和水分利用效率方面，模型模拟值和实测值的变化趋势基本一致，且PM模拟值和实测值间均较CK表现出显著性差异。【结论】AquaCrop模型能够较好地模拟覆膜条件下冬小麦生长发育过程，可以用于覆膜条件下作物生产力的模拟和预测，为AquaCrop模型的推广应用提供了可靠的数据支持。

AquaCrop模型；覆膜；土壤水分；产量；水分利用效率

0 引言

【研究意义】冬小麦是中国主要粮食作物之一，中国小麦总产量和播种面积分别占粮食作物的 21.3% 和 25.1%，而冬小麦则占小麦总产量的90.5%和面积的 86.9%[1]，因此确保冬小麦生产可持续发展对保障中国粮食安全具有十分重要的意义。关中地区降雨季节分布不均匀，主要集中在7—10月份，而冬小麦的生育期为每年的10月中旬到次年的6月[2]，期间降雨量较少，干旱缺水已经成为限制冬小麦产量的关键因素。如何有效地保蓄水分，增强作物抗旱能力，提高作物水分利用效率和产量已经成为保证该区小麦高产稳产的关键。地膜覆盖技术自1978年引入中国以来，在旱区得到大面积推广。覆膜可以有效地抑制土壤水分蒸发，改善作物供水条件，减少水分亏缺，增加土层温度[3]，并提高作物产量和水分利用效率[4]。国内外学者对模型已经展开大量应用研究，通过模型模拟不同农田水分状况下各个指标的变化，为农田水分优化和灌溉制度的设计提供参考依据和技术支持。因此，研究覆膜条件下Aquacrop模型对冬小麦的生长发育的影响对提高旱区作物产量具有重要意义。【前人研究进展】AquaCrop模型是由世界粮农组织（Food and Agriculture Organization，FAO）研究提出的关于作物生长的水分驱动模型。AquaCrop模型数据输入简洁、用户界面友好、数据输出准确[5]，是一款面向用户的简单、直观、便捷的作物生产力模拟、预测工具[6]，具有很好的发展应用前景。大量研究表明，AquaCrop模型相对于其他复杂模型比如CropSyst[7]，WOFOST[8]，CEREC-Rice[9]等具有相似的甚至更好的模拟效果。国内外学者对 AquaCrop模型开展了大量应用研究，主要集中在参数调整和模型验证方面，并证明AquaCrop模型在不同地区对于玉米[6,10]、小麦[11-13]、土豆[14]、番茄[15]等作物的生长有很好的模拟精度，可以辅助农田管理决策。【本研究切入点】国外对于 AquaCrop模型的研究较多[7-10, 12-15]，而中国对于AquaCrop模型的研究尚属起步阶段，主要集中在灌溉水平和施氮水平的参数优化和模型验证[6, 11, 16-17]。模型中自带覆盖模块，但是 AquaCrop模型设计之初对农业地表覆盖和耕作技术等考虑不足，AquaCrop模型中不具备土壤温度模块，而覆膜对土壤温度的影响较大，所以对覆膜栽培作物的模型模拟和验证很少，仅杨宁等[18]对AquaCrop 模型覆膜栽培玉米水分利用过程进行了研究，基于 AquaCrop模型覆膜条件下冬小麦农田水分利用过程的模拟与验证尚属首例。【拟解决的关键问题】本研究以关中地区雨养条件下冬小麦地膜覆盖为研究对象，对模型参数进行本地化，并验证 AquaCrop模型在水分利用和作物生长及产量方面的适用性，以便为AquaCrop模型的改进、应用推广、作物生产力预测、农业资源的综合管理提供一些基础资料。

1 材料与方法

1.1 试验区概况

试验地设在陕西省杨凌西北农林科技大学教育部旱区农业水土工程重点实验室灌溉试验站（108°24′E，34°20′N，海拔521 m）。该地区为关中旱作类型区，属于暖温带半湿润大陆性季风气候区，四季分明，多年平均气温13—15℃，多年平均降雨量在630 mm左右，降雨量年际分布不均，主要集中在7—10月份，无霜期为213 d左右。供试土壤为中壤土（砂粒﹕粉粒﹕黏粒质量比8﹕73﹕19）。1 m 土层的平均田间持水率为 23%—25%，凋萎含水率为8.5%（以上均为质量含水率），土壤 pH平均值为8.2，平均容重为1.37 g·cm-3。该站地下水埋深 5 m以下，其向上补给量可忽略不计。耕层土壤养分含量为有机碳 8.14 g·kg-1，全氮、全磷和全钾分别为0.95、0.83 和20.42 g·kg-1，速效磷、速效钾分别为21、290 mg·kg-1。

1.2 试验设计

本试验设置不覆盖（CK）、白色地膜覆盖（PM）两种处理措施，采用随机区组试验设计，每个处理重复3次，共6个小区。小区方向为东西走向，面积为10 m2（5 m×2 m），每个小区之间有宽0.5 m的保护行。试验取样和水分测定主要集中在小区中间部位，减少边界影响。覆盖处理播前，先进行人工翻耕、整地、施肥，试验所用地膜为普通聚乙烯塑料薄膜，覆膜宽度为30 cm，膜厚0.01 mm。

冬小麦供试品种为小偃22，播种日期分别为2013 年10月9日、2014年10月15日、2015年10月13日，采用条播种植，深度为5—6 cm，行距30 cm，播种量为187.5 kg·hm-2。分别于2014年6月5日，2015 年6月6日和2016年6月5日收获。各个小区统一施肥水平：基肥为 120 kg N·hm-2（尿素）和 100 kg P2O5·hm-2（过磷酸钙）；返青前统一追肥30 kg N·hm-2（尿素）。生育期内不灌溉。

1.3 观测指标

1.3.1 土壤水分和盐分 本试验使用 TRIME-TDR （TRIME-PICO-IPH TDR，德国IMKO公司）测定深层土壤含水率和盐分，测定土层最大深度为260 cm，TRIME管安装在各小区中间部位，自地表以下1 m内每10 cm为一个测量深度间隔，1 m之后每20 cm为一个测量深度间隔，本研究土壤剖面深度选为 120 cm。小麦种植前和收获时，采用打土钻的方法校正TRIME-TDR仪器。

土壤贮水量的计算公式如下：

式中，H为土壤贮水量，mm；Qi为第i层土壤体积含水率，%；hi为第i层土壤厚度，cm；n为测土壤体积含水率时的层序。

玉米生育期耗水量采用土壤水分平衡公式计算：

式中，ET为作物耗水量，mm；△W为收获期与播种期土壤贮水量的变化（按120 cm计），mm；P为作物生育期内的降水量，mm；I为生育期灌溉量，mm，旱作为0；K为时段内地下水的补给量，mm。由于试验田地下水埋深在5 m以下，可视K为0；降水入渗深度不超过2 m，可视深层渗漏为0。

1.3.2 作物生长指标 生物量的测定：在冬小麦各生育期内，每个小区选取生长状况良好、具有代表性的5株小麦，从茎基部切断，获得完整的地上部，105℃杀青30 min，在75℃恒温条件下烘干至质量恒定后称重，记录冬小麦的干物质重量。叶面积指数（LAI）：采用冠层分析仪（SunScan，英国Delta T仪器公司）测定，每个小区在不同位置测定3次取平均值。

通过LAI计算冬小麦冠层覆盖度[11]：

1.3.3 产量 冬小麦成熟后，在各个小区取1 m2的小麦植株样测产，麦穗人工脱粒，105℃杀青30 min，在75℃恒温条件下烘干至质量恒定后称重，换算成单位面积产量（kg·hm-2）。

1.3.4 水分利用效率（WUE） 水分利用效率（WUE）计算公式如下：

式中，WUE为水分利用效率，kg·hm-2·mm-1；Y为单位面积小麦产量，kg·hm-2；ET为耗水量，mm。

1.4 AquaCrop模型数据库的建立

模型输入数据包括：气象数据、作物参数数据、田间管理数据、土壤参数、模型运行初始条件和田间试验数据。

气象数据：由陕西省杨凌西北农林科技大学教育部旱区农业水土工程重点实验室灌溉试验站提供。ET0采用FAO提供的参考作物蒸发蒸腾量计算软件通过最高气温、最低气温、空气相对湿度、风速和日照时数计算。根据降雨量、ET0和温度数据建立模型气象数据文件（图1）。

作物参数数据：根据试验区冬小麦实际生长发育期输入AquaCrop模型生成作物参数数据库文件。用2014—2015年试验区冬小麦观测数据进行模型的参数调试，用2013—2014年和2015—2016年的观测数据进行模型的验证，最终确定 AquaCrop模型作物参数如表1所示（保守参数参照FAO提供的AquaCrop模型参考手册中参数的设定。由耕作方式不同导致覆膜与裸地处理的生育期有两套参数[18-20]，GDD先根据生育期设定，然后在作物生长模块选择 Growing degree-days，模型自动将各阶段的生育期天数转换为积温，再采用试错法调试。水分生产力在模型模拟中根据作物类型的特定设置区间和计算值综合考虑设定，后采用试错法进行调试。参考收获指数首先根据实际计算值进行设定，后在模型模拟中进行调试）。

表1 AquaCrop模型中的作物参数Table 1 Crop parameters in AquaCrop model

模型运行初始条件：根据试验测定数据输入土层深度、土壤水分含量和土壤盐分，生成模型运行初始条件，如表2所示。

田间试验数据：根据试验数据输入冠层覆盖度、生物量和贮水量，生成田间试验数据文件。

表2 不同土层深度下土壤水分含量和土壤盐分初始值Table 2 The initial values of soil moisture and soil salinity at different depths

1.5 模型的验证与评价方法

通过比较覆膜与裸地冬小麦冠层覆盖度、生物量和土壤贮水量等参数的模拟值和实测值，来评价模型模拟的效果。验证指标为相对误差（relative error，RE）、均方根误差（root mean square error，RMSE）、标准均方根误差（normalized root mean square error，NRMSE）和决定系数（goodness-of-fit，R2）。式中，和n分别为观测值、模拟值、观测值的平均值，模拟值的平均值和实测值与观测值的个数。RE越小表示模拟精度越高；RMSE越小表明模拟精度越高；NRMSE值越接近0表示模拟精度越高，当NRMSE小于10%，模拟效果很好；当NRMSE处于10%和20%之间，表明模拟效果好；当NRMSE处于20%和30%之间，表明模拟效果合理；当NRMSE大于30%，表明模拟效果差[21]。R2越接近于1表明模型的解释程度越高。

2 结果

2.1 冬小麦冠层覆盖度模拟结果

图1 2013—2016年小麦季降雨量、最高气温、最低气温和参考作物蒸发蒸腾量Fig. 1 Rainfall, maximum temperature, minimum temperature, and reference evapotranspiration in winter wheat seasons from 2013 to 2016

图2 2013—2016年覆膜与裸地处理下冬小麦冠层覆盖度的模拟结果Fig. 2 Simulated values of winter wheat canopy cover under the plastic mulching and non-mulching treatments from 2013-2016

叶面积指数是研究冬小麦生物量和产量形成的重要指标，模型中用冠层覆盖度代替叶面积指数可以较为清晰的表现出作物的生长状况，2013—2016年3季冬小麦的冠层覆盖度的发育动态与模型的评价指标如图2所示，在苗期和越冬期，作物生长缓慢，叶面积指数较小（因覆盖小，冠层分析仪无法准确测定），冠层覆盖度相对增长缓慢，拔节期以后叶面积指数增大，冠层覆盖度迅速增加，并于抽雄期达到最大值，之后随着生育期的变化部分叶片开始衰老，冠层覆盖度也随之下降。从图中可以看出，PM和CK在冬小麦生育期内模拟值与实测值的总体变化趋势基本一致。PM较CK生育期提前，增长较快，且最大冠层覆盖度明显高于CK，CK的模拟值与实测值的RMSE 为2.1%—8.1%，而PM的模拟值与实测值的RMSE为2.4%—2.7%，冠层覆盖度方面的模拟PM较CK模拟精度较高。总体来讲，2014—2015年（调参）相关系数R2均大于0.95，2013—2014和2015—2016年（验证）相关系数R2均大于0.85。2013—2016年模拟值和实测值之间的相关关系表现为 2.1%＜RMSE＜ 8.1%，2.9%＜NRMSE＜11.9%，可见，模型在冠层覆盖度方面的模拟精度较高。因此，AquaCrop模型可以较好模拟冬小麦冠层覆盖度的动态变化过程。

在教授的过程中，教师通过解决问题来解释数学思想，数学应用方法解决了数学问题.老师反映思维活动，学生掌握了以自己的逻辑思维生活的能力，发现数学的乐趣，找到解决问题的最佳解决方案.教学问题解决过程通过分析向学生讲授课程的方法.通过不断的学习，总结和理解数学思维方法.在高中数学教育过程中，教师需要通过教授学生的知识和技能，有效地提高学生的学习效率，并将数学思想无形地嵌入学生中.

2.2 冬小麦生物量模拟结果

地上部生物量是决定冬小麦产量的重要因素之一，图3比较了3年试验中地膜覆盖和裸地条件下冬小麦生育期内地上部生物量模拟值和实测值随时间的变化。不同处理冬小麦地上部生物量模拟值和实测值变化趋势基本一致，但是在冬小麦生育期的苗期和越冬期，模型的模拟值都普遍低于实测值，可能是冬小麦生长发育阶段（尤其是越冬期）部分小麦苗死亡，根据出苗率计算生物量，导致生物量计算值偏大，而模型在苗期和越冬期冬小麦生长速率较慢，导致模拟值普遍低于实测值。CK的模拟值与实测值的 RMSE 为0.814—1.663 t·hm-2，而PM的模拟值与实测值的RMSE为1.487—1.933 t·hm-2，在地上部生物量的模拟方面CK较PM模拟精度较高。总体来说，PM和CK模拟值和实测值相关系数R2均大于0.95，很接近于1，说明模型的解释程度很好。2013—2014年PM和CK 的NRMSE均小于20%，说明模型的模拟好；2014—2015年CK的NRMSE大于30%，说明模型模拟的不好，PM的NRMSE为22.3%；2015—2016年PM和 CK的20%＜NRMSE＜30%，说明模拟值可以合理解释实测值的趋势变化。因此，AquaCrop模型能够较好模拟地上部生物量的动态变化规律。

2.3 冬小麦土壤贮水量模拟结果

图4 2013—2016年覆膜与裸地处理下1.2 m土层贮水量的模拟结果Fig. 4 Simulated values of winter soil water content in 1.2 m depth under the plastic mulching and non-mulching treatments from 2013-2016

2013—2016年冬小麦PM和CK的1.2 m土壤贮水量的模拟值和实测值的相关关系如图4所示，作物生育期内未设置灌水，降雨导致了模拟值和实测值中土壤贮水量的增加。CK模拟值与实测值的相关系数均大于0.85，PM模拟值和实测值的相关系数均大于0.75，CK的平均相关系数R2比PM高，2015—2016 年CK的相关系数较低，这是由于人为因素在越冬期1月9日进行了灌水，实测的贮水量出现了增加的趋势，模型模拟灌溉参数没有发生改变，使得模拟值偏低，从而CK的相关系数R2相对于PM较低。在模型中设置2016年1月9日灌水30 mm，增加灌水参数后，实测值与模拟值的R2为0.95，RMSE为11.4，NRMSE为3.6。整体而言，PM的模拟值和实测值的模拟精度较 CK差，原因可能是模型中覆盖模块的设计并不是很完善，只是考虑了土壤蒸发，对降雨入渗没有进行设定，本文通过改变径流曲线可能不是十分完善，而且模块中没有考虑地膜覆盖材料随时间的变化，覆盖材料的完整程度对土壤蒸发和降雨入渗影响较大，导致 PM模拟值与实测值的差异较大。统计数据显示，PM和CK的RMSE在9.2至17.6之间，NRMSE＜5.5%，模拟值与实测值之间的误差较小，AquaCrop 模型可以较准确的模拟冬小麦PM和CK土壤水分的动态变化。

2.4 产量和水分利用效率

耗水量、生育期内降雨量和作物生长发育状况密切相关，同一年中冬小麦生育期内降雨量不同是由PM 和 CK的生育期天数和降雨的时空分布不均导致。表2给出了3年试验PM和CK产量、耗水量和水分利用效率实测值和模拟值的相对误差（RE）。试验数据表明，各处理产量实测值为 3 728.9—7 400.9 kg·hm-2，模拟值为4 098—7 174 kg·hm-2，二者的RE为-4.4%—9.0%。PM和CK实测产量呈显著性差异，3年实测值和模拟值的平均值，PM分别较CK增加 40.5%和 40.3%；2014—2015年和 2015—2016年耗水量RE为5.4%—12.9%，模拟值普遍高于实测值。类似的，PM和CK水分利用效率实测值为18.0—29.5 kg·m-3，模拟值为16.0—26.8 kg·m-3，二者的RE为-10.4%—-1.5%。2014—2015年和2015—2016 年PM和CK实测水分利用效率的平均值分别为18.1 kg·m-3和27.9 kg·m-3，模拟水分利用效率的平均值为17.3 kg·m-3和25.5 kg·m-3，模拟值略低于实测值。2014 —2015年PM水分利用效率的实测值和模拟值分别较 CK显著提高44.5%和 44.3%（P＜0.05），2015 —2016年PM水分利用效率的实测值和模拟值分别较CK显著提高63.9%和50.6%（P＜0.05）。总体来讲，AquaCrop模型能够较好的模拟作物产量和水分利用效率。

表3 2013—2016年地膜覆盖与裸地冬小麦产量和水分利用效率的模拟值和实测值Table 3 Measured and simulated values for plastic mulching and non-mulching winter wheat yield and water use efficiency from 2013-2016

3 讨论

3.1 参数的设定

作物参数的确定对于作物模型的应用起着重要作用[22]。AquaCrop模型中本身提供了典型作物的一系列生理参数，但是对不同的品种、种植区域、田间管理措施，需要利用实测数据率定模型参数，将模型本地化[23]。大量研究表明，覆膜可以明显增加表层地温[3,24-25]，明显降低土壤蒸发并促进作物蒸腾[18]，还可以提高作物产量和水分利用效率[26]。但AquaCrop模型是关于作物生长的水分驱动模型，在模型参数输入中仅有气温参数模块，而没有地温模块，对于模型中自带的地膜覆盖模块，也只考虑了地膜对土壤蒸发的影响，而忽略了地膜对土壤温度的影响，进而忽略了地膜的增温效应对作物生长发育的影响。基于考虑地膜的增温效应对作物生育期天数的影响，将地膜和裸地涉及到的地温增加的参数根据记录的物候期设定为两套参数值，并率定模型参数，得到地膜和裸地的部分作物参数（表1）。

3.2 冠层覆盖度和生物量的模拟和验证

本研究表明AquaCrop模型可以模拟覆膜条件下冬小麦的生长发育过程，并具有较高的模拟精度，这与杨宁等[18]研究覆膜玉米的结果类似。本研究对冠层覆盖度的模拟结果较好，决定系数R2均达到0.85以上，均方根误差RMSE值在2.1%—8.1%之间，这与以往的研究结果类似。比如李会等[27]研究的AquaCrop模型对不同灌溉水平处理冠层覆盖度的模拟效果非常好，决定系数在0.95以上，RMSE在2% —6%之间，而滕晓伟等[11]的研究结果为R2和RMSE分别为0.464和8.0%。EBRAHIM[28]研究结果为0.66 ＜R2＜0.98，9%＜RMSE＜45%。本研究结果表明AquaCrop模型在覆膜条件下冠层覆盖度的模拟各项指标达到了较为理想的结果，AquaCrop模型可以较好的模拟冬小麦冠层覆盖度的动态变化过程。但是试验中冠层覆盖度的测定值较少，以后试验中应增加试验数据测定频率。

地上部生物量模拟（图3）较高的相关系数（R2＞0.95）和较低的均方根误差（0.814 t·hm-2＜RMSE ＜1.933 t·hm-2）说明AquaCrop模型在覆膜条件下能够较好的模拟生物量，但是模拟结果明显高于付驰等[23]（0.036 t·hm-2＜RMSE＜0.190 t·hm-2）和张铁楠等[29]（RMSE为0.0959 t·hm-2）对春小麦地上部生物量的模拟。与滕晓伟等[11]（RMSE为1.622 t·hm-2）、HSIAO等[30]（0.58 t·hm-2＜RMSE＜6.18 t·hm-2）和MEBANE等[31]（RMSE为1.33 t·hm-2和1.53 t·hm-2）研究AquaCrop模型在地上部生物量上的结果类似。

本研究对于 AquaCrop模型在地上部生物量的研究结果显示AquaCrop模型通常高估地上部生物量（图3），但是在拔节期后的生育期内出现前期高估后期低估地上部生物量的现象（图3-a1，a2，b2），这种研究结果与KATERJI等[14]和AHMADI等[32]的研究结果一致。原因可能是冠层覆盖度很大，作物根系较为发达，后期水分亏缺会影响作物生物量的形成（图3-b1），本文中2014—2015年和2015—2016年相对于2013—2014年降雨量减少，抽雄期后作物耗水量（作物蒸腾和土壤蒸发）增加，而水分供给不足，易影响作物生物量的形成，从而造成生物量降低，而2014—2015年和2015—2016年覆膜能够蓄水保墒，可减小水分供给不足对作物产量的影响，这在2014—2015年尤为明显。而且Aqua Crop模型缺乏更复杂的生理子模块，所以不能很好地解释水分胁迫对地上部生物量积累及光合产物向籽粒转移的影响[33]。

综上所述，国内外对于冠层覆盖度和地上部生物量研究结果大体一致，但模拟值和实测值的大小关系不能得到一致结论，原因可能与测量数据的质量、数据处理的方法、校准模型参数的准确性以及一些外界因素如灌溉条件、作物生长环境、耕作条件等因素影响有关。模型本身需要改进的问题同时也影响模型模拟的准确性，AquaCrop模型本身简化了作物衰老过程，并没有考虑作物生长过程中的复杂性[32]，VANUYTRECHT等[34]、HSIAO等[30]和ZELEKE等[35]研究表明AquaCrop模型在模拟冠层衰老阶段不够准确。

3.3 土壤贮水量的模拟和验证

AquaCrop模型作为水分驱动模型，其对产量和生物量的模拟性能取决于精确的土壤水分动态模拟[32]。本研究结果表明 AquaCrop模型能够较好地模拟覆膜和裸地土壤水分贮水量的变化，相关系数达到0.75以上，且均方根误差为 9.2—17.6 mm。模拟精度高于MONTOYA等[36]和李子忠等[16]对土壤贮水量的模拟。AquaCrop模型模拟值较实测值偏高（图 4），这与HSIAO等[30]，FARAHANI等[37]的研究结果一致，但是IGBAL等[38]和MKHABELA等[39]对于小麦贮水量和 MEBANE等[31]对玉米的贮水量的模拟值却低于实测值。对此，国内外学者并不能得到一致结论。原因可能是试验地土壤的结构不同，而且模型并不能考虑土壤结构的变异性，因此，相同的水分数据在模型模拟中有很多的不确定性[35]。

图 2—4中可以看出覆膜的模拟精度低于裸地处理，原因可能是模型设计之初对农业地面覆盖模块的设计考虑不足，模型中覆膜模块只是考虑覆膜对土壤蒸发的影响[18]，未考虑对降雨入渗的影响（本文中通过改变CN曲线来区别地膜覆盖和裸地处理对降雨径流的影响），而且覆盖模块中覆盖度等参数是固定值，实际中覆盖度（覆盖材料的完整度）会随着时间的变化而发生改变。

3.4 产量和水分利用效率的模拟和验证

应用AquaCrop模型模拟冬小麦生产力形成过程，主要考虑作物产量和水分利用效率两个方面。本文在产量和水分利用效率方面模型模拟覆膜和裸地冬小麦相对误差较小，说明 AquaCrop模型能够较为准确的模拟覆膜冬小麦产量和水分利用形成过程。这与以往的研究结果一致[18]。且实测覆膜与裸地处理产量和水分利用效率均成显著性差异，这与模拟值中覆膜与裸地的差异基本一致。但是水分利用效率的模拟值均低于实测值，原因为作物生育期内耗水量的模拟值均高于实测值，这是生育后期模型没有考虑作物生长的复杂性，从而简化作物衰老的过程，使得后期土壤贮水量模拟不够准确。

本文仅应用一个作物品种在一个研究区域对Aquacrop模型的覆膜模块进行模拟和验证，在今后的工作中，还需要针对不同的作物品种和不同的研究区域进行模拟和验证，从而完善Aquacrop模型在覆膜模块中的应用。而且模型本身在设计之初对覆膜模块的考虑不足，未考虑覆膜增温效应，而且覆膜模块只考虑了覆膜对土壤蒸发的影响，关于地膜对降雨入渗的阻滞作用未涉及，这需要在以后的模型改进中加以完善。同时模型简化了作物衰老，没有考虑作物生长的复杂性，而且忽略了覆盖材料的变化，应该提高并细化时间参数和作物参数，从而提高模拟精度。在以后试验中要增加田间实测数据频率，充分掌握各类资料后，完善模型在覆膜冬小麦生长上的应用，以期为作物生产力模拟和预测以及农艺措施的优化提供一些理论依据和技术支持。

4 结论

利用3年冬小麦资料对AquaCrop模型在覆膜条件下进行模拟和验证，根据实测数据对模型参数进行率定，并模拟了覆膜和裸地2种覆盖方式对旱作冬小麦冠层覆盖度、地上部生物量、土壤贮水量的影响，计算分析了2013—2016年冬小麦的产量和水分利用效率。结果表明 AquaCrop模型在关中地区对覆膜冬小麦的生长发育和产量形成，具有较高的模拟精度，可以应用 AquaCrop模型对地膜覆盖的生产力进行简单的预测和初步的评估。

[1] 孙爽, 杨晓光, 李克南, 赵锦, 叶清, 解文娟, 董朝阳, 刘欢. 中国冬小麦需水量时空特征分析. 农业工程学报, 2013, 29(15): 72-82. SUN S, YANG X G, LI K N, ZHAO J, YE Q, XIE W J, DONG C Y, LIU H. Analysis of spatial and temporal characteristics of water requirement of winter wheat in China. Transactions of the Chinese Society of Agricultural Engineering, 2013, 29(15): 72-82. (in Chinese )

[2] 姚宁, 宋利兵, 刘健, 冯浩, 吴淑芳, 何建强. 不同生长阶段水分胁迫对旱区冬小麦生长发育和产量的影响. 中国农业科学, 2015, 48(12): 2379-2389. YAO N, SONG L B, LIU J, FENG H, WU S F, HE J Q. Effects of water stress at different growth stages on the development and yields of winter wheat in arid region. Scientia Agricultura Sinica, 2015, 48(12): 2379-2389. (in Chinese)

[3] ZHOU L M, LI F M, JIN S L, SONG Y J. How to ridges and the furrow mulched with plastic film affect soil water, soil temperature and yield of maize on the semiarid Loess Plateau of China. Field Crops Research, 2009, 113(1): 41-47.

[4] 李巧珍, 李玉中, 郭家选, 刘晓英, 徐春英. 覆膜集雨与限量补灌对土壤水分及冬小麦产量的影响. 农业工程学报, 2010, 26(2): 25-30. LI Q Z, LI Y Z, GUO J X, LIU X Y, XU C Y. Effects of field rainwater harvesting by plastic mulch and complement irrigation on soil water and yield of winter wheat. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(2): 25-30. (in Chinese)

[5] 朱秀芳, 李宜展, 潘耀忠, 史培军. AquaCrop作物模型研究和应用进展. 中国农学通报, 2014, 30(8): 270-278. ZHU X F, LI Y Z, PAN Y Z, SHI P J. A review on the research and application of Aqua Crop model. Chinese Agricultural Science Bulletin, 2014, 30(8): 270-278. (in Chinese)

[6] 项艳. AquaCrop 模型在华北地区夏玉米生产中的应用研究[D]. 泰安: 山东农业大学, 2009. XIANG Y. AquaCrop model application of summer maize planting in North China[D]. Tai’an: Shandong Agricultural University, 2009. (in Chinese)

[7] SAAB M T A, TODOROVIC M, ALBRIZIO R. Comparing AquaCrop and CropSyst models in simulating barley growth and yield under different water and nitrogen regimes. Does calibration year influence the performance of crop growth models? Agricultural Water Management, 2015, 147(1): 21-33.

[8] TODOROVIC M, ALBRIZIO R, ZIROTIC L, SAAB M T A, STOCKLE C, STEDUTO P. Assessment of AquaCrop, Cropsyst, and WOFOST models in the simulation of sunflower growth under different water regimes. Agronomy Journal, 2009, 101(3): 509-521.

[9] AMIRI E, REZAEI M, REZAEI E E, BANNAYAN M. Evaluation of Ceres-rice, AquaCrop and Oryza 2000 models in simulation of rice yield response to different irrigation and nitrogen management strategies. Journal of Plant Nutrition, 2014, 37(11): 1749-1769.

[10] PAREDES P, MELO-ABREU J P D, ALVES I, PEREIRA L S. Assessing the performance of the FAO AquaCrop model to estimate maize yields and water use under full and deficit irrigation with focus on model parameterization. Agricultural Water Management, 2014, 144(3): 81-97.

[11] 滕晓伟, 董燕生, 沈家晓, 孟鲁闽, 冯海宽. AquaCrop 模型对旱区冬小麦抗旱灌溉的模拟研究. 中国农业科学, 2015, 48(20): 4100-4110. TENG X W, DONG Y S, SHEN J X, MENG L M, FENG H K. Winter wheat irrigation simulation in arid area based on AquaCrop model. Scientia Agricultura Sinica, 2015, 48(20): 4100-4110. (in Chinese )

[12] KUMAR P, SARANGI A, SINGH D K, PARIHAR S S. Evaluation of AquaCrop model in predicting wheat yield and water productivity under irrigated saline regimes. Irrigation & Drainage, 2014, 63(4): 474-487.

[13] ZHANG W, LIU W, XUE Q, CHEN J, HAN X. Evaluation of the AquaCrop model for simulating yield response of winter wheat to water on the southern Loess Plateau of China. Water Science & Technology A Journal of the International Association on Water Pollution Research, 2013, 68(4): 821-828.

[14] KATERJI N, CAMPI P, MASTRORILLI M. Productivity,evapotranspiration, and water use efficiency of corn and tomato crops simulated by AquaCrop under contrasting water stress conditions in the Mediterranean region. Agricultural Water Management, 2013, 130(4): 14-26.

[15] DOMÍNGUEZ A, TARJUELO J M, JUAN J A D, LÓPEZ-MATA E, BREIDY J, KARAM F. Deficit irrigation under water stress and salinity conditions: The mopeco-salt mode. Agricultural Water Management, 2011, 98(8): 1451-1461.

[16] 李子忠, 徐洋, 卢宪菊, 胡克林, 江丽华, 徐钰. AquaCrop模型在大葱生物量和土壤贮水量模拟中的应用和验证. 中国农业大学学报, 2011, 16(4): 59-66. LI Z Z, XU Y, LU X J, HU K L, JIANG L H, XU Y. Evaluation of the Aqua Crop model for simulating biomass for Chinese green onion and soil water storage. Journal of China Agricultural University, 2011, 16(4): 59-66. (in Chinese)

[17] 李玥, 牛俊义, 郭丽琢, 高珍妮, 孙小花. AquaCrop 模型在西北胡麻生物量及产量模拟中的应用和验证. 中国生态农业学报, 2014, 22(1): 93-103. LI Y, NIU J Y, GUO L Z, GAO Z N, SUN X H. Application and validation of AquaCrop model in simulating biomass and yield of oil flax in Northwest China. Chinese Journal of Eco-Agriculture, 2014, 22(1): 93-103. (in Chinese)

[18] 杨宁, 孙占祥, 张立桢, 郑家明, 冯良山, 李开宇, 张哲, 冯晨. 基于改进AquaCrop 模型的覆膜栽培玉米水分利用过程模拟与验证. ,农业工程学报, 2015, 31(1): 122-132. YANG N, SUN Z X, ZHANG L Z, ZHENG J M, FENG L S, LI K Y, ZHANG Z, FENG C. Simulation of water use process by film mulched cultivated maize based on improved AquaCrop model and its verification. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(1): 122-132. (in Chinese)

[19] 王虎全, 韩思明, 唐拴虎, 李岗, 谢惠民. 渭北旱原冬小麦全程地膜覆盖超高产栽培技术研究. 干旱地区农业研究, 1998, 16(1): 24-30. WANG H Q, HAN S M, TANG S H, LI G, XIE H M. Super high yield cultivation technique of winter wheat film mulched for whole growing process in Weibei dryland. Agricultural Research in the Arid Areas, 1998, 16(1): 24-30. (in Chinese)

[20] 方彦杰, 黄高宝, 李玲玲, 汪佳. 旱地全膜双垄沟播玉米生长发育动态及产量形成规律研究. 干旱地区农业研究, 2010, 28(4): 128-134. FANG Y J, HUANG G B, LI L L, WANG J. Yield and growth dynamics of rainfed maize in the system of completely mulched alternating narrow and wide ridges with furrow planting. Agricultural Research in the Arid Areas, 2010, 28(4): 128-134. (in Chinese)

[21] JAMIESON P D, PORTER J R, WILSON D R. A test of computer simulation model ARC-WHEAT1 on wheat crops grown in New Zealand. Field Crops Research, 1991, 27(4): 337-350.

[22] 张艳红, 马永良, 廖树华. CERES-maize 模拟模型中品种参数优化方法研究. 中国农业大学学报, 2004, 9(4): 24-29. ZHANG Y H, MA Y L, LIAO S H. Method of optimizing maize variety parameters in the CERES-Maize simulation model. Journal of China Agricultural University, 2004, 9(4): 24-29. (in Chinese)

[23] 付驰, 李双双, 李晶, 王泳超, 芦玉双, 许为政, 魏湜. AquaCrop作物模型在松嫩平原春麦区的校正和验证. 灌溉排水学报, 2012, 31(5): 99-102. FU C, LI S S, LI J, WANG Y C, LU Y S, XU W Z, WEI S. Calibration and validation of AquaCrop model in spring wheat region of Songnen Plain. Journal of Irrigation and Drainage, 2012, 31(5): 99-102. (in Chinese)

[24] BU L D, LIU J L, ZHU L, LUO S S, CHEN X P, LI S Q, ROBERT L H, ZHAO Y. The effects of mulching on maize growth, yield and water use in a semi-arid region. Agricultural Water Management, 2013, 123(10): 71-78.

[25] 吴从林, 黄介生, 沈荣开. 地膜覆盖在冬小麦全生育期内增温保墒作用的试验研究. 中国农村水利水电, 2001(8): 7-9. WU C L, HUANG J S, SHEN R K. Experimental research on effects of mulch on soil’s temperature and moisture content during the whole growth period of winter wheat. China Rural Water and Hydropower, 2001(8): 7-9. (in Chinese )

[26] 王敏, 王海霞, 韩清芳, 李荣, 张睿, 贾志宽, 杨宝平. 不同材料覆盖的土壤水温效应及对玉米生长的影响. 作物学报, 2011, 37(7) : 1249-1258．WANG M, WANG H X, HAN Q F, LI R, ZHANG R, JIA Z K, YANG B P. Effects of different mulching materials on soil water, temperature, and corn growth. Acta Agronomica Sinica, 2011, 37(7): 1249-1258. (in Chinese)

[27] 李会, 刘钰, 蔡甲冰, 毛晓敏. AquaCrop模型的适用性及应用初探.灌溉排水学报, 2011, 30(3): 28-33. LI H, LIU Y, CAI J B, MAO X M. The applicability and application of AquaCrop model. Journal of Irrigation and Drainage, 2011, 30(3): 28-33. (in Chinese)

[28] EBRAHIM A. Calicration and testing of the AquaCrop model for rice under water and nitrogen management. Communications in Soil Science and Plant Analysis, 2016, 47(3): 387-403.

[29] 张铁楠, 付驰, 李晶, 顾万荣, 许为政, 芦玉双, 魏湜. 基于寒地春小麦 AquaCrop与 WOFOST模型适应性验证分析. 作物杂志,2013(3): 121-126. ZHANG T N, FU C, LI J, GU W R, XU W Z, LU Y S, WEI S. The adaptability test analysis of AquaCrop and WOFOST model based on the cold spring wheat. Crops, 2013(3): 121-126. (in Chinese)

[30] HSIAO T C, LEE H, STEDUTO P, ROJASLARA B, RAES D, FERERES E. Aquacrop—the FAO crop model to simulate yield response to water: III. Parameterization and testing for maize. Agronomy Journal, 2009, 101(3): 448-459.

[31] MEBANE V J, DAY R L, HAMLETT J M, WASTON J E, ROTH G W. Validating the FAO AquaCrop model for rain-fed maize in Pennsylvani. Agronomy Journal, 2013, 105(2): 419-427.

[32] AHMADI S H, MOSALLAEEPOUR E, KAMGAR-HAGHIGHI A A, SEPASKHAH A R. Modeling maize yield and soil water content with Aquacrop under full and deficit irrigation managements. Water Resources Management, 2015, 29(8): 2837-2853.

[33] 张万红, 刘文兆, 王芸. 基于水分驱动的 AquaCrop模型及其研究进展. 干旱地区农业研究, 2014(4): 96-101. ZHANG W H, LIU W Z, WANG Y. AquaCrop model based on water-driven principle and its research progress. Agricultural Research in the Arid Areas, 2014(4): 96-101. (in Chinese)

[34] VANUYTRECHT E, RAES D, WILLEMS P. Global sensitivity analysis of yield output from the water productivity model. Environmental Modelling & Software, 2014, 51(1): 232-332.

[35] ZELEKE K T, LUCKETT D, COELEY R. Calibration and testing of the FAO AquaCrop model for canola. Agronomy Journal, 2011, 103(6): 1610-1618.

[36] MONTOYA F, CAMARGO D, ORTEGA J F, CÓRCOLES J, DOMÍNGUEZ A. Evaluation of Aquacrop model for a potato crop under different irrigation conditions. Agricultural Water Management, 2016, 164(2): 267-280.

[37] FARAHANI H J, IZZI G, OWEIS T Y. Parameterization and evaluation of the AquaCrop model for full and deficit irrigated cotton. Agronomy Journal, 2009, 101(3): 469-476.

[38] IGBAL M A, SHEN Y, STRICEVIC R, PEI H, SUN H, AMIRI E. Evaluation of the FAO AquaCrop model for winter wheat on the North China Plain under deficit irrigation from field experiment to regional yield simulation. Agricultural Water Management, 2014, 135(2): 61-72.

[39] MKHABELA M S, BULLOCK P R. Performance of the FAO AquaCrop model for wheat grain and soil moisture simulation in western Canada. Agricultural Water Management, 2012, 110(7): 16-24.

（责任编辑 杨鑫浩）

Evaluation Analysis of AquaCrop Model in Modeling Winter Wheat Growing Development and Soil Moisture Under Plastic Mulching

LIU Xia1,3, DING DianYuan2,3, ZHANG HaoJie1,3, CHU XiaoSheng1,3, YU Kun1,3, FENG Hao1,3,4
(1Institute of Water-Saving Agriculture in Arid Areas of China, Northwest A&F University, Yangling 712100, Shaanxi;2School of Hydraulic Energy and Power Engineering, Yangzhou University, Yangzhou 225009, Jiangsu;3College of Water Resources and Architectural Engineering, Northwest A&F University, Yangling 712100, Shaanxi;4Institute of Water and Soil Conservation, Chinese Academy of Sciences and Ministry of Water Resource, Yangling 712100, Shaanxi)

【Objective】This study was conducted to evaluate the accuracy of the AquaCrop model in the simulation of winter wheat growth, soil moisture, yield, and water use efficiency under plastic mulching, which can provide a theoretical basis and scientific method for the calibration of the AquaCrop model under plastic mulching.【Method】The experiment was conducted atYangling, Shaanxi, from 2013 to 2016 including the flat planting under plastic mulching (PM) and a control treatment without mulching (CK). The AquaCrop model was calibrated using the experiment data in 2014-2015 and was validated using the data in 2013-2014 and 2015-2016. 【Result】The determination coefficient (R2) of the simulated and the measured canopy cover was between 0.86 and 0.99. The root mean square error (RMSE) of the simulated and the measured canopy cover was between 2.1% and 8.1%, indicating AquaCrop model a good simulation for canopy cover. The R2of the simulated and the measured biomass was greater than 0.95. Meanwhile, the RMSE of the simulated and the measured biomass was between 0.814 and 1.933 t·hm-2. The R2of the simulated and the measured soil water content in CK was greater than 0.85, and the R2of the simulated and the measured soil water content in PM was greater than 0.75. The RMSE of the simulated and the measured soil water content in both CK and PM was between 9.2 and 17.6 mm. The normalized root mean square error (NRMSE) of the simulated and the measured soil water content in both CK and PM was lower than 5.5. Moreover, the RE of the simulated and the measured yield in both CK and PM was from -4.4% to 9.0%. The simulated and the measured yield in PM increased 40.5% and 40.3% compared to that in CK, respectively, and there was a significant difference between CK and PM. The RE of the simulated and the measured water use efficiency in both CK and PM was from -10.4% to -1.5%. The simulated and the measured water use efficiency in PM increased by 54.1% and 47.5% compared to that in CK, respectively, and there was also a significant difference between CK and PM. The results showed that the simulated and measured values had the similar trend in wheat canopy coverage, biomass, yield, and water use efficiency. Those indicated a good performance of the AquaCrop model in modeling plastic mulching treatment. 【Conclusion】 The AquaCrop model can be used to model the winter wheat growing development and productivity under plastic mulching. This study provides a scientific method to calibrate the AquaCrop model and a good data support for the application and development of the AquaCrop model.

AquaCrop model; plastic mulching; soil water; yield; water use efficiency

2016-09-02；接受日期：2017-03-08

国家高技术研究计划（“863”计划）（2013AA102904）、高等学校学科创新引智计划（“111”计划）（B12007）

联系方式：刘匣，E-mail：liux18792966995@163.com。通信作者冯浩，E-mail：nercwsi@vip.sina.com