Loading [MathJax]/jax/output/HTML-CSS/jax.js
贵州霜冻气候变化特征
张波1,2, 于飞1, 吴战平1,*, 胡家敏1
1.贵州省山地环境气候研究所,贵州 贵阳 550002
2.贵州省山地气候与资源重点实验室,贵州 贵阳 550002
*通信作者,吴战平, E-mail: wuzp1019@sohu.com

作者简介:张波(1985—),男,山东临沂人,硕士,工程师,从事农业气象研究。E-mail: nj0622@126.com

摘要

为合理利用气候资源指导农业生产,利用贵州省1961—2018年84个气象站点地面0 cm日最低气温资料,采用线性趋势、M-K突变检验和空间分析等方法,分析了贵州地区初霜冻日、终霜冻日和无霜冻期的变化特征。结果表明:近58 a贵州地区平均初霜冻日、终霜冻日分别为12月2日和2月23日,平均无霜冻期为291 d,分别以2.34 d·(10a)-1、2.95 d·(10a)-1和5.76 d·(10a)-1的速率呈明显推迟、提前和延长的变化趋势;空间上初霜冻日呈现西部早、东部晚,北部早、南部晚的变化特征,69%的站点呈显著推迟的变化趋势;终霜冻日呈现南部早、西北部和东北部晚的变化特征;无霜冻期整体呈延长趋势,82%的站点通过显著性水平检验;海拔高度是影响贵州地区初霜冻日、终霜冻日和无霜冻期的关键因子,无霜冻期随纬度和海拔高度升高呈缩短趋势。初霜冻日推迟、终霜冻日提前和无霜冻期的延长,为贵州农业生产提供了充足的热量资源,利于粮食生产。

关键词: 初霜冻日; 终霜冻日; 无霜冻期; 贵州; 气候特征
中图分类号:S162 文献标志码:A 文章编号:1004-1524(2020)04-0685-11
Climate characteristics of frost in Guizhou
ZHANG Bo1,2, YU Fei1, WU Zhanping1,*, HU Jiamin1
1. Guizhou Institute of Mountainous Environment and Climate, Guiyang 550002, China
2. Guizhou Key Laboratory of Mountainous Climate and Resources, Guiyang 550002, China
Abstract

The distribution characteristics of frost in Guizhou was analyzed, which would provide scientific basis for rational utilization of climate resources and supervise agriculture.The first frost day, the last frost day and the frost free period in Guizhou were analyzed based on the daily minimum temperature of 0 cm of 84 meteorological stations in Guizhou Province from 1961 to 2018 by using linear trend, M-K mutation test and spatial analysis. The results showed that in recent 58 years, the average first frost day and the last frost day in Guizhou was December 2 and February 23, the average frost free period was 291 d, with the rates at 2.34 d·(10a)-1, 2.95 d·(10a)-1 and 5.76 d·(10a)-1, showed significant trends of delay, advance and extension ( P<0.01), respectively. In the spatial distribution, the first frost day showed the characteristics of earlier in the west, later in the east, earlier in the north and later in the south, 69% of the sites showed a trend of significant delay. The last frost day showed the characteristics of earlier in the south, later in the northwest and northeast. The frost free period was generally prolonged, and 82% of the stations passed the significance level test. Altitude was the key factor affecting the first frost, the last frost and the frost free period in Guizhou, the frost free period decreased with the increase of latitude and altitude. The first frost event occured later, the last frost event occured eailer, and the frost-free period becomed longer, which would provide sufficient heat resources and more favorable agricultural production in Guizhou.

Keyword: the first frost day; the last frost day; frost-free preiod; Guizhou; climatic characteristics

霜冻是农业生产中最常见的气象灾害之一, 每年给农业生产造成较大损失, 严重影响着农业的稳定发展。在农业生产中通常采用无霜冻期的长短来表征作物受霜冻灾害影响的概率, 无霜冻期短, 农作物生长周期短, 农作物遭受霜冻害的概率就越高; 而无霜冻期的长短是由当年的初霜冻日和终霜冻日决定。因此, 了解和掌握初霜冻日、终霜冻日和无霜冻期的变化趋势对减少低温灾害对农业生产的影响具有重要的科学意义[1, 2, 3, 4]

国内外关于初霜冻日、终霜冻日和无霜冻期已开展一系列相关研究。Erlat等[5]以日最低气温≤ 0 ℃为霜冻指标, 分析了1950— 2010年土耳其地区霜冻日变化情况, 结果表明, 土耳其的霜冻日均呈显著减少趋势。Terando等[6]针对北美地区霜冻日的研究结果表明, 1981— 2010年北美中部部分平原区霜冻日数呈增加趋势, 而其余地区则表现为显著减少趋势; 王国复等[7]、许艳等[8]、叶殿秀等[9]和宁晓菊等[10]的研究结果表明, 中国大部分地区初霜日推迟、终霜日提前, 无霜期延长, 主要与日最低气温和地面最低地温的升高有关; 李芬等[11]利用山西逐日最低地温资料, 分析了山西终霜冻日的时空分布特征; 任景全等[12]研究表明, 吉林省初霜日显著推迟, 终霜日显著提前, 且无霜期显著延长; 刘静等[13]分析了秦岭南北无霜期时空变化特征, 结果表明, 秦岭南北无霜期随纬度增加而减少, 无霜期呈显著延长趋势且发生增多突变; 周晓宇等[14]分析了东北地区初终霜日和无霜期的气候变化特征, 表明东北地区自南向北随纬度的升高和海拔的升高, 初霜日提前、终霜日推后、无霜期缩短; 张磊等[15]利用逐日最低气温资料, 研究了宁夏全区与不同区域的无霜期基本特征和变化趋势, 发现宁夏平均无霜期以4.7 d· (10a)-1的速率延长, 不同区域的无霜期均呈延长趋势; 马尚谦等[16]利用逐日地面0 cm最低气温资料, 分析了淮河流域异常初、终霜日时空变化, 研究发现, 气候态的转型促进了初霜日的推迟和终霜日的提前。针对贵州霜冻研究也有相关报道, 郜红娟等[17]利用1960— 2013年逐日最低气温格点数据, 分析了贵州霜冻日变化特征, 但未对初霜冻日和终霜冻日的气候特征做定量分析。

本文在前人研究的基础上, 以地面0 cm日最低气温≤ 0 ℃作为霜冻日指标, 系统分析贵州省初霜冻日、终霜冻日和无霜冻期的气候特征, 以期为合理利用气候资源和农业防灾减灾提供科学依据。也可充分发挥气象服务在脱贫攻坚中“ 趋利避害、减灾增效” 的作用, 促进贵州省对低温霜冻敏感的特色产业的健康可持续发展[18]

1 材料与方法
1.1 数据来源

选取1961— 2018年贵州省84个地面观测站的地面0 cm日最低气温数据, 气象数据来自贵州省气象信息中心。贵阳站点7个, 遵义站点14个, 安顺站点6个, 毕节站点8个, 铜仁站点10个, 六盘水站点3个, 黔东南站点16个, 黔南站点12个, 黔西南站点8个。

1.2 方法

以地面0 cm日最低气温≤ 0 ℃作为霜冻日指标, 将本年春(秋)季地面0 cm日最低气温≤ 0 ℃的最后1天(第1天)分别定义为终(初)霜冻日, 终、初霜冻日之间的间隔日数为无霜冻期。为了方便统计和描述, 将初霜冻日、终霜冻日均采用自1月1日起的日序列来表示, 即1月1日的日序为1, 1月2日的日序为2, 2月1日的日序为32, 以此类推。同时将霜冻日按照分级标准, 分成不同等级:较多年平均初霜冻日提前1~5 d为正常初霜冻日, 较多年平均初霜冻日提前6~10 d为偏早初霜冻日, 较多年平均初霜冻日提前≥ 11 d为特早初霜冻日; 较多年平均终霜冻日推迟1~5 d为正常终霜冻日, 较多年平均终霜冻日推迟6~10 d为偏晚终霜冻日, 较多年平均终霜冻日推迟≥ 11 d为特晚终霜冻日。

霜冻时间变化特征采用气候倾向率[19]、Mann-Kendell[20]方法分析; 利用ArcGIS技术分析初霜冻日、终霜冻日和无霜冻期的空间分布特征。

气候倾向率反映气候因子的变化趋势, 由一元线性趋势方程表示, 如公式1所示,

xt=att+a0 (t=1, 2, 3, …, n)。(1)

式中:xt为气象要素的拟合值, ata0为回归系数, t为时间序列, 10at称作气候倾向率, 表示气象要素每10 a的变化趋势。

Mann-Kendall是一种非参数趋势检验法, 监测气候要素长期变化趋势, 是目前比较常用的趋势诊断方法。在Mann-Kendall检验中, 原假设H0为时间序列数据(x1~xn), 是n个独立的、随机变量同分布的样本; 假设H1是双边检验, 对于所有的k, jn, 且kj, xk~xj的分布是不相同的, 检验的统计变量S计算如下式:

S=k=1n-1j=k+1nSgn(xj-xk)。(2)

其中, S为正态分布; xj为时间序列的第j个数据值; n为数据样本的长度; Sgn是符号函数, 标准化统计量Z由下式计算得出:

{Z=(S-1)/n(n-1)(2n+5)/18S>0Z=0S=0Z=(S-1)/n(n-1)(2n+5)/18S<0

Z为正值表示增加趋势, 负值表示减少趋势。当检验序列是否发生突变时, 常通过下式进行计算检验统计量Sk:

Sk=i=1kji-1aij(k=2, 3,4,,n)(3)

其中, aij= {1Xi>Xj0XiXj(1ji),

定义:

UFk= [Sk-E(Sk)]Var(Sk)(k=1, 2, 3, …, n)。 (4)

式中:E (Sk)=k (k+1)/4, Var (Sk)=k (k-1)(2k+5)/72, UFk为标准正态分布。若UFk值大于0, 则表明序列呈上升趋势; 小于0则表明呈下降趋势; 当它们超过显著性水平对应的临界值时, 表明上升或下降趋势显著。

2 结果与分析
2.1 时间变化特征

2.1.1 初霜冻日、终霜冻日和无霜冻期时间变化特征

图1所示为贵州省1961— 2018年初霜冻日、终霜冻日和无霜冻期时间变化特征, 由图1-a可以看出, 近58 a贵州地区多年平均初霜冻日为12月2日, 最早初霜冻日为11月12日, 发生在1992年, 最晚初霜冻日为12月16日, 发生在1991年, 相差34 d; 初霜冻日随时间变化呈显著推迟的变化趋势, 推迟速率为2.34 d· (10a)-1。从年代际变化趋势来看, 20世纪60、70年代初霜冻日分别为11月30日和11月28日, 80和90年代初霜冻日为12月1日, 均早于多年平均值; 进入21世纪初霜冻日推迟明显, 21世纪00年代初霜冻日为12月6日, 21世纪10年代初霜冻日为12月12日。

图1 初霜冻日、终霜冻日和无霜冻期时间变化特征Fig.1 Variation characteristics of the first frost, the last frost and frost-free period

从终霜冻日时间变化特征(图1-b)可以看出, 贵州地区多年平均终霜冻日为2月23日, 最早终霜冻日为1月27日, 发生在1973年, 最晚终霜冻日为3月17日, 发生在1969年, 相差49 d。近58 a终霜冻日随时间变化呈显著提前趋势, 变化率为2.95 d· (10a)-1, 终霜冻日平均提前了17.11 d。从年代际来看, 终霜冻日呈现延迟-提前-延迟-提前的波动变化趋势, 具体为20世纪60年代终霜冻日为3月4日, 晚于多年平均值; 70、80年代终霜冻日为2月24日, 与平均终霜冻日基本持平; 90年代终霜冻日呈提前态势, 终霜冻日为2月18日; 进入21世纪的00年代终霜冻日又回到多年平均终霜冻日水平, 到21世纪10年代终霜冻日为2月15日, 又呈提前的变化趋势。

无霜冻期时间变化趋势(图1-c)表明, 多年平均无霜冻期为291 d, 最小无霜冻期为253 d, 发生在1969年, 最大无霜冻期为325 d, 发生在2015年, 相差73 d。无霜冻期随时间变化呈显著延长的变化趋势, 速率为5.76 d· (10a)-1, 平均延长了33.41 d。从年代际来看, 20世纪60、70和80年代无霜冻期天数少于多年平均值, 分别为278、285、286 d; 90年代到21世纪10年代无霜冻期天数均多于多年平均值。

2.1.2 初霜冻日、终霜冻日和无霜冻期突变特征

图2所示为贵州省1961— 2018年初霜冻日、终霜冻日和无霜冻期的突变特征。从初霜冻日M-K突变检验结果(图2-a)可以看出, 贵州地区初霜冻日在1994年以前呈上下波动, 1994年以后呈递增的变化趋势, 并在2005年以后UF统计量超过显著性水平0.05的临界线, 表明初霜冻呈显著推迟趋势。根据UF和UB曲线在2条临界线之间的交点, 确定2004年为初霜冻日的突变年。

图2 初霜冻日(a)、终霜冻日(b)和无霜冻期(c)突变检验特征Fig.2 Mann-Kendall test of the first frost day(a), the last frost day(b) and frost-free period (c)

终霜冻日突变检验结果(图2-b)显示, 终霜冻日在1974年以前呈上下波动趋势, 1974年以后呈下降趋势, 并在1981年超过显著性水平0.05的临界线, 1981至1990年呈小幅度的上下波动, 1990年后UF统计量呈下降趋势。根据UF和UB曲线在2条临界线之间的交点, 确定1980年为终霜冻日的突变年。

无霜冻期突变检验结果(图2-c)显示, 无霜冻期在1979年之前UF统计量呈上下波动, 之后呈递增趋势, 并在1998年后UF统计量超过显著性水平0.05的临界线, 表明无霜冻期呈显著延长趋势。根据UF和UB曲线在2条临界线之间的交点, 确定1995年为无霜冻期的突变年。

2.2 空间变化特征

2.2.1 初霜冻日、终霜冻日和无霜冻期空间分布

由图3可以看出, 初霜冻日出现在11月上旬至12月下旬, 最早初霜冻日为11月7日, 发生在威宁, 最晚初霜冻日为12月20日, 发生在罗甸, 两地初霜冻日天数相差44 d。在空间分布上呈现东部晚、西部早, 南部晚、北部早的变化特征, 初霜冻日较早的区域主要分布在毕节市大部、六盘水市北部和遵义市局地等地区, 初霜冻日在11月7日— 11月15日; 初霜冻日较晚的地区主要分布在黔西南州望谟、册亨、南北盘江流域、黔南州罗甸, 以及遵义赤水等低热河谷区, 初霜冻日在12月13日— 12月20日; 其余大部分地区初霜冻日多发生在11月下旬— 12月中旬。终霜冻日空间上与初霜冻日呈相反的分布趋势, 西部晚、东部早, 北部晚、南部早, 终霜冻日期出现在1月下旬— 4月上旬, 平均终霜冻日为2月23日。最早终霜冻日为1月26日, 发生在望谟, 最晚终霜冻日为4月7日, 发生在威宁。

图3 初霜冻日、终霜冻日和无霜冻期空间分布Fig.3 The spatial distribution of the first frost day, the last frost day and frost-free period

无霜冻期空间分布特征如图3所示, 无霜冻期空间分布特征和初霜冻日空间分布基本一致, 自东向西、自南向北呈递减趋势, 无霜冻期较短的地区主要分布在毕节市大部、六盘水市北部和遵义市局地等地区, 无霜冻期在215~242 d; 无霜冻期较长的地区主要分布在黔西南州望谟、册亨、南北盘江流域、黔南州罗甸, 以及遵义赤水等低热河谷区, 长度在325~351 d。最短无霜冻期发生在威宁, 为215 d, 最长无霜冻期发生在赤水, 达351 d。

2.2.2 初霜冻日、终霜冻日和无霜冻期气候倾向率空间分布

从图4可以看出, 近58 a贵州省初霜冻日气候倾向率变化范围为-0.64~6.54 d· (10a)-1, 平均气候倾向率为2.56 d· (10a)-1, 全省除了仁怀和从江等地的初霜冻日气候倾向率呈递减趋势, 其余97%站点气候倾向率呈递增趋势。贵州地区初霜冻日整体呈推迟的变化趋势, 其中有20%的站点通过0.05显著性检验, 49%的站点通过了0.01显著性检验; 空间上, 赤水、黔西、修文、开阳、镇宁、织金、施秉和平坝等地初霜冻日气候倾向率较高, 变化范围为4.01~6.54 d· (10a)-1

图4 初霜冻日、终霜冻日和无霜冻期气候倾向率空间分布特征
▲, P< 0.01; +, P< 0.05.
Fig.4 Spatial distribution of trends of the first frost day, last frost day and frost-free period

终霜冻日气候倾向率变化范围在-7.81~0.92 d· (10a)-1, 平均气候倾向率为-3.29 d· (10a)-1, 全省除了仁怀和岑巩等地的终霜冻日气候倾向率呈递增趋势, 其余站点气候倾向率呈递减趋势。贵州地区终霜冻日整体呈提前的变化趋势, 其中有32%的站点通过0.05显著性检验, 37%的站点通过了0.01显著性检验。空间上终霜冻日气候倾向率较高的区域主要分布在三都、思南、兴仁、晴隆、玉屏、清镇、平坝和荔波等地, 变化范围为-7.82~-5.01 d· (10a)-1(图4)。

无霜冻期气候倾向率变化范围在-5.18~13.82 d· (10a)-1, 平均气候倾向率为5.73 d· (10a)-1, 全省除了仁怀、岑巩、白云和江口等地的无霜冻期气候倾向率呈递减趋势, 其余95%站点气候倾向率呈递增趋势。无霜冻期整体呈延长的变化趋势, 82%的站点通过显著性检验, 其中13%的站点通过0.05显著性检验, 69%的站点通过了0.01显著性检验。空间上无霜冻期气候倾向率较高的区域主要分布在乌当、清镇、兴仁、榕江、三都、思南、黔西、普安和修文等地, 变化范围在9.02~13.82 d· (10a)-1(图4)。

2.2.3 平均初霜冻日、终霜冻日和无霜冻期临界线空间分布特征

图5所示为不同准气候周期平均初霜冻日、终霜冻日和无霜冻期临界线空间分布特征。从平均初霜冻日(12月2日)临界线空间位置(图5-a)可以看出, 不同准气候周期贵州地区平均初霜冻日临界线空间分布差异显著, 随着准气候周期的变化, 平均初霜冻日临界线整体上呈向北偏移的趋势。1961— 1990年平均初霜冻日临界线主要分布在贵州南部和东南部的安龙、贞丰、长顺、平塘、三都、剑河、榕江和从江一线; 1971— 2000年平均初霜冻日临界线较1961— 1990年变化不大, 走势基本一致, 仅黔东南部分区域略有北偏的趋势; 1981— 2010年, 平均初霜冻日临界线向北偏移的趋势明显, 主要分布在兴仁、关岭、镇宁、长顺、惠水、麻江、黄平、余庆、凤冈和正安一线; 1991— 2018年间平均初霜冻日临界线整体向北偏移, 主要分布在盘州、普定、安顺、清镇、花溪、贵阳、瓮安和金沙一线。

图5 不同准气候周期平均初霜冻日(a)、平均终霜冻日(b)和平均无霜冻期(c)的临界线空间变化Fig.5 Spatial distribution of trends of the first frost day (a), the last frost day (b) and frost-free period (c)

平均终霜冻日(2月23日)临界线空间位置变化特征如图5-b所示, 随着准气候周期的变化, 平均终霜冻日临界线整体呈向北、向西偏移的趋势。1961— 1990年平均终霜冻日临界线主要分布在安龙、关岭、平塘、榕江和从江一线; 1971— 2000年平均终霜冻日临界线主要分布在晴隆、镇宁、惠水、三都、剑河和锦屏一线; 1981— 2010年与1991— 2010年平均终霜冻日临界线走势基本一致, 主要分布在晴隆、普定、安顺、龙里、黄平、余庆和德江一线, 其中1991— 2010年间平均终霜冻日在铜仁地区临界线呈向西偏移趋势。

平均无霜冻期(291 d)临界线空间位置变化特征如图5-c所示, 空间分布上和平均初霜冻日临界线变化趋势基本一致。1961— 1990年平均无霜冻期主要分布在兴义、关岭、平塘、三都、剑河和锦屏一线; 1971— 2000和1981— 2010年平均无霜冻期临界线分布基本一致, 主要分布在兴仁、普定、惠水、贵阳、都匀、雷山、余庆、正安一线, 与1961— 1990年相比, 平均无霜冻期基本覆盖黔西南州、黔南州、黔东南州和铜仁市; 1991— 2018年平均无霜冻期临界线向北和向西偏移趋势明显, 主要分布在盘州、修文、麻江、福泉、息烽和金沙一线。

2.2.4 不同等级霜冻的概率分布

从图6可以看出:贵州地区偏早初霜冻发生概率在1%~19%, 平均概率为7%; 空间上整体呈现北部高, 南部低的分布特征; 偏早初霜冻发生概率较高的区域主要分布在东北部的印江、绥阳、岑巩、凯里、台江和西北部的威宁等地, 发生概率在12%~19%; 发生概率较低的区域分布在西南部望谟、册亨和兴义, 以及南部边缘的罗甸、平塘、荔波和从江等地, 发生概率在4%左右。

图6 偏早初霜冻和特早初霜冻概率分布特征Fig.6 Spatial distribution of probability of the earlier first frost and the earliest first frost

如图6所示, 特早初霜冻的发生概率在2%~29%, 平均概率为17%, 空间上呈现西北部高, 西南部边缘和南部边缘低的特征; 特早初霜冻发生概率较高的区域主要分布在毕节市的大部、六盘水市北部等区域, 以及中东部的瓮安、福泉和麻江等地, 发生概率在23%~29%; 发生概率较低的区域分布在西南部望谟、册亨、南部边缘的罗甸, 以及北部的赤水等地, 概率在9%左右。

从图7可以看出, 贵州地区偏晚终霜冻发生概率在2%~22%, 平均概率为8%。空间上, 偏晚终霜冻发生概率较高的区域主要分布在北部的播州、德江、松桃和东南部雷山、贵定与黎平等地, 发生概率在15%~22%; 发生概率较低的区域分布在西南部、南部边缘、东南部从江, 以及西北部威宁等地, 发生概率在6%左右。

图7 偏晚终霜冻和特晚终霜冻概率分布特征Fig.7 Spatial distribution of probability of the later last frost and the latest last frost

如图7所示, 特晚终霜冻的发生概率在7%~38%, 平均概率为25%。空间上, 特晚终霜冻发生概率较高的区域主要分布在西部大部、黔东南州北部、铜仁市西部和遵义市东部等地, 发生概率在27%~38%; 发生概率较低的区域分布在西南部望谟、册亨和南部边缘的罗甸, 以及北部的赤水等地, 概率在13%左右。

2.3 初霜冻日、终霜冻日和无霜冻期影响因子分析

利用初霜冻日、终霜冻日和无霜冻期分别与纬度和海拔高度, 以及气温建立相关性, 结果如表1所示。纬度和海拔高度在不同程度上影响着贵州地区初霜冻日、终霜冻日和无霜冻期, 纬度和海拔高度与初霜冻日、无霜冻期呈负相关, 与终霜冻日呈正相关, 具体表现为, 初霜冻日、终霜冻日和无霜冻期与海拔高度的相关系数分别为-0.74、0.63和-0.64, 相关系数均通过0.01的信度检验; 与纬度的相关系数分别为-0.11、0.31和-0.26, 表明海拔高度对贵州地区初霜冻日、终霜冻日和无霜冻期的影响更显著。

表1 初霜冻日、终霜冻日和无霜冻期的影响因子 Table 1 Influence factors of the first frost day, the last frost day and frost-free period

与气温的相关性分析表明, 初霜冻日和无霜冻期与气温呈显著正相关, 终霜冻日与气温呈显著负相关, 相关系数均通过了0.01的信度检验。秋季平均最低气温是影响初霜冻日的关键因子, 春季平均最低气温是影响终霜冻日和无霜冻期的关键因子。

3 结论与讨论

本研究基于贵州省1961— 2018年84个气象站点地面0 cm日最低气温资料, 以地面0 cm日最低气温≤ 0 ℃作为霜冻指标, 研究贵州省初霜冻日、终霜冻日和无霜冻期的变化特征, 得出如下结论:(1) 1961— 2018年, 贵州地区平均初霜冻日为12月2日, 平均终霜冻日为2月23日, 平均无霜冻期为291 d, 随时间变化初霜冻日、终霜冻日和无霜冻期分别呈推迟、提前和延长的变化趋势, 并且存在明显突变年。(2) 在空间分布上, 初霜冻日和无霜冻期空间分布基本一致, 自东向西、自南向北初霜冻日、无霜冻期呈递减的变化趋势, 初霜冻日较早和无霜冻期较短的地区主要分布在毕节市大部、六盘水市北部和遵义市局地等地区, 终霜冻日空间上呈西部晚、东部早, 北部晚、南部早趋势。(3) 1961— 2018年, 贵州地区偏早初霜冻平均发生概率为7%, 空间上整体呈现北部高、南部低的分布特征; 特早初霜冻的平均发生概率为17%, 空间上呈现西北部高, 西南部边缘和南部边缘低的分布特征; 偏晚终霜冻平均发生概率为8%, 特晚终霜冻日的平均发生概率为25%, 空间上特晚终霜冻发生概率较高的区域主要为西部大部、黔东南州北部、铜仁市西部和遵义市东部等地。(4)海拔高度是影响贵州地区初霜冻日、终霜冻日和无霜冻期的关键因子; 秋季平均最低气温是影响初霜冻日的关键因子, 春季平均最低气温是影响终霜冻日和无霜冻期的关键因子。

在全球气候变暖的背景下, 贵州地区初霜冻日、终霜冻日和无霜冻期分别呈推迟、提前和延长的变化趋势, 这与许艳等[8]对中国西南地区霜期的分析结果基本一致, 但初霜冻日推迟幅度、终霜冻日提前的幅度和无霜期延长幅度明显大于山西[11]、陕西[22]、宁夏[15]和东北 [12,14]等地区。宁晓菊等[10]研究得出, 全国初霜日推迟幅度、终霜日提前幅度和无霜期延长幅度北方大于南方、东部大于西部, 本研究结果与此不一致。贵州受自身地理位置的影响, 海拔高度是影响初霜冻日、终霜冻日和无霜冻期的关键因子, 其结果与北方高纬地区的纬度是影响初霜冻日、终霜冻日和无霜冻期关键因子[14]不一致, 这是由于贵州气候的地域性差异表现为水平距离近但地形起伏较大的山区, 气温随海拔高度的升高而降低, 立体气候特征明显。

此外, 本研究仅对贵州地区初霜冻日、终霜冻日和无霜冻期气候变化趋势进行了定量分析, 未对霜冻危害程度进行研究; 初霜冻日的推后、终霜冻日的提前和无霜期的延长, 为贵州农业生产提供了充足的热量资源, 但有可能使农业产业结构发生调整。因此, 霜冻灾害影响评估及其对农业产业发展调整方面有待进一步研究。

(责任编辑 侯春晓)

参考文献
[1] 明惠青, 唐亚平, 孙婧, . 近50年辽宁无霜期积温时空演变特征[J]. 干旱地区农业研究, 2011, 29(2): 276-280.
MING H Q, TANG Y P, SUN J, et al. Analysis of spatio-temporal evolvement characteristics of frostless period accumulated temperature in Liaoning in recent 50 years[J]. Agricultural Research in the Arid Areas, 2011, 29(2): 276-280. (in Chinese with English abstract) [本文引用:1]
[2] 戴君虎, 王焕炯, 葛全胜. 近50年中国温带季风区植物花期春季霜冻风险变化[J]. 地理学报, 2013, 68(5): 593-601.
DAI J H, WANG H J, GE Q S. Changes of spring frost risks during the flowering period of woody plants in temperate monsoon area of China over the past 50 years[J]. Acta Geographica Sinica, 2013, 68(5): 593-601. (in Chinese with English abstract) [本文引用:1]
[3] 中华人民共和国气象行业标准. QX/T 88-2008作物霜冻害等级[S]. 北京: 中国气象出版社, 2008. [本文引用:1]
[4] 韩荣青, 李维京, 艾婉秀, . 中国北方初霜冻日期变化及其对农业的影响[J]. 地理学报, 2010, 65(5): 525-532.
HAN R Q, LI W J, AI W X, et al. The climatic variability and influence of first frost dates in Northern China[J]. Acta Geographica Sinica, 2010, 65(5): 525-532. (in Chinese with English abstract) [本文引用:1]
[5] ERLAT E, TÜRKEŞ M. Analysis of observed variability and trends in numbers of frost days in Turkey for the period 1950-2010[J]. International Journal of Climatology, 2012, 32(12): 1889-1898. [本文引用:1]
[6] TERANDO A, EASTERLING W E, KELLER K, et al. Observed and modeled twentieth-century spatial and temporal patterns of selected agro-climate indices in north America[J]. Journal of Climate, 2012, 25(2): 473-490. [本文引用:1]
[7] 王国复, 许艳, 朱燕君, . 近50年我国霜期的时空分布及变化趋势分析[J]. 气象, 2009, 35(7): 61-67.
WANG G F, XU Y, ZHU Y J, et al. The spatial and temporal characteristics and trend of frost days over China for the past 50 years[J]. Meteorological Monthly, 2009, 35(7): 61-67. (in Chinese with English abstract) [本文引用:1]
[8] 许艳, 王国复, 王盘兴. 近50 a中国霜期的变化特征分析[J]. 气象科学, 2009, 29(4): 4427-4433.
XU Y, WANG G F, WANG P X. Climatic change of frost in China in recent 50 a[J]. Scientia Meteorologica Sinica, 2009, 29(4): 4427-4433. (in Chinese with English abstract) [本文引用:2]
[9] 叶殿秀, 张勇. 1961—2007年我国霜冻变化特征[J]. 应用气象学报, 2008, 19(6): 661-665.
YE D X, ZHANG Y. Characteristics of frost changes from 1961 to 2007 over China[J]. Journal of Applied Meteorological Science, 2008, 19(6): 661-665. (in Chinese with English abstract) [本文引用:1]
[10] 宁晓菊, 张丽君, 杨群涛, . 1951年以来中国无霜期的变化趋势[J]. 地理学报, 2015, 70(11): 1811-1822.
NING X J, ZHANG L J, YANG Q T, et al. Trends in the frost-free period in China from 1951 to 2012[J]. Acta Geographica Sinica, 2015, 70(11): 1811-1822. (in Chinese with English abstract) [本文引用:2]
[11] 李芬, 张建新, 武永利, . 近50年山西终霜冻的时空分布及其影响因素[J]. 地理学报, 2013, 68(11): 1472-1480.
LI F, ZHANG J X, WU Y L, et al. Spatial and temporal distribution and its impact factors of the last frost over Shanxi Province from 1961 to 2010[J]. Acta Geographica Sinica, 2013, 68(11): 1472-1480. (in Chinese with English abstract) [本文引用:2]
[12] 任景全, 郭春明, 王丽伟, . 1961—2015年吉林省霜时空分布特征及其影响因素[J]. 气象与环境学报, 2018, 34(4): 119-125.
REN J Q, GUO C M, WANG L W, et al. Spatial and temporal distribution characteristics and impact factors of the frost over Jilin Province from 1961 to 2015[J]. Journal of Meteorology and Environment, 2018, 34(4): 119-125. (in Chinese with English abstract) [本文引用:1]
[13] 刘静, 殷淑燕. 1960—2014年秦岭南北无霜期时空变化特征及对比分析[J]. 长江流域资源与环境, 2017, 26(4): 615-623.
LIU J, YIN S Y. Spatio-temporal variations and comparative analysis of the frost-free period in the northern and southern regions of the Qinling Mountains from 1960 to 2014[J]. Resources and Environment in the Yangtze Basin, 2017, 26(4): 615-623. (in Chinese with English abstract) [本文引用:1]
[14] 周晓宇, 赵春雨, 崔妍, . 1961—2013年中国东北地区初终霜日及无霜期的气候变化特征[J]. 自然资源学报, 2017, 32(3): 494-506.
ZHOU X Y, ZHAO C Y, CUI Y, et al. Variation characteristics of the date of first frost and last frost and the frost-free period in the northeast China during 1961-2013[J]. Journal of Natural Resources, 2017, 32(3): 494-506. (in Chinese with English abstract) [本文引用:2]
[15] 张磊, 张晓煜, 李红英, . 1961—2010年宁夏无霜期变化特征[J]. 生态环境学报, 2013, 22(5): 801-805.
ZHANG L, ZHANG X Y, LI H Y, et al. Characteristics of frost-free days changes over Ningxia from 1961 to 2010[J]. Ecology and Environmental Sciences, 2013, 22(5): 801-805. (in Chinese with English abstract) [本文引用:2]
[16] 马尚谦, 张勃, 杨文义, . 1960—2015年淮河流域异常初、终霜日时空变化及其影响因子[J]. 自然资源学报, 2019, 34(2): 359-373.
MA S Q, ZHANG B, YANG W Y, et al. Spatial and temporal changes of the abnormal first and last frosts and their influencing factors in Huaihe River Basin from 1960 to 2015[J]. Journal of Natural Resources, 2019, 34(2): 359-373. (in Chinese with English abstract) [本文引用:1]
[17] 郜红娟, 蔡广鹏, 陆艺, . 1960—2013年贵州省霜冻日数时空变化特征[J]. 农业科学研究, 2017, 38(4): 16-19.
GAO H J, CAI G P, LU Y, et al. Spatial and temporal variations of frost days in Guizhou Province from 1960 to 2013[J]. Journal of Agricultural Sciences, 2017, 38(4): 16-19. (in Chinese with English abstract) [本文引用:1]
[18] 吴红梅. 精准脱贫背景下现代山地农业产业体系构建机制探索: 基于贵州“校农结合”扶贫产业实践的思考[J]. 安徽农业大学学报(社会科学版), 2018, 27(6): 16-21.
WU H M. Construction of modern mountain agricultural industry system under the background of targeted poverty alleviation: reflections based on guizhou's “university-agriculture cooperation” model of industrialization poverty alleviation[J]. Journal of Anhui Agricultural University ( Social Sciences), 2018, 27(6): 16-21. (in Chinese with English abstract) [本文引用:1]
[19] 吴洪宝, 吴蕾. 气候变率诊断和预测方法[M]. 北京: 气象出版社, 2010. [本文引用:1]
[20] 符淙斌, 王强. 气候突变的定义和检测方法[J]. 大气科学, 1992, 16(4): 482-493.
FU C B, WANG Q. The definition and detection of the abrupt climatic change[J]. Chinese Journal of Atmospheric Sciences, 1992, 16(4): 482-493. (in Chinese with English abstract) [本文引用:1]
[21] 范晓辉, 王麒翔, 王孟本. 山西省近50年无霜期变化特征研究[J]. 生态环境学报, 2010, 19(10): 2393-2397.
FAN X H, WANG Q X, WANG M B. Variations in frost-free period in Shanxi Province, China: 1959-2008[J]. Ecology and Environmental Sciences, 2010, 19(10): 2393-2397. (in Chinese with English abstract) [本文引用:1]
[22] 柏秦凤, 李星敏, 朱琳. 近50年陕西省无霜期的变化及果区霜冻风险分布[J]. 干旱区资源与环境, 2013, 27(8): 65-70.
BAI Q F, LI X M, ZHU L. The changes of the frost-free periods from 1961 to 2010 and its impact on apple industry in Shaanxi Province[J]. Journal of Arid Land Resources and Environment, 2013, 27(8): 65-70. (in Chinese with English abstract) [本文引用:1]