近年来,中国玉米单产增长速度表现出明显的放缓趋势,因此研究各驱动因子对玉米单产的影响和贡献,探寻中国玉米单产增速变缓的原因,可为促进中国玉米单产持续增长提供理论支持,为保障国家粮食安全提供有效的途径。本文首先采用集合经验模态分解方法,分析1981—2010年中国玉米单产的变化趋势,发现全国尺度上和24个省(自治区、直辖市)均出现了玉米单产增速变缓的现象。然后利用统计方法和EPIC模型两种方法分析了1981—2010年5个主要的驱动因子(气候、肥料、灌溉、种植面积和品种)对玉米单产增长的影响和贡献。研究得到的主要结论如下:(1)1981—2010年中国玉米单产总体上呈增长趋势,但是随时间变化,单产增加到一定程度出现了增速变缓现象,尤其是在1990年之后。其中31个省(自治区、直辖市)中有24个地区的玉米单产增速出现了逐步变缓现象。(2)统计方法的结果表明,肥料、灌溉、种植面积和降水显著地促进了全国玉米单产的增产,当这4个因子每增加1%时,将分别促进全国玉米单产增长0.27%、0.15%、0.08%和0.09%。而气温和云覆盖率则对全国玉米单产变化产生负面影响,温度每增加1%和辐射每减少1%,将导致全国玉米单产下降0.87%和0.35%。在1981—2010年,肥料的快速增长(148.51%)促使全国玉米单产增产40.10%,综合其他驱动因子的贡献来看,肥料是对全国玉米单产贡献最高的因子,过量施肥会导致肥料利用率显著降低,因此导致玉米单产增速变缓。在区域尺度上,肥料也是两大玉米主产区(北方春玉米区和黄淮海夏玉米区)玉米增产最重要的因子,30a里北方春玉米区和黄淮海夏玉米区肥料平均促使玉米单产分别增长了32.83%和36.77%。相较于肥料,气候变化、玉米种植面积和灌溉的变化对玉米单产增长的贡献较小。受气候变化的影响,全国玉米单产略有降低,其中温度升高导致全国玉米单产下降了2.83%,降水和辐射的减少则导致玉米单产分别下降了0.17%和1.54%。(3)EPIC模型模拟的结果表明,1981—2010年快速增长的肥料投入对玉米单产增长的贡献最高且一直呈逐步增加趋势,但由于肥料的投入量已临近边际效益递减的边缘,贡献的增长速度出现变缓,肥料利用率显著降低,导致玉米单产增速放缓。品种改良对玉米单产增长的贡献仅次于肥料,且在全国尺度上贡献增速出现加快趋势。气候变化和灌溉的变化对全国和北方春玉米区的玉米单产有正面影响,气候变化和灌溉系数的下降不利于黄淮海夏玉米区玉米增产。玉米种植面积的增长整体上对玉米增产有不利影响,相较于其他驱动因子,它的贡献最低。(4)无论是统计方法,还是作物模型方法,都认为肥料是对玉米单产增长贡献最高的驱动因子,过量施肥引起肥料对玉米单产贡献的增速变缓,继而引起玉米单产增速变缓,但统计法下肥料对玉米的贡献总是小于模拟法。对于全国尺度上玉米种植面积的增长和灌溉系数的下降,统计法认为两者分别使玉米单产增长和下降,而EPIC模拟法的结果则完全相反。关于气候变化对玉米的影响,统计方法认为气候变化导致全国和北方春玉米区玉米单产下降,而EPIC模拟结果则认为气候变化促进全国和北方春玉米区玉米单产增长,黄淮海夏玉米的玉米单产下降。EPIC模型模拟法可以研究品种更替对玉米单产的影响,统计方法则无法定量分析品种对玉米单产的影响。由于数据质量、数据种类、模型选择的不同以及简化处理问题等,统计结果还存在一定的不确定性。同样的,空间输入数据、模型的参数化和模型的校准等过程特使作物模型模拟法也存在不确定性。今后将对统计法和EPIC模型模拟的不确定性来源进行研究,通过多个统计模型和作物模型的对比选择最优的模型,更全面和综合地考察更多因子对玉米单产的影响,加强各因子影响的机理方面的研究,提高结果的精准度。