我国冬小麦和玉米生产对气候变化的响应

我国冬小麦和玉米生产对气候变化的响应第1页/共36页内容提要研究背景和目标冬小麦和玉米产量与关键气候因子时空分布特征冬小麦和玉米产量对区域气候变化的响应分析主要结论第2页/共36页1.1研究背景分析探讨我国气候变化对粮食作物产量的可能影响,可为解决我国粮食安全问题提供科学依据，对实现粮食生产的持续发展具有十分重大的现实意义；尽管这方面的研究很多，但大部分是通过复杂的模型模拟，得出了各种各样的结论，但仍然缺少一个大家能普遍接受的结论。通过过去的研究，我们发现气候变化对农业生产的影响主要表现在农业产量上。但是在全国范围内通过产量和气候主要因子直接挂钩的研究还比较少。我们的目标是想找出一种简单的产量和主要气候因子之间的关系，并用这种简单关系来探讨冬小麦和玉米对气候变化的相应。第3页/共36页1.2研究目标分析研究区域1982～1998年冬小麦、玉米产量和气候因子的时空分布特征，以此为基础，分别建立代表区域冬小麦和玉米的产量与其对应生长期内积温和降水因子的拟和关系式；建立较大范围的作物产量和气象因子关系，结合模拟的未来气候变化情景，对未来粮食作物产量进行预测，并据此分析未来区域气候变化对我国粮食安全的影响。第4页/共36页1.3研究思路1.4数据来源与分析，研究区域选择表2-4本文研究区域的选择

作物数据：1982～1998年研究区域作物产量（总产、单产）；全国小麦、玉米1949～2004的总产、单产和播种面积。（来源：中国种植业信息网农作物数据库、FAO、省市年鉴、网络）气象数据：1982～1998年与研究区域对应站点的日均温和日降水。（来源：WMO)因为本文研究的是我国冬小麦和玉米产量对气候变化的响应，结合这两种作物的栽培和所掌握的数据情况，选取研究区域如下表：第5页/共36页

2冬小麦、玉米产量与气候因子的时空分布第6页/共36页

各研究区域冬小麦单产在1997年前总体大幅度增长，1997年以后总体减产。

单产区域性差异较大，平均单产以黄淮海产区的青岛、北京和济南为较高产区，徐州次之，西安产区最低，成都产区次低。北京减产幅度最大，青岛、济南以及西安单产较为稳定。

图4｜6

研究区域冬小麦单产年际变化2.1冬小麦产量的时空分布特征第7页/共36页图4－14东北研究区域春玉米单产量年际变化

图4—15西南研究区玉米单产年际变化

2.2

玉米产量的时空分布特征

研究区玉米单产在1996~1998年之间总体达到历史单产新高，1998以后总体有下降趋势。第8页/共36页图4|16黄淮海平原研究区夏玉米单产年际变化

单产区域性差异较大，平均单产以东北产区最高（长春、沈阳和哈尔滨）；以黄淮海产区青岛、北京和济南为较高，徐州次之；西南（成都）产区最低；西安产区次低。而且东北产区单产整体比较稳定，北京减产幅度最大，青岛、徐州、济南以及西安产量较为稳定，成都玉米单产增长最为稳定。2.2

玉米产量的时空分布特征第9页/共36页1982

~1998年冬小麦和玉米生长期内对应站点气温和降水的变化同样存在时间上和空间上的差异（图略）。在小麦和玉米生长期内，积温年际变化总体有微弱增加趋势，变暖的季节主要出现在冬季，这与北半球及中国气温变化的特点是一致的。各站点冬小麦生长期内降雨量呈递减趋势，玉米生育期的降雨量在地理位置上随着纬度的降低，有明显增加趋势。

2.3

气候因子的时空分布特征第10页/共36页3冬小麦和玉米生产对气候变化的响应产量分离方法研究区冬小麦和玉米气候单产提取以及和气候因子的拟和关系全国冬小麦和玉米产量和气候因子的关系我国冬小麦和玉米生产对气候变化的响应预测第11页/共36页

趋势单产定义

3.1

产量分离方法——

气候单产，表征气候波动对实际单产的贡献，随气候波动而增加或减少。——

实际单产；——

趋势单产，表征社会经济条件和生产力技术水平对实际单产的贡献，随科技进步和社会经济条件改善呈现上升趋势，逐渐逼近一个极限值；第12页/共36页趋势单产分离公式：5a滑动平均和5a线性调和滑动平均

3.1

产量分离方法——第j年趋势单产滑动平均值，即第j年趋势单产

；——滑动平均系数；——滑动步长年间第（j-k+1）年的实际单产；P——滑动平均时间步长。jj-p+1pj-k+1第13页/共36页3.1

产量分离方法5a调和滑动平均系数

满足下式

——滑动平均步长p区间内，离历史年份j越近年份的单产信息对j年理论趋势单产越大

。第14页/共36页3.1

产量分离方法全国小麦全国玉米第15页/共36页对分离出的气候单产

与对应站点冬小麦生育期内积温∑T与降雨总量∑P之间的关系进行二次多项式回归拟合分析。

式中：C为常数项，a、b、c、d、e为气候因子系数。3.2研究区冬小麦和玉米气候单产与气候因子的拟和第16页/共36页徐州冬小麦气侯单产两种方法的提取效果均比较显著，二次多项式拟合效果良好；济南冬小麦气侯单产经过5a线性调和滑动平均处理提取后的回归拟合效果显著；所有区域的气候单产拟合效果总体不好。注：*为5a滑动平均处理；**为5a线性调和滑动平均处理；样本数n=13，即1986~1998年时间序列。

表5—31986~1998年冬小麦研究区气候单产二次多项式拟合分析表

3.2研究区冬小麦和玉米气候单产与气候因子的拟和第17页/共36页原因分析在气象单产拟合年限内，中国冬小麦产量已经基本处于历史的高产水平，生产力技术水平的大幅度提高使得趋势单产成为实际单产的主要贡献因子，几乎贡献了95%以上，所以各站点的气候单产贡献整体微弱，直接表现为趋势单产和气候因子的相关系数很小。气候变化的单产响应区域差异很大，雨热季节分布不一，灌溉条件状况各异等，造成站点气候单产对实际单产贡献不一。样本空间小（n=16），而且经过滑动处理，发生样本委缩（n=13），不能很好地体现气候时间趋势性，因而波动性气候单产提取效果欠佳，在某些站点尤其明显。由于以上各因素，各区域冬小麦气候单产的提取方法不具有普适性。理论上，可以通过改善样本空间状况(实际增加样本数，或者通过信息膨胀处理增加样本数，如非线性信息扩散方法等)、样本品质（提高样本数据的统计精确性）以及拟合方法等实现比较理想的趋势单产提取和模拟。3.2研究区冬小麦和玉米气候单产与气候因子的拟和第18页/共36页

济南和青岛玉米气侯单产两种方法的提取效果均比较显著，二次多项式拟合效果良好；西安玉米气侯单产经过5a滑动平均处理提取后的回归拟合效果显著；其它区域的气候单产拟合效果总体不好。具体原因与小麦气候单产拟合效果欠佳的原因比较相似。

表5—91986~1998年研究区玉米气象单产二次多项式拟合分析表

注：*为5a滑动平均处理；**为5a线性调和滑动平均处理；样本数n=13，即1986~1998年序列。

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

所选区域玉米实际单产数据系列和对应站点生长期内积温和降雨总量记录系列具有良好的二次多项式拟合关系。西安拟合效果稍差。可以利用这些拟合二次多项关系式进行未来以上站点局域气候变化情景下的玉米单产预测。

第21页/共36页3.4较大区域冬小麦和玉米单产与气候因子的拟和表5—71983~1998年研究区冬小麦平均单产二次多项式拟合分析表

总体上，所选冬小麦研究区平均实际单产和气候单产与气象数据具有良好的二次多项式拟合关系，可以用于所选冬小麦产区的区域气候波动情景下的单产预测。第22页/共36页3.4较大区域冬小麦和玉米单产与气候因子的拟和表5—131982~1998年研究区玉米平均单产二次多项式拟合分析表

东北春玉米产区和黄淮海夏玉米产区气候单产与历史记录具有良好的二次多项式拟合关系，可以用于所选玉米产区的区域气候波动情景下的单产预测。第23页/共36页3.5冬小麦和玉米生产对气候变化的响应预测GCM模型模拟的2010到2050年气候变化情景

未来到2050年中国地区由于温室效应气温可能上升1.4℃，降水总的趋势也是增加的，但增加幅度不大，到2050年可能只增加4.2%；降水增加幅度小于气温升高幅度。第24页/共36页3.5冬小麦和玉米生产对气候变化的响应预测第25页/共36页3.5冬小麦和玉米生产对气候变化的响应预测注：*东北春玉米生长期从5月1日至9月10日，生长期按133d计算；黄淮海夏玉米生长期从6月10日至9月30日，生长期按112d计算；西南秋玉米生长期从8月1日至11月10日，生长期按101d计算。第26页/共36页3.5冬小麦和玉米生产对气候变化的响应预测表5—16基于GCM模拟2010到2050年气候变化情景的中国冬小麦产量响应中国冬小麦单产对未来区域气候波动，即气温升高，雨量增加的响应行为总体为增产趋势，但是这种增长是具体的，具有明显的时空分布特征，并且也是复杂的，受非气候因子波动的综合影响。第27页/共36页中国玉米单产对未来区域气候波动，即气温升高，雨量增加的响应行为总体为减产效应，明显削弱了生产技术水平提高对实际单产的贡献，总体表现为微弱的增长趋势。3.5冬小麦和玉米生产对气候变化的响应预测表5—17GCM模拟2010到2050年气候变化情景的中国玉米产量响应

第28页/共36页3.6冬小麦和玉米生产安全分析区域气候波动下的中国冬小麦产量分析

如以所选区域平均处理的拟合二项式为例，2010年、2020年、2030年、2040年、2050年气候单产（式5-7）预测值分别为250kg/ha、121kg/ha、-16kg/ha、-145kg/ha、-433kg/ha，对应时间全国的气候单产yw（式5-8）预测值分别为169kg/ha、188kg/ha、206kg/ha、223kg/ha、256kg/ha。而对应时间的研究区平均处理的实际单产（式5-9）预测值分别为4909kg/ha、5117kg/ha、5284kg/ha、5416kg/ha、5666kg/ha，表现为稳定的增长趋势。第29页/共36页3.6冬小麦和玉米生产安全分析

冬小麦生产安全分析若根据式5-9的预测实际单产，未来50年内如果中国小麦种植面积保持1995年~2004年的平均水平2.66005×107ha，即使考虑未来人口增长，则2010年、2020年、2030年、2040年、2050年人均小麦粮食占有量由2004年的70.2kg增长到94.6kg、98.6kg、101.9kg、104.4kg、109.2kg。如果未来50年总产保持1995年~2004年的平均水平，则2010年、2020年、2030年、2040年、2050年节省小麦种植耕地面积比例分别达21.75%、24.92%、27.30%、29.07%、32.20%。在未来50年模拟情景气候变化条件下，即使考虑未来50年中国人口的增长，中国小麦产量总体保持稳定增长，人均小麦可供应量仍保持稳定提高。总体看来，随着生产技术水平的稳定提高和种植面积的保证，未来中国小麦产量的气候波动响应不会对中国粮食生产安全产生负效应。第30页/共36页3.6冬小麦和玉米生产安全分析

区域气候波动下的中国玉米产量分析如以全国平均处理的拟合二项式为例，2010年、2020年、2030年、2040年、2050年气候单产预测值yw（式5-10）分别为-48.0kg/ha、-173.4kg/ha、-270.8kg/ha、-34.3kg/ha、-469.9kg/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，为中低增产区，各预测区表现为不同程度的稳定的增长趋势。

第31页/共36页3.6冬小麦和玉米生产安全分析玉米生产安全分析若根据式5-12的预测实际单产，未来50年内如果中国小麦种植面积保持1995年~2004年的平均水平2.436801×107ha，即使考虑未来人口增长，则2010年、2020年、2030年、2040年、2050年人均小麦粮食占有量由2004年的99.5kg增长到109.5kg、110.6kg、111.9kg、113.2kg、116.1kg。如果未来50年总产保持1995年~2004年的平均水平，则2010年、2020年、2030年、2040年、2050年节省玉米种植面积比例分别达21.07%、21.90%、22.82%、23.68%、25.58%。虽然未来50年模拟情景气候变化对中国玉米产量为明显的减产效应，而且即使考虑未来50年中国人口的增长，中国玉米产量总体仍然保持稳定增长，人均玉米可供应量仍保持小幅度稳定提高。总体看来，随着生产技术水平的稳定提高和种植面积的保证，未来中国玉米产量的气候波动响应在未来人口高峰期间，不会对中国粮食生产安全产生负效应第32页/共36页

主要结论自1982年以来，研究区域