滁州市初霜日、終霜日變化特征

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

關鍵詞：霜；趨勢檢驗；突變；周期分析

中圖分類號：P423.4 文獻標志碼：A 文章編號：1674-7909（2024）13-121-4

DOI：10.19345/j.cnki.1674-7909.2024.13.027

0 引言

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

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

1 資料與方法

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

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

2 結果與分析

2.1 初霜日變化特征

2.1.1 初霜日變化趨勢

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

2.1.2 初霜日突變分析

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

2.1.3 初霜日周期變化

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

2.2 終霜日變化特征

2.2.1 終霜日變化趨勢

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

2.2.2 終霜日突變分析

對1959—2018年滁州市平均終霜日的人工觀測數據進行Mann-Kendall突變檢驗，結果如圖3所示。滁州市平均終霜日突變年份為1996年，突變前（1959—1995年）平均終霜日為4月1日，突變后（1996—2018年）平均終霜日為3月25日，提前7 d。

2.2.3 終霜日周期變化

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

3 結論

對滁州市初霜日、終霜日變化特征進行統計分析，得出如下結論。

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

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

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

參考文獻：

［1］中國氣象局.地面氣象觀測規范［M］.氣象出版社：北京，2003：21-22.

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

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

［4］寧曉菊，張麗君，楊群濤，等.1951年以來中國無霜期的變化趨勢［J］.地理學報，2015，70（11）：1811-1822.

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

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

［7］羅靜，費杜秋.1961—2017年陜西省初終霜日時空變化特征分析［J］.安徽農業科學，2019，47（20）：228-234，254.

［8］姬興杰，李鳳秀，王紀軍.1971—2010年河南省霜期的時空分布特征及其對氣溫的響應［J］.氣象與環境學報，2015，31（1）：67-75.

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

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

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

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

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

作者簡介：張鑫童（1987—），男，碩士，工程師，研究方向：應用氣象。