張志高， 劉晴， 劉慧穎， 郭超凡， 蔡茂堂， 尹紀(jì)媛

1. 安陽(yáng)師范學(xué)院 資源環(huán)境與旅游學(xué)院,河南 安陽(yáng) 455000； 2. 中國(guó)地質(zhì)科學(xué)院 地質(zhì)力學(xué)研究所,北京100081

霜凍常出現(xiàn)在春、 秋兩季農(nóng)作物生長(zhǎng)期,是由寒潮和強(qiáng)冷空氣活動(dòng)引致土壤或植物表面溫度降至0 ℃以下而造成作物受害的現(xiàn)象,是中高緯度地區(qū)影響農(nóng)業(yè)生產(chǎn)的主要災(zāi)害之一[1-2].目前,全球正在經(jīng)歷以變暖為主要特征的氣候變化,全球地表溫度正以前所未有的速度上升,2011-2020年全球平均表面溫度相比1850-1900年上升1.09 ℃[3].全球變暖將改變氣溫、 輻射以及降水等要素的地理分布,在一定程度上對(duì)初、 終霜日和霜期長(zhǎng)短產(chǎn)生影響,進(jìn)而影響農(nóng)業(yè)生產(chǎn)[4-5].因此,近年來(lái)在全球氣候變暖背景下,霜凍的時(shí)空演變特征及其影響因素與災(zāi)害評(píng)估受到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注[6-9].Frich等[10]研究發(fā)現(xiàn)20世紀(jì)中后期全球年霜凍日數(shù)不斷減少.Crimp等[11]研究發(fā)現(xiàn)澳大利亞南部地區(qū)霜凍時(shí)間偏晚,東南部地區(qū)的霜凍頻率呈增加趨勢(shì).寧曉菊等[12]分析了1951年以來(lái)中國(guó)無(wú)霜期的變化趨勢(shì),發(fā)現(xiàn)中國(guó)大部分區(qū)域呈初霜日推后、 終霜日提前無(wú)霜期延長(zhǎng)的趨勢(shì),在變化幅度方面北方大于南方、 東部大于西部.張?chǎng)蔚萚13]研究發(fā)現(xiàn)北疆初霜日推遲,終霜日提前,霜期縮短,速率分別為2.2,-1.7,-3.9 d/10a,初、 終霜日和霜期的主周期均為2a左右.任景全等[14]發(fā)現(xiàn)吉林省初霜日顯著推遲,終霜日顯著提前且無(wú)霜期顯著延長(zhǎng),其主周期均為3a左右.馬尚謙等[15]發(fā)現(xiàn)淮河流域呈初、 終霜日和無(wú)霜期分別以2.15,-2.49,4.38 d/10a的速率呈推遲、 提前和延長(zhǎng)的趨勢(shì),且在在20世紀(jì)90年代的變化最為顯著.在霜凍變化影響因素方面,王國(guó)復(fù)等[16]認(rèn)為中國(guó)大部地區(qū)霜期縮短,初霜日推遲,終霜日提前主要與日最低氣溫和0 cm最低地溫的不斷升高有關(guān)； 馬彬等[17]認(rèn)為極渦與副熱帶高壓的位置、 強(qiáng)度、 面積變化驅(qū)動(dòng)了中國(guó)農(nóng)業(yè)區(qū)霜凍多年氣候狀態(tài)的轉(zhuǎn)變； 馬尚謙等[18]認(rèn)為淮河流域氣候的轉(zhuǎn)型促進(jìn)了初霜日的推遲和終霜日的提前； 陳少勇等[19]探討了西北地區(qū)影響初霜凍發(fā)生的環(huán)流特征,認(rèn)為與副熱帶高壓正相關(guān),與極渦負(fù)相關(guān),初霜日推遲是氣候變暖的結(jié)果； Liu等[20]研究發(fā)現(xiàn)氣候變暖減少了霜凍事件,但植物生長(zhǎng)季節(jié)的延長(zhǎng)可能在生長(zhǎng)季誘發(fā)更高的霜凍頻率.此外還有學(xué)者分析了霜凍災(zāi)害對(duì)農(nóng)作物的影響,代立芹等[21]認(rèn)為河北省冬小麥霜凍害以輕度為主,霜凍導(dǎo)致冬小麥氣象產(chǎn)量下降、 減產(chǎn)率增加； 高文波等[22]構(gòu)建了西南茶區(qū)茶樹(shù)春梢晚霜凍害概率和強(qiáng)度指數(shù),發(fā)現(xiàn)西南茶區(qū)霜凍終日和初日均呈顯著提前趨勢(shì),茶樹(shù)晚霜凍害危險(xiǎn)性呈下降趨勢(shì)； 呂佳佳等[23]采用數(shù)理統(tǒng)計(jì)和墨西哥帽小波法分析了寒地大豆霜凍的時(shí)空變化特征,發(fā)現(xiàn)霜凍對(duì)大豆產(chǎn)量的影響為霜凍等級(jí)及出現(xiàn)日期的協(xié)同作用．

河南省位于黃河中下游,平原面積廣闊,屬北亞熱帶向暖溫帶過(guò)渡的大陸性季風(fēng)氣候,氣候分明,雨熱同期,適宜多種農(nóng)作物生長(zhǎng),被譽(yù)為“中原糧倉(cāng)”,在我國(guó)農(nóng)業(yè)生產(chǎn)中占有舉足輕重的地位.然而由于河南省地處南北氣候過(guò)渡帶和山區(qū)到平原的過(guò)渡帶,對(duì)氣候變化的適應(yīng)能力較弱,加之農(nóng)業(yè)基礎(chǔ)條件薄弱,受霜凍等氣象災(zāi)害影響較為嚴(yán)重.因此,本研究基于1960-2019年逐日氣象觀測(cè)數(shù)據(jù),運(yùn)用線性傾向估計(jì)、 Mann-Kendall突變檢驗(yàn)、 相關(guān)分析及R/S分析等方法對(duì)河南省初霜、 終霜及霜期的時(shí)空演變進(jìn)行研究,并探討其與地理因子的關(guān)系,以期為河南省霜凍災(zāi)害的預(yù)報(bào)與防御及農(nóng)業(yè)可持續(xù)發(fā)展提供科學(xué)依據(jù)和參考．

1 資料與方法

1.1 研究區(qū)概況

河南省位于中國(guó)中東部、 黃河中下游,地理位置為31°23′-36°22′N,110°21′-116°39′E.全省地勢(shì)西高東低,太行山脈、 伏牛山脈、 桐柏山脈、 大別山脈沿省界呈半環(huán)形分布,中、 東部為華北平原,西南部為南陽(yáng)盆地.全省總面積16.7×104km2,其中平原和盆地占55.7%,山地和丘陵占44.3%,農(nóng)業(yè)用地占74.25%.河南省氣候自南向北由亞熱帶向暖溫帶氣候逐漸過(guò)渡,氣候溫和,四季分明,年平均氣溫在12～16 ℃,年降水量約為500～1 400 mm,無(wú)霜期為190～230 d,日照時(shí)數(shù)1 740～2 310 h.河南省地處南北氣候過(guò)渡帶,同時(shí)自東向西由平原向丘陵山地氣候過(guò)渡,霜凍等氣象災(zāi)害較為頻繁．

審圖號(hào)： 豫S(2019年)013號(hào)．圖1 河南省氣象代表站點(diǎn)分布

1.2 數(shù)據(jù)來(lái)源

河南省17個(gè)氣象站點(diǎn)1960-2019年逐日地面0 cm最低溫度數(shù)據(jù)來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http： //data.cma.cn)中國(guó)地面氣候資料日值數(shù)據(jù)集.以地面0cm最低溫度≤0 ℃作為霜凍指標(biāo),參考前人研究[24],中國(guó)北方初霜日一般在8月1日后,因此以8月1日作為統(tǒng)計(jì)起始點(diǎn),將8月1日定義為1,用日序表示初、 終霜日,初霜日為8月1日以后第一次出現(xiàn)地面最低溫度≤0 ℃的日期,終霜日為最后一次地面最低溫度≤0 ℃的日期,霜凍終日的翌日至初日的前一天之間的日數(shù)為無(wú)霜期.利用算術(shù)平均法建立1960-2019年河南省初霜期、 終霜期和無(wú)霜期序列.研究區(qū)域及代表站點(diǎn)分布如圖1所示．

1.3 R/S分析

重標(biāo)極差分析法(R/S)最初由英國(guó)科學(xué)家赫斯特提出,用于分析尼羅河水庫(kù)流量和儲(chǔ)蓄能力[25],后被廣泛用于氣象學(xué)領(lǐng)域的趨勢(shì)預(yù)測(cè).R/S分析法用氣候要素時(shí)間序列計(jì)算出Hurst指數(shù)(H值)來(lái)揭示氣候要素的趨勢(shì),如果0

此外,采用線性傾向估計(jì)法分析河南省霜凍的時(shí)間變化特征； 利用相關(guān)分析和線性回歸法分析霜凍日期與地理因子的關(guān)系,變量間相關(guān)性強(qiáng)度依據(jù)相關(guān)系數(shù)絕對(duì)值R大小分為5個(gè)等級(jí),分別為極強(qiáng)相關(guān)(0.8

2 結(jié)果與分析

2.1 河南省霜凍日期時(shí)間變化特征

近60a來(lái)河南省初霜平均日期為11月8日,其中1981年初霜日最早,為10月24日.2011年初霜日最晚,為11月25日(圖2a).進(jìn)行線性擬合可知,河南省初霜日期顯著推遲,傾向率為2.66 d/10a.1960-2019年河南省終霜平均日期為3月28日,其中2007年終霜日最早,為 3月12日.1961年終霜日最晚,為4月14日(圖2b),60a來(lái)終霜日以-2.88 d/10a的傾向率顯著提前.1960-2019年河南省無(wú)霜期平均為224 d,其中1962年最短,為202 d,2016年最長(zhǎng),為249 d(圖2c),60a來(lái)無(wú)霜期顯著延長(zhǎng),傾向率為5.54 d/10a．

年代際變化方面,20世紀(jì)70年代初霜日最早,平均為11月3日,之后逐漸推遲,2000年左右距平轉(zhuǎn)正,21世紀(jì)10年代初霜日最晚,平均為11月18日.20世紀(jì)60年代終霜日最晚,為4月3日,之后逐年代提前,2000年以后提前至3月21日.20世紀(jì)70年代無(wú)霜期最短為214 d,之后逐漸延長(zhǎng),2000年左右距平轉(zhuǎn)正,21世紀(jì)10年代最長(zhǎng),為241 d,相比20世紀(jì)70年代延長(zhǎng)27 d(表1)．

圖2 1960-2019年河南省初霜日、 終霜日和無(wú)霜期年際變化

表1 河南省初霜日、 終霜日和無(wú)霜期年代際變化

2.2 河南省霜凍日期空間分布特征

2.2.1 霜凍的平均空間分布特征

河南省各地初霜日出現(xiàn)在10月末至11月中上旬,空間差異明顯,大致與緯線平行,呈“北早南晚”的空間分布格局,其中河南省西南部地區(qū)初霜日在11月14日以后； 中部地區(qū)在11月3-10日,北部安陽(yáng)站最早,為10月30日(圖3a).近60a河南省各地平均終霜日自南向北逐漸推遲,南部地區(qū)終霜日多數(shù)出現(xiàn)在3月中下旬； 北部地區(qū)終霜日多數(shù)出現(xiàn)在3月下旬和4月上旬； 中部地區(qū)終霜日多數(shù)出現(xiàn)在3月26日-4月2日(圖3b).近60a河南省各地平均無(wú)霜期為206～244 d,大致與緯線平行,自南向北呈逐漸減少趨勢(shì),西南部地區(qū)無(wú)霜期在237 d以上； 豫北地區(qū)無(wú)霜期普遍小于213 d,中部地區(qū)無(wú)霜期多在213～229 d(圖3c)．

審圖號(hào)： 豫S(2019年)013號(hào)．圖3 1960-2019年河南省初霜日、 終霜日和無(wú)霜期的空間分布

2.2.2 霜凍日期變化趨勢(shì)空間分布特征

河南省17個(gè)站點(diǎn)初霜日期傾向率均為正值,均呈推遲趨勢(shì),全省初霜日氣候傾向率在0.18～3.60 d/10a之間,其中西南部地區(qū)推遲幅度較小,多在2.53 d/10a以下； 中部和北部地區(qū)推遲幅度較大,多在2.68 d/10a以上(圖4a).近60a全省17個(gè)站點(diǎn)終霜日傾向率均為負(fù)值,終霜日均呈提前趨勢(shì),全省終霜日氣候傾向率為-0.39～-0.14 d/10a,西南部提前幅度較大,多在-0.35 d/10a以上(圖4b).60a來(lái)河南省無(wú)霜期顯著延長(zhǎng),氣候傾向率為1.70～6.74 d/10a,東南部無(wú)霜期延長(zhǎng)幅度較小,在4.04 d/10a以下,西南部為5.69～6.24 d/10a,中部地區(qū)無(wú)霜期延長(zhǎng)幅度較大,在5.79 d/10a以上(圖4c)．

審圖號(hào)： 豫S(2019年)013號(hào)．圖4 1960-2019年河南省初霜日、 終霜日和無(wú)霜期趨勢(shì)系數(shù)的空間分布

2.3 霜凍日期突變特征

近60a來(lái)河南省初霜日正序列曲線UF整體呈上升趨勢(shì),反序列UB曲線呈下降趨勢(shì),二者相交于2004年,且在α=0.05的信度線內(nèi),確定初霜日的突變點(diǎn)為2004年,突變前平均初霜日期為11月5日,突變后為11月18日,推遲了13 d(圖5a).終霜日正序列曲線UF呈下降趨勢(shì),反序列曲線UB呈上升趨勢(shì),UF和UB曲線相交于1998年,且在α=0.05的信度線內(nèi),確定終霜日的突變點(diǎn)為1998年,突變后平均終霜日期為3月22日,較突變前的4月2日提前了11 d(圖5b).無(wú)霜期正序列曲線UF呈下降趨勢(shì),反序列UB曲線呈上升趨勢(shì),UF和UB曲線相交于1998年,確定終霜日的突變點(diǎn)為1998年,突變后無(wú)霜期平均為238 d,較突變前的217 d延長(zhǎng)了21 d(圖5c)．

2.4 霜凍日期的未來(lái)趨勢(shì)變化

根據(jù)R/S分析,河南省初霜日、 終霜日和無(wú)霜期的H指數(shù)平均值分別為0.88,0.91,0.96(表2),均明顯大于0.5,說(shuō)明序列的趨勢(shì)與原序列呈正相關(guān),未來(lái)河南省霜凍日期變化與過(guò)去60a的變化趨勢(shì)保持一致,即初霜日期持續(xù)推遲,終霜日期持續(xù)提前,無(wú)霜日持續(xù)延長(zhǎng)．

圖5 1960-2019年河南省初霜日、 終霜日和無(wú)霜期的M-K檢驗(yàn)

表2 河南省初霜日、 終霜日和無(wú)霜期的R/S分析

2.5 霜凍日期與地理因子的關(guān)系

河南省初霜日、 終霜日和無(wú)霜期與海拔和經(jīng)緯度的相關(guān)關(guān)系和回歸系數(shù)如表3所示.初霜日與緯度呈極顯著負(fù)相關(guān)關(guān)系,緯度每增加1°,初霜日提前4.72 d,初霜日與海拔和經(jīng)度呈極弱負(fù)相關(guān)關(guān)系,海拔每上升10 m,經(jīng)度每增加1°,初霜日分別提前0.01 d和0.17 d.終霜日與緯度呈極顯著正相關(guān)關(guān)系,緯度每增加1°,終霜日推遲5.06 d,終霜日與海拔和經(jīng)度呈極弱正相關(guān)關(guān)系,海拔每上升10 m,經(jīng)度每增加1°,終霜日分別推遲0.01 d和0.43 d.無(wú)霜期與緯度呈極顯著負(fù)相關(guān)關(guān)系,緯度每增加1°,無(wú)霜期縮短9.7 d.無(wú)霜期與海拔和經(jīng)度呈極弱負(fù)相關(guān)關(guān)系,海拔每上升10 m,經(jīng)度每增加1°,無(wú)霜期分別縮短0.05 d和4.58 d．

表3 河南省初霜日、 終霜日和無(wú)霜期與地理因子的關(guān)系

3 結(jié)論與討論

3.1 結(jié)論

1960-2019年河南省初、 終霜平均日期分別為11月8日和3月28日,無(wú)霜期平均為224 d,近60a來(lái)河南省初霜日以2.66 d/10a的傾向率顯著推遲,終霜日顯著提前,傾向率為-2.88 d/10a,無(wú)霜期顯著延長(zhǎng),傾向率為5.54 d/10a.初霜日突變點(diǎn)為2004年,突變后初霜日推遲13 d,終霜日和無(wú)霜期于1998年發(fā)生突變,突變后終霜日提前11 d,無(wú)霜期延長(zhǎng)21 d．

近60a河南省各地初霜日自南向北逐漸提前,終霜日自南向北逐漸推遲,多數(shù)地區(qū)終霜日出現(xiàn)在3月下旬至4月上旬,全省無(wú)霜期在206～244 d之間,自南向北逐漸減?。?/p>

河南省初霜日和無(wú)霜期均與緯度呈極顯著負(fù)相關(guān)關(guān)系,終霜日和緯度呈極顯著正相關(guān)關(guān)系,河南省初、 終霜日和無(wú)霜期主要受緯度影響．

3.2 討論

在全球變暖背景下,近60a河南省初霜日推遲,終霜日提前,無(wú)霜期顯著延長(zhǎng),這與中國(guó)大部分地區(qū)研究結(jié)果一致[12，27],這將延長(zhǎng)作物生長(zhǎng)季[28],利于作物適宜種植區(qū)域由較低海拔向較高海拔區(qū)域推進(jìn)擴(kuò)展,然而生長(zhǎng)季延長(zhǎng)又增加了作物在霜凍中的暴露,因此,在未來(lái)全球變暖背景下,各地還需因地制宜地采取調(diào)整農(nóng)作物種植制度等措施,以減少霜凍危害．

本研究探討了霜凍與地理因子的關(guān)系,認(rèn)為河南省初、 終霜日和無(wú)霜期主要受緯度影響,限于數(shù)據(jù)的可獲取性及篇幅所限,環(huán)流指數(shù)和城市化等因素對(duì)霜凍的影響有待進(jìn)一步研究．