一 主要粮食作物种植区农业气候资源变化

本研究所用数据来自中国气象局气象数据共享数据网，包括1961～2010年601个地面气象台站（剔除数据缺失的站点）的逐日平均气温、最高气温、最低气温、气温日较差、降水量和日照时数资料（见图1）。作物生育期资料来自《中国农业物候图集》（张福春，1987）（见表1）。

图1 研究区气象台站的分布

表1 中国小麦、玉米和水稻的生育期

二 农业气候资源变化对主要粮食作物产量的影响

（一）资料来源

气象资料来自中国气象局气象数据共享数据网，包括1961～2010年601个地面气象台站（剔除数据缺失的站点）的逐日平均气温、气温日较差和降水量资料；作物生育期资料来自《中国农业物候图集》。1961～2008年作物的省级产量资料来自《新中国农业60年统计资料》^[1]，并通过中华人民共和国农业部种植业管理司网站（http：//www.zzys.moa.gov.cn/）收集整理了2009年和2010年的作物产量资料，包括1961～2010年省级冬小麦、玉米、水稻（单季稻和双季稻）的种植面积、单产和总产资料。研究采用的统计分析方法包括气候倾向率（魏凤英，2007）、一阶差分（Nicholls，1997；Lobell et al.，2008）和多元线性回归（魏凤英，2007）等。

（二）研究方法

（1）气候变量和产量的一阶差分（年际变化）计算如下：

Y（k）=X（k+1）-X（k）

式中，Y（k）是一阶差分值，X（k）是年气候变量或产量，k是序列。

（2）为阐明气候变量一阶差分值和产量一阶差分值之间的相关关系，建立平均气温、气温日较差、降水量与产量的多元线性回归方程如下