山西省無霜期時空變化特征研究

張蕾 韓小英 蘆艷珍 高洋 楊三維

摘要 依據近38年山西109個氣象站逐日地面最低溫度資料，計算得到初霜日、終霜日和無霜期，通過正交函數分析模型和小波分析等方法擬合時空變化，并對分布特征進行分析。結果表明：無霜期距平主要模態呈東北向西南、山區向盆地增加;時間序列表現為一定的年代際變化，無霜期呈延長趨勢，線性傾向率為0.8 d/10年;未來一段時期內無霜期小幅變長，伴隨終霜日提前和初霜日延后。本研究旨在為未來農業生產布局、農業自然災害應對提供決策參考。

關鍵詞 山西省;無霜期;時空變化;特征分析

中圖分類號：S161 文獻標識碼：A 文章編號：2095–3305（2020）07–0–03

DOI：10.19383/j.cnki.nyzhyj.2020.07.037

Study on the Temporal-Spatial variations of Frost free Period in Shanxi from 1979 to 2016

ZHANG Lei et al （College of Agricultural Economics and Management of Shanxi Agricultural university， Taiyuan， Shanxi 030006）

Abstract Based on the daily minimum ground surface temperature datasets of 109 meteorological stations in Shanxi Province during the period 1979-2016， the temporal-spatial distribution characteristics of the first frost day（FFD）， the last frost day（LFD） and the frost free period （FFP） are analyzed by using Empirical Orthogonal Function （EOF） and wavelet analysis. The results show that： ①The regional distribution of FFP departure is increasing from northeast to southwest and from mountain to basin; the time evolution trend is， the FFP exhibits slight decadal variation and extension trend， with increase rate of 0.8d/10a; ②The oscillation in FFP will become longer essentially in the future， along with FFD delay and LFD advance. This paper aims at providing a reference for the distribution of agricultural production and avoiding agricultural natural disasters in Shanxi in the future.

Key words Shanxi Province; Frost free period; Temporal-spatial distribution characteristics; Characteristic analysis

無霜期是一年里終霜日起至初霜日止，期間不再有霜出現的持續日數，其長短與作物生長期密不可分，反映當地熱量資源狀況，直接關乎農作物的種植熟制、作物布局和品種選擇。不同地域無霜期長短不同，年際變化較大，通常用當地無霜期多年平均值來表示。

據聯合國政府間氣候變化專門委員會（IPCC）公布的第5次氣候評估報告，1880—2012年全球地表平均溫度約升高0.85℃[1]。2009年哥本哈根聯合國氣候大會將全球氣候變化提上國際議案。2019年聯合國秘書長古特雷斯發出警告：全球變暖即將步入臨界點。20世紀80年代以來，我國學者采用不同資料和研究方法，深入研究我國氣溫變化。寧曉菊等[2]分析初、終霜日及無霜期變化特征，指出我國80%以上區域呈初霜日推后、終霜日提前和無霜期延長的趨勢，無霜期隨緯度或海拔增加而減少，年際波動隨緯度增加或隨海拔降低而減少;許艷等[3]認為，中國大部分地區霜期有縮短趨勢;李艷麗等[4]總結了山西霜凍災害的防治對策和措施;李芬等[5]發現山西省平均終霜日以波動為主，中西部和南部提前幅度較大，西北部和中東部推后幅度較大;王正旺等[6]分析山西東南部霜凍特征認為，無霜期在氣候變暖背景下有延長趨勢。

山西地形復雜多樣，地勢東北高西南低，多數地區海拔在1 000～2 000 m，晝夜溫差較大，霜凍災害多發。隨著全球變暖，各地初、終霜日和無霜期出現新變化，缺乏最新時段研究成果。

利用1979—2016年間山西109個地面氣象觀測站的逐日地面最低溫度數據，運用經驗正交函數分析（EOF）和Morlet小波分析方法，研究山西近38年來無霜期的時空變化，為優化農業產業空間布局、種植結構調整和提升有機旱作農業發展建設水平等提供決策參考。

1 資料與方法

1.1 資料來源與處理

以中國地面氣候資料日值數據集（V3.0）1979—2016年間山西省109個氣象站（圖1）逐日地面最低溫度和各站地理信息資料為基礎，通過數據核校、修正，建立標準數據庫，確保所用資料的系統、完整和可靠性。以地面最低溫度≤0℃為霜日標準，對各站點初、終霜日和無霜期日數進行統計分析。將日期序列統一轉換為數字序列，賦值1～365（閏年為1～366），如1月1日賦值為1，1月2日賦值為2，…，依此類推，對新的初、終霜日期和無霜期序列進行EOF分析，通過主要模態探究其時空分布特征;對所有站點逐年數據平均值進行小波分析，研究其振蕩周期與時間演變規律，預測未來氣候的演變趨勢，并分析對農業生產的影響。

1.2 研究方法

用EOF分析山西無霜期的主要結構特征，1979—2016年各站無霜期場作為時間-空間函數，即：

i=1，2，…，m j=1，2，…n

式中，m為109個氣象站，n為38個時間點。將X分解為時間函數Z和空間函數V兩部分，進行氣象向量場的自然正交展開：

X=VZ（2）

（3）

其中：

Morlet小波分析法用于獲得信息的頻率特征，以109站平均初、終霜日及無霜期序列作為Xn（n=1，…，N）的連續小波變換，定義為小波函數ψ0尺度化及xn的卷積：

式中，*表示共軛復數，N是時間系列個數（38年），用于小波函數標準化。

通過小波分析，得到時間系列振幅和相位的信息，具有非正交性。

式中t為時間，ω0是無量綱頻率。

2 結論與分析

2.1 初霜日、終霜日及無霜期空間分布

經統計得1979—2016年山西省初霜日、終霜日情況（表1）。

圖2a表明，山西終霜日由西南向東北部逐漸推后：東西部山地丘陵區、西南部及中部盆地平原區（忻定平原、太原盆地等）基本在4月底前結束霜期，太行山、呂梁山以及北部較晚，部分地區延續到5月底之后;初霜（圖2b）由東北向西南依次出現，東北部、太行山和呂梁山海拔較高地區最早進入霜期，在9月底之前，西南部丘陵區、中部盆地區及東西兩側平原較晚，西南部的運城市、臨猗縣、永濟市等地10月27日后才進入霜期。山西省無霜期整體呈現為東北部和高海拔地區較短、西南部及盆地區較長（圖2c），呂梁山、五臺山等高海拔地區無霜期不足80 d，運城地區則長達200 d以上。

2.2 無霜期距平模態時空分布

無霜期距平EOF第一模態（圖3）解釋了近38年山西省無霜期變化的57.45%。無霜期呈一致延長趨勢，其空間分布為：南部變化梯度大于北部，西部整體增速大于東部（臨縣、柳林縣、石樓縣一帶無霜期延長23 d以上，東部的陽泉市、平定縣、昔陽縣等地增加21 d）。結合時間序列（PC1），20世紀末前后，無霜期一致縮短變為一致延長，線性擬合表明，近38年無霜期增長速率為0.8 d/10年，相關度達90%。無霜期延長側面反映了氣候變暖，伴隨初霜日期推后及終霜日期提前。

EOF2（圖略）解釋方差為8.05%，空間分布呈帶狀，表現為南部、北部縮短，中部延長的“﹣+﹣”三極子型。PC2反映一定的10～15年內振蕩的年際變化。EOF3（圖略）空間分布表現為“南減北增”的蹺蹺板型，以臨縣—太原—盂縣為界，以北區域無霜期增加，以南地區無霜期縮短，PC3以不等幅年際振蕩為主。

2.3 初、終霜日和無霜期序列的突變分析及時間演變規律

根據山西無霜期時間序列的小波分析結果（圖4c）：近38年來無霜期在約14年尺度上存在明顯正負變化波動，為主要周期;在5～8年尺度上變化較為明顯，為次要周期;2000年以來，無霜期出現大于30年尺度的明顯正波動，與初霜日推遲有關。綜上：分別以1979—1991年、1991—2004年和2004—2016年為周期，無霜期均由負距平轉為正距平，依此推測，2016年后的3～5年，全省無霜期仍為增加趨勢，在2020年左右開始新的負距平周期。

3 討論

通過對山西省無霜期、初霜日期和終霜日期的時空變化特征分析，結論如下：（1）全省無霜期分布呈東北向西南、山區向平川增加的趨勢。即東北部、太行和呂梁山區等較短，西南及中部盆地較長。（2）無霜期主要的空間模態為西南部增長快于東北部，其大值中心位于中部偏西的臨縣、柳林和石樓縣等地，無霜期在20世紀末之前小幅縮短，之后延長。（3）小波分析表明，無霜期存在明顯的正負波動，主周期約為14年，預計未來幾年山西省初霜日較常年推遲，終霜日提前，無霜期將小幅變長。

理論上講，無霜期延長（初霜日延后、終霜日提前）可減少春秋霜凍對作物生長的危害，使作物生長期變長，復種指數提升，晚熟作物品種面積增加，作物品種布局和選擇范圍擴大、熟制增加，總體上對農業生產有利，但具體種植結構的調整與優化還要視耕作區種植制度、品種特性和作物生長發育的節率與氣候變化的協調與否來決定。例如，近年來局部地區早霜凍天氣和倒春寒天氣時有發生，一定程度上增加了作物品種布局的復雜性和不確定性，建議結合氣候變化規律和當地生產實際，因地制宜合理布局農業生產，科學調整耕作模式。

參考文獻

[1] Pachauri K，Meyer A.Climate Change 2014 Synthesis Report[J].Environmental Policy Collection， 2014，27（2）：408.

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

[3] 許艷，王國復，王盤興.近50年中國霜期的變化特征分析[J].氣象科學，2009，29（4）：427–433.

[4] 李艷麗，王迎春，孫忠富.山西省霜凍災害現狀及其防御對策分析[J].中國農業資源與區劃，2006（1）：57–59.

[5] 李芬，張建新，武永利，等.近50年山西終霜凍的時空分布及其影響因素[J].地理學報，2013，68（11）：1472–1480.

[6] 王正旺，劉新瑩，武志躍，等.山西東南部氣候變暖背景下的霜凍特征分析[J].干旱地區農業研究，2014，32（5）：234–238.

責任編輯：黃艷飛