免费文献传递   相关文献

Simulation of maize production under climate change scenario in Northeast China

气候变化情景下东北地区玉米产量变化模拟



全 文 :中国生态农业学报 2008 年 11 月 第 16 卷 第 6 期
Chinese Journal of Eco唱Agriculture,Nov.2008,16(6):1448唱1452   
DOI: 10.3724/SP.J.1011.2008.01448
 
气候变化情景下东北地区玉米产量变化模拟倡 倡
张建平1,2 赵艳霞3 王春乙3 杨晓光1 何 勇4
(1.中国农业大学资源与环境学院 北京 100193; 2.重庆市气象科学研究所 重庆 401147;
3.中国气象科学研究院 北京 100081; 4.国家气候中心 北京 100081)
摘 要 WOFOST作物模型在东北地区玉米适应性验证的基础上,结合气候模型 BCC唱T63 输出的未来 60 年
(2011 ~2070 年)气候情景资料,模拟分析了未来气候变化情景下我国东北地区玉米生育期和产量变化情况。
模拟结果显示:未来气候变化情景下,玉米生育期将缩短,其中,中熟玉米平均缩短 3.4 d,晚熟玉米平均缩短
1.1 d;玉米产量将相应下降,中熟玉米平均减产 3.5%,晚熟玉米平均减产 2.1%。
关键词 气候变化 玉米 生育期 产量 模拟研究 WOFOST作物模型
中图分类号:P467/S165  文献标识码:A  文章编号:1671-3990(2008)06-1448-05
Simulation of maize production under climate change scenario in Northeast China
ZHANG Jian唱Ping1,2, ZHAO Yan唱Xia3, WANG Chun唱Yi3, YANG Xiao唱Guang1, HE Yong4
(1.College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193, China;
2.Chongqing Institute of Meteorological Sciences, Chongqing 401147, China; 3.Chinese Academy of
Meteorological Sciences, Beijing 100081, China; 4.National Climate Center, Beijing 100081, China)
Abstract Based on adaptability test of WOFOST crop model for maize production in Northeast China and using a 60 -year cli唱
mate scenario data (2011 ~2070) produced by BCC -T63 model, maize growth and yield change umder climate change scenar唱
io in Northeast China was simulated and analyzed.The simulation results show a shortened maize growth cycle.The growth cycle
of medium唱mature maize shortens by 3.4 days and that of late唱mature maize shortens by 1.1 days.Yield of maize also decreases
by different degrees; medium唱mature maize decreases by 3.5% while late唱mature maize decreases by 2.1%.
Key words Climate change, Maize, Growth cycle, Yield, Simulation, WOFOST crop model
(Received Feb.13, 2008; accepted June 2, 2008)
  东北地区作为我国最大的粮食生产基地之一,
一直以来都备受诸多学者和专家的关注。 金之庆
等 [1,2]评价了全球气候变化对东北玉米生产造成的
可能影响,并初步提出相应的粮食生产对策;尚宗
波 [3]分析了全球气候变化对沈阳地区春玉米生产
的可能影响;吴金栋等[ 4]和马树庆[ 5]针对东北的粮
食生产也作了相关的研究。 所有这些研究成果对
东北地区的农业生产都有不同程度的指导意义。
农业生产很大程度上受制于环境条件,特别是
气候条件的变化。 有关研究指出[ 6] ,在未来 100 年
内,中国平均地表温度还将上升 1.1 ~5.2 ℃,而东
北地区作为我国未来百年增温幅度最大的地区之
一,到本世纪末,增温幅度分别达到 6.1 ℃(A2 排
放情景)和 4.5 ℃(B2 排放情景),可见,东北地区
未来的农业生产面临着严峻的气候变化的挑战。
鉴于此,本文主要以东北玉米为主要研究对象,探
讨分析了未来 60 年(2011 ~2070 年)玉米生育期和
产量的变化趋势,以便为相关部门提供决策依据。
1 研究区概况及所需资料
研究区为东经 120°~135°、北纬 40°~50°的东
北区域,主要以黑龙江、吉林和辽宁 3 省为研究对
象,地跨寒温带、中温带和暖温带,属温带大陆性季
风气候。 东北地区北临北半球冬季的寒极———东
西伯利亚,冬季强大的冷空气南下,使之成为同纬
度各地中最寒冷的地区,与同纬度的其他地区相
比,温度一般低 15℃左右。 夏季受低纬度海洋湿热
气流影响,气温则高于同纬度各地。 因此,东北地
区年温差大大高于同纬度各地。 该区域大部分地
倡 国家科技支撑计划课题(2006BAD04B09)资助
张建平(1975 ~) ,男,博士研究生,主要从事气候变化及生态与农业气象研究。 E唱mail:jeepjohn@163.com
收稿日期:2008-02-13 接受日期:2008-06-02
第 6期 张建平等:气候变化情景下东北地区玉米产量变化模拟
区的年降水量为 400 ~700 mm,东部降水集中在
5 ~9月,西部则集中在 6 ~8 月,70% ~80%的降水
集中在 6 ~7 月。 主要农作物包括春玉米、水稻、大
豆和甜菜,是我国主要的粮食生产基地。
研究所需资料主要包括作物资料、气象资料和
土壤资料[7,8] 。 作物资料有发育期观测资料和田间
生物量试验资料,来源于农气试验站、农气观测站以
及农气表-1;气象资料包括最高气温、最低气温、日
照时数、水汽压、风速和降雨量(均为日值),收集了
东北地区 81 个基本气象站 1961 ~2004 年历年逐日
气象要素,来源于国家气象信息中心(部分气象站点
分布见图 1);土壤资料有土壤质地类型、土壤水分含
量、田间持水量、土壤容重及凋萎湿度,来源于中国土
壤图集、中国土种志以及中国土壤质地类型图。
图 1 气象站点分布图
Fig.1 Distribution chart of meteorological stations
2 作物模型简介
WOFOST作物模型[9 -13]由荷兰瓦赫宁根大学开
发,起源于世界粮食研究中心(CWFS)组织,是一个
根据气象和土壤条件模拟作物根、茎、叶和穗生物量
以及土壤水分的动态模型。 模型的主要用途包括作
物产量预测与评估、土地资源定量评价、气候变化影
响的量化评估、风险分析和年际间产量变化分析等。
2.1 模型的基本框架
WOFOST模型用于模拟 1 年生作物的生态生
理过程,主要包括光合作用、呼吸作用等,并描述这
些过程如何受环境的影响。 模型主要包括 3 大模
块,分别是作物模块、气候模块和土壤模块。 根据 3
大模块的计算可以得到潜在生产力 (没有限制水
平)和有限制水平的生产力(图 2)。
图 2 WOFOST模型的基本框架
Fig.2 Structure diagram of the WOFOST model
2.2 模型中主要子模型
2.2.1 发育子模型
玉米生长发育模型采用“积温法”,公式如下:
D r, t = TeiTSUM j    ( j =1,2) (1)
T ei =
0      (T i≤Tb)
T i -Tb   (Tb <T i <Tmax,e)    (2)
Tmax, e -Tb   (T i≥Tmax,e)
式中, D r, t 为 t时刻的发育速率(d -1), Tei 为有效温
度(℃), TSUM j( j =1,2) 为完成某一发育阶段所
需的有效积温 (℃· d),T i 为日平均气温 (℃),取
日最高气温和最低气温的平均值, Tb 为作物发育下
限温度(℃),Tmax,e为作物发育上限温度(℃)。
2.2.2 土壤水分平衡子模型
WOFOST模型中对土壤水分平衡过程的模拟
分为考虑和不考虑地下水影响两种情形。 不考虑
地下水影响时,忽略毛管上升水,只考虑降水、灌
溉、水分从地表到根区的入渗、地表到大气的蒸发、
作物蒸腾、根区水分向下的渗漏、地表径流等过程。
θt =INup +( IN low -Ta)RD Δt (3)
INup =P +Ie -E s +SS t /Δt -SR (4)
IN low =CR -Perc (5)
式中,θt 为 t时刻实际土壤含水量(cm3 · cm -3 ),INup
为上层土壤水分通量(cm· d -1),Ta 为实际蒸腾速率
(cm· d -1 ),P 为降水量 ( cm · d -1 ), Ie 为灌溉量
(cm· d -1),Es 为土壤蒸发速率(cm· d -1 ),SS t 为地
表积水(cm),SR为地表径流速率(cm· d -1 ),CR 为
9441
中国生态农业学报 2008 第 16卷
毛管上升水速率(cm· d -1 ),Perc 为渗透速率(cm·
d -1),RD为实际根长(cm),Δt为时间步长(d)。
2.3 模型适应性调整
2.3.1 作物参数的确定及区域化
采取“试错法”,根据东北地区玉米多年发育期
观测资料确定区域玉米的各阶段发育参数,玉米品
种熟性根据有效积温来划分,其他作物参数来源于
多年试验观测以及国内外玉米的研究成果。 区域
化的基本原则是在整个区域尺度上,首先进行分
区,然后以各分区的代表样点调试模型参数;基本
原理是假设在农作物种植区域尺度上,所有作物品
种都同属一种生态型;调整方案是以区域为空间尺
度,使用多点同期的面上产量资料进行调试。 表 1
列出了模型所需要的主要作物生长参数。
表 1 模型主要作物参数
Tab.1 The main crop parameters of the model
参数
Parameter
说明
Explanation
单位
Unit
取值
Value
TSUM0 播种 ~出苗需要的总热量Temperature sum from sowing to emergence ℃· d
162.2(晚熟,Late唱mature)
149.2(中熟,Medium唱mature)
TSUM1 出苗 ~抽雄需要的总热量Temperature sum from emergence to anthesis ℃· d
1 020.6(晚熟,Late唱mature)
943.1(中熟,Medium唱mature)
TSUM2 抽雄 ~成熟需要的总热量Temperature sum from anthesis to maturity ℃· d
968.8(晚熟,Late唱mature)
880.5(中熟,Medium唱mature)
TDWI 初始总干物重Initial total crop dry weight kg· hm -2 10.61
LAIEM 出苗时叶面积指数Leaf area index at emergence 0.005 9
RGRLAI 叶面积指数最大增长速率Maximum relative increase in LAI d -1 0.01
CVL 叶生长同化物转化效率Conversion efficiency of assimilate into leaves 0.72
CVO 贮存器官生长同化物转化效率Conversion efficiency of assimilate into storage organs 0.72
CVR 根生长同化物转化效率Conversion efficiency of assimilate into roots 0.72
CVS 茎生长同化物转化效率Conversion efficiency of assimilate into stems 0.69
RML 叶相对维持呼吸速率Relative maintain respiration rate of leaves kg ( CH2O) · kg
-1 · d -1 0.015
RMO 贮存器官相对维持呼吸速率Relative maintain respiration rate of storage organs kg ( CH2O) · kg
-1 · d -1 0.010
RMR 根相对维持呼吸速率Relative maintain respiration rate of roots kg ( CH2O) · kg
-1 · d -1 0.010
RMS 茎相对维持呼吸速率Relative maintain respiration rate of stems kg ( CH2O) · kg
-1 · d -1 0.010
2.3.2 生长模拟检验
为验证模型的适应性,选择东北地区有代表性
的 10 个站点进行连续 5 年的验证,由于篇幅有限,
本文只列出了部分站点的模拟验证结果(见表 2)。
从表中可以看出,玉米开花期模拟值与实测值
相差 1 ~6 d,平均相差 2.8 d;成熟期模拟值与实测
值相差 1 ~5 d,平均相差 3.3 d;产量误差范围在
9.34% ~-5.23%之间。 对模拟值与实测值进行
统计检验分析,两者无显著差异(P =0.05,n =10),
故所建模型能较好地反映东北玉米生长发育情况。
3 气候变化对东北地区玉米生长的影响
3.1 气候变化对玉米生育期的影响
根据枟气候变化国家评估报告枠研究表明:到
2070 年,东北地区年平均地表温度增温幅度为 3.8
℃(A2 情景)和 3.4 ℃(B2 情景)。 为初步分析未来
气候变化影响下我国东北地区玉米的生产情况,基
于作物模型,结合气候模型输出的气候资料,对东
北地区不同种类(晚熟和中熟)春玉米生育期进行
了 60 年(2011 ~2070 年)的数值模拟分析,模拟结
果见图 3a 和图 3b。
综合来看,2011 ~2070 年东北春玉米生长期会
有所缩短。 其中,中熟玉米生育期变化范围在
2.2 ~-8.5 d之间,平均缩短 3.4 d(图 3a);晚熟玉
0541
第 6期 张建平等:气候变化情景下东北地区玉米产量变化模拟
表 2 2002及 2003年各站点玉米发育及产量模拟验证结果
Tab.2 Verification results of maize growth and yields of different stations tested by model in 2002 and 2003
年份
Year
站号
Station
number
品种
Variety
开花期(月 -日)
Flowering period
(month唱day)
成熟期(月 -日)
Mature period
(month唱day)
生育期
Development
period ( d)
产量
Yield
( kg· hm -2 )
模拟
Simulation
实测
Observation
模拟
Simulation
实测
Observation
模拟
Simulation
实测
Observation
模拟
Simulation
实测
Observation
误差
Error
(%)
2002
50844 晚熟Late唱mature 07 -16 07 -22 09 -30 09 -27 156 153 11 371 10 400 +9.34
50742 中熟Medium唱mature 07 -28 07 -29 09 -23 09 -18 140 135 5 056 5 182 -2.43
54266 晚熟Late唱mature 07 -26 07 -27 09 -28 09 -24 152 148 7 909 7 524 +5.12
54049 中熟Medium唱mature 08 -01 07 -28 09 -22 09 -26 149 153 8 427 8 120 +3.78
2003
50844 晚熟Late唱mature 07 -29 07 -31 09 -24 09 -27 150 153 10 121 10 600 -4.52
50742 中熟Medium唱mature 08 -07 08 -08 10 -09 10 -08 156 155 4 680 4 400 +6.36
54266 晚熟Late唱mature 07 -30 08 -03 09 -29 09 -24 148 143 5 415 5 714 -5.23
54049 中熟Medium唱mature 07 -27 07 -30 09 -19 09 -20 139 140 8 468 8 500 -0.38
米生育期变化范围在 9.0 ~-7.8 d 之间,平均缩短
1.1 d(图 3b)。
3.2 气候变化对玉米产量的影响
产量是农作物生长发育状况的最终表现,本研
究利用气候模型输出的未来 60 年(2011 ~2070 年)
的气候数据,结合作物模型,就未来气候变化对东
北玉米(晚熟和中熟)产量变化的影响进行了分析
与评估。 模拟的产量变化趋势见图 3c 和图 3d。
从模拟结果可以看出,2011 ~2070 年东北地区
玉米产量总体呈下降趋势,但下降幅度不大。 中熟
玉米产量变化范围在 7.4% ~-11.4%,平均减产
3.5% (图 3c ); 晚 熟 玉 米 产 量 变 化 范 围 在
8.1% ~-10.0%,平均减产 2.1%(图 3d)。
4 结论和讨论
在气候变化大背景下,未来 60 年(2011 ~2070
年)东北玉米生长期可能会受到影响,绝大部分地
区玉米生长期可能要缩短。 气候变化导致的最终
影响是农业产量发生变化。 根据模拟结果,在今后
60 年东北玉米产量整体呈下降趋势。 中熟玉米平
均减产 3.3%,晚熟玉米平均减产 2.7%。 在一定范
围内,通过适当调整栽培措施,可减少或减缓气候
变化对农业造成减产的影响。
低温冷害一直是东北地区农作物主要的农业
气象灾害。 因此,气候变暖一定程度有利于该地区
农业的发展,但随着温度的不断上升,利会趋于向
害的方向发展,故应采取相应的农业适应对策与措
施来应对气候变暖这种大势所趋的现象。 一是调
整玉米播种期。 其他条件(玉米品种类型、生产水
平以及作物、气候和土壤参数)不变的情况下,适当
调整玉米播种期(在原有基础上提前播种 10 d)。
对东北地区进行了未来 10 年的模拟分析。 结果发
现,中晚熟玉米减产均有不同程度回升,中熟玉米
减产可回升 2.2 个百分点,晚熟玉米减产可回升不
到 1 个百分点(均为 10 年产量的平均值)。 二是改
变玉米品种类型。 马树庆[ 5]研究指出,农作物生长
季气温每上升 1 ℃,≥10 ℃的持续期将延长 7 ~8
d,积温增加 200 ℃· d 以上,农业气候带就会向北
偏移。 因此可考虑把目前东北地区主要以中晚熟
玉米种植全部改为晚熟玉米,经调整熟制后,各地
玉米所需生育期会延长,玉米产量有所增加,减产
幅度大大降低,由原来的 3.1%降为 1.8%。
模拟结果表明,通过改变玉米品种类型,其产
量会相应增加,减产会有一定程度的回升。 这也
正是几十年来东北地区非但没有受到气候变暖而
导致粮食减产的影响,反而经各地普遍采取适应
气候变暖的对策后,玉米单产成倍提高的原因
所在。
当然调整品种结构的作用也不是无限的,如气
候进一步变暖超出了品种结构调整的适应能力时,
应考虑调整种植制度。 如温度升高到对生长期直
接产生危害时,还应考虑改变作物种类。 另外,还
可考虑应用先进技术和增加农业投入的方法。
参考文献
[1] 金之庆,葛道阔,石春林, 等.东北平原适应全球气候变
1541
中国生态农业学报 2008 第 16卷
图 3 东北中熟与晚熟玉米生育期(a、b)和产量(c、d)变化
Fig.3 Changes of lifetime and yield of medium唱mature maize (a,c) and late唱mature maize(b,d) in Northeast China
  化的若干粮食生产对策的模拟研究[ J].作物学报,
2002,28(1):21 -31
[2] 金之庆,葛道阔,郑喜莲,等.评价全球气候变化对我国
玉米生产的可能影响 [ J].作物学报,1996,22 (5):
513 -524 
[3] 尚宗波.全球气候变化对沈阳地区春玉米生长的可能影
响[J].植物学报,2000,42(3):300 -305
[4] 吴金栋,王石立,张建敏.未来气候变化对中国东北地区
水热条件影响的数值模拟研究[ J].资源科学,2000,
22(6):36 -42
[5] 马树庆.气候变化对东北区粮食产量的影响及其适应性
对策[J].气象学报,1996, 54(4): 484 -492
[6] 枟气候变化国家评估报告枠编写委员会.气候变化国家
评估报告[M].北京:科学出版社,2007
[7] 中国农业年鉴编辑委员会.中国农业年鉴[M].北京:中
国农业出版社,2001 ~2005
[8] 全国土壤普查办公室.中国土种志(1 ~6 卷) [M].北
京:农业出版社,1993
[9] 张黎.基于遥感信息的水分胁迫条件下华北冬小麦生长
模拟研究[D].北京:中国气象科学研究院,2005
[10] 张建平,赵艳霞,王春乙.气候变化对南方双季稻发育
和产量影响 [ J].气候变化研究进展,2005,1 (4):
151 -156 
[11] 张建平,赵艳霞,王春乙.气候变化对我国华北地区冬
小麦发育和产量影响 [ J].应用生态学报, 2006,
17(7):1179 -1184 
[12] 张建平,李永华,高阳华,等.未来气候变化对重庆地区
冬小麦产量的影响[ J].中国农业气象,2007,28(3):
268 -270
[13] 刘布春.应用于低温冷害预报的东北玉米区域动力模
型的研究[D].北京:中国农业大学,2003
2541