摘 要：根據(jù)1959—2018年滁州市及所轄6個(gè)縣（市）共7個(gè)國(guó)家氣象觀測(cè)站地面氣象觀測(cè)資料，運(yùn)用Mann-Kendall趨勢(shì)檢驗(yàn)、非參數(shù)突變檢驗(yàn)及Morlet小波分析等方法，探究滁州市平均初霜日、終霜日的變化趨勢(shì)、突變特征及其周期。結(jié)果表明：滁州市平均初霜日的推遲幅度為2.22 d/10 a，平均終霜日的提前幅度為1.86 d/10 a；平均初霜日存在明顯的突變點(diǎn)，突變年份為2003年，平均終霜日突變年份為1996年；在28 a特征時(shí)間尺度上，平均初霜日變化周期約為20 a，約經(jīng)歷了3個(gè)提前—推遲趨勢(shì)轉(zhuǎn)換期；在28 a特征時(shí)間尺度上，平均終霜日變化周期約為20 a；在12 a特征時(shí)間尺度上，平均終霜日變化周期約為10 a。

關(guān)鍵詞：霜；趨勢(shì)檢驗(yàn)；突變；周期分析

中圖分類號(hào)：P423.4 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1674-7909（2024）13-121-4

DOI：10.19345/j.cnki.1674-7909.2024.13.027

0 引言

霜是指水汽在地面和近地面物體上凝華而成的白色松脆的冰晶，或是由露凍結(jié)而成的冰珠［1］，易在晴朗微風(fēng)的夜間生成。氣象學(xué)上一般將秋季第一次出現(xiàn)的霜稱為“初霜”，春季最后一次出現(xiàn)的霜稱為“終霜”。在全球氣候變暖的背景下，研究初、終霜出現(xiàn)日期的變化，對(duì)合理安排農(nóng)事活動(dòng)、降低農(nóng)業(yè)生產(chǎn)風(fēng)險(xiǎn)具有重要意義［2］。馬尚謙等［3］研究發(fā)現(xiàn)，近56 a淮河流域平均初霜日以2.15 d/10 a的速率呈推遲趨勢(shì)，終霜日以2.49 d/10 a的速率呈提前趨勢(shì)；寧曉菊等［4］根據(jù)1951年以來國(guó)內(nèi)824個(gè)氣象站點(diǎn)的資料分析發(fā)現(xiàn)，中國(guó)80%以上區(qū)域呈現(xiàn)初霜日推后、終霜日提前的趨勢(shì)；潘淑坤等［5］、李輯等［6］、羅靜等［7］、姬興杰等［8］分別對(duì)新疆、遼寧、陜西、河南等地區(qū)的霜期變化進(jìn)行研究，結(jié)果均與上述趨勢(shì)相符。

滁州市位于皖東江淮之間，地勢(shì)西高東低，屬亞熱帶濕潤(rùn)季風(fēng)氣候區(qū)，是國(guó)家大型商品糧生產(chǎn)基地。滁州市現(xiàn)轄天長(zhǎng)、明光2市，來安、全椒、定遠(yuǎn)、鳳陽(yáng)4縣，以及瑯琊、南譙2區(qū)。全市共建有國(guó)家氣象觀測(cè)站7個(gè)，分別位于市本級(jí)瑯琊區(qū)和所轄6個(gè)縣（市）。目前，針對(duì)滁州市初霜日、終霜日變化特征的研究尚屬空白，且現(xiàn)有研究［9-11］多以日最低氣溫或地表最低溫度判斷霜日，而除溫度外，相對(duì)濕度和風(fēng)速也是霜形成的重要條件。以1959—2018年滁州市7個(gè)國(guó)家氣象觀測(cè)站霜的人工觀測(cè)資料為基礎(chǔ)，對(duì)當(dāng)?shù)爻跛?、終霜日的時(shí)空變化特征進(jìn)行分析，為合理利用農(nóng)業(yè)氣候資源和有效防御霜凍提供一定的科學(xué)依據(jù)。

1 資料與方法

1959—2018年滁州市7個(gè)國(guó)家氣象觀測(cè)站霜的人工觀測(cè)數(shù)據(jù)（有霜日記錄1，無(wú)霜日記錄0）來源于全國(guó)綜合氣象信息共享平臺(tái)（CIMISS）。經(jīng)篩選，剔除無(wú)霜日，確定初霜日、終霜日。根據(jù)原始數(shù)據(jù)，分別以年度總天數(shù)中的第幾天表示初霜日和終霜日，重新建立各站初霜日、終霜日的數(shù)據(jù)序列［5］。例如，1959年3月17日為終霜日，10月30日為初霜日，則該年度終霜日和初霜日分別以75和302計(jì)。以新建立的初霜日、終霜日的數(shù)據(jù)序列為基礎(chǔ)，進(jìn)而對(duì)所有數(shù)據(jù)序列的變化特征進(jìn)行分析。

Mann-Kendall（簡(jiǎn)稱M-K）檢驗(yàn)法是時(shí)間序列數(shù)據(jù)趨勢(shì)檢驗(yàn)中使用廣泛的非參數(shù)檢驗(yàn)方法，適用于水文、氣象等非正態(tài)分布的數(shù)據(jù)分析，計(jì)算簡(jiǎn)便［12］。研究利用該方法分析滁州市初霜日、終霜日的時(shí)間變化特征。小波分析也稱多分辨率分析，可以反映系統(tǒng)在不同周期中的變化趨勢(shì)，并能對(duì)系統(tǒng)未來發(fā)展趨勢(shì)進(jìn)行定性估計(jì)［13］。研究采用Morlet小波分析法探究滁州市初霜日、終霜日的周期變化特征。M-K檢驗(yàn)與Morlet小波分析均使用Matlab（R2017a）軟件完成。

2 結(jié)果與分析

2.1 初霜日變化特征

2.1.1 初霜日變化趨勢(shì)

對(duì)1959—2018年滁州市平均初霜日的人工觀測(cè)數(shù)據(jù)進(jìn)行Mann-Kendall趨勢(shì)檢驗(yàn)，結(jié)果見表1（Z為正值表示增加趨勢(shì)，負(fù)值表示減少趨勢(shì)；Z的絕對(duì)值在大于等于1.28、1.64、2.32時(shí)，分別表示通過了置信度90%、95%、99%的顯著性檢驗(yàn)）［12］。由表1可知，滁州市平均初霜日的變化趨勢(shì)通過了置信度99%的顯著性檢驗(yàn)，存在顯著的推遲趨勢(shì)（Z為正值），平均初霜日的推遲幅度為2.22 d/10 a。

2.1.2 初霜日突變分析

對(duì)1959—2018年滁州市平均初霜日的人工觀測(cè)數(shù)據(jù)進(jìn)行Mann-Kendall突變檢驗(yàn)，給定顯著性水平α=0.05，即μ0.05=±1.96，UF統(tǒng)計(jì)量用于衡量數(shù)據(jù)序列中上升和下降極值的數(shù)量之差，UB統(tǒng)計(jì)量用于衡量極值之間的時(shí)間間隔，結(jié)果如圖1所示。滁州市平均初霜日自1996年起呈持續(xù)推遲趨勢(shì)，并于2005年達(dá)到顯著水平（超過顯著性水平0.05臨界線），根據(jù)UF和UB曲線交點(diǎn)的位置確定突變點(diǎn)為2003年，突變前（1959—2002年）平均初霜日為11月5日，突變后（2003—2018年）平均初霜日為11月13日，推遲8 d。

2.1.3 初霜日周期變化

由上述分析可知，滁州市平均初霜日存在顯著的變化趨勢(shì)和突變點(diǎn)，選擇 Matlab 小波工具箱中的Morlet 復(fù)小波函數(shù)對(duì)平均初霜日進(jìn)行周期分析，結(jié)果如圖2所示。圖2（a）可以反映初霜日時(shí)間序列在不同時(shí)間尺度上的周期變化，實(shí)線表示小波系數(shù)實(shí)部值為正，虛線為負(fù)，對(duì)應(yīng)初霜日處于推遲或提前期。由圖2（a）可知，滁州市平均初霜日變化具有明顯的周期性特征，主要表現(xiàn)為25～30 a尺度上出現(xiàn)的推遲—提前交替的準(zhǔn)3次振蕩，在整個(gè)時(shí)間尺度上出現(xiàn)4個(gè)偏多中心和3個(gè)偏少中心，且這一時(shí)間尺度的周期變化在整個(gè)分析時(shí)段表現(xiàn)得非常穩(wěn)定，具有全域性。圖2（b）可以反映初霜日時(shí)間序列的波動(dòng)能量隨尺度的分布情況，以確定初霜日變化過程中存在的主周期。小波方差圖2（b）中只存在一個(gè)明顯的峰值，對(duì)應(yīng)28 a的時(shí)間尺度，說明28 a左右（時(shí)間尺度）的周期振蕩最強(qiáng)，為初霜日變化的主周期。根據(jù)小波方差檢驗(yàn)的結(jié)果，繪制平均初霜日變化過程的主周期趨勢(shì)圖［見圖2（c）］，在28 a特征時(shí)間尺度上，初霜日變化的平均周期約為20 a，大約經(jīng)歷了3個(gè)提前—推遲趨勢(shì)轉(zhuǎn)換期。

2.2 終霜日變化特征

2.2.1 終霜日變化趨勢(shì)

對(duì)1959—2018年滁州市平均終霜日的人工觀測(cè)數(shù)據(jù)進(jìn)行Mann-Kendall趨勢(shì)檢驗(yàn)，結(jié)果見表2。由表2可知，滁州市平均終霜日的變化趨勢(shì)通過了置信度99%的顯著性檢驗(yàn)，表現(xiàn)出顯著的提前趨勢(shì)（Z為負(fù)值），平均終霜日的提前幅度為1.86 d/10 a。

2.2.2 終霜日突變分析

對(duì)1959—2018年滁州市平均終霜日的人工觀測(cè)數(shù)據(jù)進(jìn)行Mann-Kendall突變檢驗(yàn)，結(jié)果如圖3所示。滁州市平均終霜日突變年份為1996年，突變前（1959—1995年）平均終霜日為4月1日，突變后（1996—2018年）平均終霜日為3月25日，提前7 d。

2.2.3 終霜日周期變化

由上述分析可知，滁州市平均終霜日存在顯著的變化趨勢(shì)和突變點(diǎn)，選擇Matlab 小波工具箱中的Morlet 復(fù)小波函數(shù)對(duì)全市平均終霜日進(jìn)行周期分析，結(jié)果如圖4所示。由圖4（a）可知，滁州市平均終霜日變化具有明顯的周期性特征，具體表現(xiàn)為24～32 a尺度上出現(xiàn)的推遲—提前交替的準(zhǔn)3次振蕩，在9～18 a尺度上出現(xiàn)的準(zhǔn)7次振蕩，在3～9 a尺度上出現(xiàn)的準(zhǔn)14次振蕩。小波方差圖4（b）中存在3個(gè)較為明顯的峰值，其中最大峰值對(duì)應(yīng)28 a的時(shí)間尺度，說明28 a左右（時(shí)間尺度）的周期振蕩最強(qiáng)，為終霜日變化的第一主周期，12 a、6 a的時(shí)間尺度依次為第二、第三主周期。根據(jù)小波方差檢驗(yàn)的結(jié)果，繪制平均終霜日變化過程的第一、第二主周期趨勢(shì)圖［見圖4（c）］，在28 a特征時(shí)間尺度上，終霜日變化的平均周期約為20 a，在12 a特征時(shí)間尺度上，終霜日變化的平均周期約為10 a。

3 結(jié)論

對(duì)滁州市初霜日、終霜日變化特征進(jìn)行統(tǒng)計(jì)分析，得出如下結(jié)論。

①滁州市平均初霜日存在顯著的推遲趨勢(shì)，推遲幅度為2.22 d/10 a；平均終霜日表現(xiàn)出顯著的提前趨勢(shì)，提前幅度為1.86 d/10 a。

②平均初霜日存在明顯的突變點(diǎn)，突變年份為2003年，突變后平均初霜日為11月13日，推遲8 d；平均終霜日也存在明顯的突變點(diǎn)，突變年份為1996年，突變后平均終霜日為3月25日，提前7 d。

③平均初霜日主要表現(xiàn)為25～30 a時(shí)間尺度上的年代際周期，在28 a特征時(shí)間尺度上，初霜日變化的平均周期約為20 a，大約經(jīng)歷了3個(gè)提前—推遲趨勢(shì)轉(zhuǎn)換期；平均終霜日具體表現(xiàn)為24～32 a、9～18 a及3～9 a時(shí)間尺度上的年代際周期，在28 a特征時(shí)間尺度上，終霜日變化的平均周期約為20 a，在12 a特征時(shí)間尺度上，終霜日變化的平均周期約為10 a。

參考文獻(xiàn)：

［1］中國(guó)氣象局.地面氣象觀測(cè)規(guī)范［M］.氣象出版社：北京，2003：21-22.

［2］王媛媛，張勃.1971年至2010年隴東地區(qū)霜凍與積溫變化特征［J］.資源科學(xué)，2012，34（11）：2181-2188.

［3］馬尚謙，張勃，唐敏，等.1960—2015年淮河流域初終霜日時(shí)空變化分析［J］.中國(guó)農(nóng)業(yè)氣象，2018，39（7）：468-478.

［4］寧曉菊，張麗君，楊群濤，等.1951年以來中國(guó)無(wú)霜期的變化趨勢(shì)［J］.地理學(xué)報(bào)，2015，70（11）：1811-1822.

［5］潘淑坤，張明軍，汪寶龍，等.1960—2011年新疆初終霜日及無(wú)霜期的變化特征［J］.干旱區(qū)研究，2013，30（4）：735-742.

［6］李輯，嚴(yán)曉瑜，王穎.遼寧省近50年霜的氣候變化特征［J］.氣象，2010，36（11）：38-45.

［7］羅靜，費(fèi)杜秋.1961—2017年陜西省初終霜日時(shí)空變化特征分析［J］.安徽農(nóng)業(yè)科學(xué)，2019，47（20）：228-234，254.

［8］姬興杰，李鳳秀，王紀(jì)軍.1971—2010年河南省霜期的時(shí)空分布特征及其對(duì)氣溫的響應(yīng)［J］.氣象與環(huán)境學(xué)報(bào)，2015，31（1）：67-75.

［9］郜紅娟，蔡廣鵬，陸藝，等.1960—2013年貴州省霜凍日數(shù)時(shí)空變化特征［J］.農(nóng)業(yè)科學(xué)研究，2017，38（4）：16-19，42.

［10］周曉宇，趙春雨，崔妍，等.1961—2013年中國(guó)東北地區(qū)初終霜日及無(wú)霜期的氣候變化特征［J］.自然資源學(xué)報(bào)，2017，32（3）：494-506.

［11］拉巴次仁，索朗加措，白瑪.1981—2010年西藏霜凍日數(shù)的變化特征［J］.地理學(xué)報(bào)，2014，69（5）：690-696.

［12］張盛霖，鄧高燕，黃勇奇.Mann-Kendall檢驗(yàn)法在Excel中的實(shí)現(xiàn)與應(yīng)用［EB/OL］.北京：中國(guó)科技論文在線［2014-06-27］.https：//www.paper.edu.cn/releasepaper/content/201406-448.

［13］王朋朋，程杰，袁峰，等.時(shí)間序列模型分析威海市近50年氣溫變化特征［J］.南水北調(diào)與水利科技，2015，13（6）：1051-1055.

作者簡(jiǎn)介：張?chǎng)瓮?987—），男，碩士，工程師，研究方向：應(yīng)用氣象。