王萬(wàn)滿 多杰才旦

摘要 利用青海省共和縣氣象局觀測(cè)的1961—2017年逐日最低氣溫資料，采用線性傾向估計(jì)、絕對(duì)變率、保證率等分析方法，對(duì)共和地區(qū)≤0℃霜凍初日、終日、無(wú)霜期變化特征進(jìn)行分析，為有效利用共和地區(qū)農(nóng)業(yè)氣候資源提供科學(xué)參考依據(jù)，對(duì)該地區(qū)的農(nóng)業(yè)生產(chǎn)起到積極有效的指導(dǎo)和促進(jìn)作用。

關(guān)鍵詞 霜凍;氣候特征;對(duì)比分析;共和

中圖分類號(hào)：P468.02 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：2095-3305（2020）02-025-04

DOI： 10.19383/j.cnki.nyzhyj.2020.02.012

Comparative Analysis of Frost Climate Characteristics and Historical Characteristics in Qinghai in the Past 10 Years

WANG Wan-man? et al（Meteorological Bureau of Hainan Tibetan Autonomous Prefecture，Gonghe，Qinghai 813099）

Abstract Based on the daily minimum temperature data from 1961 to 2017 observed by the Meteorological Bureau of Gonghe County，the change characteristics of the first day，the last day，and the frost free period of the frost in Gonghe district at ≤0℃ were analyzed by using the linear tendency estimation，absolute variability and assurance rate，which provided a scientific reference for the effective utilization of the agricultural climate resources in Gonghe district and played a positive and effective role in guiding and promoting agricultural production in the region.

Key words? ?Frost;Climatic characteristics;Comparative analysis;Gonghe

霜凍是指植物生長(zhǎng)季節(jié)里因土壤表面和植株體溫度降低至0℃或0℃以下而引起植物凍害，是一種常見的農(nóng)業(yè)氣象災(zāi)害，歷來(lái)就是農(nóng)業(yè)科學(xué)、氣象科學(xué)及其他相關(guān)科學(xué)研究的重點(diǎn)問題[1-2]。作物剛開始生長(zhǎng)階段或成熟階段出現(xiàn)霜凍，都會(huì)對(duì)農(nóng)作物造成嚴(yán)重危害。為了有效地預(yù)防霜凍災(zāi)害，許多學(xué)者對(duì)霜凍的發(fā)生規(guī)律、形成機(jī)理、預(yù)防措施及對(duì)農(nóng)業(yè)生產(chǎn)的影響等方面進(jìn)行了研究[3-5]。筆者對(duì)共和地區(qū)≤0℃霜凍初日、終日、無(wú)霜期變化特征進(jìn)行分析，為有效利用共和地區(qū)農(nóng)業(yè)氣候資源提供科學(xué)參考依據(jù)，對(duì)該地區(qū)的農(nóng)業(yè)生產(chǎn)起到積極有效的指導(dǎo)和促進(jìn)作用。

1 資料與方法

1.1 資料來(lái)源

選取青海省共和縣氣象局觀測(cè)的1961—2017年共57年的逐日最低氣溫?cái)?shù)據(jù)，計(jì)算霜凍初終日期間的無(wú)霜期日數(shù)和同時(shí)期的積溫。采用1981—2010年的30年平均值為多年平均值。把日最低氣溫≤0℃作為霜凍標(biāo)準(zhǔn)[6]，終霜日定義為前半年最后一次出現(xiàn)日最低氣溫≤0℃的日期，初霜日定義為后半年首次出現(xiàn)日最低氣溫≤0℃的日期。無(wú)霜期是一年中終霜后至初霜前的一整段時(shí)間[7]，在此期間內(nèi)沒有霜凍出現(xiàn)。

1.2 研究方法

1.2.1 霜凍初日、終日的時(shí)間序列 采用Julian日換算方法[8]，即將霜凍初日、終日出現(xiàn)日期轉(zhuǎn)換為距離1月1日的實(shí)際日數(shù)，得到霜凍初日、終日的時(shí)間序列。

1.2.2 線性傾向估計(jì)法 運(yùn)用線性傾向估計(jì)法[9]對(duì)霜凍初日、終日、無(wú)霜期的變化規(guī)律進(jìn)行分析。

1.2.3 80%保證率 保證率是指某一界限以上或以下要素各值出現(xiàn)概率的總和。利用如下方法計(jì)算80%保證率[10-11]。

式中，P80%是80%保證率;S為均方差;為序列的均值;C為變異系數(shù)，可由下式算出。

1.2.4 絕對(duì)變率 用以下式子計(jì)算初終霜日及無(wú)霜期的絕對(duì)變率。

應(yīng)用DPS7.05、Excel 2003軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)處理、分析及繪圖。

2 結(jié)果與分析

2.1 共和地區(qū)霜凍的氣候變化趨勢(shì)

2.1.1 霜凍終日變化趨勢(shì) 由圖1可見，1961—2017年近57年高寒地區(qū)共和地區(qū)霜凍終日日序呈現(xiàn)出波動(dòng)提前的趨勢(shì)，以每10年3.6 d的速率提前，霜凍終日日序與年份間的線性相關(guān)系數(shù)為0.478，通過99.99%的顯著性檢驗(yàn)。近57年≤0℃霜凍終日日序提前了20.6 d，提前趨勢(shì)極顯著。由表1可知，共和地區(qū)≤0℃霜凍終日歷年平均日序?yàn)?29.9 d（5月10日），最早日序?yàn)?06.0 d（1992年4月17日），最晚日序?yàn)?66.0 d（1969年6月16日），最早與最晚日序差值60 d?！?℃霜凍終日日序絕對(duì)變率為9.7，變異系數(shù)為0.095，80%保證率在5月3日。

由≤0℃霜凍終日年代平均值（表2）可以看出，20世紀(jì)60、70、80年代共和地區(qū)≤0℃霜凍終日落后于歷年平均日期，20世紀(jì)90年代開始≤0℃霜凍終日呈提前趨勢(shì)。

2.1.2 霜凍初日變化趨勢(shì) 從圖2可以看出，近57年來(lái)共和地區(qū)≤0℃霜凍初日日序呈現(xiàn)出波動(dòng)延遲的趨勢(shì)，其氣候傾向率為4.0 d/10年，≤0℃霜凍初日日序與年份間的線性相關(guān)系數(shù)為0.687，通過99.99%的顯著性檢驗(yàn)。近57年來(lái)≤0℃霜凍初日日序推遲了22.8 d，推遲趨勢(shì)極顯著。從表1可見，共和地區(qū)≤0℃霜凍歷年平均日序?yàn)?72.6 d（9月30日），初霜日發(fā)生日序最早為241.0 d（1962年8月30日），最晚為286.0 d（2000年10月14日），最早與最晚日序差值為45 d?！?℃霜凍初日日序的絕對(duì)變率為7.9，變異系數(shù)為0.036，80%保證率在9月19日。由初霜日年代平均值（表2）可以看出，20世紀(jì)60、70、80年代≤0℃霜凍初日早于歷年平均日期，20世紀(jì)90年代開始≤0℃終日則呈推后趨勢(shì)。

上述結(jié)果表明，共和地區(qū)初霜日的絕對(duì)變率要小于終霜日的絕對(duì)變率，表明初霜日較終霜日穩(wěn)定，年際差異相對(duì)較小。從最小值和最大值的差值來(lái)看，終霜日為60 d，初霜日為45 d，說明共和地區(qū)的終霜日年際差異較大，穩(wěn)定性差，對(duì)農(nóng)作物（馬鈴薯、碗豆等）和蔬菜的苗期以及果樹的花期造成不利影響。

2.1.3 無(wú)霜期的變化趨勢(shì) 農(nóng)作物的生長(zhǎng)期也與無(wú)霜期有著密切關(guān)系，無(wú)霜期長(zhǎng)，生長(zhǎng)期也長(zhǎng)。由圖3可見，近57年來(lái)無(wú)霜期日數(shù)呈現(xiàn)出波動(dòng)延長(zhǎng)的趨勢(shì)，其氣候傾向率為7.6 d/10年，無(wú)霜期日數(shù)與年份間的線性相關(guān)系數(shù)為0.719，通過99.99%顯著性水平檢驗(yàn)，即近57年來(lái)無(wú)霜期顯著延長(zhǎng)了43.4 d。由表1可知，共和地區(qū)無(wú)霜期日數(shù)的歷年平均值為141.6 d，絕對(duì)變率為14.1，變異系數(shù)為0.129，80%保證率在121 d。由無(wú)霜期年代平均值（表2）可以看出，20世紀(jì)60—80年代無(wú)霜期日數(shù)短于歷年平均日數(shù)，其中20世紀(jì)80年代最短，只有113.5 d;20世紀(jì)90年代以來(lái)無(wú)霜期日數(shù)長(zhǎng)于歷年平均日數(shù)，其中2011—2017年最長(zhǎng)，為153.7 d?？梢?，隨著年份的增加，共和地區(qū)霜凍的終日提前、初日延遲、無(wú)霜期延長(zhǎng)，這與許多學(xué)者[12-14]的研究結(jié)論相一致，說明隨著全球氣候變暖，不同地區(qū)霜凍的年際變化規(guī)律具有相似性。

2.2 近10年霜凍氣候變化特征

2.2.1 近10年霜凍終日變化趨勢(shì) 從圖4可以看出，近10年來(lái)共和地區(qū)≤0℃霜凍終日呈現(xiàn)出波動(dòng)推遲的趨勢(shì)，以每年1.2 d的速率推遲，≤0℃霜凍終日與年份間的線性相關(guān)系數(shù)為0.429，未通過0.10顯著性水平檢驗(yàn)。2008—2011年≤0℃霜凍終日呈逐年提前趨勢(shì)，2012—2017年≤0℃霜凍終日在波動(dòng)中呈現(xiàn)推遲趨勢(shì)，其中2014年≤0℃霜凍終日出現(xiàn)在5月19日，2015年出現(xiàn)在4月26日，兩者相差23 d，年際間變化幅度大。

2.2.2 近10年霜凍初日變化趨勢(shì) 從圖5可以看出，近10年來(lái)共和地區(qū)≤0℃霜凍初日呈現(xiàn)出波動(dòng)推遲的趨勢(shì)，以每年0.62 d的速率推遲，≤0℃霜凍初日與年份間的線性相關(guān)系數(shù)為0.708，通過0.05顯著性水平檢驗(yàn)。2008—2010年≤0℃霜凍初日呈逐年提前趨勢(shì)，2011—2017年≤0℃霜凍初日在波動(dòng)中呈現(xiàn)推遲趨勢(shì)，2017年出現(xiàn)在10月10日，推遲趨勢(shì)顯著。

2.2.3 近10年無(wú)霜期變化趨勢(shì) 由圖6可見，近10年來(lái)共和地區(qū)≤0℃無(wú)霜期呈現(xiàn)出波動(dòng)中縮短的趨勢(shì)，以每年0.56 d的速率縮短，10年間共縮短了5.6 d?！?℃無(wú)霜期與年份間的線性相關(guān)系數(shù)為0.199，未通過0.10顯著性水平檢驗(yàn)，說明近10年無(wú)霜期縮短趨勢(shì)不顯著。2008—2012年≤0℃無(wú)霜期呈逐年延長(zhǎng)趨勢(shì)，2012—2013年≤0℃無(wú)霜期呈急劇縮短趨勢(shì)，2015—2017年無(wú)霜期在波動(dòng)中呈現(xiàn)延長(zhǎng)趨勢(shì)。

2.3 初、終霜日的概率分布

從表3可以看出，終霜日較氣候平均值偏早的概率中，以提早1～5 d的概率最大，為52.0%，平均每2年一遇;其次是提早6～10 d和11～15 d的概率均為20.0%，平均每5年一遇;提早16～20 d和21～25 d的概率均為4.0%，平均每25年一遇。

初霜日較氣候平均值偏早的概率中，以提早1～5 d和6～10 d的概率最大，均為29.0%，平均每3年一遇;其次是提早11～15 d的概率，為22.6%，平均每4年一遇;提早16～20 d的概率為12.9%，平均每7年一遇;提早21～25 d的概率為6.5%，平均每15年一遇。

終霜日較氣候平均值偏晚的概率中，以偏晚1～5 d的概率最大，為37.0%，平均每2.7年一遇;其次是偏晚6～10 d的概率，為29.6%，平均每3.4年一遇;偏晚11～15 d的概率為11.1%，平均每9年一遇;偏晚21～25 d的概率為14.8%，平均每6.8年一遇;偏晚26～30 d的概率為7.4%，平均每14年一遇。

初霜日較氣候平均值偏晚的概率中，以偏晚1～5 d的概率最大，為54.2%，平均每2年一遇;其次是偏晚6～10 d的概率，為41.7%，平均每2.4年一遇;偏晚11～15 d的概率為4.2%，平均每24年一遇。

2.4 無(wú)霜期保證率

從氣候資源利用的角度評(píng)價(jià)，無(wú)霜期的變化范圍越大，表明熱量資源的穩(wěn)定性越差，可利用的程度越低;從防御災(zāi)害的角度分析，無(wú)霜期變化越大，農(nóng)業(yè)生產(chǎn)上躲避霜凍的難度越大，遭受的危害性越高。

用分組法計(jì)算保證率，將57年無(wú)霜期分為11組，統(tǒng)計(jì)各組出現(xiàn)的次數(shù)，計(jì)算各組出現(xiàn)的頻率，然后各組的頻率依次累加，其累計(jì)頻率就是無(wú)霜期保證率。

從圖7可以看出，共和地區(qū)平均無(wú)霜期為141.6 d，近57年無(wú)霜期長(zhǎng)于平均值的年份占40%，而短于平均值的年份占60%，長(zhǎng)于平均值的年份比短于平均值的年份少20%。無(wú)霜期的年際變化范圍較寬，最長(zhǎng)和最短的差值（極差）為75 d，接近平均無(wú)霜期的一半。在這個(gè)范圍內(nèi)，集中在平均值前后10 d的有25年，占43.9%;集中在平均值前后15 d的有35年，占61.4%;集中于平均值前后20 d的有39年，占68.4%。說明歷年無(wú)霜期的值不是很分散。保證率為95%、90%、80%的無(wú)霜期，分別比平均無(wú)霜期短33、19、11 d。

3 結(jié)論

（1）共和地區(qū)近57年來(lái)≤0℃霜凍的終日歷年平均日期為5月10日，初日歷年平均日期為9月30日?！?℃霜凍的無(wú)霜期日數(shù)的歷年平均值為141.6 d。

（2）近57年來(lái)共和地區(qū)≤0℃霜凍的終日以每10年3.6 d的速率呈極顯著的提前趨勢(shì)，初日以每10年4.0 d的速率呈顯著的后延趨勢(shì)，終日提前的幅度小于初日后延的幅度;無(wú)霜期以每10年7.6 d的速率呈現(xiàn)極顯著延長(zhǎng)趨勢(shì)。初、終霜日的絕對(duì)變率小于無(wú)霜期的絕對(duì)變率。霜凍終日80%保證率在5月3日，霜凍初日80%保證率在9月19日，無(wú)霜期80%保證率在121 d。近10年霜凍終日呈不顯著推遲趨勢(shì)，霜凍初日則呈顯著推遲趨勢(shì)，由于受到霜凍終日和初日均呈現(xiàn)推遲的影響，無(wú)霜期呈不顯著縮短趨勢(shì)。

（3）≤0℃霜凍的終霜日較氣候平均值偏早的概率中，以提早1～5 d的概率最大，為52.0%，平均每2年一遇;初霜日較氣候平均值偏早的概率中，以提早1～5 d和6～10 d的概率最大，均為29.0%，平均每3年一遇;終霜日較氣候平均值偏晚的概率中，以偏晚1～5 d的概率最大，為37.0%，平均每2.7年一遇;其次是偏晚6～10 d的概率為29.6%，平均每3.4年一遇;初霜日較氣候平均值偏晚的概率中，以偏晚1～5 d的概率最大，為54.2%，平均每2年一遇;其次是偏晚6～10 d的概率為41.7%，平均每2.4年一遇。

參考文獻(xiàn)

[1] 項(xiàng)陽(yáng)，張慶奎，尉傳陽(yáng).阜陽(yáng)市霜凍特征[J].中國(guó)農(nóng)學(xué)通報(bào)，2016，32（11）：156-160.

[2] 丁文魁，楊曉玲，殷玉春，等.武威市初、終霜日氣候特征[J].氣象科技，2016，44（4）：647-652.

[3] 陳芳，汪青春，殷萬(wàn)秀.青海省近45年霜凍變化特征及其對(duì)主要作物的影響[J].氣象科技，2009，37（1）：35-41.

[4] 嚴(yán)應(yīng)存，周秉榮，陳國(guó)茜.1961—2010年青海省霜凍災(zāi)害變化特征及風(fēng)險(xiǎn)區(qū)劃[J].氣象科技，2015，45（5）：986-991.

[5] 劉義花，周強(qiáng)，魯延榮，等.湟水河流域近50年來(lái)農(nóng)業(yè)氣候資源變化[J].中國(guó)農(nóng)學(xué)通報(bào)，2016，32（12）：163-170.

[6] 趙恒和，郭連云，趙年武，等.共和盆地界限溫度初終日和積溫對(duì)氣候變化趨勢(shì)的響應(yīng)[J].水土保持研究，2012，19（4）：207-211.

[7] 張玉書，班顯秀，紀(jì)瑞鵬，等.遼寧省氣候資源分析[J].氣象科技，2004，32（1）：39-43.

[8] 竺可楨，宛敏渭.物候?qū)W（增訂版）[M].北京：科學(xué)出版社，1980.

[9] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京：氣象出版社，2007：37-59.

[10] 吳洪寶，吳蕾.氣候變率診斷和預(yù)測(cè)方法[M].北京：氣象出版社，2005：208-244.

[11] 符淙斌，王強(qiáng).氣候突變的定義和檢測(cè)方法[J].大氣科學(xué)，1992，16（4）：482-493.

[12] 馬國(guó)柱.中國(guó)北方地區(qū)霜日的變化與區(qū)域增暖相互關(guān)系[J].地理學(xué)報(bào)，2003，589（S1）：31-37.

[13] 張山清，普宗朝，李景林，等.氣候變暖背景下新疆無(wú)霜凍期時(shí)空變化分析[J].資源科學(xué)，2013，35（9）：1908-1916.

[14] 張霞，錢錦霞.氣候變暖背景下太原市霜凍發(fā)生特征及其對(duì)農(nóng)業(yè)的影響[J].中國(guó)農(nóng)業(yè)氣象，2010，31（1）：111-114.

責(zé)任編輯：鄭丹丹