中国玉米生育期变化及其影响因子

中国玉米生育期变化及其影响因子研究

翟治芬

1,2

，胡 玮，严昌荣，刘 勤，刘 爽

1111

（1中国农业科学院农业环境与可持续发展研究所，北京100081；2农业部规划设计研究院，北京100125）

摘要：【目的】在全球气候变化背景下，分析中国玉米播种期和成熟期的变动情况以及气候资源变化特征。【方法】在收集整理全国2 414个县的玉米生育期数据的基础上，绘制了1970s和2000s中国玉米的播种期与收获期分布图；在整理全国618个气象站点1971—2010年气象资料的基础上，绘制了1970s时段和2000s时段中国年均温度、降雨和太阳辐射量空间分布图。以农业种植一级区为基本单位，建立不同区域农业气候资源变化与玉米生育期变化的回归方程，并将PRECIS模型中B2情景数据代入方程组预测2030s中国玉米的生育期。【结果】与1970s时段相比，2000s时段东北大豆春麦甜菜区的玉米播种期基本保持不变；其它各农业种植一级区的玉米播种期均提前约1—15 d；除东北大豆春麦甜菜区和北部高原小杂粮甜菜区春玉米的成熟期平均推迟了11 d和3 d，2000s时段其它玉米种植区域的成熟期平均提前3—12 d。2000s时段云贵高原稻玉米烟草区的玉米生育期缩短约5 d，黄淮海棉麦油烟果区、华南双季稻热带作物甘蔗区和西北绿洲麦棉甜菜葡萄区的玉米生育期基本保持不变；其它各区域玉米生育期均有所延长。与2000s时段相比，B2情景下，2030s东北大豆春麦甜菜区的春玉米播种期将推迟2—5 d，其它各农业种植一级区的玉米播种期将提前2—19 d；东北大豆春麦甜菜区、北部高原小杂粮甜菜区和华南双季稻热带作物甘蔗区的玉米成熟期将推迟4—15 d，黄淮海棉麦油烟果区、长江中下游稻棉油桑茶区、川陕盆地稻玉米薯类柑橘桑区和云贵高原稻玉米烟草区的玉米成熟期将提前2—12 d，南方丘陵双季稻茶柑橘区和西北绿洲麦棉甜菜葡萄区的玉米成熟期则基本保持不变。2030s时段黄淮海棉麦油烟果区和云贵高原稻玉米烟草区的玉米生育期则将缩短3—6 d，其它区域的玉米生育期将延长2—15 d。【结论】中国气候正朝着增温、变干和低辐射的方向发展。受温度、降雨和太阳辐射量变化的影响，中国不同农业种植区域内玉米生育期变动明显，其中除东北大豆春麦甜菜区外玉米播种期以提前为主，玉米成熟期的变动则较为复杂，玉米的生育期则以延长为主。

关键词：生育期；播种期；成熟期；温度；降雨；太阳辐射量

Change of Maize Growth Period and Its Impact Factor in China

ZHAI Zhi-fen1,2, HU Wei1, YAN Chang-rong1, LIU Qin1, LIU Shuang1

(1Institute of Environmental and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences，Beijing 100081;

2

Chinese Academy of Agricultural Engineering, Beijing 100125)

Abstract: 【Objective】 In the background of global climate change, the sowing time, maturity time and growth period of maize, and the changes of climate resource in China were analyzed.【Method】After collection of maize growth period data of 2414 counties and meteorological data of 618 meteorological stations from 1971 to 2010, by using the Kriging method, the distribution of maize sowing time and maturity time of 1970s and 2000s were drawn, and the annual temperature, precipitation and solar radiation were calculated in those two periods. Taking the crops area as the basic unit, the regression equation between agro-climatic resources and the change of maize sowing time and maturity time were built in different areas. Put B2 scenario data from PRECIS model into the equations, the maize growth period were predicted in 2030s. 【Result】 Compared with 1970s, the sowing time in 2000s was unchanged in Northeast soybean, spring wheat and beet area; while in other agricultural planting areas were advanced about 1 to 15 收稿日期：2012-05-15；接受日期：2012-08-28 基金项目：“十二五”国家科技支撑计划课题（2012BAD09B01）、国家“973”项目（2012CB955904）、农业部公益性行业科研专项（200903009） 联系方式：翟治芬，E-mail：zhaizhifen0821@163.com。通信作者严昌荣，E-mail：yancr@ieda.org.cn

4588 中 国 农 业 科 学 45卷

days. And the maturity time was delayed 11 days and 3 days in Northeast soybean, spring wheat, beet area and Northern plateau small grains, beet area, while other areas advanced 3 to 12 days. The maize growth period shortened about five days in Yunnan-Guizhou Plateau rice, corn, tobacco area, which remained unchanged in cotton, wheat, oil, tobacco, fruit area, South double cropping of rice, tropical crops, sugarcane area and Northwest Oasis wheat, cotton, beet, grape area, while in other regions has extended longer. In 2030s of B2 scenarios, the maize sowing time will delay 2 to 5 days in Northeast soybean, spring wheat, beet area, while in other agricultural planting areas may be advanced 2 to 19 days comparing to 2000s. The maturity time will delay 4-15d in Northeast soybean, spring wheat, beet area, Northern plateau small grains, beet area and South double cropping of rice, tropical crops, sugarcane area, which will advance 2 to 12 days in cotton, wheat, oil, tobacco, fruit area，Yangtze River rice, cotton, oil, mulberry, tea area，Sichuan and Shaanxi Basin rice, corn, potato, citrus, mulberry area and Yunnan-Guizhou Plateau rice, corn, tobacco area, and it will essentially unchanged in South Hills double cropping of rice, tea, citrus area and Northwest Oasis wheat, cotton, beet, grape area. Comparing to 2000s, the maize growing period in 2030s will reduce 3 to 6 days in cotton, wheat, oil, tobacco, fruit area and Yunnan-Guizhou Plateau rice, corn, tobacco area, and in other regions it will extend 2 to 15 days. 【Conclusion】The climate is changing to warming, dry and low radiation. By impact of temperature, precipitation and solar radiation, there exists an objective change of maize growth stage in different agricultural planting areas, and the sowing time will advance in most part of China, the maturity time will become more complex, and the growing period of maize will extend in most part of China.

Key words: growth period; sowing time; maturity time; temperature; precipitation; solar radiation

0 引言

【研究意义】气候变化是近百年来最显著的地球环境变化之一，农业生产受其影响明显而广泛[1-4]。粮食作物对当前及未来气候适宜性问题引起人们的普遍关注[5-7]。其中，作物的生长和发育进程作为反映作物潜在产量的重要指标，对指导不同地区的种植生产具有十分重要的意义。【前人研究进展】许多学者对作物生长发育节奏、作物种植界限北移、可种植海拔高度上限抬升以及产量增减等受气候变化的影响进行了在玉米杂交种中，大量的研究工作：Evans[8]研究表明，籽粒生长时期长短与产量高低关系密切。IPCC综合报告认为，在北半球，特别是在高纬度地区，生长期在过去的40年中每10年大约变长1—4 d[9]。Menzel等[10]利用1959—1993年国际物候花园里观察到的资料得出，由于温度升高，欧洲的生长期增长了10.8 d。Myneni等[11]利用1981—1991的卫星资料，估计北半球的生长期提前了（8±3）d，延迟了（4±2）d。国内，徐铭志等[12]研究表明在近40年中，气候生长期在全国范围平均增加了6.6 d。薛昌颖等[13]研究表明，水稻生育期的响应表现为播种和移栽日期显著提前和移栽-抽穗期长度显著延长，导致生育期延长。张明捷等[14]认为由于小麦生育期内冬季气候明显变暖，其生长期和灌浆期延长。张建立[15]研究表明，豫南地区夏玉米生产中气象因子的限制作用由大到小依次为降海拔高度、水、光照和温度。Dowswell等[16]研究表明，纬度与温度对生育期的影响具有并列重要的作用。同

时，一些学者认为生育期长短与纬度和海拔高度呈正向对应关系[17-18]。【本研究切入点】作物生育期的变化是作物本身生理过程和环境条件综合作用的结果。由于作物生长影响因子的多样性及相互间关系的复杂性，在评估作物气候适宜性时，各因子选取及其关系的表达非常困难。故以往温度、降水等单项气候因子的适宜性评价报道较多，综合气候适宜性变化研究成果比较少见。而单因子气候适宜性评价模型只考虑某一因子在其余影响因子变化不大的条件下的作用，不能反映作物生长受气候因子影响的真实情况。综合评价各气候因子对玉米生长的影响程度，极有必要。同时由于中国幅员辽阔，不同地区的玉米生育期存在较大差异，调查所有地区玉米生育期的数据，需要消耗大量的人力物力，大范围内分析玉米生育期特点的研究较少。【拟解决的关键问题】本文分析了1971—1980和2001—2010时段内中国玉米播种期和成熟期及气候因子（温度、降雨和太阳辐射量）的变化特点；并以农业种植一级区为基本单元，运用回归分析的方法建立年均温度、年均降雨和年均太阳辐射量变化以及经度、纬度和海拔与玉米播种期和成熟期变化量的回归方程，预测B2情景下2031—2040时段中国玉米播种期和成熟期的分布情况。

1 研究方法

1.1 数据来源

本文采用的1971—1980时段（1970s）玉米生育期数据为日期型数据，主要参考《中国主要农作物气

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4589

[19]

候资源图集》并进行校准；2001—2010时段（2000s）中国的可持续发展模式基本相符[20]。为了保持数据的一致性，利用最短距离算法提取距离618个气象站点最近的气象格点，再用克里格插值法对情景数据进行空间插值，然后运用ArcGIS10.0软件的多点提取工具从插值结果中提取气象站点在2030s内的气候资源数据。气象资料分析要素包括温度、降水、日照等，各时段的数据由日值数据经统计分析求得年均值。

由于中国幅员辽阔，不同区域气候、土壤、地形差异较大，导致作物种植制度不同，因此在分析作物生育期时需结合农业种植区划。本文主要参考中国农业种植分区[21]，将中国分为10个农业种植大区：东北大豆春麦甜菜区、北部高原小杂粮甜菜区、黄淮海棉麦油烟果区、长江中下游稻棉油桑茶区、南方丘陵双季稻茶柑橘区、华南双季稻热带作物甘蔗区、川陕盆地稻玉米薯类柑橘桑区、云贵高原稻玉米烟草区、西北绿洲麦棉甜菜葡萄区和青藏高原青稞小麦油菜区，具体分区结果见图1。

玉米生育期数据来源于实地调研。1970s和2000s两个时段的玉米生育期数据获取地点、方法及标准一致；排除数据缺失和数据异常情况，实际进行统计整理的县共有2 414个。由于得到的是日期型数据，而进行数据分析的数据时须是数值型数据，所以将所有样本点的日期型数据转化为日序。

采用的1970s和2000s的气候资源数据来源于国家气象局“气象数据共享网”，考虑到资料年代的统一性、代表性和稳定性，筛选的气象站点为618个，其分布见图1。采用的2031—2040时段（2030s）的气候资源数据为PERCIS（Providing Regional Climate for Impact Study）模型B2情景下中国区订正后的栅格数据，由中国农业科学院许吟隆研究员的科研团队提供。B2情景为SRES情景中区域可持续发展情景。在该情景下，世界体现出区域化倾向，世界人口缓慢但持续地增长，温室气体排放维持中-低水平，该情景与

01-10分别代表东北大豆春麦甜菜区、北部高原小杂粮甜菜区、黄淮海棉麦油烟果区、长江中下游稻棉油桑茶区、南方丘陵双季稻茶柑橘区、华南双季稻热带作物甘蔗区、川陕盆地稻玉米薯类柑橘桑区、云贵高原稻玉米烟草区、西北绿洲麦棉甜菜葡萄区和青藏高原青稞小麦油菜区

Note: 01-10 express: Northeast soybean, spring wheat and beet area; Northern plateau small grains, beet area; Huanghuaihai cotton, wheat, oil, smoke, fruit area; Yangtze river rice, cotton, oil, mulberry, tea area; South hills double cropping of rice, tea, citrus area; South double cropping of rice, tropical crops, sugarcane area; Sichuan and Shaanxi Basin rice, corn, potato, citrus, mulberry area; Yunnan-Guizhou Plateau rice, corn, tobacco area; Northwest oasis wheat, cotton, beet, grape area and Qinghai-Tibet Plateau highland barley wheat, rapeseed area, respectively

图1 中国农业种植区分布图

Fig. 1 The agricultural planting areas of China

4590 中 国 农 业 科 学 45卷

1.2 研究方法

1.2.1 太阳辐射量计算 太阳总辐射量是任一特定时段内（如日、月、年）水平面上太阳总辐射强度的累计值（MJ·m-2）。其太阳辐射量的日值计算公式为[22]：

Qn

TI0

2(0sinsincoscossin0) （1）

式中，Qn为日太阳辐射总量（MJ·m-2·d-1）；T为时间周期（24×60 min·d-1）；I0为太阳常数为0.0820（MJ·m-2·min-1）；ρ为日地距离系数，无量纲；φ为

地理纬度（rad）；δ为太阳赤纬（rad）；ω0为日落时角（rad）。

1.2.2 克里格插值法 克里格插值法是以空间自相关为基础，利用区域化变量的原始数据和变异函数的结构特点，对未知点的区域化变量进行线性无偏最优估计的一种插值方法。本文采用普通克里格法[23]进行插值计算（Ordinary Kriging, OK），其插值公式为:

n

Z

i

Z(xi

) （2）

i1

式中，Z为未知点x的预测值；Z(xi)为测量点x的测量值；λi是测量值对预测值影响大小的系数，由克里格方程组求得，而克里格方程组主要影响因素是根据原始数据的分布建立的协方差函数或变异函数模型，详细内容见文献[24]。普通克里格法的优化包括邻域搜捕优化和半变异函数理论模型、块金方差、基台值、步长等方面。

1.2.3 回归分析 本文采用SAS软件分别对玉米的播种期变化日数和成熟期变化日数与气象站点的年均降雨量、年均温度和年均太阳辐射量变化量及站点经度、纬度和海拔高度进行多元线性回归[25]，其一般公式如下：

Y=aX1+bX2+cX3+dX4+eX5+eX6+ε （3） 式中，Y表示生育期要素变化值（播种期或成熟期日序数）；X1为经度，X2为纬度，X3为海拔高度，X4为年均降雨量变化量，X5为年均温度变化量，X6为年均太阳辐射量变化量；a、b、c、d、e、f分别代表各要素的相关系数；ε为误差项系数。

另外，本文利用SAS软件中逐步回归的方法筛选显著影响玉米的播种期变化和成熟期变化的气象因子。

2 结果

2.1 1970s与2000s玉米生育期变化分析

2.1.1 玉米播种期变化 图2为1970s和2000s时段玉米播种期日序等值线分布图。根据实地调查结果，春玉米和夏玉米的种植分界基本以黄淮海地区北界区为分界，即图2中蓝色与绿色线的交界区域，界区以北为春玉米种植区，界区以南为夏玉米种植区。从图2可以看出，以春玉米与夏玉米的种植分界区为界限，玉米播种期分别向北推迟、向南提前。春玉米种植区播种期最早出现在4月10日，最晚出现在5月11日，时间跨度长达30 d左右；夏玉米种植区播种期最早出现在5月10日，最晚出现在6月21日，时间跨度长达45 d左右。

在1970s时段，西北绿洲麦棉甜菜葡萄区内春玉米播种期出现在4月11日—5月1日，东北大豆春麦甜菜区春玉米和北部高原小杂粮甜菜区内春玉米播种期出现在4月21日—5月11日；黄淮海棉麦油烟果区和长江中下游稻棉油桑茶区夏玉米播种期出现在5月22日—6月21日，南方丘陵双季稻茶柑橘区和华南双季稻热带作物甘蔗区夏玉米的播种期最早出现在5月6日—5月23日。特别的，在秦岭大巴山区与云贵高原的东北部地区由于地势较高，夏玉米的播种期出现在5月30日—6月21日。与1970s时段相比，2000s时段东北大豆春麦甜菜区的玉米播种期基本保持不变，个别地区存在推迟1—2 d的现象；其它各农业种植一级区的玉米播种期均有不同程度的提前，其中北部高原小杂粮甜菜区和西北绿洲麦棉甜菜葡萄区的春玉米播种期提前时间较短，约1—3 d；川陕盆地稻玉米薯类柑橘桑区和长江中下游稻棉油桑茶区的夏玉米播种期提前时间较长，约10—15 d；其它夏玉米种植区的播种期约提前3—8 d。

2.1.2 玉米成熟期变化 通过分析1970s与2000s时段年玉米成熟期日序分布，从图3可以看出，以春玉米与夏玉米的种植分界区为界限，玉米成熟期分别向北向南提前。春玉米种植区成熟期最早出现在9月1日，最晚出现在9月27日，时间跨度约1个月；夏玉米种植区成熟期最早出现在8月21日，最晚出现在10月2日，时间跨度长达1.5个月。

在1970s时段，西北绿洲麦棉甜菜葡萄区内春玉米成熟期出现在8月28日—9月20日，东北大豆春麦甜菜区春玉米和北部高原小杂粮甜菜区内春玉米成熟期出现在9月12日—9月25日；黄淮海棉麦油烟果区玉米成熟期出现在9月22日—10月2日；长江中下游稻棉油桑茶区、南方丘陵双季稻茶柑橘区和华南双季稻热带作物甘蔗区夏玉米的成熟期最早出现在

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4591

图2 1970s和2000s玉米播种期分布图

Fig. 2 Distribution of maize sowing time in 1970s and 2000s

4592 中 国 农 业 科 学 45卷

8月21日—9月1日；在秦岭大巴山区与云贵高原稻玉米烟草区，夏玉米的成熟期出现在9月11日—9月21日。与1970s时段相比，2000s时段东北大豆春麦甜菜区和北部高原小杂粮甜菜区春玉米的成熟期分别平均推迟了11 d和3 d。其它区域的玉米成熟期均有不同程度的提前，其中，川陕盆地稻玉米薯类柑橘桑区和云贵高原稻玉米烟草区的成熟期提前时间较长，约10—15 d；其它区域的成熟期约提前3—8 d。

由于玉米播种期和成熟期的变动，与1970s时段相比，2000s时段黄淮海棉麦油烟果区、华南双季稻热带作物甘蔗区和西北绿洲麦棉甜菜葡萄区的玉米生育期基本保持不变；南方丘陵双季稻茶柑橘区和川陕盆地稻玉米薯类柑橘桑区的玉米生育期延长了2—3 d，北部高原小杂粮甜菜区、长江中下游稻棉油桑茶区延长了4—6 d，东北大豆春麦甜菜区平均延长了11 d；而云贵高原稻玉米烟草区的玉米生育期则平均缩短了5 d左右。

2.2 1970s与2000s气候资源变化分析

2.2.1 年均温度 基于国家气象局“气象数据共享网”1971—1980和2001—2010时段的气候资源数据，运用ArcGIS10.0软件分别绘制两个时段中国年均温度分布图（图4）。由图4看出，中国年均温度分布整体呈南高北低梯级分布的特点，青藏高原青稞小麦油菜区和东北大豆春麦甜菜区平均温度最低，平均温度为2.79—4.95℃，北部高原小杂粮甜菜区与西北绿洲麦棉甜菜葡萄区次之，平均温度为5.78—8.24℃，华南双季稻热带作物甘蔗区温度最高，平均温度为21.48—22.27℃，其它地区平均温度为12.59—18.28℃。即在1970s时段，大兴安岭地区及青海西南与西藏的东北部地区的平均温度为-5.00—-0.31℃，不能满足玉米种植所需的热量资源；黑龙江、吉林、青海和西藏的大部分地区、内蒙古北部地区以及新疆的东北地区平均温度为0.36—5.18℃，华北地区（包括陕西省）的北部、吉林省、宁夏回族自治区、甘肃和新疆大部、内蒙古的西部、四川省西北部以及西藏的东南部地区的平均温度为5.20—10.47℃，华北地区和四川省中南部、陕西的南部、云南省的部分地区以及南疆的阿克苏地区平均温度为10.71—15.66℃，长江中下游地区包括重庆、湖北、湖南、江西、浙江、福建及广西的北部和云南的大部分地区平均温度为15.17—20.08℃，广西的中南部、广东的大部和海南省的平均温度均在20℃以上，西沙地区达27.18℃。

2000s时段，中国年均温度的基本分布情况与1970s相似，但与1970s相比，温度条带北移0.2°左右，全国年均温度增加1.03℃左右。具体的，川陕盆地稻玉米薯类柑橘桑区、云贵高原稻玉米烟草区和华南双季稻热带作物甘蔗区增温幅度较低，平均增加0.66—0.79℃，黄淮海棉麦油烟果区和南方丘陵双季稻茶柑橘区次之，平均增加0.66—0.79℃，其它地区平均增温幅度在1.13—1.17℃。由以上分析可以看出，在全球气候变化的背景，中国的年均气温在逐步上升。

2.2.2 年均降雨量 中国的玉米种植区内70%以上的地区以雨养种植为主，因此降水量是玉米生长重要指标之一。通过分析可知，中国年均降雨量从西北到东南呈条带状增加的趋势（图5），1970s和2000s时段中国平均年降雨量，分别为800.49和794.28 mm。 在1970s时段，西北绿洲麦棉甜菜葡萄区年均降雨量仅128.88 mm，远小于玉米生长所需的水分，因此该区域的玉米种植以灌溉为主；北部高原小杂粮甜菜区、东北大豆春麦甜菜区和黄淮海棉麦油烟果区的年均降雨量分别为405.33、563.12和 685.25 mm，基本可满足玉米生长的需要，遇到干旱年份或时段进行补充灌溉；其它区域的年均降雨量均大于1 000 mm，完全可以保证玉米生长期内的水分需求（青藏高原青稞小麦油菜区非玉米种植区，不做讨论）。

与1970s时段相比，

2000s时段中国年平均降水量减少了6.21 mm，其中：西北绿洲麦棉甜菜葡萄区年均降雨量将增加至152.19 mm，但仍不能满足玉米生长所需；北部高原小杂粮甜菜区的年均降雨量减少26.50 mm，使该地区的水资源更为短缺；东北大豆春麦甜菜区和黄淮海棉麦油烟果区的年均降雨量基本保持不变；其它区域的年平均降雨量均有所减少，但对玉米种植不会产生较大影响。

2.2.3 年均太阳辐射量 太阳辐射是地面的主要能量来源，是天气和气候形成及变化的基础，其分布决定了区域气候特点和农业生产布局，对维持粮食稳产、增产有重要的作用[26]。通过对不同站点的年均太阳辐射量的计算可以看出，1970s和2000s时段中国年均太阳辐射量，分别为5 166.84和5 064.87 MJ·m-2（图6）。1970s时段，云贵高原稻玉米烟草区的东北部地区和川陕盆地稻玉米薯类柑橘桑区的年均太阳辐射量最低，平均值介于3 143.86—3 998.92 MJ·m-2；南方丘陵双季稻茶柑橘区和长江中下游稻棉油桑茶区平均太阳辐射量的平均值分别为4 452.12和4 673.25 MJ·m-2；

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4593

图3 1970s与2000s玉米成熟期分布图

Fig. 3 Distribution of maize maturity time in 1970s and 2000s

4594 中 国 农 业 科 学 45卷

图4 1970s与2000s年均温度分布图

Fig. 4 Distribution of average annual temperature during 1970s and 2000s in China

东北大豆春麦甜菜区、北部高原小杂粮甜菜区、黄淮海棉麦油烟果区、华南双季稻热带作物甘蔗区及云贵高原稻玉米烟草区的大部分地区年均太阳辐射量的平西北绿洲麦棉甜均值介于4 593.46—5 542.05 MJ·m-2；

菜葡萄区年均太阳辐射量的平均值介于5 569.47—

6 183.55 MJ·m-2，而在青藏高原青稞小麦油菜区最大值达7 957.77 MJ·m-2。

由图6可以看出，2000s时段，中国年均太阳辐射量＞7 000 MJ·m-2的范围缩小，＜5 000 MJ·m-2的范围扩大至整个长江中下游地区，全国年太阳辐射量平

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4595

图5 1970s与2000s中国平均降雨量分布图

Fig. 5 Distribution of average annual precipitation during 1970s and 2000s in China

均减少了102.06 MJ·m-2。其中，黄淮海棉麦油烟果区长江中下的减少量最大，平均减少了312.06 MJ·m-2；游稻棉油桑茶区次之，平均减少了129.59 MJ·m-2；

但云贵高原稻玉米烟草区和华南双季稻热带作物甘蔗区的年太阳辐射量略有增加，分别增加了2.89%和1.65%。

4596 中 国 农 业 科 学 45卷

图6 1970s与2000s年均太阳辐射量分布图

Fig. 6 Distribution of average annual solar radiation during 1970s and 2000s in China

2.3 气候资源与玉米生育期变化影响分析

以2000s相对于1970s各气象站点的年均气温、年均降雨量和年均太阳辐射量的变化量以及站点的经度、纬度和海拔高度为自变量，建立不同农业种植一级区内显著性<0.05的玉米播种期和成熟期变化量的拟合方程（表1）。玉米播种期变化量的拟合方程中，

华南双季稻热带作物甘蔗区和黄淮海棉麦油烟果区的决定系数较高，分别是0.7596、0.7028；长江中下游稻棉油桑茶区、西北绿洲麦棉甜菜葡萄区和南方丘陵双季稻茶柑橘区的决定系数次之；其它4个区域的决定系数则较低。玉米成熟期变化量的拟合方程中，东北大豆春麦甜菜区和北部高原小杂粮甜菜区的决定系

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4597

数最高，分别是0.9075和0.8434；南方丘陵双季稻茶柑橘区、长江中下游稻棉油桑茶区和西北绿洲麦棉甜菜葡萄区次之。

利用逐步回归的方法，筛选到各农业种植一级区内显著性＜0.01的影响玉米播种期和成熟期变化的气象因子（表）。由表可以看出，温度、降雨和太阳辐射量均是影响中国玉米种植的气象因子，且在温度和降雨基本满足玉米生长需求时，太阳辐射量的影响效果更为明显。具体的，在东北大豆春麦甜菜区，温度对于玉米成熟具有制约作用；在北部高原小杂粮甜菜区，光、温、水、热均对玉米播种具有显著影响；在黄淮海棉麦油烟果区，降雨和太阳辐射量对玉米种植同样具有制约作用；在西北绿洲麦棉甜菜葡萄区，温

表 玉米播种期与成熟期变化拟合表

度是影响玉米播种的主要气象因子。在南方丘陵双季稻茶柑橘区，降雨是影响玉米播种的主要气象因子；在华南双季稻热带作物甘蔗区，温度将对玉米的成熟产生显著的影响；在川陕盆地稻玉米薯类柑橘桑区和云贵高原稻玉米烟草区，降雨对玉米种植具有制约作用。

2.4 2030s玉米生育期预测

2.4.1 2030s气候资源变化分析 基于PERCIS模型B2情景下中国区校订后的气象栅格数据，采用的克里格插值法及ArcGIS中软件多点提取工具，提取2030s时段618个气象站点上的B2情景数据，将之与2000s时段的气候数据作对比，绘制得到2030s中国气候资源变化分布图（图7）。

Table The fitting table of the change of maize sowing time and maturity time

种植区域名称 Planting areas name 东北大豆春麦甜菜区 and beet area

北部高原小杂粮甜菜区 Northern plateau small grains, beet area

黄淮海棉麦油烟果区 smoke, fruit area

长江中下游稻棉油桑茶区 Yangtze River rice, cotton, oil, mulberry, tea area 南方丘陵双季稻茶柑橘区 South hills double cropping of rice, tea, citrus area

华南双季稻热带作物甘蔗区 South double cropping of rice, tropical crops, sugarcane area

播种期回归方程

Regression equation of sowing time

Y1=-12.7764+0.0934X1-0.0213X20.0132X6

Y2=4.8793-0.1042X1+0.0687X2 +0.0004X3-0.0137X4+1.069X5- 0.0031X6

Y3=-107.4959+0.2922X1+1.9708+0.0124X6

Y4=22.3903-0.5093X1+0.9505X2-0.0053X3+0.0027X4-1.319X5+ 0.0112X6

Y5=80.8232-0.6018X1-0.6006X2 +0.0013X3+0.0061X4-1.0266X5-0.0173X6

Y6=45.8109-0.2621X1-0.9708X2 +0.0141X3+0.122X4-3.893X5+ 0.0134X6

0.3514

X6

0.7596

X6

0.5379

X4 X6

0.6081

X6

0.7028

X4 X6

0.435

X4 X5 X6

0.3182

Northeast soybean, spring wheat +0.0082X3-0.0151X4+1.496X5+

决定系数主要影响因子成熟期回归方程

R

2

决定系数 主要影响因子

R2 0.9075

The main impact factorX5, X6

The main impact factor X6

Regression equation of sowing time

Z1=-123.708+0.3973X1+1.777X2+0.0159X3+0.0083X4+2.0993X5+ 0.0168X6

Z2=-43.0154-0.4431X1+2.3589X2+0.0012X3-0.006X4+0.2158X5+ 0.0058X6

Z3=101.2042-0.8788X1-0.03X2- 0.0044X3-0.0026X4-0.37X5+ 0.0041X6

Z4=56.9861-0.8888X1+1.3783X2-0.0044X3+0.0016X4-0.1462X5+ 0.0107X6

Z5=157.1496-1.1555X1-0.9794X2-0.0016X3+0.0031X4+0.6555X5- 0.0145X6

Z6=-98.2998+0.08927X1+0.0503X2+0.0091X3-0.0038X4-5.3383X5 -0.0251X6

Z7=66.4064-0.6617X1-0.1974X2-0.0024X3-0.00162X4+3.4153X5+ 0.0054X6

0.8434 X6

0.4827 X4, X6

Huanghuaihai cotton, wheat, oil, X2+0.0003X3+0.2556X4+1.3866X5

0.5777 X6

0.649 X6

0.4981 X5, X6

川陕盆地稻玉米薯类柑橘桑区 Y7=149.7043-0.7703X1-2.5712X2Sichuan and Shaanxi Basin rice, +0.0009X3+0.0008X4-1.2098X5-corn, potato, citrus, mulberry area 0.0155X6 云贵高原稻玉米烟草区 Yunnan-Guizhou Plateau rice, corn, tobacco area 西北绿洲麦棉甜菜葡萄区 Northwest oasis wheat, cotton, beet, grape area

Y8=123.6642-1.6169X1+1.5828 X2-0.0014X3+0.0531X4-1.1443X5+0.008X6

Y9=18.8168+0.1688X1-0.8673X2-0.0023X3+0.0091X4+1.0743X5- 0.0018X6

0.3918 X4

0.3488X4 Z8=-55.9515-0.6012X1+4.1049X2-0.0009X3+0.0149X4-0.6063X5- 0.0101X6

0.3842 X4, X6

0.5915X5 Z9=-49.865+0.02167X1+0.5760X2-0.0003X3+0.0136X4+2.3299X5- 0.0115X6

0.5331 X6

X1、X2、X3、X4、X5、X6分别表示气象站点经度、纬度、海拔高度、年均降雨量变化量、年均温度变化量和年均太阳辐射量变化量，Y表示玉米播种

期变化量， Z表示玉米成熟期变化量

X1, X2, X3, X4, X5, X6 means longitude, latitude, altitude of the meteorological station, the change of average annual precipitation, the change of average annual temperature and the change of average annual solar radiation, respectively. Y represents the change of maize sowing time, while Z represents the change of maize maturity time

4598 中 国 农 业 科 学 45卷

图7 B2情景下2030s中国气候资源变化分布图

Fig. 7 Change of Chinese climate resource during 2030s in B2 scenario

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4599

与2000s相比，在B2情景下2030s中国温度升高的趋势趋于平缓，全国年平均温度可能增加1.07℃左右，且不同区域温度变动情况不同，可能出现的变动范围为-5.86—7.46℃。由图7-a可以看出， 东北大豆春麦甜菜区的大部分地区、北部高原小杂粮甜菜区的大部、黄淮海棉麦油烟果区、长江中下游稻棉油桑茶区大部分地区、川陕盆地稻玉米薯类柑橘桑区的西部、云贵高原稻玉米烟草区、华南双季稻热带作物甘蔗区以及西北绿洲麦棉甜菜葡萄区的新疆中东部地区均存在不同程度的增温，温度变化范围为1.54—7.46℃；东北大豆春麦甜菜区的大兴安岭地区和长白山地区、北部高原小杂粮甜菜区的东部地区和长江中下游稻棉油桑茶区的东部存在一定程度的降温，温度变化范围在-0.64—0.93℃；青藏高原青稞小麦油菜区和西北绿洲麦棉甜菜葡萄区的东部地区的降温幅度较大，温度变化范围在-5.86—-0.76℃。

2030s中国平均年降水量比2000s将平均减少34.44 mm，减少率为4.34%。由图7-b可以看出，东北大豆春麦甜菜区和黄淮海棉麦油烟果区的东部、长江中下游稻棉油桑茶区、南方丘陵双季稻茶柑橘区、华南双季稻热带作物甘蔗区、云贵高原稻玉米烟草区的南部、青藏高原青稞小麦油菜区东部的降雨均表现为减少的趋势，降雨变率范围在-34.63—-1.1%；西北绿洲麦棉甜菜葡萄区内甘肃省的中部和新疆的北部降雨减少趋势更为显著，降雨变率范围在-65.01—-32.06%；东北大豆春麦甜菜区的和黄淮海棉麦油烟果区的东部、北部高原小杂粮甜菜区的大部、川陕盆地稻玉米薯类柑橘桑区以及云贵高原稻玉米烟草区的北部的降雨略有增加，降雨变率范围在0.14%—27.57%；青藏高原青稞小麦油菜区的中西部地区、西北绿洲麦棉甜菜葡萄区内新疆的南部和内蒙古的西部地区以及北部高原小杂粮甜菜区的北部降雨的增加趋势较为明显，降雨变率范围在64.98%—255.56%，然而由于这些地区降雨量较低，降雨增加的相对值仅为70.59—114.02 mm，远不能满足农业生产的需求。可见，在B2情景下未来中国的水资源将更加短缺。

B2情景下，2030s中国年均太阳辐射量预计比2000s时段平均减少168.91 MJ·m-2，减少率为1.95%。由图7-c可以看出中国大部分区域的太阳辐射量的变率范围在-10%—0%，特别提到的是，青藏高原青稞小麦油菜区和华南双季稻热带作物甘蔗区的年太阳辐射量减少较为显著，可能平均减少11.87%和8.69%；东北大豆春麦甜菜区、黄淮海棉麦油烟果区和云贵高原

稻玉米烟草区内的部分区域年均太阳辐射量可能略有增加，平均增幅在0.07%—11.47%；川陕盆地稻玉米薯类柑橘桑区增幅较大，增幅范围为14.75%—27.01%。

2.4.2 2030s玉米生育期预测 将B2情景下2030s相对于2000s各气象站点的年均气温、年均降雨量和年均太阳辐射量的变化量代入表1的回归方程，得到该情景下2030s时段中国玉米播种期和收获期变动量，并绘制玉米生育期预测分布图（图8）。

从图8-a可以看出，以春玉米与夏玉米的种植分界区为界限，2030s东北和西北地区的春玉米播种期逐渐将向北推迟，黄淮海地区的夏玉米播种期逐渐向南推迟，南方各地的夏玉米播种期则表现为逐渐向南提前。春玉米种植区播种期最早出现在4月1日，最晚出现在5月15日，时间跨度长达45 d左右；夏玉米种植区播种期最早出现在4月18日，最晚出现在6月19日，时间跨度长达60 d左右。具体的，西北绿洲麦棉甜菜葡萄区内春玉米播种期将出现在4月1日—4月25日，东北大豆春麦甜菜区春玉米和北部高原小杂粮甜菜区内春玉米播种期将出现在4月20日—5月15日；黄淮海棉麦油烟果区的夏玉米播种期将出现在5月19日—7月5日，其它区域的夏玉米的播种期将出现在4月18日—5月15日。与2000s相比，2030s玉米播种期变化较大：东北大豆春麦甜菜区的春玉米播种期将推迟2—5 d；其它各农业种植一级区的玉米播种期则表现为提前，长江中下游稻棉油桑茶区、川陕盆地稻玉米薯类柑橘桑区和华南双季稻热带作物甘蔗区的夏玉米播种期提前时间最长，约9—19 d；北部高原小杂粮甜菜区和西北绿洲麦棉甜菜葡萄区的春玉米播种期提前时间最短，约2—4 d左右；其它区域的播种期约提前5—10 d左右。

由2030s玉米成熟期预测结果可以看出，以春玉米与夏玉米的种植分界区为界限，玉米成熟期分别向北推迟，向南提前（图8-b）。在B2情景下，春玉米种植区成熟期最早出现在8月20日，最晚出现在10月17日，时间跨度约60 d左右；夏玉米种植区成熟期最早出现在8月8日，最晚出现在9月27日，时间跨度亦为50 d左右。具体的，西北绿洲麦棉甜菜葡萄区内春玉米成熟期出现在8月20日—9月27日，东北大豆春麦甜菜区春玉米和北部高原小杂粮甜菜区内春玉米成熟期出现在9月17日—10月17日；黄淮海棉麦油烟果区的玉米成熟期出现在9月9日—9月27日；其它各夏玉米种植区的成熟期出现在8月8日—9

4600 中 国 农 业 科 学 45卷

图8 B2情景下2030s玉米播种期和成熟期预测分布图

Fig. 8 Distribution of maize sowing time and maturity time during 2030s in B2 scenario

月7日。与2000s相比，2030s时段东北大豆春麦甜菜区的玉米成熟期将推迟4—15 d，北部高原小杂粮甜菜区和华南双季稻热带作物甘蔗区将平均推迟4 d；川陕

盆地稻玉米薯类柑橘桑区和云贵高原稻玉米烟草区的玉米成熟期将提前6—12 d；黄淮海棉麦油烟果区和长江中下游稻棉油桑茶区将提前2—8 d；南方丘陵双季

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4601

稻茶柑橘区和西北绿洲麦棉甜菜葡萄区的玉米成熟期则基本保持不变。

与2000s时段相比，B2情景下2030s西北绿洲麦棉甜菜葡萄区和南方丘陵双季稻茶柑橘区生育期将延长2—4 d，北部高原小杂粮甜菜区、长江中下游稻棉油桑茶区将延长5—6 d；华南双季稻热带作物甘蔗区、川陕盆地稻玉米薯类柑橘桑区的玉米和东北大豆春麦甜菜区将延长11—15 d；而黄淮海棉麦油烟果区和云贵高原稻玉米烟草区的玉米生育期则将缩短3—6 d。

3 讨论

全球气候变化背景下，中国表现为暖干化的趋势，同时这一暖干化特征有向南和向东扩展的趋势[27-28]，本文也得到了相同的结果，中国年平均气温条带北移，温度上升效应十分显著，年降水量和年太阳辐射量呈减少的趋势。

本文依据《中国农业区划图集》中农业种植区域的划分，在本研究组已有的主要农作物生育期研究工作的基础上，综合考虑气候资源对玉米播种期和成熟期的影响，分析了2000s相比1970s和B2情景下2030s相比2000s中国玉米播种期和成熟期空间位移特征，及由此带来的玉米生育期变化。研究结果表明，受气候变化影响，中国玉米的生育期以延长为主。这一结论与前人结果[9-12]一致。本文的研究结果提高了玉米生长发育与气象的关系的分析精度，并对选育适宜的区域性栽培品种、预测玉米产量等具有辅助作用。

玉米生育期的变化是一个十分复杂的问题，本文在讨论玉米播种期和成熟期变化时，只是基于气候要素指标，并没有考虑土壤、社会经济和品种更替等因素的影响，同时增加多重考虑因素后玉米生育期变动与其影响因素间的相关关系的表达形式也有待于进一步研究。此外，2030s玉米生育期的预测仅引用了B2情景，存在很大的不确定性。以上不足有待进一步完善。

全球气候变化背景下，中国表现为暖干化的趋势，即年平均气温条带北移，温度上升效应十分显著，年降水量和年太阳辐射量呈减少的趋势，同时这一暖干化特征有向南和向东扩展的趋势[27-28]。本文依据《中国农业区划图集》中农业种植区域的划分，在本研究组已有的主要农作物生育期研究工作的基础上，综合考虑气候资源对玉米播种期和成熟期的影响，分析了2000s相比1970s和B2情景下2030s相比2000s中国玉米播种期和成熟期空间位移特征，及由此带来的玉

米生育期变化。研究结果表明，受气候变化影响，中国玉米的生育期以延长为主。这一结论与前人结果[9-12]一致。

本研究一方面可以提高玉米生长发育与气象关系的分析精度，如在农业气象灾害的研究中，生育期可以作为其分析的起止日期，从而得到各地与玉米生育进程密切相关的气象条件；另一方面对选育适宜当地栽培的品种、预测玉米产量等研究有很大的辅助作用，如通过分析品种自身的生长周期特点与各地的生育期相关性，得出该品种适宜推广的种植区域。但是，玉米生育期的变化是一个十分复杂的问题，本文在讨论玉米播种期和成熟期变化时，只是基于气候要素指标，并没有考虑土壤、社会经济和品种更替等因素的影响。对于2030s玉米生育期的预测仅引用了B2情景，存在很大的不确定性。以上不足有待进一步完善。

4 结论

（1）与1970s相比，2000s中国年均温度增加1.03℃左右，间隔为5℃的温度条带北移0.2°左右。与2000s相比，在B2情景下2030s中国温度升高的趋势趋于平缓，全国年平均温度可能增加1.07℃左右。2000s中国的年平均降水量比1970s时段减少了6.21 mm；B2情景下，2030s的年均降雨量将比2000s平均减少34.44 mm。与1970s相比，2000s时段中国平均太阳辐射量大于7 000 MJ·m-2的范围缩小，

低于5 000 MJ·m-2的范围扩大至整个长江中下游地区。

B2情景下，2030s中国年均太阳辐射量预计比2000s时段平均减少168.91 MJ·m-2。可见，全球气候变化背景下，中国气候正朝着增温、变干和低辐射的方向发展。 （2）以春玉米与夏玉米的种植分界区为界限，玉米播种期分别向北推迟、向南提前，玉米成熟期分别向北向南提前。与1970s时段相比，2000s时段东北大豆春麦甜菜区的玉米播种期基本保持不变；其它各农业种植一级区的玉米播种期均有不同程度的提前；2000s时段除东北大豆春麦甜菜区和北部高原小杂粮甜菜区春玉米的成熟期有所推迟外，其它玉米种植区域的成熟期平均提前3—12 d。由于玉米播种期和成熟期的变动，2000s时段云贵高原稻玉米烟草区的玉米生育期比1970s时段缩短约5 d，黄淮海棉麦油烟果区、华南双季稻热带作物甘蔗区和西北绿洲麦棉甜菜葡萄区的玉米生育期基本保持不变；其它各区域玉米生育期延长均有所延长。

（3）B2情景下以春玉米与夏玉米的种植分界区为

4602 中 国 农 业 科 学 45卷

界限，2030s东北和西北地区的春玉米播种期逐渐将向北推迟，黄淮海地区的夏玉米播种期逐渐向南推迟，南方各地的夏玉米播种期则表现为逐渐向南提前；玉米成熟期则以分界区为界限分别向北推迟，向南提前。与2000s相比，除东北大豆春麦甜菜区外，2030s时段玉米播种期均有所提前；玉米成熟期变化情况更为复杂。2030s黄淮海棉麦油烟果区和云贵高原稻玉米烟草区的玉米生育期则将缩短3—6 d，其它各区域的玉米生育期均表现为延长趋势。

References

[1] 林而达. 气候变化的危险水平与可持续发展的适应能力建设. 气

候变化研究进展, 2005, 1(2): 76-79.

Lin E D. Damage level of climate change and building the adaptive capacity for sustainable development. Advances in Climate Change Research, 2005, 1(2): 76-79. (in Chinese)

[2] 缪启龙, 丁园圆, 王 勇. 气候变暖对中国亚热带北界位置的影响.

地理研究, 2009, 28(3): 634-642.

Miu Q L, Ding Y Y, Wang Y. The effect of climate warming on the northern boundary of Chinese subtropical. Geographical Research, 2009, 28(3): 634-642. (in Chinese)

[3] 巢清尘, 胡国权, 赵宗慈.气候变化的风险、挑战与决策. 气候变化

研究进展, 2009, 5(4): 249.

Chao Q C, Hu G Q, Zhao Z C. Climate change risk, challenges and decision. Advances in Climate Change Research, 2009, 5(4): 249. (in Chinese)

[4] 李 爽, 王 羊, 李双成. 中国近30年气候要素时空变化特征. 地理

研究, 2009, 28(6): 1593-1605.

Li S, Wang Y, Li S C. Spatio-temporal patterns of climatic change in China in recent 30 years. Geographical Research, 2009, 28(6): 1593-1605. (in Chinese)

[5] 崔读昌.气候变暖对我国农业生产的影响与对策. 中国农业气象,

1992, 13(2): 16-20.

Cui D C. The possible influence of climate getting warmer on agricultural production in China. China Journal of Agro-meteorology, 1992, 13(2): 16-20. (in Chinese)

[6] 赵俊芳, 杨晓光, 刘志娟. 气候变暖对东北三省春玉米严重低温冷

害及种植布局的影响. 生态学报, 2009, 29(12):6544-6551.

Zhao J F, Yang X G, Liu Z J. Influence of climate warming on serious low temperature and cold damage and cultivation pattern of spring maize in Northeast China. Acta Ecologica Sinica, 2009, 29(12): 6544-6551. (in Chinese)

[7] 邓振镛, 张 强, 蒲金涌, 刘德祥, 郭 慧, 王全福, 赵 鸿, 王鹤龄.

气候变暖对中国西北地区农作物种植的影响. 生态学报, 2008, 28(8): 3760-3768.

Deng Z Y, Zhang Q, Pu J Y, Liu D X, Guo H, Wang F Q, Zhao H, Wang H L. The impact of climate warming on crop planting and production in northwestern China. Acta Ecologica Sinica, 2008, 28(8): 3760-3768. (in Chinese)

[8] Evans L T. Crop Evolution, Adaptation and Yield. London: Cambridge

University Press, 1993:259.

[9] Watson R T, Albritton D L, Barker T. Climate Change 2001——

Synthesis Report. United States of America: Cambridge University Press, 2001:46.

[10] Menzel A, Fabian P. Growing season extended in Europe. Nature,

1999, 397: 659.

[11] Myneni R B, Keeling C D, Tucker C J, et al. Increased plant growth in

the northern high latitudes from 1981-1991. Nature, 1997, 386: 698-702.

[12] 徐铭志, 任国玉. 近40年中国气候生长期的变化. 应用气象学报,

2004, 15(3): 306-312.

Xu M Z, Ren G Y. Change in growing season over China: 1961-2000. Journal of Applied Meteorology Science, 2004, 15(3): 306-312. (in Chinese)

[13] 薛昌颖, 刘荣花, 吴 骞. 气候变暖对信阳地区水稻生育期的影响.

中国农业气象,2010,31(3):353-357.

Xue C Y, Liu R H, Wu Q. Effect of climate warming on rice growing stages in Xinyang. Chinese Journal of Agrometeorology, 2010, 31(3): 353-357. (in Chinese)

[14] 张明捷, 王运行, 赵桂芳. 濮阳冬小麦生育期气候变化及其对小麦

产量的影响. 中国农业气象, 2009, 30(2): 223-229.

Zhang M J, Wang Y H, Zhao G F. Climate change during winter wheat growing period and its impacts on winter wheat yield in Puyang of Henan province. Chinese Journal of Agro-meteorology, 2009, 30(2): 223-229. (in Chinese)

[15] 张建立. 气候因子对豫南夏玉米生长发育的影响. 河南农业科学,

2011, 40(1): 54-57.

Zhang J L. Effects of climatic factors on growth and development of maize in southern Henan province. Journal of Henan Agricultural Sciences, 2011, 40(1): 54-57. (in Chinese)

[16] Dowswell C R, Paliwal R L, Cantrell R P. Maize in the Third World.

USA: West view Press, 1996:35-37.

[17] 闫洪奎, 杨 镇, 徐 方. 玉米生育期和生育阶段的纬度效应研究.

中国农学通报, 2010, 26(12): 324-329.

Yan H K, Yang Z, Xu F. Effect of latitude on growth period and phase in maize. Chinese Agricultural Science Bulletin, 2010, 26(12):

22期 翟治芬等：中国玉米生育期变化及其影响因子研究 4603

324-329. (in Chinese)

[18] 陈学君, 曹广才, 贾银锁, 吴东兵, 陈 婧, 于亚雄, 李 唯, 李 杰.

玉米生育期的海拔效应研究. 中国生态农业学报, 2009, 17(3): 527-532.

Chen X J, Cao G C, Jia Y S, Wu D B, Chen J, Yu Y X, Li W, Li J. Influence of elevation on growth duration of maize. Chinese Journal of Eco-Agriculture, 2009, 17(3): 527-532. (in Chinese)

[19] 崔读昌. 中国主要农作物气候资源图集. 北京:气象出版社, 1984.

Cui D C. Climate resources atlas of main crops in China. Beijing:Meteorological Press, 1984. (in Chinese)

[20] The Intergovernmental Panel on Climate Change. Climate change

2007: The physical science basis, summary for policymakers. (2007. 03—07). http: // www. ipec. Ch/.

[21] 中国农业科学院农业区划研究所. 中国农业区划图集. 北京: 农业

出版社, 1999.

Institute of Agricultural Regional Planning, Chinese Academy of Agricultural Science. Agriculture Regional Atlas of China. Beijing: Agriculture Press, 1999. (in Chinese)

[22] 中国气象局. 太阳能资源评估方法. 中华人民共和国气象行业标

准，QX/T—2006.

China Meteorological Administration. Solar resource assessment methods. The People's Republic of Meteorological industry standards, QX/T-2006. (in Chinese)

[23] 李丽娟, 王 娟, 李海滨. 无定河流域降雨量空间变异性研究. 地

理研究, 2002, 21(4): 434-440.

Li L J, Wang J, Li H B. Analysis of spatial variability of rainfall in

Wuding river basin. Geographical Research, 2002, 21(4): 434-440. (in Chinese)

[24] 黄 英. 统计学方法与应用. 北京: 地质大学出版社, 2011.

Huang Y. Statistical Methods and Applications. Beijing: Geological University Press, 2011. (in Chinese)

[25] 张建同, 孙昌言, 王世进. 应用统计学. 北京: 清华大学出版社,

2010.

Zhang J T, Sun C Y, Wang S J. Applied Statistics. Beijing: Tsinghua University Press, 2010. (in Chinese)

[26] Pinker R T, Zhang B, Dutton E G. Do satellites detect trends in surface

solar radiation. Science, 2005, 308(5723): 850-854.

[27] 杨晓光, 刘志娟, 陈 阜. 全球气候变暖对中国种植制度可能影响

Ⅰ.气候变暖对中国种植制度北界和粮食产量的可能影响分析.中国农业科学, 2010, 43(2): 329-336.

Yang X G, Liu Z J, Chen F. The possible effects of global warming on cropping systems in ChinaⅠ.The possible effects of climatic warming on northern limits of cropping systems and crop yields in China. Scientia Agricultura Sinica, 2010, 43(2): 329-336. (in Chinese) [28] 马柱国, 黄 刚, 甘文强, 陈明林. 近代中国北方干湿变化趋势的

多时段特征. 大气科学, 2005, 29(5): 671-681.

Ma Z G, Huang G, Gan W Q, Chen M L. Multi-scale temporal characteristics of the dryness/wetness over northern China during the last century. Chinese Journal of Atmospheric Sciences, 2005, 29(5): 671-681. (in Chinese)

（责任编辑 郭银巧）