辽宁省玉米初霜冻特征及风险评估

郑 驰,杨利霞

(1.镇巴县气象局,陕西镇巴 723600；2.汉中市气象局,陕西汉中 724200)

辽宁省玉米初霜冻特征及风险评估

郑 驰1,杨利霞2

(1.镇巴县气象局,陕西镇巴 723600；2.汉中市气象局,陕西汉中 724200)

利用1984—2013年辽宁省33个气象台站0cm地温资料和1992—2011年辽宁省14个地市玉米种植分布数据,统计分析辽宁省初霜冻日期的特征。结果表明：辽宁平均初霜冻日均出现在每年10月,76%的站点平均初霜冻日期在10月15日之前,平均初霜冻日期自沿海向内陆逐渐提前；初霜冻日为10月5日的10a平均变化等值线明显北抬,表明初霜冻对辽宁玉米种植的影响范围在逐步缩小；从玉米初霜冻风险等级区划图来看,高风险区主要集中在松辽平原的阜新、昌图、开原三地,低风险集中在气候温暖、湿润的大连、长海、旅顺三地；在气候变暖的大气候背景条件下,可以适当调整种植制度,在辽宁省北部地区扩大玉米种植面积。

玉米；初霜冻；灾害；风险评估；辽宁省

玉米是辽宁地区的主要粮食作物,初霜冻是影响玉米产量、品质的主要因素。在已有的研究中,常将下半年日最低0cm地温≤0℃的第一天定义为初霜冻日,或将下半年百叶箱所观测的日最低气温≤2℃的第一天定义为初霜冻日[1-2]。韩荣青等[3]的研究表明利用气温定义的初霜冻日较初霜冻的实际观测日期偏晚的区域面积大,偏晚幅度高,而利用地温定义的初霜冻日普遍更接近于观测日期,并且指出20世纪60、70、80和90年代初霜冻发生地带北移明显,2000年以后初霜冻出现日期较前4个10年都明显偏迟。王萍[4]、张林梅等[5]的研究表明初霜冻的年代际变化还和初霜冻发生的早晚有关,初霜冻发生越早的地方年际变化越小,反之年际变化越大。

利用辽宁省玉米主要产区30a玉米生长季内的0cm地温资料,结合辽宁省玉米种植面积分布特征以及玉米产量数据对辽宁省的玉米初霜冻日特征进行研究,并进行风险评估,对科学种植、有效防御玉米初霜灾害、提高玉米产量、减少经济损失具有重要意义。

1 数据与研究方法

0cm地温气象资料采用1984—2013年辽宁省33个气象站观测数据,资料来源于中国气象科学数据 共享服务网 (http：// cdc.cma.gov.cn/home.do)；玉米播种面积、玉米单产数据和辽宁省行政面积采用1992—2011年 《辽宁统计年鉴》资料。

初霜冻只有发生在玉米成熟之前才会对玉米产量、品质造成影响,因此重点分析9—10月的0cm地温数据。将9—10月第1次日最低0cm地温≤0℃的当日定义为初霜冻日,根据QX/T 88—2008《作物霜冻害等级》并结合当地实际情况,将霜冻分为：轻霜冻、中霜冻、重霜冻、严重霜冻四个等级(表1)。

表1 霜冻等级划分

灾害风险一般定义为一定概率下灾害造成的破坏或损失。评估程序主要包括致灾因子危险性、孕灾环境暴露性、承灾体易损性和人类防灾减灾能力四个方面的评估[6-8]。致灾因子危险性是某一地区某一时间段内发生霜冻的严重程度以及发生霜冻频率的综合体现,因此将两者结合起来可以组成一种较客观地反应危险性的指标；孕灾环境暴露性是指研究风险源和承灾体的接触暴露关系；承灾体易损性是指同一级别的灾害对不同承灾体造成的损失,这里将玉米的相对单产作为承灾体易损性的评价指标；人类防灾减灾能力要考虑当地科学技术水平、不可防御风险等社会综合因素。玉米霜冻灾害风险指数基于致灾因子危险性、孕灾环境敏感性、承灾体易损性、人类防灾减灾能力的综合影响构造,因此可以作为评估玉米霜冻的灾害风险指数的综合指标。王晾晾等[6]对吉林等地玉米霜冻风险指数的研究中构造了模型

其中,R为霜冻灾害风险指数,Hb为致灾因子危险性,Vc为孕灾环境暴露性,Vd为承灾体易损性, D为人类防灾减灾能力,Wb为致灾因子对灾害风险的权重系数；Wc为物理暴露性对灾害风险指数的权重系数；Wd为灾害敏感性对灾害风险的权重系数；We为人类防灾减灾能力对灾害风险指数的权重系数。

参照文献[8],Hb、Vc、Vd计算公式：

式中,j=1,2,3,4,依次代表轻、中、重、严重霜冻；n为总年数；Dj(℃)为j等级霜冻强度下的温度中值；Gj(℃)为j等级霜冻强度出现的次数；Sp(hm2)为各地市20a玉米平均种植面积；S(km2)为各地市行政面积。G(kg)为各地市20a玉米平均单产；¯G(kg)为整个研究区域20a玉米单产平均值。

式中各因子的权重系数采用专家打分法获得,参考郭志华等[6]的研究结果确定各致灾因子对灾害风险指数的权重系数 (见表2)。由于防灾减灾能力占权重系数小,且评估过程复杂不易操作,故在评估过程中忽略此项。

表2 各致灾因子权重系数

2 结果与分析

2.1 初霜冻日特征分析

2.1.1 平均初霜冻日特征 统计1984—2013年辽宁省33个气象站平均初霜冻日(见表3),发现,辽宁平均初霜冻日均出现在每年10月,其中,平均初霜冻日最早为10月1日(新宾),最晚为10月26日(长海)。辽宁内陆北部和西部地区平均初霜冻日大多在5日前后,最南部的旅顺、大连、长海平均初霜冻日则在10月下旬,总体来说,平均初霜冻日自沿海向内陆逐渐提前。

2.1.2 初霜冻日期10a平均变化特征 辽宁省各地以种植晚熟玉米为主,10月中旬前为玉米主要收获期。统计结果显示辽宁33个站点中有76%的站点初霜冻日在10月15日之前,这段时间,尤其是10月上旬发生霜冻会减少玉米产量、降低玉米品质,造成一定经济损失。以10月5日为标准,通过分析1984—1993年、1994—2003年、2004—2013年内平均初霜冻日为10月5日的等日期线空间分布(图1),了解近30a来辽宁省初霜冻的变化趋势和影响范围。从图1可以看出,近30年来平均初霜冻日为10月5日的等日期线明显北移,2004—2013年平均初霜冻日为10月5日之前的影响范围比前20a大幅减少,仅有叶柏寿和昌图的偏北地区。

表3 1984—2013年辽宁省33个气象站平均初霜冻日

图1 辽宁省1984—2013年初霜冻日期为10月5日的10a平均变化

2.2 初霜冻风险评估

2.2.1 致灾因子危险性评估 利用公式(2)计算各站点的Hb,对Hb采用相对化处理方法进行无量纲化,使用自然断点法对无量纲化后的Hb进行分级。表4为致灾因子危险性分区阈值范围,据此制作辽宁省玉米霜冻致灾因子分区(图2,见第25页)。可以看出,在辽宁省的中南部地区由于初霜冻发生频次较低、初霜冻日期较晚、强度也较弱,因此这些地区玉米初霜冻致灾因子危险性较低；辽宁省北部、东北部以及西北部地区由于初霜冻发生频次高、初霜冻发生日期早、初霜冻强度大,造成这些地区玉米致灾因子危险性高；辽宁省玉米初霜冻致灾因子危险性整体由南向北随纬度升高而增加。

表4 致灾因子危险性分区闭值范围

2.2.2 孕灾环境暴露度评估 由于《辽宁省统计年鉴》中对于行政面积和玉米种植面积数据以市为单位进行统计,在利用公式 (3)计算辽宁各地市玉米孕灾环境暴露度值 (表5)进行玉米初霜冻风险评估时将同一地市的站点赋予相同的孕灾环境暴露度值。考虑计算方便并兼顾利于实际操作的原则,将孕灾环境暴露度分为5级 (表6)。

表5 孕灾环境暴露度

表6 孕灾环境暴露度分区闭值范围 hm2/km2

从孕灾环境暴露度来看辽宁省玉米种植并不是均一分布的,地域差异较大。辽宁省东南部的玉米种植密度较小,孕灾环境暴露度较小；辽宁省的东北部地区玉米种植密度大,孕灾环境暴露度较大,其中沈阳、铁岭为高暴露区,孕灾环境暴露度值在20hm2/km2以上；阜新、锦州为次高暴露区。

2.2.3 承灾体易损性评估 同一级别的霜冻强度对单位面积玉米造成的经济损失的差异取决于玉米的单产,单产越大经济损失越大,反之越小。因此,把玉米的相对单产作为承灾体易损性评估指标。采用自然断点法划分承灾体易损性分区阈值范围(表7)。

表7 承灾体易损性分区闭值范围

利用公式(4)计算各地市的承灾体易损性,再根据分区阈值进行分级。结果显示,达到次高及以上易损等级的主要分布在辽宁东北和东部地区；中易损区集中在辽宁省中部地区；辽宁南部和西北部地区易损性都在0.9以下,属于低易损区。由此也说明辽宁省各地玉米产量差异性较大,其中东北部的沈阳、铁岭、抚顺等地区相对单产高,是承灾体易损较大的地区。

2.2.4 初霜冻灾害风险评估与区划 为了消除量纲以便参与综合评估,采用相对化处理方法对致灾因子危险性、孕灾环境暴露性、承灾体易损性分别利用公式 (5)进行去量纲化处理

式中,Kmax为各因子的极大值,K为因子值,K*为去量纲后的值。

图2 辽宁省初霜冻致灾因子危险性分区

图3 辽宁省玉米初霜冻风险等级区划

将去量纲后的各因子带入公式(1)中计算33个县(市)的玉米初霜风险指数R。利用SPSS软件对风险指数R进行样本聚类分析,将33个县(市)分为低风险区、次低风险区、中风险区、次高风险区、高风险区共5类风险区。再使用surfer软件基于33个县(市)的风险指数R采用克里金插值法绘制辽宁省各地区玉米初霜冻风险区划图(图3)。从图3可以看出,①低风险区位于大连沿海地区,主要由于沿海地区温度相对较高,致灾因子危险度低；且沿海地区不以种植玉米为主要经济来源,种植面积小、相对单产低,孕灾环境暴露度小、承灾体易损性也小。②次低风险区主要位于辽西—辽南,虽然辽西—辽南地区处于致灾因子危险度低值地区、承灾体易损性也处于低值区,但是孕灾环境暴露度处于中暴露度区。③中风险区主要位于40°N～41°N区域附近,这些地区的致灾因子危险性、孕灾环境暴露性、承灾体易损性都处于中等水平。④次高风险区主要位于41°N～42°N区域附近,其中叶柏寿、建昌主要是由于致灾因子危险性高造成的。海城地区的致灾因子不属于高值区,但结合玉米种植分布数据来看其孕灾环境暴露性和承灾体易损性都处于高值区,因此海城地区也处于玉米初霜冻次高风险区。⑤高风险区主要位于辽宁省东北部,一方面是由于这些地方温度条件不好导致致灾因子大,另一方面这些地区玉米种植密度较大、相对产量较高导致孕灾环境暴露性和承灾体易损性也大,因此综合表现为玉米初霜冻为高风险区。

3 结论

(1)整体来看,1984—2013年辽宁省平均初霜冻日均出现在每年10月,76%的站点平均初霜冻日期在10月15日之前,且平均初霜冻日由南向北逐步提前。初霜冻日为10月5日的10 a平均变化等值线明显向北抬,表明初霜冻对辽宁玉米种植的影响范围在逐步缩小。

(2)从玉米初霜冻风险评估划分的五个等级中,阜新、昌图、开原属高风险区；辽西—辽南地区的长海、大连、旅顺属低风险区,适宜玉米的种植；其它地区风险居中,通过培育抗寒玉米品种可以大幅增加玉米种植面积,获得良好经济效益。

(3)在气候变暖的大气候背景条件下,可以适当调整种植制度,在辽宁省北部地区扩大玉米种植面积。

[1] 李星敏,朱琳,贺文丽,等.基于GIS的陕西省农业气候资源与区划 [M].陕西：陕西科学技术出版社,2013：75-76.

[2] YeDianxiu,ZhangYong.Characteristicsoffrost changesfrom1961to2007overChina[J].Journal ofAppliedMeteorologicalScience,2008,19 (6)：661-665.

[3] 韩荣青,李维京,艾婉秀,等.中国北方初霜冻日期变化及其对农业的影响 [J].地理学报, 2010,65(5)：526-532.

[4] 王萍,李廷全,闫平.黑龙江省近35年气候变化对粳稻发育期及产量的影响 [J].中国农业气象,2008,29(3)：268-271.

[5] 张林梅,李建丽,张建.阿勒泰地区霜冻变化特征分析[J].陕西气象,2014(6)：13-18.

[6] 王晾晾,连萍.东北三省玉米霜冻灾害风险评估及区划[J].自然灾害学报,2008,12(5)：56 -60.

[7] 郭志华,刘祥梅,肖文发,等.基于GIS的中国气候分区及综合评价[J].资源科学,2007,29 (6)：2-9.

[8] 唐为安,田红,杨元建,等.基于GIS的低温冷冻灾害风险区划研究[J].地理科学,2012,32 (3)：356-361.

P423.4

A

郑驰,杨利霞.辽宁省玉米初霜冻特征及风险评估[J].陕西气象,2015(2)：21-26.

1006-4354(2015)02-0021-06

2014-09-13

郑驰(1990—),男,汉族,陕西汉中人,学士,技术员,从事基层气象服务。