張家口壩上高寒區(qū)近53年霜凍變化特征及其對農(nóng)業(yè)的影響

吳偉光 郝瑛 劉星燕

摘要 利用1965—2017年張家口壩上高寒區(qū)康保、沽源、尚義和張北4個(gè)氣象站逐日最低氣溫資料，運(yùn)用一元線性回歸、Mann Kendall檢驗(yàn)法對該地區(qū)初（終）霜凍日和無霜期變化特征進(jìn)行了分析。結(jié)果表明，張家口壩上高寒區(qū)有平均初霜凍日推遲、平均終霜凍日提早、無霜期延長的趨勢，但是還存在一定的不穩(wěn)定性，平均終霜日1987—2012年幾乎以等振幅波動性變化；平均終霜凍日進(jìn)入21世紀(jì)后出現(xiàn)突變現(xiàn)象，突變幾率明顯增大。

關(guān)鍵詞 張家口壩上高寒區(qū)；初（終）霜凍日；變化特征；M K檢驗(yàn)

中圖分類號：S425 文獻(xiàn)標(biāo)識碼：A 文章編號：2095-3305（2018）04-045-03

DOI： 10.19383/j.cnki.nyzhyj.2018.04.018

Abstract Based on daily minimum temperature data of 4 meteorological stations including Kangbao， Guyuan， Shangyi and Zhangbei from 1965 to 2017 in high cold region of Bashang in Zhangjiakou， the climate change characteristics of the first and last frost date and frost free period were analyzed by using one dimensional linear regression method and Mann Kendall test method. The results showed that the average first frost date had a delaying trend and the situation was quite the reverse for the average last frost in Bashang. Moreover， the frost free period was extended. The first frost date had significantly delayed， but there was still some instability， the average last frost date fluctuated almost by equal amplitude from 1987 to 2012. The average last frost date had a sudden change， and the probability of mutation increased significantly in the 21st century.

Key words High cold region of Bashang in Zhangjiakou； First and last frost date； Change characteristics； Mann Kendall test

霜凍發(fā)生在春秋季節(jié)，是一種較為常見的農(nóng)業(yè)氣象災(zāi)害。農(nóng)業(yè)氣象學(xué)上，霜凍是指日平均氣溫高于0℃的溫暖時(shí)期，地面和作物體表面的溫度迅速降低到0℃以下，使作物遭受凍害或死亡的現(xiàn)象，是由于低溫造成的作物生理現(xiàn)象。在全球增暖[1-3]的氣候背景下，國內(nèi)外有不少學(xué)者對各地的霜凍發(fā)生規(guī)律進(jìn)行了研究。研究表明，美國、加拿大和北歐的霜凍日數(shù)都有減少趨勢[4-6]，國內(nèi)有很多學(xué)者也得出了類似的結(jié)論[7-9]。陳乾金等[10-11]發(fā)現(xiàn)，霜凍的異常變化特征還存在明顯的地區(qū)差異。

張家口壩上高寒區(qū)是以農(nóng)業(yè)為主的貧困地區(qū)[12]，地處40°N以北，海拔在1 200～1 600 m，是霜凍的多發(fā)、頻發(fā)嚴(yán)重區(qū)。作為該地一種重要的農(nóng)業(yè)氣象災(zāi)害，霜凍是影響農(nóng)作物單產(chǎn)水平及總產(chǎn)量波動的重要因素之一[13]，但目前對張家口壩上高寒區(qū)霜凍的研究文獻(xiàn)較少，且大都是局部地區(qū)或針對個(gè)別農(nóng)作物的霜凍防御研究[14-15]。筆者旨在分析氣候變暖背景下，壩上高寒區(qū)4個(gè)代表站霜凍的變化特征和突變特征，以期為霜凍的預(yù)報(bào)服務(wù)、選取合適生育期品種等工作提供理論參考。

1 資料與方法

1.1 資料來源

考慮數(shù)據(jù)的時(shí)間序列長度、完整性和高寒區(qū)的代表性，選取壩上康保、沽源、尚義、張北4個(gè)高寒地區(qū)觀測站的逐日0 cm氣溫、逐年初霜凍日和逐年終霜凍日資料。時(shí)間序列為1965—2017年，初霜凍日為9月1日起0 cm地溫≤0℃的第1天；終霜凍日為次年3月份起0 cm地溫≤0℃的最后1天。資料來源為張家口市氣象局資料室?？当！⒐猎?、尚義、張北4個(gè)縣站1965—2017年氣溫觀測資料，其中張北為國家基準(zhǔn)站，其余3站為國家基本站。

1.2 研究方法

采用一元線性回歸方程擬合計(jì)算初（終）霜凍日的傾向率。方程為y（t）=a0+a1t，y為初、終霜凍日擬和值，t是年序列，t=1，2，3，…，53（年），其線性趨勢傾向系數(shù)由最小二乘法求得。dy（t）/dt=a1，a1×10年為（終）霜凍日的傾向率，單位為d/10年。

用Mann Kendall（M K）方法檢驗(yàn)初（終）霜凍日和無霜期的突變點(diǎn)。

2 結(jié)果與分析

2.1 初、終霜凍日的氣候統(tǒng)計(jì)特征及年際變化特征

2.1.1 初、終霜凍日和無霜期的氣候統(tǒng)計(jì)特征 由表1可以看出，近53年張家口市壩上高寒區(qū)平均初霜凍日在9月上、中旬，平均終霜凍日在5月下旬或6月上旬。最早初霜凍日為8月15日（康保），最晚初霜凍日為10月3日（張北）；最早終霜凍日為4月28日（康保和張北），最晚終霜凍日為6月26日（康保和尚義）。各站最早、最晚初霜凍日的差值沽源最小，為28 d，其他3站最早、最晚終霜凍日的差值均在40 d以上。從線性趨勢可以看出，康保、沽源、尚義和張北均有初霜凍推遲、終霜凍結(jié)束早、無霜期延長的趨勢，其中尚義更為顯著，初霜凍日平均每10年推遲2.7 d，終霜凍日平均每10年提早4.4 d。

2.1.2 初、終霜凍日的年際變化特征 圖1為張家口市壩上4縣平均初、終霜凍日年際變化及一元線性回歸和四階多項(xiàng)式擬合曲線。由圖1a可知，張家口市壩上4縣的平均初霜凍日在振動中逐漸推遲，平均每10年推遲0.9 d。1965—1983年、2000—2014年初霜凍出現(xiàn)的時(shí)間波動較大，1972年康保、沽源、張北初霜凍均出現(xiàn)在8月下旬，2012年張北和尚義8月22日出現(xiàn)初霜凍，分別創(chuàng)造和接近歷史極值；1984—1999年初霜凍出現(xiàn)的時(shí)間相對比較平穩(wěn)。

由圖1b可知，近53年終霜凍日在波動變化中逐漸提前，平均每10年提前0.7 d。1977—1987年終霜凍日波動較大，雖然平均終霜凍日在振動中提前，但是終霜凍日出現(xiàn)的偏晚概率增大了。

2.2 初、終霜凍日的突變性分析

從初霜凍日UF曲線（圖2a）可以看出，初霜凍日只是階段性超過了0.01顯著性水平臨界線，未出現(xiàn)明顯超過0.01顯著性水平臨界線的點(diǎn)，說明初霜凍日推遲并不是太明顯。1973—1982年UF和UB 2條曲線出現(xiàn)多個(gè)交點(diǎn)，判斷初霜凍日在此階段可能發(fā)生突變；終霜凍日UF和UB 2條曲線于1973年出現(xiàn)交點(diǎn)，可以判定初霜凍日突變可能發(fā)生在1973年前后。UF曲線階段性超出了0.01顯著性水平臨界線，未出現(xiàn)明顯超過0.01顯著性水平臨界線的點(diǎn)，說明終霜凍日提前不是很顯著（圖2b）。

2.3 霜凍對農(nóng)業(yè)生產(chǎn)的影響

霜凍對農(nóng)作物的損傷，是由于突然降溫使農(nóng)作物體內(nèi)的水分結(jié)冰，特別是植物細(xì)胞內(nèi)體積增大導(dǎo)致細(xì)胞機(jī)能受到破壞而形成的[16]。張家口壩上地區(qū)初霜凍早、終霜凍晚、無霜期短，主要農(nóng)作物有小麥、莜麥、馬鈴薯、青玉米、甜菜、蔬菜和豆類及油料類雜糧。張家口壩上地區(qū)的平均終霜凍日為5月中旬，對馬鈴薯和錯(cuò)季蔬菜影響較大，其他多數(shù)作物播種期一般在5月中下旬，影響不是很明顯，極端天氣除外。

初霜凍對作物的影響和危害非常大，雖然初霜凍日有推遲的趨勢，但自20世紀(jì)90年代起，突發(fā)提前的頻率在逐漸加大，極端天氣出現(xiàn)的概率越來越高。如2012年8月22—24日，張家口市壩上沽源、張北、康保、尚義、塞北、察北和壩下崇禮、赤城、懷安、宣化縣相繼出現(xiàn)大幅降溫天氣，致使以上10個(gè)縣區(qū)的63鄉(xiāng)鎮(zhèn)、829個(gè)行政村遭受了自1961年有氣象記載以來50年一遇的特大低溫冷凍災(zāi)害，即將收獲的青玉米、馬鈴薯、架豆、雜豆、谷黍等農(nóng)作物不同程度減產(chǎn)，甚至大面積絕收，受災(zāi)人口112 520戶、315 103人，受災(zāi)面積96 650.68 hm2，造成農(nóng)業(yè)損失44 647.94萬元。

據(jù)資料統(tǒng)計(jì)，近53年來張家口壩上地區(qū)的初霜凍日70%的年份出現(xiàn)在8月下旬至9月上中旬，這期間正是糧食作物的成熟期，在此期間初霜凍出現(xiàn)的越早對作物的危害越嚴(yán)重。

3 結(jié)論與討論

（1）張家口壩上4縣的平均初霜凍日在9月上、中旬，平均終霜凍日在5月下旬或6月上旬，并且有平均初霜凍日推遲、平均終霜凍日提早、無霜期延長的趨勢，這與國內(nèi)外學(xué)者分析的大部分地區(qū)霜凍日數(shù)有減少的趨勢一致。

（2）張家口壩上4縣平均初霜凍日1965—1983年、2000—2014年初霜凍出現(xiàn)的時(shí)間波動較大，1984—1999年相對比較平穩(wěn)。20世紀(jì)80—90年代，初霜凍日有顯著推遲的趨勢，但是還存在一定的不穩(wěn)定性。平均終霜日1974—1985年出現(xiàn)的時(shí)間相對較早，在5月上、中旬；1987年以后幾乎以等振幅波動性變化。

（3）康保、沽源、張北和尚義4縣平均初霜凍日20世紀(jì)70—80年代突變相對較多，平均終霜凍日進(jìn)入21世紀(jì)后突變幾率明顯增大。

（4）初霜凍日推遲、終霜凍日提早、無霜期延長可以為農(nóng)作物提供更多的熱量資源和生長時(shí)間，對農(nóng)作物的生長和發(fā)育很有利，但是霜凍是冷空氣突然入侵強(qiáng)降溫造成的，是一種時(shí)間尺度很短的農(nóng)業(yè)氣象災(zāi)害，它的發(fā)生可以與長時(shí)間尺度的氣候變化趨勢相反。在氣候變暖的大背景下，出現(xiàn)初霜凍日突發(fā)提前的頻率在逐漸增大。霜凍作為一種嚴(yán)重的氣象災(zāi)害，一旦出現(xiàn)，對農(nóng)作物將造成無法挽救的損傷，直接影響到農(nóng)民的收成。研究霜凍出現(xiàn)的本地指標(biāo)以及環(huán)流形式是未來比較緊要的工作。

參考文獻(xiàn)

[1] IPCC. Climate change 2007： The physical science basis. Contribution of working group 1 to the fourth assessment report of the intergovernmental panel on climate change[R]. Cambridge， United Kingdom and New York， NY USA： Cambridge University Press，2007：996.

[2] 任國玉，初子瑩，周雅清.中國氣溫變化研究最新進(jìn)展[J].氣候與環(huán)境研究，2005，10（4）：701-706.

[3] 王紹武，蔡靜寧，朱錦紅，等.中國氣候變化的研究[J].氣候與環(huán)境研究，2002，7（2）：137-145.

[4] EASTERLING D R.Recent changes in the frost days and the frost free season in the United States[J]. Bulletin of the American Mathematical Society，2002，83（9）：1327-1332.

[5] BONSAL B R，ZHANG X，VINCENT L A，et al. Characteristics of daily and extreme temperature over Canada[J]. Journal of Climate，2001，14（5）：1959-1976.

[6] HEINO R， BRáZDIL E. F RLAND H，et al. Progress in the study of climate extremes in Northern and central Europe[J].Climate Change，1999，42（1）：15l-181.

[7] 葉殿秀，張勇.1961—2007年我國霜凍變化特征[J].應(yīng)用氣象學(xué)報(bào)，2008，19（6）：661-664.

[8] 李鳳云，王玉山，吳澤新，等.魯西北霜凍日期變化特征及防御措施[J].蘭州大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，46（S1）：126-128.

[9] 杜春英，李帥，郭建平， 等.東北三省玉米初霜凍時(shí)空分布特征及預(yù)報(bào)方法探討[J].安徽農(nóng)業(yè)科學(xué)，2010，38（1）：237-240.

[10] 陳乾金，張水山.華北異常初（終）霜凍氣候特征的研究[J].自然災(zāi)害學(xué)報(bào)，1995，4（3）：33-39.

[11] 陳乾金，夏洪星，張水山.我國江淮流域近40年異常初終霜凍的分析[J].應(yīng)用氣象學(xué)報(bào)，1995，6（1）：50-56.

[12] 劉愛梅，李景宇，苗志成.蔬菜作物需水與節(jié)水灌溉氣象服務(wù)[J].氣象科技，2010，38（3）：373-376.

[13] 劉明亮，陳百明.我國近期糧食生產(chǎn)的波動及其與農(nóng)業(yè)自然災(zāi)害發(fā)生狀況的相關(guān)分析[J].災(zāi)害學(xué)，2000（12）：78-85.

[14] 黃山江，劉星燕，孫躍飛，等.一次強(qiáng)降溫天氣對杏扁的凍害分析及防御思考[J].湖北農(nóng)業(yè)科學(xué)，2012，51（20）：4503-4506.

[15] 劉愛梅，苗志成，劉星燕，等.張家口氣候資源與特色農(nóng)業(yè)[J].河北氣象，2004，23（4）：26-27.

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

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