1-例冬小麦和玉米生产对气候变化课件

1、研究背景和目标研究背景和目标冬小麦和玉米产量与关键气候因子时空分布特征冬小麦和玉米产量与关键气候因子时空分布特征冬小麦和玉米产量对区域气候变化的响应分析冬小麦和玉米产量对区域气候变化的响应分析主要结论主要结论1.1 研究背景 分析探讨我国气候变化对粮食作物产量的可分析探讨我国气候变化对粮食作物产量的可能影响能影响,可为解决我国粮食安全问题提供科学依可为解决我国粮食安全问题提供科学依据，对实现粮食生产的持续发展具有十分重大的据，对实现粮食生产的持续发展具有十分重大的现实意义；现实意义；农业产量气象预报农业产量气象预报按预报时效：按预报时效：年景预报、趋势预报和定量预报三种。年景预报、趋势预报和定

2、量预报三种。 a a、年景预报。根据对未来天气气候的长期预测年景预报。根据对未来天气气候的长期预测结果，结合作物的生长发育过程及其对气象条件的需结果，结合作物的生长发育过程及其对气象条件的需求，评述气象条件对作物生长和产量形成的利弊影响求，评述气象条件对作物生长和产量形成的利弊影响而编制的农业气象预报。一般在作物播种前后发布，而编制的农业气象预报。一般在作物播种前后发布，它对指导种植计划有非常重要的意义。它对指导种植计划有非常重要的意义。 b b、趋势预报。根据作物播种前后至收获前两个月趋势预报。根据作物播种前后至收获前两个月的气象条件及其对作物生长发育的利弊影响和对未来的气象条件及其对作物生

3、长发育的利弊影响和对未来天气气候的预测及影响评价而对作物产量丰歉趋势进天气气候的预测及影响评价而对作物产量丰歉趋势进行预测的一种农业气象预报。一般在行预测的一种农业气象预报。一般在收获前两个月收获前两个月发发布，对指导农业生产和田间管理有一定的现实意义。布，对指导农业生产和田间管理有一定的现实意义。 c c、定量预报。根据作物播种前后至收获前一个月定量预报。根据作物播种前后至收获前一个月的气象条件及其对作物生长发育的利弊影响和对未的气象条件及其对作物生长发育的利弊影响和对未来天气气候的预测及影响评价而对作物产量进行预来天气气候的预测及影响评价而对作物产量进行预测、预报的一种农业气象预报。一般在

4、测、预报的一种农业气象预报。一般在收获前一个收获前一个月月发布，对指导粮食储运、流通和消费有一定的参发布，对指导粮食储运、流通和消费有一定的参考价值。考价值。 按预报对象：按预报对象：粮食作物产量预报、经济作物产量预粮食作物产量预报、经济作物产量预报和名特优小宗农产品产量预报报和名特优小宗农产品产量预报三种。三种。 a、粮食作物产量预报。粮食作物产量预报在我粮食作物产量预报。粮食作物产量预报在我国气象部门开展得比较广泛，从国家到省、地、县国气象部门开展得比较广泛，从国家到省、地、县基本都已投入业务运行。其内容主要包括粮食总产基本都已投入业务运行。其内容主要包括粮食总产量预报、小麦产量预报、玉米

5、产量预报、双季早稻量预报、小麦产量预报、玉米产量预报、双季早稻产量预报、中（一季）稻产量预报、双季晚稻产量产量预报、中（一季）稻产量预报、双季晚稻产量预报、大豆产量预报等。预报、大豆产量预报等。 b、经济作物产量预报。经济作物产量预报主要包经济作物产量预报。经济作物产量预报主要包括棉花产量预报、油菜产量预报和花生产量预报等。括棉花产量预报、油菜产量预报和花生产量预报等。目前我国多数省、区气象部门都开展此项业务，目前我国多数省、区气象部门都开展此项业务， c、名特优小宗农产品产量预报。主要包括绿豆、名特优小宗农产品产量预报。主要包括绿豆、红小豆、板栗、金丝小枣等农产品产量预报。此项红小豆、板栗、

6、金丝小枣等农产品产量预报。此项预报是国内个别省、区、县气象部门根据用户的特预报是国内个别省、区、县气象部门根据用户的特殊要求而开展的非常规业务预报。殊要求而开展的非常规业务预报。 1.2 研究目标v分析分析研究区域研究区域1982198219981998年年冬小麦、玉米冬小麦、玉米产量和气候因子的时空产量和气候因子的时空分布特征，分布特征，以此为基础，以此为基础，分别建立代表区域冬小麦和玉米的产量分别建立代表区域冬小麦和玉米的产量与其对应生长期内积温和降水因子的拟和关系式；与其对应生长期内积温和降水因子的拟和关系式；v建立较大范围的作物产量和气象因子关系，结合模拟的未来气候建立较大范围的作物产

7、量和气象因子关系，结合模拟的未来气候变化情景，对未来粮食作物产量进行预测，并据此分析未来区域变化情景，对未来粮食作物产量进行预测，并据此分析未来区域气候变化对我国粮食安全的影响气候变化对我国粮食安全的影响。 1.3 研究思路 1.4数据来源与分析，研究区域选择表表2-4 2-4 本文研究区域的选择本文研究区域的选择 v作物数据：作物数据：1981982 219981998年研究区域作物产量（总产、单产）；全国小麦、年研究区域作物产量（总产、单产）；全国小麦、玉米玉米1941949 9 2004 2004的总产、单产和播种面积。（来源：中国种植业信息网的总产、单产和播种面积。（来源：中国种植业信

8、息网农作物数据库、农作物数据库、FAOFAO、省市年鉴、网络）省市年鉴、网络）v气象数据：气象数据：1981982 219981998年与研究区域对应站点的日均温和日降水。（来源：年与研究区域对应站点的日均温和日降水。（来源：WMO)WMO)v因为本文研究的是我国冬小麦和玉米产量对气候变化的响应，结合这两种因为本文研究的是我国冬小麦和玉米产量对气候变化的响应，结合这两种作物的栽培和所掌握的数据情况，选取研究区域如下表：作物的栽培和所掌握的数据情况，选取研究区域如下表：北部冬小麦区 黄淮海平原冬小麦区 西南冬小麦区 北京 济南 青岛 徐州 西安 成都 东北春玉米区 黄淮海夏玉米区 西南玉米区 哈

9、尔滨 长春 沈阳 北京 济南 青岛 徐州 西安 成都 2 冬小麦、玉米产量与气候因子的时空分布冬小麦、玉米产量与气候因子的时空分布 北京5.9229.453.003.504.004.505.005.506.006.50-10-5051015单产变化率青岛6.29414.943.003.504.004.505.005.506.006.50-15-10-505101520济南5.52412.923.003.504.004.505.005.506.00单产/103kg.ha.-1-15-10-5051015变化率/%徐州5.41415.613.003.504.004.505.005.506.00-4

10、0-30-20-1001020西安4.61742.732.002.503.003.504.004.505.00-30-20-1001020304050成都4.819108.092.002.503.003.504.004.505.0019821983198419851986198719881989199019911992199319941995199619971998199920002001200220032004年份-100-50050100150n 各研究区域冬小麦单产在1997年前总体大幅度增长，1997年以后总体减产。n 单产区域性差异较大，平均单产以黄淮海产区的青岛、北京和济南为较高产

11、区，徐州次之，西安产区最低，成都产区次低。北京减产幅度最大，青岛、济南以及西安单产较为稳定。 图图46 研研究究区区域域冬冬小小麦麦单单产产年年际际变变化化2.1冬小麦产量的时空分布特征图图414 东北研究区域春玉米单产量年际变化东北研究区域春玉米单产量年际变化 成都4.634327.750.501.001.502.002.503.003.504.004.505.0019821983198419851986198719881989199019911992199319941995199619971998199920002001200220032004年份单产/103kg ha-1-100-500

12、50100150200250300350变化率/%单产变化率图图415 西南研究区玉米单产年际变化西南研究区玉米单产年际变化 2.2 玉米产量的时空分布特征n 研究区玉米单产在19961998年之间总体达到历史单产新高，1998以后总体有下降趋势。 图图4|16 黄黄淮淮海海平平原原研研究究区区夏夏玉玉米米单单产产年年际际变变化化 n 单产区域性差异较大，平均单产以东北产区最高（长春、沈阳和哈尔滨）；以黄淮海产区青岛、北京和济南为较高，徐州次之；西南（成都）产区最低；西安产区次低。而且东北产区单产整体比较稳定，北京减产幅度最大，青岛、徐州、济南以及西安产量较为稳定，成都玉米单产增长最为稳定。2

13、.2 玉米产量的时空分布特征n 1982 1998年冬小麦和玉米生长期内对应站点气温和降水的变化同样存在时间上和空间上的差异。n 在小麦和玉米生长期内，积温年际变化总体有微弱增加趋势，变暖的季节主要出现在冬季，这与北半球及中国气温变化的特点是一致的。n 各站点冬小麦生长期内降雨量呈递减趋势，玉米生育期的降雨量在地理位置上随着纬度的降低，有明显增加趋势。 2.3 气候因子的时空分布特征3 冬小麦和玉米生产对气候变化的响应冬小麦和玉米生产对气候变化的响应v 产量分离方法v研究区冬小麦和玉米气候单产提取以及和气候因子的拟合关系v全国冬小麦和玉米产量和气候因子的关系v我国冬小麦和玉米生产对气候变化的响

14、应预测n 趋势单产定义趋势单产定义 iy( )iy t( )iy w3.1 产量分离方法 气候单产，表征气候波动对实际单产的贡献，随气候波动气候单产，表征气候波动对实际单产的贡献，随气候波动而增加或减少而增加或减少 。 实际单产；实际单产； 趋势单产，表征社会经济条件和生产力技术水平对实际单产趋势单产，表征社会经济条件和生产力技术水平对实际单产的贡献，随科技进步和社会经济条件改善呈现上升趋势，逐渐逼的贡献，随科技进步和社会经济条件改善呈现上升趋势，逐渐逼近一个极限值；近一个极限值；n趋势单产分离公式：5a滑滑动动平均和平均和5a线线性性调调和滑和滑动动平均平均 111( )pjjkijkiky

15、tcy( )jiy t( )jiyt1j kc 1j kiy 3.1 产量分离方法第第j 年趋势单产滑动平均值，即第年趋势单产滑动平均值，即第j 年趋势单产年趋势单产 ；滑动平均系数；滑动平均系数；滑动步长年间第（滑动步长年间第（j-k+1）年的实际单产；年的实际单产；P 滑动平均时间步长。滑动平均时间步长。jj-p+1pj-k+111(1)pj kkcpkk 11111pjkkjkcc1j kc 1jkc3.1 产量分离方法n 5a调和调和滑滑动动平均系数平均系数 满满足下式足下式 滑动平均步长滑动平均步长p区间内，离历史年份区间内，离历史年份j越近年份的单产信息越近年份的单产信息对对j年理

16、论趋势单产越大年理论趋势单产越大 。- 0 . 3- 0 . 2- 0 . 100 . 10 . 20 . 30 . 40 . 50 . 6195319561959196219651968197119741977198019831986198919921995199820012004年 份气候单产/103kg ha-1滑 动调 和3.1 产量分离方法全国小麦全国玉米n 对分离出的气候单产对分离出的气候单产 与对应站点冬小麦生育期内积温与对应站点冬小麦生育期内积温T与降雨与降雨总量总量P之间的关系进行二次多项式回归拟合分析。之间的关系进行二次多项式回归拟合分析。 式中式中：C为常数项，为常数项，

17、a、b、c、d、e为气候因子系数。为气候因子系数。wiy22()()wiyCaTbPcTPdTeP 3.2 研究区冬小麦和玉米气候单产与气候因子的拟合研究区冬小麦和玉米气候单产与气候因子的拟合n 徐州冬小麦气侯单产两种方法的提取效果均比较显著，二次多项式拟合效果徐州冬小麦气侯单产两种方法的提取效果均比较显著，二次多项式拟合效果良好；良好；济南济南冬小麦气侯单产经过冬小麦气侯单产经过5a线性调和滑动平均处理提取后的回归拟合线性调和滑动平均处理提取后的回归拟合效果显著；所有区域的气候单产拟合效果总体不好。效果显著；所有区域的气候单产拟合效果总体不好。注：*为5a滑动平均处理；*为5a线性调和滑动平

18、均处理；样本数n = 13，即19861998年时间序列。 表表53 19861998年冬小麦研究区气候单产二次多项式拟合分析表年冬小麦研究区气候单产二次多项式拟合分析表 3.2 研究区冬小麦和玉米气候单产与气候因子的拟合研究区冬小麦和玉米气候单产与气候因子的拟合n 原因分析原因分析在气象单产拟合年限内，中国冬小麦产量已经基本处于历史的高产水平，生产力技在气象单产拟合年限内，中国冬小麦产量已经基本处于历史的高产水平，生产力技术水平的大幅度提高使得趋势单产成为实际单产的主要贡献因子，几乎贡献了术水平的大幅度提高使得趋势单产成为实际单产的主要贡献因子，几乎贡献了95%以上，所以各站点的气候单产贡献

19、整体微弱，直接表现为趋势单产和气候因子的相以上，所以各站点的气候单产贡献整体微弱，直接表现为趋势单产和气候因子的相关系数很小。关系数很小。气候变化的单产响应区域差异很大，雨热季节分布不一，灌溉条件状况各异等，造气候变化的单产响应区域差异很大，雨热季节分布不一，灌溉条件状况各异等，造成站点气候单产对实际单产贡献不一。成站点气候单产对实际单产贡献不一。样本空间小（样本空间小（n=16n=16），），而且经过滑动处理，发生样本萎缩（而且经过滑动处理，发生样本萎缩（n=13n=13），），不能很好地体不能很好地体现气候时间趋势性，因而波动性气候单产提取效果欠佳，在某些站点尤其明显。现气候时间趋势性，因

20、而波动性气候单产提取效果欠佳，在某些站点尤其明显。由于以上各因素，各区域冬小麦气候单产的提取方法不具有普适性。理论上，可以由于以上各因素，各区域冬小麦气候单产的提取方法不具有普适性。理论上，可以通过改善样本空间状况通过改善样本空间状况( (实际增加样本数，或者通过信息膨胀处理增加样本数，如非实际增加样本数，或者通过信息膨胀处理增加样本数，如非线性信息扩散方法等线性信息扩散方法等) )、样本品质（提高样本数据的统计精确性）以及拟合方法等实、样本品质（提高样本数据的统计精确性）以及拟合方法等实现比较理想的趋势单产提取和模拟。现比较理想的趋势单产提取和模拟。3.2 研究区冬小麦和玉米气候单产与气候因

21、子的拟合研究区冬小麦和玉米气候单产与气候因子的拟合n 济南和青岛玉米气侯单产两种方法的提取效果均比较显著，二次多项式拟合效果良好；西安玉米气侯单产经过5a滑动平均处理提取后的回归拟合效果显著；其它区域的气候单产拟合效果总体不好。具体原因与小麦气候单产拟合效果欠佳的原因比较相似。 表表59 19861998年研究区玉米气象单产二次多项式拟合分析表年研究区玉米气象单产二次多项式拟合分析表 注：*为5a滑动平均处理；*为5a线性调和滑动平均处理；样本数n = 13，即19861998年序列。 3.3 研究区冬小麦和玉米实际单产与气候因子的拟合研究区冬小麦和玉米实际单产与气候因子的拟合表55 1983

22、1998年冬小麦研究区实际单产二次多项式拟合分析表n 所选冬小麦区实际单产数据系列和对应生长期内积温和降雨总量记录系列具有良好的二次多项式拟合关系。徐州拟合效果稍差。可以利用这些拟合二次多项关系式进行未来以上站点局域气候变化情景下的冬小麦单产预测。 注：样本数n = 16，即19831998年时间序列。 3.3 研究区冬小麦和玉米实际单产与气候因子的拟合研究区冬小麦和玉米实际单产与气候因子的拟合表511 19831998年研究区玉米实际单产二次多项式拟合分析表注：样本数n = 17，即19821998年序列。n 所选区域玉米实际单产数据系列和对应站点生长期内积温和降雨总量记录系列具有良好的二次

23、多项式拟合关系。西安拟合效果稍差。可以利用这些拟合二次多项关系式进行未来以上站点局域气候变化情景下的玉米单产预测。 3.4 较大区域冬小麦和玉米单产与气候因子的拟合大区域冬小麦和玉米单产与气候因子的拟合表57 19831998年研究区冬小麦平均单产二次多项式拟合分析表n 总体上，所选冬小麦研究区平均实际单产和气候单产与气象数据具有良好的二次多项式拟合关系，可以用于所选冬小麦产区的区域气候波动情景下的单产预测。 3.4 较大区域冬小麦和玉米单产与气候因子的拟合大区域冬小麦和玉米单产与气候因子的拟合表513 19821998年研究区玉米平均单产二次多项式拟合分析表n 东北春玉米产区和黄淮海夏玉米产

24、区气候单产与历史记录具有良好的二次多项式拟合关系，可以用于所选玉米产区的区域气候波动情景下的单产预测。 3.5 冬小麦和玉米生产对气候变化的响应预测温度变化（） 降水变化（） 年份 冬 春 夏 秋 年平均 冬 春 夏 秋 年平均 2010 0.38 0.35 0.34 0.35 0.35 1.4 1.3 1.2 0.9 1.1 2020 0.69 0.63 0.62 0.64 0.65 2.6 2.3 2.1 1.6 1.9 2030 0.93 0.86 0.84 0.87 0.88 3.5 3.2 2.8 2.2 2.6 2040 1.12 1.04 1.02 1.05 1.06 4.2 3

25、.8 3.4 2.5 3.2 2050 1.48 1.37 1.34 1.38 1.40 5.5 5.1 4.5 3.5 4.2 GCM模型模拟的模型模拟的2010到到2050年气候变化情景年气候变化情景 n 未来到2050年中国地区由于温室效应气温可能上升1.4，降水总的趋势也是增加的，但增加幅度不大，到2050年可能只增加4.2%；降水增加幅度小于气温升高幅度。3.5 冬小麦和玉米生产对气候变化的响应预测3.5 冬小麦和玉米生产对气候变化的响应预测注：*东北春玉米生长期从5月1日至9月10日，生长期按133d计算；黄淮海夏玉米生长期从6月10日至9月30日，生长期按112d计算；西南秋玉米

26、生长期从8月1日至11月10日，生长期按101d计算。3.5 冬小麦和玉米生产对气候变化的响应预测表516 基于GCM模拟2010到2050年气候变化情景的中国冬小麦产量响应n 中国冬小麦单产对未来区域气候波动，即气温升高，雨量增加的响应行为总体为增产趋势，但是这种增长是具体的，具有明显的时空分布特征，并且也是复杂的，受非气候因子波动的综合影响。 n中国玉米单产对未来区域气候波动，即气温升高，雨量增加的响应行为总体为减产效应，明显削弱了生产技术水平提高对实际单产的贡献，总体表现为微弱的增长趋势。3.5 冬小麦和玉米生产对气候变化的响应预测表517 GCM模拟2010到2050年气候变化情景的中

27、国玉米产量响应 3.6 冬小麦和玉米生产安全分析n区域气候波动下的中国冬小麦产量分析 如以所选区域平均处理的拟合二项式为例，2010年、2020年、2030年、2040年、2050年气候单产（式5-7）预测值分别为250kg/ha、121 kg/ha、-16 kg/ha、-145 kg/ha、-433 kg/ha，对应时间全国的气候单产yw（式5-8）预测值分别为169kg/ha、188 kg/ha、206 kg/ha、223 kg/ha、256kg/ha。而对应时间的研究区平均处理的实际单产（式5-9）预测值分别为4909kg/ha、5117 kg/ha、5284kg/ha、5416kg/h

28、a、5666kg/ha，表现为稳定的增长趋势。3.6 冬小麦和玉米生产安全分析n 冬小麦生产安全分析 若根据式5-9的预测实际单产，未来50年内如果中国小麦种植面积保持1995年 2004年的平均水平2.66005107 ha，即使考虑未来人口增长，则2010年、2020年、2030年、2040年、2050年人均小麦粮食占有量由2004年的70.2kg增长到94.6 kg、98.6 kg、101.9 kg、104.4 kg、109.2 kg。 如果未来50年总产保持1995年2004年的平均水平，则2010年、2020年、2030年、2040年、2050年节省小麦种植耕地面积比例分别达21.7

29、5%、24.92%、27.30%、29.07%、32.20%。 在未来50年模拟情景气候变化条件下，即使考虑未来50年中国人口的增长，中国小麦产量总体保持稳定增长，人均小麦可供应量仍保持稳定提高。 总体看来，随着生产技术水平的稳定提高和种植面积的保证，未来中国小麦产量的气候波动响应不会对中国粮食生产安全产生负效应。3.6 冬小麦和玉米生产安全分析n 区域气候波动下的中国玉米产量分析 如以全国平均处理的拟合二项式为例，2010年、2020年、2030年、2040年、2050年气候单产预测值yw（式5-10）分别为-48.0kg/ha、-173.4 kg/ha、-270.8 kg/ha、-34.3

30、kg/ha、-469.9 kg/ha。而对应时间的东北春玉米产区站点平均处理的实际单产（式5-12）预测值分别为6200kg/ha、6265kg/ha、6339kg/ha、6411kg/ha、6574kg/ha，为高增产区；而对应时间的黄淮海夏玉米产区站点平均处理的实际单产（式5-13）预测值分别为5343kg/ha、5426kg/ha、5499kg/ha、5554kg/ha、5671kg/ha，为中低增产区，各预测区表现为不同程度的稳定的增长趋势。 3.6 冬小麦和玉米生产安全分析n 玉米生产安全分析 若根据式5-12的预测实际单产，未来50年内如果中国小麦种植面积保持1995年 2004年的平均水平2.436801107ha ，即使考虑未来人口增长，则2010年、2020年、2030年、2040年、2050年人均小麦粮食占有量由2004年的99.5kg增长到109