关键词:
作物模型
品种基因型参数
参数率定
抽样方法
区域差异性
摘要:
马铃薯是具有高营养价值的优质农作物,在保障粮食安全中发挥重要作用,我国是全球马铃薯种植第一大国,面对供需不平衡和气候变化的潜在不利影响,需要采取提高产量和适应环境的对应战略。作物生长模型服务于农业现代化,可以在作物生长发育模拟、农田精准管理、气候演变分析等领域进行精准、大范围监测和预测,但作物模型研究多始于田块尺度,在实际模拟中需要考虑不同环境条件下参数的异质性以打破空间应用限制,尤其是作物基因型参数的变异,需要通过田间试验进行参数本地化工作,也称参数率定。田间试验费时费力,率定参数所需的样本难以获得,所以如何通过率定基因型参数来进一步提升模型的精度成为DSSAT马铃薯模型应用中的热点问题,也对指导农业生产、保障粮食安全、实现农业现代化具有重要意义。本研究旨在为解决作物模型基因型参数本地化率定缺乏样本的问题,开展适合不同样本量进行基因型参数率定的研究,选取东北地区十个试验站点,探索用不同的抽样方法对DSSAT-SUBSTOR模型的马铃薯基因型参数进行率定,基于此分析同一站点不同抽样方法下基因型参数率定结果的差异、同一抽样方法不同站点基因型参数率定结果的差异和不同站点不同抽样方法基因型参数率定结果的差异。主要得到以下结论:第一,改变抽样方法可以提升基因型参数率定的准确性。同一站点不同抽样方法基因型参数率定结果显示三种抽样方法的验证精度从高到低排序为交叉验证法>自助法>留出法,即单站点小样本的参数率定工作中交叉验证法为最优抽样方法。不同站点不同抽样方法基因型参数率定结果显示三种抽样方法的验证精度从高到低排序为留出法>自助法>交叉验证法,即多站点下随机抽样进行参数率定表现最优的是留出法。第二,马铃薯基因型参数分布呈现区域差异性。同一抽样方法下DSSAT-SUBSTOR模型的5个基因型参数率定结果在不同站点均存在差异,其中块茎生长胁迫系数(PD)、光周期系数(P2)、块茎开始生长的临界温度上限(TC)均呈现出0.01水平的显著性差异,说明基因型参数对产量的响应变化与环境有关。具体对比站点间参数率定的差异从大到小排序为光周期系数(P2)>块茎开始生长的临界温度上限(TC)>块茎生长胁迫系数(PD)>叶面积扩展率(G2)>潜在块茎生长率(G3),说明在本研究选取的东北作区影响马铃薯生长发育产生最大差异的因素是太阳辐射和日长。第三,基因型参数对产量的响应受环境影响。同一抽样方法下,对比用单站点和多站点率定的最优参数组进行模拟的结果,单站点率定的最优参数组进行模拟的结果优于多站点,即单一环境下率定的基因型参数比不同环境间率定的基因型参数更接近基因型参数的真值,说明环境间的产量变化往往大于品种参数变化带来的产量变化。因此同一作物同站点的不同品种是否可以用同一组品种基因型参数以简化DSSAT模型实际应用进行参数本地化的率定工作,未来可进一步研究与讨论。第四,样本量不同,适用的抽样方法也不同。同一抽样方法下,对比用单站点和多站点率定的最优参数组进行模拟的结果,多站点的马铃薯基因型参数率定结果中三种抽样方法间的差异较小,单站点的马铃薯基因型参数率定结果中交叉验证法的参数率定表现最优。说明样本量不同,适用的抽样方法也不同,且样本数量越多,抽样方法引起的结果差异越小,小样本适合用交叉验证法,大样本受抽样方法的影响较小,可根据实际需求选取。