1、收稿日期院2023-02-15作者简介院 王子一 渊1986-冤袁 男袁 吉林白城人袁 硕士学位袁 工程师袁 主要从事气候及农业气象相关研究袁遥基金项目院四平市气象局科研项目野四平市霜冻特征研究渊202202冤冶资助遥1引言霜冻是一种冻害袁当出现霜冻时袁植物附近的土壤和空气温度会在非常短的时间内降至 0益或以下袁直接导致植物内部的水分冻结袁破坏植物细胞的结构形态袁发生不可逆的损伤袁从而使其死亡1遥在农业上这种情况会直接导致整个区域减产或毁种遥 除了海南省之外全国几乎所有省份都会发生霜冻灾害遥 霜冻主要分为早霜冻和晚霜冻两种袁早霜冻一般出现在每年的秋季袁第一次出现的霜冻就叫做早霜冻或者初霜冻曰
2、而晚霜冻一般发生在每年的春季袁 最后一次出现的霜冻就叫做晚霜冻或者终霜冻遥 头年的初霜冻第一天到来年的终霜冻最后一天之间的日数称为霜期袁 而年初的终霜冻最后一天到秋天的初霜冻第一天之间的日数则称为无霜期2遥一般东北地区的霜冻成因主要与冷空气入侵导致的辐射冷却同时平流层冷却降温等因素有密切关系3遥四平市作为吉林省重要的粮食产区之一袁其粮食产量对吉林省的整体粮食产量有重要影响遥但每年冷空气影响四平市的时间较早袁 而且影响时间较长袁特别是每年初春极易发生霜冻袁春霜冻使农作物发芽及幼苗受到影响袁 从而造成农作物产量减产袁降低农民收入袁影响本地农业经济遥 由此可见袁 分析研究四平市的初霜冻变化规律及其发
3、生的气象学条件袁做好防灾减灾预报服务工作袁对于防御四平市农业气象灾害有重要意义遥2资料与方法2.1资料本文采用 1986要2020 年袁四平市城区及三个所辖县渊市冤的初尧终霜冻日期确定霜期日数遥将所得到的日期按照儒历日定义日期袁 并由此建立所选站点初渊终冤霜冻日的数据序列4遥2.2方法2.2.1Mann-Kendall 突变分析法Mann-Kendall 突变分析是以趋势分析为基础发展的监测序列变化趋势突变的方法袁 这是一种非参数统计检验方法遥 非参数检验也称为非分布检验袁该方法的优点是袁该样本是不需要遵守一个特定的分布袁并且它不受一些异常值影响袁更适合用于类型变量和顺序变量袁计算也相对较为简
4、单遥可用在长时间监测气象要素在数列的变化趋势和确定突变点袁且监测精度较高遥其优点是不需要服从一定的分布袁也不受少数异常值的干扰5-7遥2.2.2小波分析法小波分析是用一簇具有震荡性尧 能够迅速衰减到零的小波函数系来表示或逼近某一信号或函1986要2020 年四平市霜冻变化特征分析王子一1安红善2岳振铎1巩宪伟1渊1.四平市气象局袁吉林四平136000曰2.长春市气象局袁吉林长春130012冤摘要院 利用四平市 4 个观测站 1986要2020 年的初 渊终冤 霜冻日期资料袁 通过 Morlet 小波分析尧Mann-Kendall 突变分析等方法袁对四平市初渊终冤霜冻及霜期的演变规律尧变化特征和
5、四平市霜冻年际尧年代际变化及其周期特征进行了分析遥 研究表明院四平市近 35 a 的年霜期日数整体于 2000 年之后呈现出下降趋势曰2003要2020 年表现出霜期缩短的趋势袁初霜日整体呈现推后趋势袁终霜日的整体变化趋势却呈现提前的趋势曰因四平市区城市化程度高于其他 3 个县渊市冤袁城市下垫面减少袁平均温度高于梨树县尧伊通县和双辽市袁所以气候倾向率高于其他 3 个县渊市冤遥关键词院Morlet 小波分析曰霜冻曰霜期曰特征第 30 卷第 2 期2023 年 6 月气象灾害防御METEOROLOGICAL DISASTER PREVENTIONVol.30 No.2June.2023第 2 期图
6、 1四平渊a冤尧双辽渊b冤尧梨树渊c冤尧伊通渊d冤霜期日数周期变化数袁 小波方差图可用来确定信号中不同种尺度扰动的相对强度和存在的主要时间尺度袁即主周期遥3结果分析3.1霜期日数变化特征3.1.1平均霜期日数的变化特征为了进一步分析霜期变化特征以及分布特点袁根据四平市霜期日数做逐年特征变化分析袁可以看出袁 近 35a 四平市总体霜期日数波动比较剧烈袁1998 年到 2002 年变化最为显著袁 从 2002 年后总体呈现逐年减少的趋势遥根据四平市霜期日数变化曲线可以看出袁四平市整体呈现霜期缩短的趋势袁四平尧梨树尧伊通和双辽霜期日数的显著性检验相关系数分别为-0.773尧-0.373尧-0.523
7、 和-0.726袁 均在 0.01 显著性水平之内遥 年霜冻日数每 10a 分别减少16.3d尧17.5d尧13.1d 和 7.8d袁2003 年之前大致呈现出正距平趋势袁之后呈现出负距平趋势遥 虽然在整个 35a 略有波动袁但整体霜期日数下降袁可以从侧面反映出 1986要2020 年冬季温度在 2000年之后呈上升趋势袁这与 2000 年之后全球变暖袁气温升高有关8-9遥3.1.2霜期日数的周期特征将四平市 4 个台站近 35a 的年霜期日数作为被检测信号袁实现一维连续小波变换遥从图 1 可以看出袁 四平霜期日数存在一个25a 的变化周期袁周期性较强曰双辽霜期日数存在一个 17a 左右的主变
8、化周期袁周期强度相对较强袁同时在 1997 年之后还存在一个 7a 的变化周期袁相对于之前的 17a 的周期强度较弱曰 梨树存在一个 17a 的主周期袁 周期在 1996要2004 年相对较强曰伊通年霜期日数在 1986要2020 年不同尺度上的周期振荡呈现一个 25a 的周期变化袁 且在1986要2013 年较为明显袁1990要2004 年存在一个8a 左右的周期遥3.2初尧终霜日的变化特征根据四平市 1986要2020 年的气象数据袁 按照儒历日数排列袁将下半年地面 0cm 日最低温度臆0益的第一天定义为初霜日袁 可得到四平市初霜日变化曲线渊图 2冤遥 按照同样的方法袁将上半年最后一次出现
9、地面 0cm 日最低温度臆0益的那天定义为终霜日袁可得到四平市终霜日变化曲线渊图 3冤遥四平市整体初霜日呈现推后趋势袁 四平尧梨abcd353025201510519851990199520002005201020152020年份3025201510519851990199520002005201020152020年份353025201510519851990199520002005201020152020年份3025201510519851990199520002005201020152020年份王子一袁等院1986要2020 年四平市霜冻变化特征分析35气象灾害防御第 30 卷图 2四平市
10、初霜日变化图 3四平市终霜日变化树尧 双辽和伊通初霜日显著性检验相关系数分别为 0.449尧0.429尧0.191 和 0.520袁 只有双辽初霜日未通过显著性检验袁表明其推后趋势在 0.05 显著性水平上是不显著的遥 1998 年和 2012 年 4 站同时出现了两次整体的大幅度正距平遥 1998 年四平和梨树初霜日最为靠后袁 均为当年的第 292 天曰2012 年梨树初霜日最为靠后袁 为当年的第 305天遥 四平市 1986要2020 年最早初霜日出现在2002 年的双辽袁为当年的第 253 天曰特别值得注意的是袁在 2011 年四平市 4 站初霜日基本都在当年的第 260 和第 261
11、天出现袁初霜日基本重叠遥四平市终霜日的整体变化趋势呈现提前的趋势袁四平尧梨树尧双辽和伊通终霜日显著性检验相关系数分别为-0.772尧-0.275尧-0.557 和-0.680袁梨树的终霜日未通过显著性检验袁 表明其提前趋势在 0.05 显著性水平上是不显著的遥1990要1992 年四平市 4 站终霜日整体呈现正距平袁 表明这个期间的终霜日相对靠后袁其中 1991 年四平和伊通的终霜日最为靠后袁 分别为当年的第 157 和第 159天曰而在 1999要2001 年袁四平市 4 站终霜日整体呈现显著的负距平袁表明这段时间里袁四平市整体终霜日提前袁其中梨树提前幅度最大袁为 2001 年的第 53 天
12、遥1986要2020 年袁四平尧梨树尧伊通和双辽的平均霜期日数分别为 210.66d尧201.83d尧216.20d 和211.89d袁平均霜期日数最少的是梨树袁最多的是伊通遥 四平尧梨树尧伊通和双辽的平均初霜日分别为第 271.57 天尧 第 272.69 天尧 第 268.49 天和第270.17 天袁平均初霜日出现最早的是伊通袁最晚的是梨树遥 四平尧梨树尧伊通和双辽的平均终霜日分别为第 114.69 天尧第 107.97 天尧第 118.80 天和第113.86 天袁平均终霜日出现最早的是梨树袁最晚的是伊通遥四平和双辽的平均霜期日数以及平均初尧终霜日都非常接近遥3.3霜期日数的突变特征对
13、四平市 4 站的霜期日数进行 M-K 突变检验袁通过检验阈值在 95%的水平下袁Ub 和 Uf 曲线相交年份表明四平和双辽年霜冻日数在2000 年发生了一次突变袁梨树年霜冻日数在 1989 年发生了一次突变袁 伊通年霜冻日数在 2000 年和 2004年分别发生了一次突变遥 四平尧梨树和伊通 Uf 曲线分别于 2002 年和 2007 年开始完全突破显著性检验的临界袁 说明此后的时间这 3 个地区的霜期日数呈现出明显逐年缩短的趋势曰双辽 Uf 曲线仅在 2006 年和 2016 年突破临界值袁 之后虽整体霜期日数缩短袁但波动较大遥3.4霜期日数的周期特征3.4.1初霜日小波分析将四平市 4 个
14、台站近 35 年的年初霜日作为被检测信号袁实现一维连续小波变换袁得到 4 个县渊市冤初霜日的小波分析图渊图 4冤遥通过图 4 的小波分析可以看出袁 四平存在一个 20a 左右的主周期变化袁 在这个 20a 周期变化中 袁1986要1989 年 尧1996要2003 年 尧2010要2016年呈负位相袁 代表初霜日相对提前曰1990要1995年尧2004要2009 年尧2017要2020 年呈正位相袁代表初霜日相对推后遥 双辽存在一个 20a 左右的主周期变化袁同时在 1991要2014 年双辽初霜日还存在一个比较明显的 8a 左右的副周期遥 主周期中1986要1990 年尧1997要2004
15、年尧2008要2015 年呈现负位相袁表明双辽初霜日相对提前曰在 1991要1996 年尧2005要2007 年尧2016要2020 年呈现正位相袁表明双辽初霜日相对推后遥 梨树呈现一个 10a左右的周期变化遥 在这个 10a 的周期变化中袁31030029028027026025019851990199520002005201020152020四平梨树双辽伊通四平梨树双辽伊通19851990199520002005201020152020年份170150130110907050年份36第 2 期图 4四平渊a冤尧双辽渊b冤尧梨树渊c冤尧伊通渊d冤初霜日周期变化1986要1989 年尧1993
16、要1996 年尧1999要2001 年尧2005要2008 年尧2012要2016 年呈正位相袁 表示初霜日相对推后曰1990要1992 年尧1997要1998 年尧2002要2004 年尧2009要2011 年呈负位相袁 代表初霜日相对提前遥 伊通初霜日主要呈现一个 20a 左右的主周期变化且周期明显袁在 1989要2009 年还存在一个 10a 左右的副周期遥 主周期中 1986要1989 年尧1997要2002 年尧2011要2015 年呈现负位相袁表明伊通初霜日相对提前袁在 1990要1996 年尧2003要2010 年尧2016要2020 年呈正位相袁 表明伊通初霜日相对推后遥3.4.2终霜日小波分析将四平市 4 个台站近 35 年的年终霜日作为被检测信号袁实现一维连续小波变换袁得到 4 站终霜日的小波分析图渊图 5冤遥通过图 5 小波分析可以看出袁 四平终霜日存在一个 10a 左右的主周期变化遥 在这个 10a 周期变化中袁1986要1987 年尧1993要1995 年尧1998要2001 年尧2005要2007 年 尧2012要2015 年 呈 正 位相袁 说明终霜日相对