气候变化对河南省小麦和玉米气候资源利用效率的影响
2019-11-22徐延红李树岩
徐延红,李树岩
(1.中国气象局 河南省农业气象保障与应用技术重点实验室,河南 郑州 450003;2. 洛阳市气象局,河南 洛阳 471000;3.河南省气象科学研究所,河南 郑州 450003)
气候资源是农业生产发展的重要环境条件,特别是光、热、水等重要的农业气候资源,直接关系到农业生产[1-2]。如何科学有效地利用农业气候资源,按照各类农作物生长所需,实现光、热、水等气候因子的最佳组合,提高农业生产效率,最终实现农业气候资源的可持续利用,是农业气候资源研究的最终目标[3]。多年来,针对作物对气候资源的利用进行了多方面的研究,作物品种涉及广泛,包含了玉米、水稻、小麦、高粱、谷子、大豆和棉花等,形成了多种农业气候资源的利用效率及程度评价的方法和模式[4-13],为开展河南地区气候资源利用效率评估提供了技术支撑。李夏子[5],王晓煜[6]等从产量潜力损失率的角度分析作物对气候资源的利用状况。徐玲玲[7],Guo J P[8]等利用气候适宜度理论对一个地区的农业气候资源配置和利用情况进行评价。崔读昌[9],高涛[10],辛吉武[11],钱锦霞[12],宋梦美[13],李淑娅[14]等利用作物干物质或产量与光热水资源的比值,定量评价气候资源利用情况,该方法能够直观反映气候资源对产量的贡献,对合理评价生态气候资源具有重要意义。然而,当前研究多数是分析单一气象要素利用率,较少分析作物对不同气候资源综合利用率的时空分布。
河南省地处黄淮海平原腹地,光热资源丰富,是中国主要的小麦-玉米轮作种植区。其小麦播种面积约占全国小麦面积的25%,总产量占全国总产的24%左右,每年向市场提供的商品小麦占全国的25%~30%,均居全国第一[15]。玉米作为河南省的第二大作物,面积和总产量约占全国的10%[16]。小麦、玉米产量关系到国家粮食安全和人民温饱问题[17]。河南省地形复杂,横跨第二、三级地貌台阶,包括山地、丘陵、平原和盆地[18],自然资源地理分布差异较大,加之作物管理水平和技术等方面的限制,使得该区域的农业发展不平衡,粮食作物气候资源利用率区域间分布不均,蕴藏着巨大的潜力。受气候变化影响,近50a河南省平均气温呈明显增加趋势,降水量年际间波动增大,时空分配不均,局地干旱、高温、低温冷害等极端气候事件加剧,农业生产的不稳定性增加[19-20],气候变化对充分合理利用气候资源、保障粮食安全提出了新的要求。因此,通过对河南省农业气候资源利用效率的综合分析,研究其时空分布的变化特征是指导粮食生产和发展地区经济的需要。本研究选用河南省29个农业气象观测站1981—2016年气象数据和冬小麦、夏玉米产量数据,分析冬小麦、夏玉米生育期内气候资源变化特征,研究粮食生产的光能、热量和降水资源利用效率的现状和变化规律,对于提高农业气候资源利用效率和粮食产量以及适应未来气候变化具有重要指导意义。
1 资料与方法
1.1 数据来源
气象数据来源于河南省29个农业气象观测站,为1981—2016年的逐日气象资料,包括平均气温、最高气温、最低气温、日照时数、降水量,辐射资料由日照时数估算得到[21]。作物产量数据来源于河南省农业统计年鉴。
1.2 研究方法
1.2.1 光能利用效率(Light Utilization Efficiency,LUE) 作物干物重的 90%以上来自光合作用,籽粒产量主要取决于光能利用率的高低,即光合产物中贮存的能量占有效辐射能或占太阳总辐射能的百分比。
(1)
式中,LUE为光能利用率(kg·MJ-1),Y为粮食作物产量(kg·hm-2), ∑Q为作物生长季内的太阳总辐射(MJ·hm-2),本研究中小麦生长季为10月—翌年5月,玉米生长季为6—9月。
1.2.2 热量利用效率 (Heat Utilization Efficiency,HUE) 农业气象学的热量条件以积温表示,≥0℃积温表示全年生长期的热量条件,≥10℃的积温表示喜温作物生长的热量条件。
(2)
式中,HUE为热量资源利用效率(kg·℃-1·d-1·hm-2),∑T为作物生长季积温(℃·d)。小麦生长季∑T为≥0℃的积温,玉米生长季∑T为≥10℃的积温。
1.2.3 降水量利用效率 (Precipitation Utilization Efficiency,PUE) 降水量利用效率反映降水对粮食产量形成的贡献,即每生产 1kg粮食所消耗的自然降水量,本研究忽略了灌溉、地下水供给,并非严格意义上的水分利用效率。
(3)
式中,PUE为降水资源利用效率(kg·mm-1·hm-2),∑R为作物生长季内总降水量(mm)。
1.2.4 气候资源利用效率综合指数(Utilization Efficiency Comprehensive Index,CUE*)
(1)各因子对气候资源利用效率的贡献率
对式(1)两边取对数,得到
(4)
根据研究区域内每个站点每年作物产量、太阳辐射以及光能利用效率,通过一元线性回归,分别得到近36 a(1981—2016年)ln(Y)、ln(1/Q)和ln(LUE)的变化趋势及相应的一次线性回归方程,其斜率分别为k1、k2、k3。辐射变化对光能利用效率变化的贡献率为[22]
(5)
(6)
式中,i为站点序号,n为研究区域内总站点数,CQLi表示某站点辐射变化对LUE变化的贡献率,CQL为研究区域内所有站点辐射变化对LUE变化的平均贡献率,CQL值越大表示辐射变化对LUE变化的贡献越大。
同样,对式(2)、(3)取对数,用上述方法统计积温、降水量变化对HUE、PUE变化的贡献,用CTH、CRP分别表示积温变化对HUE变化的平均贡献率、降水量变化对PUE变化的平均贡献率。
(2)气候资源利用效率综合指数
对光能、热量和降水量利用效率的综合评价采用综合利用效率指数。首先对LUE、HUE和PUE作离差标准化处理,使原始值映射到[0~1]之间,记为LUE*、HUE*、PUE*,具体计算公式如下式(7)
(7)
式中,Zi为某要素的标准化值,Xi为原始要素值,Xmax为某要素序列的最大值,Xmin为某要素序列的最小值。
综合光能、热量和降水量利用效率,定义粮食作物对农业气候资源的综合利用效率指数CUE*为:
CUE*=wlLUE*+whHUE*+wpPUE*
(8)
(9)
(10)
(11)
式中,wl、wh、wp分别为光能、热量和降水量利用效率的权重系数,根据光、热、水资源变化对气候资源利用效率的贡献率来确定。
2 结果与分析
2.1 气候资源变化特征
2.1.1 太阳总辐射 如图1所示,除卢氏站点外,近36 a全省太阳总辐射均呈不同程度的减少趋势,小麦生长季平均每10 a减少56 MJ·m-2,玉米生长季平均每10 a减少62 MJ·m-2,豫东、豫北地区受人类活动影响,大气能见度降低,辐射资源下降较快;豫西地区多山区,所受影响相对较少,辐射资源下降较慢。
空间分布来看,豫北地区雨日数相对较少,大气透明度高,辐射较强,小麦生长季辐射在3 175~3 215 MJ·m-2,玉米生长季大部分地区辐射在2 250 MJ·m-2以上。豫南、豫东地区阴雨日数较多,大气透明度降低,太阳辐射较少,小麦生长季辐射为3 115 ~3 155 MJ·m-2,玉米生长季辐射在2 230~2 250 MJ·m-2。豫西、豫西南受地形影响,辐射资源亦较低,小麦生长季辐射资源多在3 055~3 095 MJ·m-2,玉米生长季辐射资源多在2 200~2 230 MJ·m-2,其中内乡站点辐射资源为全省最低。
2.1.2 积温 在全球变暖的背景下,河南省作物生长季的热量条件呈不同程度的上升趋势(图2)。小麦生长季≥0℃的积温倾向率为52~151 ℃·d·10a-1,其中豫西北、驻马店和南阳地区增幅较大,增速在100 ℃·d·10a-1以上。玉米生长季,≥10℃的积温增加速率集中在7~47 ℃·d·10a-1,伊川、郑州站点增加幅度最为明显,气候倾向率达87 ℃·d·10a-1。
由于地形、地貌以及大气环流条件的不同,热量资源空间分布差异较大,表现为北部少于南部,丘陵山区少于平原地区。小麦生长季≥0℃积温在1 850~2 450℃·d,豫南、豫中地区热量资源丰富,积温在2 250℃·d以上,豫西、豫北地区积温低于2 150℃·d。玉米生长季≥10℃积温在2 600~3 040℃·d范围内变化,豫西、豫北北部热量资源相对匮乏,积温低于2 920℃·d,信阳地区高达3 040℃·d,其他大部地区集中在2 960~3 000℃·d。
2.1.3 降水量 如图3所示,小麦生长季降水量变化气候倾向率为-1 mm·10a-1,玉米生长季降水量变化气候倾向率为-3 mm·10a-1,均未通过0.05的信度检验(P>0.05)。降水量变化趋势空间分布不一致,小麦生长季,豫北、豫中以及信阳地区降水量呈减少趋势,东部的商丘、周口、南部的驻马店、南阳东部降水量增加明显,气候倾向率为6~15 mm·10a-1。玉米生长季,豫西、豫东南降水量呈减少趋势,豫北、豫东降水量有所增加,其中伊川、三门峡、汝州、巩义减少速率在20 mm·10a-1以上,为全省降水量减少最明显的区域,新乡、商丘增加速率在20 mm·10a-1以上,为全省降水增加较明显的地区。
图1 小麦、玉米生长季太阳总辐射的分布Fig.1 Distribution of total solar radiation during wheat and maize growing seasons
图2 小麦、玉米生长季积温的分布Fig.2 Distribution of accumulated temperature during wheat and maize growing season
图3 小麦、玉米生长季降水量的分布Fig.3 Distribution of precipitation during wheat and maize growing seasons
受季风气候影响,加之地理条件复杂,河南省降水量水平及垂直差异显著。小麦生长季降水总量在150~440 mm,由北向南明显递增,开封-郑州-洛阳-三门峡连线以北降水量在230 mm以下,最大降水中心位于信阳,降水量达440 mm。玉米生长季降水总量在360~600 mm,豫西北降水量低于400 mm,为我省降水量的低值区,豫南地区降水量较充沛,达520 mm以上。
2.2 气候资源利用效率的时空分布
2.2.1 光能利用效率 全省光能利用效率如图4所示,小麦生长季,在辐射和产量变化的共同作用下,除三门峡外,其他28个站点小麦LUE呈显著增加趋势(P<0.05)。焦作-许昌-驻马店连线以东地区,LUE增幅较大,气候倾向率在0.35~0.65 kg·MJ-1·10a-1。玉米生长季,除三门峡站点LUE不显著增加外,其余站点LUE均呈显著增加趋势(P<0.05),增加速率为0.12~0.83 kg·MJ-1·10a-1。
河南省小麦平均光能利用率为1.48 kg·MJ-1,而玉米达到2.27 kg·MJ-1,超过世界陆地植物平均光能利用率1.69 kg·MJ-1的水平[21]。从分布图上可以看出(图4),区域间差异也较大,表现为豫东北、豫东平原>豫中部>豫西丘陵旱地、淮南。豫东北、豫东平原光能利用率最高,小麦达1.77 kg·MJ-1以上,玉米达2.50 kg·MJ-1;豫中、豫西北、豫东南居中,小麦LUE在1.50~1.68 kg·MJ-1,玉米在2.10~2.40 kg·MJ-1;豫西、豫西南小麦LUE在1.32 kg·MJ-1以下,玉米不超过1.90 kg·MJ-1。小麦平均光能利用率大部分集中在1.32~1.68 kg·MJ-1范围内,玉米集中分布在2.00~2.40 kg·MJ-1。
2.2.2 热量利用效率 如图5所示,小麦生长季,除三门峡、巩义、郑州、林州外,其余站点HUE显著增加,气候倾向率为0~0.69 kg·℃-1·d-1·hm-2·10a-1。玉米生长季,除卢氏站点外,其余站点HUE均显著增加,气候倾向率为0.04~0.46 kg·℃-1·d-1·hm-2·10a-1。
空间上,热量利用效率并不随热量资源的多寡而呈现规律性的变化(图5)。小麦生长季,热量利用效率在1.30~2.50 kg·℃-1·d-1·hm-2,呈现北高南低、东高西低的特点。淮河以南虽热量丰富,但小麦产量不高,热量资源利用效率低,信阳HUE为1.30 kg·℃-1·d-1·hm-2。旱区由于缺水,热量资源利用效率低于灌区。豫西、豫西南热量资源利用效率在1.70 kg·℃-1·d-1·hm-2以下,豫东商丘、开封、周口,以及豫东北濮阳、安阳、新乡HUE较高,达2.3 kg·℃-1·d-1·hm-2。玉米生长季,HUE在1.32~1.78 kg·℃-1·d-1·hm-2范围,呈西南向东北递增的特征。豫东北、豫东热量利用效率较高,达1.72 kg·℃-1·d-1·hm-2以上。中部、西北部热量资源较为丰富,但常常因为水分限制,热量资源不能充分利用,HUE为1.48 ~1.64 kg·℃-1·d-1·hm-2。豫西、豫西南地区HUE为全省最低,数值在1.32~1.40 kg·℃-1·d-1·hm-2。
2.2.3 降水量利用效率 如图6所示,小麦生长季,全省PUE气候倾向率在0.5~10.1 kg·mm-1·hm-2·10a-1之间,79%的站点增加趋势显著,其中封丘、濮阳站点变化最为明显,PUE增速达7.5 kg·mm-1·hm-2·10a-1。玉米生长季,PUE气候倾向率在0.5 ~3.6 kg·mm-1·hm-2·10a-1之间,卢氏、郑州站点虽产量显著增加,但降水量变化不显著,两者作用导致PUE变化不显著,其余站点均显著增加(P<0.05)。经统计,河南省小麦、玉米平均降水量利用效率分别高达20.6、11.8 kg·mm-1·hm-2,均高于全国平均水平[23]。
从空间分布上看,降水量利用效率的高低和降水量密切相关,表现为由北部向南部逐渐减小的趋势。淮河以南降水量丰富,相对小麦需水量而言,为水分盈余区。由于阴雨多、光照不足等小麦生产限制因子,该区小麦PUE在8 kg·mm-1·hm-2以下,为降水量利用效率最低的区域。豫西南浅山丘陵地区,小麦PUE在13 kg·mm-1·hm-2以下,玉米PUE在10.5 kg·mm-1·hm-2以下。中部地区小麦PUE集中在18~28 kg·mm-1·hm-2,玉米在10.5~13.5 kg·mm-1·hm-2。豫北为半湿润易旱气候区,水分亏缺量大,但由于灌溉条件较好,气候资源协同作用,降水量利用效率较高,小麦PUE超过33 kg·mm-1·hm-2,玉米达14.5 kg·mm-1·hm-2。
2.3 气候资源利用效率综合指数及评估
2.3.1 气候资源利用效率的贡献率 计算1981—2016年各气象因子变化对气候资源利用效率变化的影响程度,结果见表1。由表中可见,小麦季研究区域内总辐射变化对光能利用效率升高的平均贡献率CQL为10%,玉米季辐射变化的贡献率CQL为15%。小麦季温度变化对热量利用效率的贡献率CTH为30%,玉米CTH为6%,温度变化对小麦气候资源利用效率变化的影响程度高于玉米,这是由于玉米季温度升高的幅度较小麦季不明显,且玉米生长的积温需求比小麦高得多,热量增加对玉米生产的影响相对小麦较小。小麦季降水量变化对降水量利用效率的贡献率CRP为20%,低于玉米生长季CRP(30%)。不同气象要素变化对相应资源利用率变化的影响程度不同,温度变化对小麦HUE增加的贡献率最大,降水次之,辐射最小;玉米季贡献率依次为降水>辐射>积温。贡献率除了受光、热、水资源自身增加或减少的速率影响外,还与农业生产对气候要素变化的敏感性相关,如玉米生长季豫西、豫北北部干旱地区对水资源的变化更为敏感,表现为该地区降水量变化的贡献率高于光、热要素;小麦生长季豫北地区温度变化的贡献率又高于光 、水资源的贡献率。
表1 光、温、水因子对气候资源利用效率变化的贡献率
图4 小麦、玉米生长季光能利用效率分布图Fig.4 Distribution of light utilization efficiency during wheat and maize growing seasons
图5 小麦、玉米生长季热量利用效率分布图Fig.5 Distribution of heat utilization efficiency during wheat and maize growing seasons
图6 小麦、玉米生长季降水量利用效率分布图Fig.6 Distribution of precipitation utilization efficiency during wheat and maize growing seasons
2.3.2 气候资源利用效率综合指数 根据各气象因子的贡献率及小麦、玉米季的气候资源利用效率,计算气候资源利用效率综合指数,如图7所示,小麦生长季的气候资源利用效率综合指数呈增加趋势,气候倾向率为0~0.15·10a-1,除三门峡、巩义、郑州、林州外,其余站点均通过0.05显著性检验(P<0.05),60%站点利用效率综合指数CUE*变化集中在0.05~0.10·10a-1。玉米生长季,研究区域内所有站点CUE*均呈显著增加趋势(P<0.05),气候倾向率为0.03~0.09·10a-1,其中开封、汝州、襄城、伊川、周口、永城站点增加趋势最为明显,气候倾向率为0.07~0.09·10a-1。
从空间分布上看,河南省小麦对光、温、水气候资源的利用效率综合指数分布在0.24~0.50,呈西南向东北递增的趋势,濮阳、安阳东部、鹤壁、新乡综合利用效率较高,为0.50。三门峡、洛阳、南阳、信阳综合利用效率较低,为0.24~0.31。中部地区CUE*集中在0.35~0.47。这一趋势与作物生长季内积温利用效率的空间分布特征一致。玉米气候资源利用效率综合指数为0.21~0.31,与降水量利用效率的相关性最高,呈北高南低、东高西低的特征,高值区分布在豫北焦作、新乡、鹤壁、安阳,低值区分布在三门峡、洛阳、南阳、驻马店。未来分别提高小麦热量利用效率和玉米降水量利用效率有助于获得较高的气候资源利用效率综合指数。总体上,小麦气候资源综合利用效率高于玉米。
图7 小麦、玉米生长季气候资源利用效率综合指数分布图Fig.7 Distribution of comprehensive utilization efficiency index during wheat and maize growing seasons
3 讨 论
1)不同地区不同作物的气候资源利用效率不同,玉米积温、降水利用效率以及综合指数低于小麦,而光能利用效率高于小麦,从区域上来看,豫东北光、热、水资源利用率普遍高于豫西南,这与光照、热量、水分等条件、不同作物本身的生物学特性以及当地的生产力水平等多因素相关。高涛等[24]研究认为内蒙古粮食作物平均热量利用效率≥0℃达0 .74 kg·℃-1·d-1·hm-2,≥10℃达0 .85 kg·℃-1·d-1·hm-2,降水利用率为6.47 kg·mm-1·hm-2。钱锦霞等[12]认为山西省玉米平均光能利用率为1.46 kg·MJ-1,热量利用效率为 1.37 kg·℃-1·d-1·hm-2,降水量利用效率为12.28 kg·mm-1·hm-2。王晓煜等[6]对气候变化背景下东北三省主要粮食作物的气候资源利用效率研究表明,玉米对气候资源利用效率高于高粱、水稻、谷子、春小麦,其光能利用率达1.5 kg·MJ-1。与前人研究相比,河南地区小麦、玉米的光、热、水资源利用效率明显为高,但区域间差异更明显。
2)本文在研究小麦、玉米光、热、水单项利用率时空变化特征的基础上,分析辐射、积温、降水量变化对LUE、HUE、PUE变化的贡献,并用指数化处理方法计算了光、热、水综合利用效率,与前人研究相比[9-14],从气候角度探讨了资源利用率变化原因,反映气候变化对资源利用率的影响程度,为河南省作物合理利用光、热、水资源,挖掘高产潜力,粮食作物的合理布局以及适应未来气候变化提供理论指导。但是本文只考虑了自然降水影响,而灌溉是河南省大部分地区农作物生产的重要水分来源,因此计算的水分利用效率与实际情况有一定的差距,有待今后进一步完善。
3)气候变暖对河南省夏玉米-冬小麦种植模式产量的影响既是机遇又是挑战。温度升高使夏玉米-冬小麦轮作系统茬口时间出现一定时间的延长[25],可以选择适合早播、生育期长的夏玉米品种,调整玉米、小麦的播期和种植结构,提高气候资源利用效率。目前,土地资源、水资源、养分资源等减少,土壤质量退化等问题制约着农业的可持续发展[26]。从农业气候资源开发利用来看,设施农业可以人为控制农业生物生长发育所需要的生态条件,复种指数高,可大幅提高资源利用效率,实现资源替代。另一方面,河南省地形复杂,山地丘陵在40%左右,这部分土地气候资源丰富多样,为农林牧综合发展、特色农业提供了优越条件,因此,在深度挖掘耕地农业气候资源潜力的同时,可注重非耕地农业气候资源的开发利用。
4 结 论
全球气候变化背景下,河南省1981—2016年太阳总辐射呈逐渐减少的变化趋势,小麦生长季平均每10 a减少56 MJ·m-2,玉米生长季平均每10 a减少62 MJ·m-2,空间分布不均。热量资源表现为全省一致的增加趋势,小麦生长季≥0℃的积温气候倾向率为52~151℃·d·10a-1,玉米生长季≥10℃的积温气候倾向率为7~87℃·d·10a-1,豫南、豫中热量资源丰富,豫西、豫北地区相对匮乏。降水量总体变化趋势不显著,空间上由北向南递增。
1981—2016年除三门峡站点外,河南省光能利用效率呈显著增加趋势(P<0.05),小麦生长季LUE在1.14~1.77 kg·MJ-1,玉米生长季LUE在1.8~2.5 kg·MJ-1,空间分布表现为豫东、豫东北>豫中、豫西北、豫东南>豫西、豫西南。小麦、玉米生长季HUE分别为1.3~2.5 kg·℃-1·d-1·hm-2和1.32~1.78 kg·℃-1·d-1·hm-2,呈北高南低、东高西低的特点。小麦季PUE为8~43 kg·mm-1·hm-2,玉米季PUE为9.5~15.1 kg·mm-1·hm-2,存在一定的纬向差异,豫东北、豫东是降水资源利用效率的高值区,三门峡、洛阳、南阳、信阳地区较低。
在光、热、水资源共同作用下,河南省气候资源利用效率综合指数呈增加趋势,增加速率分别为0~0.15 10a-1(小麦)、0.03~0.09 10a-1(玉米)。从区域上看,小麦利用效率综合指数分布在0.24~0.50,温度变化对HUE增加的贡献率最大,提高热量利用率有助于增加小麦气候资源利用效率综合指数。玉米生长季利用效率综合指数分布在0.21~0.31,降水量变化对资源利用率变化的影响大于积温和辐射,提高降水资源利用率有利于增加玉米气候资源的综合指数。豫北的光、热、水综合利用效率为全省的高值区,西部和南部地区为全省的低值区,农业气候资源利用率的潜力较大。