羅靜 費(fèi)杜秋

摘要?選用1961年1月1日—2017年12月31日陜西省34個(gè)地面氣象觀測(cè)站的地面日最低溫度，基于不同基準(zhǔn)氣候期，分析陜西省初終霜日時(shí)空變化特征。結(jié)果表明，陜西省全省初霜日期在1999年發(fā)生突變，終霜日期未發(fā)生突變;近57年陜西省全省及各分區(qū)初霜日期均呈推后趨勢(shì)，終霜日期均呈提前趨勢(shì)。陜北平均初霜日期和終霜日期推后和提前趨勢(shì)最顯著，其中陜北的靖邊推后趨勢(shì)和提前趨勢(shì)分別達(dá)3.0和4.0 d/10 a;近57年陜西省各地初霜日期和終霜日期由北到南呈顯著的帶狀分布。初霜日期出現(xiàn)最早和終霜日期出現(xiàn)最遲的均是吳旗—靖邊—橫山—榆林—神木一線，初霜日期出現(xiàn)最遲和終霜日期出現(xiàn)最早的均是陜南中部地區(qū)、關(guān)中太白地區(qū)、漢中寧強(qiáng)地區(qū)。

關(guān)鍵詞?初終霜日;時(shí)空變化特征;陜西省

中圖分類(lèi)號(hào)?P467文獻(xiàn)標(biāo)識(shí)碼?A

文章編號(hào)?0517-6611（2019）20-0228-07

doi：10.3969/j.issn.0517-6611.2019.20.062

開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Analysis of Temporal and Spatial Changes of First and Last Frost Day in Shaanxi Province from 1961 to 2017

LUO Jing，F(xiàn)EI Du?qiu?（Qinghai Engineering Consulting Center，Xining，Qinghai 810001）

Abstract?The ground?day minimum temperature data of 34 surface meteorological observatories in Shaanxi Province from January 1，1961 to December 31，2017 were selected.Based on different reference climatic periods，the temporal and spatial variation characteristics of the first and last frost days in Shaanxi Province were analyzed.The results showed that the first frost date in Shaanxi Province was abruptly changed in 1999，and there was no mutation in the last frost date.The first frost date of the province and all districts in Shaanxi Province showed a post?pushing trend，and the last frost date both were in an early trend.The average first frost date and last frost date in northern Shaanxi were the most significant.Among them，the post?pushing trend and the advance trend of Jingbian in northern Shaanxi reached 3.0 and 4.0 d/10 a respectively.Nearly 57 years，the first frost date and the last frost date of Shaanxi Province showed a significant band distribution from north to south.The earliest date of the first frost and the latest date of the late frost were the Wuqi?Jingbian?HengshanYulin?Shenmu line.The latest date of the first frost date and the earliest date of last frost date were the central part of southern Shaanxi，Guanzhong Taibai area and Hanzhong Ningqiang area.

Key words?First and last frost date;Temporal and spatial changes;Shaanxi Province

近年來(lái)，中緯度地區(qū)氣候暖干化是全球氣候變化最顯著的特征，是備受矚目的全球熱點(diǎn)環(huán)境問(wèn)題之一[1-2]。IPCC第五次評(píng)估報(bào)告指出，全球地表持續(xù)升溫，1880—2012年全球平均溫度已升溫0.85 ℃，其中區(qū)域間升溫幅度在0.65～1.06 ℃[3]?，F(xiàn)在以溫度升高為主要特征的氣候變化以及并將繼續(xù)對(duì)自然生態(tài)系統(tǒng)和人類(lèi)社會(huì)經(jīng)濟(jì)產(chǎn)生重大影響，越來(lái)越引起社會(huì)和公眾的關(guān)注。

霜凍作為一種用溫度來(lái)表征的農(nóng)業(yè)氣象災(zāi)害，與溫度變化緊密聯(lián)系[4]。初霜常使晚熟的大田作物莖葉提前凋萎，籽粒成熟不飽滿(mǎn)或品質(zhì)變壞，甚至不能成熟。終霜常使大田作物不能正常生長(zhǎng)發(fā)育，對(duì)作物的苗期影響極大。作為衡量溫度變化的重要指標(biāo)，分析初、終霜日期時(shí)空變化趨勢(shì)有利于減少極端低溫對(duì)農(nóng)業(yè)生產(chǎn)帶來(lái)的危害，也有利于因地制宜地選擇作物品種、及時(shí)調(diào)整種植制度等。目前已有較多研究在初、終霜日出現(xiàn)早晚、無(wú)霜期長(zhǎng)短等方面取得重要進(jìn)展[5-9]。周曉宇等[5]利用東北地區(qū)1961—2014年日最低氣溫及平均氣溫資料，分析了其近54年初霜日、終霜日及無(wú)霜期的氣候變化特征。陳少勇等[7]利用中國(guó)西北地區(qū)135個(gè)測(cè)站1961—2009年逐日地面最低氣溫和NCEP/NCAR資料，采用線性趨勢(shì)分析、Mann-Kendall、子波分析、功率譜分析和合成分析等方法，分析了近49年中國(guó)西北地區(qū)初霜凍的氣候變化特征。馬尚謙等[9]基于淮河流域 1960—2015年61個(gè)氣象站點(diǎn)地面 0cm 日最低氣溫資料，采用線性?xún)A向估計(jì)、反距離加權(quán)、Mann-Kendall 突變檢驗(yàn)、滑動(dòng)T檢驗(yàn)方法，分析近 56年淮河流域初霜日、終霜日和無(wú)霜期的時(shí)空變化特征及突變年份。

目前，針對(duì)陜西省的霜凍研究主要集中在果區(qū)霜凍風(fēng)險(xiǎn)區(qū)劃方面[10-11]，將初、終霜日與無(wú)霜期的變化結(jié)合起來(lái)分析三者的變化趨勢(shì)及突變的研究尚未開(kāi)展。陜西省地處中國(guó)南北氣候過(guò)渡帶，屬暖溫帶半濕潤(rùn)季風(fēng)氣候區(qū)，是氣候變化的脆弱區(qū)與敏感區(qū)，同時(shí)該地區(qū)也是中國(guó)冬小麥、蘋(píng)果、獼猴桃、水蜜桃等的主產(chǎn)區(qū)，而霜凍災(zāi)害是對(duì)其威脅最大的氣象災(zāi)害?；诖?，該研究以地面 0 cm 日最低溫度作為基礎(chǔ)數(shù)據(jù)，參照《作物霜凍害等級(jí)》標(biāo)準(zhǔn)和中國(guó)科學(xué)技術(shù)藍(lán)皮書(shū)第5號(hào)《氣候》[12]，基于不同基準(zhǔn)氣候期和不同年代，探索陜西省初、終霜日、無(wú)霜期的變化規(guī)律，旨在通過(guò)系統(tǒng)分析該區(qū)域近57年來(lái)的霜日變化特征，為合理利用農(nóng)業(yè)氣候資源和有效防御霜凍提供一定的科學(xué)依據(jù)。

1?資料與方法

1.1?研究區(qū)概況

陜西省位于中國(guó)西北部，地處31°42′～39°35′N(xiāo)、105°29′～111°15′E（圖1）。境內(nèi)山塬起伏，河川縱橫，南北長(zhǎng)，東西窄。以北山和秦嶺為界，全省分為陜北高原、關(guān)中平原和陜南山地三大地貌區(qū)[13]，其中陜北高原面積最大，約占全省總面積的40%。遠(yuǎn)離水汽源地的特殊地理位置，使其形成顯著的大陸性季風(fēng)氣候。全省年平均降水量約653 mm，年平均氣溫約11.6℃，極端最高氣溫45.2 ℃，極端最低氣溫-20.6 ℃。蒸發(fā)量北大南小，年平均水面蒸發(fā)量1 200～2 100 mm，無(wú)霜期218 d左右。河川年平均徑流量420.2×108 m3。全省耕地總面積480萬(wàn)hm2，林地963萬(wàn)hm2，草地318萬(wàn)hm2。種植的農(nóng)作物主要有小麥、玉米、大豆、水稻，經(jīng)濟(jì)作物主要有獼猴桃、蘋(píng)果、蔬菜等。

1.2?資料選取?選用1961年1月1日—2017年12月31日陜西省34個(gè)地面氣象觀測(cè)站的地面日最低溫度，數(shù)據(jù)來(lái)自中國(guó)國(guó)家氣象局氣象信息中心的中國(guó)地面氣候資料日值數(shù)據(jù) V3.0[14]。其中，10個(gè)臺(tái)站位于陜北高原，11個(gè)臺(tái)站位于關(guān)中平原，13個(gè)臺(tái)站位于陜南山地。所有數(shù)據(jù)經(jīng)過(guò)嚴(yán)格的質(zhì)量控制，主要包括極值控制、缺測(cè)站點(diǎn)時(shí)間一致性檢驗(yàn)[2]。對(duì)個(gè)別臺(tái)站的缺測(cè)數(shù)據(jù)進(jìn)行插補(bǔ)處理，主要采用的是均值替換法。首先對(duì)于明顯錯(cuò)誤數(shù)據(jù)進(jìn)行檢測(cè)，確定錯(cuò)誤數(shù)據(jù)或者斷點(diǎn)，剔除或用多年平均值取代錯(cuò)誤或斷點(diǎn)數(shù)據(jù)。經(jīng)過(guò)訂正處理后的34個(gè)臺(tái)站資料具有較好的連續(xù)性。

1.3?霜日指標(biāo)計(jì)算

霜凍是指在春秋轉(zhuǎn)換季節(jié)土壤表面和植物表面溫度下降到（0 ℃以下）足以使植物遭受傷害甚至死亡的一種農(nóng)業(yè)氣象災(zāi)害[15]。根據(jù)韓榮青等[6]的分析，北方地區(qū)初霜日期基本在 8 月 1 日后，因此，以 8 月 1 日作為統(tǒng)計(jì)起始點(diǎn)，將 8 月 1 日定義為 ?用日序表示初、終霜日。初霜日期為 8 月 1 日以后第一次出現(xiàn)地面0 cm 日最低溫度≤0 ℃的日期，終霜日期為 5 月 31日之前最后一次出現(xiàn)地面0 cm日最低溫度≤0 ℃的日期，無(wú)霜期為一年中終霜日至初霜日之間的日數(shù)。

1.4?M-K突變檢驗(yàn)

Mann-Kendall檢驗(yàn)方法是一種廣泛應(yīng)用于徑流、降水、氣溫、水文序列等趨勢(shì)分析和顯著性檢驗(yàn)的非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法[16-17]，其優(yōu)點(diǎn)是既不需要時(shí)序遵循特定分布，也不受少數(shù)異常值的干擾[18]，而且可以明確突變開(kāi)始的時(shí)間和突變區(qū)域[19-20]。

對(duì)于n個(gè)樣本量的時(shí)間序列x，構(gòu)造一秩序列：

Sk=ki=1ri?（k=2，3，…，n）（1）

ri=1xj>xi0xj≤xi?（i=1，2，…，j） （2）

式中，秩序列Sk是第j時(shí)刻數(shù)值大于i時(shí)刻數(shù)值個(gè)數(shù)的累計(jì)數(shù)ri的累計(jì)值。在時(shí)間序列隨機(jī)獨(dú)立的假定下，定義統(tǒng)計(jì)量：

UFk=[Sk-E（Sk）]/Var（Sk）?（k=1，2，…，n） （3）

式中，UF1=0;E（Sk）、Var（Sk）分別是累計(jì)值Sk的均值和方差。在x?x2，…，xn相互獨(dú)立且具有相同連續(xù)分布時(shí)，它們可由下式算出：

E（Sk）=n（n-1）4

Var（Sk）=n（n-1）（2n+5）72?（k=1，2，…，n）（4）

對(duì)時(shí)間序列x=（x?x2，…，xn），當(dāng)n增加時(shí)（n>10），UFk很快收斂于標(biāo)準(zhǔn)正態(tài)分布。當(dāng)原假設(shè)為該序列無(wú)趨勢(shì)時(shí)，一般采用雙邊趨勢(shì)檢驗(yàn)，給定顯著性水平α，在正態(tài)分布表中查出臨界值Uα/2，若|UFk︱Uα/2，拒絕原假設(shè)，認(rèn)為趨勢(shì)顯著。

按時(shí)間序列x逆序xn，xn-?…，x?再重復(fù)上述過(guò)程，計(jì)算UBk，讓UBk =-UFk （k=n，n-?…，1），UB1=0。若UFk或UBk的值大于0，表明序列呈上升趨勢(shì)，反之序列呈下降趨勢(shì)。當(dāng)UFk或UBk的值超過(guò)臨界值時(shí)，表明上升或下降趨勢(shì)顯著。超過(guò)臨界線的范圍確定為出現(xiàn)突變的時(shí)間區(qū)域。若UFk或UBk 2條曲線出現(xiàn)交點(diǎn)，并且交點(diǎn)在臨界線之間，則交點(diǎn)對(duì)應(yīng)的時(shí)刻便是突變開(kāi)始的時(shí)間。

2?結(jié)果與分析

2.1?陜西省初霜日時(shí)空變化特征

2.1.1?初霜日突變檢測(cè)。

由Mann-Kendall突變檢驗(yàn)以及滑動(dòng)T檢驗(yàn)對(duì)初霜日期進(jìn)行的突變分析如圖2和圖3所示。陜西省全省初霜日期UF和UB曲線相交于1999、2001、2003年（r=0.05）（圖2a），而根據(jù)圖 3a 可知，陜西省全省初霜日期自1971年以來(lái)可能存在突變點(diǎn)的年份為 1998、1999年;綜合分析認(rèn)定，陜西省全省初霜日期在1999年發(fā)生突變。陜北地區(qū)初霜日期UF和UB曲線相交于2000年（r=0.05）（圖2b），而根據(jù)圖 3b 可知，陜北地區(qū)初霜日期自1971年以來(lái)可能存在突變點(diǎn)的年份為 1998、1999、2000、2001、2003、2005年;綜合分析認(rèn)定，陜北地區(qū)初霜日期在2000年發(fā)生突變。關(guān)中地區(qū)初霜日期UF和UB曲線相交于2001、2003年（r=0.05）（圖2c），而根據(jù)圖 3c 可知，關(guān)中地區(qū)初霜日期自1971年以來(lái)不存在可能的突變點(diǎn)年份;綜合分析認(rèn)定，關(guān)中地區(qū)初霜日期未發(fā)生突變。陜南地區(qū)初霜日期UF和UB曲線相交于2004年（r=0.05）（圖2d），而根據(jù)圖 3d 可知，陜南地區(qū)初霜日期自1971年以來(lái)可能存在突變點(diǎn)的年份為 1998、2002、2003、2004年;綜合分析認(rèn)定，陜南地區(qū)初霜日期在2004年發(fā)生突變。

2.1.2?初霜日時(shí)間變化特征。

由陜西省全省及各分區(qū)1961—2017年初霜日期的年際及年代際變化特征（圖4）可知，總體來(lái)看，近57年陜西省全省及各分區(qū)初霜日期均呈推后趨勢(shì)，陜北平均初霜日期最早且推后趨勢(shì)最顯著。1981年陜西省全省及關(guān)中、陜南地區(qū)初霜日期均歷史最早。近57年陜西省全省初霜日期以1.8 d/10 a的速率呈顯著推后趨勢(shì);全省平均初霜日期最早的發(fā)生在1981年10月10日，最遲的發(fā)生在2006年11月10日;從年代際上來(lái)看，1961—2000年全省初霜日期基本持平，2000年之后呈快速推后狀態(tài)（圖4a）。陜北地區(qū)初霜日期以1.9 d/10 a的速率呈顯著推后趨勢(shì);陜北地區(qū)終霜日期最早的發(fā)生在1972年9月22日，最遲的發(fā)生在2006年10月24日;從年代際上來(lái)看，1961—2000年陜北地區(qū)初霜日期基本維持在10月3—6日，2000年之后呈快速推后狀態(tài)，初霜日期推后至10月14日（圖4b）。關(guān)中地區(qū)這初霜日期以1.7 d/10 a的速率呈推后趨勢(shì);關(guān)中地區(qū)初霜日期最早的發(fā)生在1981年10月8日，最遲的發(fā)生在1975年11月10日;從年代際上來(lái)看，關(guān)中地區(qū)初霜日期變化特征與陜北地區(qū)基本一致，1961—2000年關(guān)中地區(qū)初霜日期基本維持在10月22—25日，2000年之后呈快速推后狀態(tài)，初霜日期推后至11月3日（圖4c）。陜南地區(qū)初霜日期以1.7 d/10 a的速率呈推后趨勢(shì);陜南地區(qū)初霜日期最早的發(fā)生在1981年10月18日，最遲的發(fā)生在2011年12月1日;從年代際上來(lái)看，1961—2000年陜南地區(qū)初霜日期基本維持在11月6—10日，2000年之后呈快速推后狀態(tài)，初霜日期推后至11月16日（圖4d）。

2.1.3?初霜日空間變化特征。從陜西省初霜日期空間分布（圖5a）可以看出，近57年陜西省各地初霜日期由北到南呈顯著帶狀分布特征。初霜日期出現(xiàn)最早的主要是陜北地區(qū)，其中吳旗—靖邊—橫山—榆林—神木一線的初霜日期在9月29日—10月9日，延安中東部地區(qū)、陜南的留壩地區(qū)初霜日期在10月9—19日;初霜日期出現(xiàn)最遲的主要是陜南地區(qū)，其中陜南中部地區(qū)、關(guān)中的太白地區(qū)、漢中的寧強(qiáng)地區(qū)初霜日期在11月9—29日，比最早出現(xiàn)初霜日期晚了將近60 d;省內(nèi)其余地區(qū)的初霜日期在10月19日—11月9日。

從陜西省初霜日期年變化率（圖5b）可以看出，近57年陜西省除石泉地區(qū)外，全省各地初霜日期均呈推后趨勢(shì)。推后趨勢(shì)最顯著的是陜北的靖邊和陜南的鎮(zhèn)安地區(qū)，推后趨勢(shì)達(dá)3.0 d/10 a以上，此外陜北西北部地區(qū)、陜南中部地區(qū)、關(guān)中東部地區(qū)的推后趨勢(shì)也達(dá)2.0 d/10 a以上。推后趨勢(shì)最不顯著的是陜北的綏德、關(guān)中的耀縣、陜南的寧強(qiáng)地區(qū)，推后趨勢(shì)不足1.0 d/10 a。省內(nèi)其余地區(qū)的推后趨勢(shì)基本在1.0～2.0 d/10 a。

2.2?陜西省終霜日時(shí)空變化特征

2.2.1?終霜日突變檢測(cè)。

由Mann-Kendall突變檢驗(yàn)以及滑動(dòng)T檢驗(yàn)對(duì)終霜日期進(jìn)行的突變分析如圖6和圖7所示。陜西省全省終霜日期UF和UB曲線相交于2006年（r=0.05）（圖6a），而根據(jù)圖 7a 可知，陜西省全省終霜日期自1971年以來(lái)，可能存在突變點(diǎn)的年份為 1979、1980、1983、1996年;綜合分析認(rèn)定，陜西省全省終霜日期未發(fā)生突變。陜北地區(qū)終霜日期UF和UB曲線相交于2004年（r=0.05）（圖6b），而根據(jù)圖 7b 可知，陜北地區(qū)終霜日期自1971年以來(lái)可能存在突變點(diǎn)的年份為 1980—1984、1997、2000—2002、2004、2008年;綜合分析認(rèn)定，陜北地區(qū)全省終霜日期在2004年發(fā)生突變。關(guān)中地區(qū)終霜日期UF和UB曲線相交于2000年（r=0.05）（圖6c），而根據(jù)圖 7c 可知，關(guān)中地區(qū)終霜日期自1971年以來(lái)可能存在突變點(diǎn)的年份為 1996年;綜合分析認(rèn)定，關(guān)中地區(qū)終霜日期未發(fā)生突變。陜南地區(qū)終霜日期UF和UB曲線相交于1996、1998年（r=0.05）（圖6d），而根據(jù)圖 7d 可知，陜南地區(qū)終霜日期自1971年以來(lái)不存在可能的突變點(diǎn)年份;綜合分析認(rèn)定，陜南地區(qū)全省終霜日期未突變。

2.2.2?終霜日時(shí)間變化特征。由1961—2017年陜西省全省及各分區(qū)終霜日期的年際及年代際變化特征（圖8）可知，總體來(lái)看，近57年陜西省全省及各分區(qū)終霜日期均呈提前趨勢(shì)，陜北終霜日期最遲且提前趨勢(shì)最顯著，關(guān)中次之。1961年陜西省全省及關(guān)中、陜南地區(qū)終霜日期均歷史最遲。近57年陜西省全省終霜日期以2.2 d/10 a的速率呈顯著提早趨勢(shì);全省平均終霜日期最遲的發(fā)生在1961年4月28日，最早的發(fā)生在2008年3月29日;從年代際上來(lái)看，1961—2000年全省終霜日期基本持平，2000年之后呈快速提早狀態(tài)。陜北地區(qū)終霜日期以2.5 d/10 a的速率呈顯著提早趨勢(shì);陜北地區(qū)終霜日期最遲的發(fā)生在1967年5月18日，最早的發(fā)生在2010年4月17日;從年代際上來(lái)看，除20世紀(jì)60年代外，整體上呈階梯狀分布，20世紀(jì)70年代比80、90年代遲8 d左右，20世紀(jì)80、90年代比2000年以后遲8 d左右。關(guān)中地區(qū)終霜日期以2.2 d/10 a的速率呈顯著提早趨勢(shì);關(guān)中地區(qū)終霜日期最遲的發(fā)生在1961年4月27日，最早的發(fā)生在2013年3月22日;從年代際上來(lái)看，除20世紀(jì)60年代外，整體上呈階梯狀分布，20世紀(jì)70、80年代比90年代遲4 d左右， 20世紀(jì)90年代比2000年以后遲4 d左右。陜南地區(qū)終霜日期以1.9 d/10 a的速率呈提早趨勢(shì);陜南地區(qū)終霜日期最遲的發(fā)生在1961年4月11日，最早的發(fā)生在1996年3月6日;從年代際上來(lái)看，整體上呈階梯狀分布，20世紀(jì)60年代比70、80、90年代遲5 d左右， 20世紀(jì)90年代比2000年以后遲5 d左右。

2.2.3?終霜日空間變化特征。

從陜西省終霜日期空間分布（圖9a）可以看出，近57年陜西省各地終霜日期與初霜日期分布特征一致，由北到南呈顯著帶狀分布特征。終霜日期出現(xiàn)最遲的主要是陜北地區(qū)，其中吳旗—靖邊—橫山—榆林一線的終霜日期在5月3—13日，延安地區(qū)、榆林東部地區(qū)、陜南的留壩地區(qū)終霜日期在4月23日—5月3日。終霜日期出現(xiàn)最早的主要是陜南地區(qū)，終霜日期基本在4月3日以前，其中陜南中部地區(qū)、關(guān)中太白地區(qū)、漢中寧強(qiáng)地區(qū)終霜日期在3月14—24日，比最遲出現(xiàn)終霜日期晚了近30 d。省內(nèi)其余地區(qū)的終霜日期在4月3—23日。

從陜西省終霜日期年變化率（圖9b）可以看出，近57年陜西省全省各地終霜日期均呈提前趨勢(shì)。提前趨勢(shì)最顯著的是陜北的靖邊和陜南的石泉地區(qū)，提前趨勢(shì)達(dá)4.0 d/10 a以上，此外榆林中部地區(qū)、陜南鎮(zhèn)安地區(qū)、關(guān)中的長(zhǎng)武、武功、蒲城地區(qū)的提前趨勢(shì)也達(dá)3.0 d/10 a以上。提前趨勢(shì)最不顯著的是陜北的定邊和吳旗、關(guān)中的咸陽(yáng)地區(qū)、陜南的漢中、商州、鎮(zhèn)坪地區(qū)，提前趨勢(shì)不足1.0 d/10 a。省內(nèi)其余地區(qū)的提前趨勢(shì)基本在1.0～3.0 d/10 a。

3?結(jié)論

（1）近57年陜西省全省初霜日期在1999年發(fā)生突變。其中，陜北地區(qū)和陜南地區(qū)初霜日期分別在2000和2004年發(fā)生突變，關(guān)中地區(qū)初霜日期未發(fā)生突變。陜西省全省、關(guān)中地區(qū)和陜南地區(qū)終霜日期均未發(fā)生突變，但陜北地區(qū)終霜日期在2004年發(fā)生突變。

（2）近57年陜西省全省及各分區(qū)初霜日期均呈推后趨勢(shì)，陜北平均初霜日期最早且推后趨勢(shì)最顯著。其中，陜西省全省初霜日期以1.8 d/10 a的速率呈顯著推后趨勢(shì)，陜北地區(qū)、關(guān)中地區(qū)和陜南地區(qū)初霜日期分別以1.9、1.7和1.7 d/10 a的速率呈推后趨勢(shì)，推后趨勢(shì)最顯著的是陜北的靖邊和陜南的鎮(zhèn)安地區(qū)，推后趨勢(shì)達(dá)3.0 d/10 a以上。近57年陜西省全省及各分區(qū)終霜日期均呈提前趨勢(shì)，陜北終霜日期最遲且提前趨勢(shì)最顯著。其中，陜西省全省終霜日期以2.2 d/10 a的速率呈顯著提早趨勢(shì)，陜北地區(qū)、關(guān)中地區(qū)和陜南地區(qū)終霜日期分別以2.5、2.2和1.9 d/10 a的速率呈提早趨勢(shì)，提前趨勢(shì)最顯著的是陜北的靖邊和陜南的石泉地區(qū)，提前趨勢(shì)達(dá)4.0 d/10 a以上。

（3）近57年陜西省各地初霜日期由北到南呈顯著帶狀分布特征。初霜日期出現(xiàn)最早的是吳旗—靖邊—橫山—榆林—神木一線，出現(xiàn)最遲的是陜南中部地區(qū)、關(guān)中太白地區(qū)、漢中寧強(qiáng)地區(qū)，比最早出現(xiàn)初霜日期晚了近60 d。陜西省各地終霜日期與初霜日期分布特征一致，由北到南呈顯著帶狀分布特征。終霜日期出現(xiàn)最遲的是吳旗—靖邊—橫山—榆林一線，出現(xiàn)最早的是陜南中部地區(qū)、關(guān)中太白地區(qū)、漢中寧強(qiáng)地區(qū)，比最遲出現(xiàn)終霜日期晚了近45 d。

參考文獻(xiàn)

[1] IPCC.Emissions scenarios：A special report of working group III of the intergovernmental panels on climate change[M].Cambridge：Cambridge Press，2000.

[2] 周丹，張勃，任培貴，等.基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的陜西省近50 a干旱特征分析[J].自然資源學(xué)報(bào)，2014，29（4）：677-688.

[3] 沈永平，王國(guó)亞.IPCC第一工作組第五次評(píng)估報(bào)告對(duì)全球氣候變化認(rèn)知的最新科學(xué)要點(diǎn)[J].冰川凍土， 2013，35（5）：1068-1076.

[4] 《大氣科學(xué)辭典》編委會(huì).大氣科學(xué)辭典[M].北京：氣象出版社，1994.

[5] 周曉宇，趙春雨，崔妍，等.1961—2013年中國(guó)東北地區(qū)初終霜日及無(wú)霜期的氣候變化特征[J].自然資源學(xué)報(bào)，2017，32（3）：494-506.

[6] 韓榮青，李維京，艾婉秀，等.中國(guó)北方初霜凍日期變化及其對(duì)農(nóng)業(yè)的影響[J].地理學(xué)報(bào)，2010，65（5）：525-532.

[7] 陳少勇，鄭延祥，樓望萍，等.中國(guó)西北地區(qū)初霜凍的氣候變化特征[J].資源科學(xué)，2013，35（1）：165-172.

[8] TERANDO A，EASTERLING W E，KELLER K，et al.Observed and modeled twentieth?century spatial and temporal patterns of selected agro-climate indices in North America[J].Journal of climate，2012，25（2）：473-490.

[9] 馬尚謙，張勃，唐敏，等.1960—2015 年淮河流域初終霜日時(shí)空變化分析[J].中國(guó)農(nóng)業(yè)氣象，2018，39（7）：468-478.

[10] 王景紅，劉璐，高峰，等.陜西省富士系蘋(píng)果晚霜凍指標(biāo)試驗(yàn)綜合方法研究[J].中國(guó)農(nóng)學(xué)通報(bào)，2014，30（25）：218-221.

[11] 柏秦鳳，李星敏，朱琳.近50年陜西省無(wú)霜期的變化及果區(qū)霜凍風(fēng)險(xiǎn)分布[J].干旱區(qū)資源與環(huán)境，2013，39（7）：65-70.

[12] 國(guó)家科學(xué)技術(shù)委員會(huì).中國(guó)科學(xué)技術(shù)藍(lán)皮書(shū) 第5號(hào) 氣候[M].北京：科學(xué)技術(shù)文獻(xiàn)出版社，1990：25-26.

[13] 史興民，陳玲俠，張瑞.陜西省近 50 年地貌災(zāi)害的時(shí)空分布[J].干旱區(qū)研究，2011，28（2）：349-354.

[14] 國(guó)家氣象信息中心.中國(guó)地面氣候資料日值數(shù)據(jù)（V3.0）[EB/OL].[2019-04-05].http//data.cma.cn/data/cdc-detail /dataCode/SURF_CLI_CHN_MUL_DAY_V3.0.html.

[15] 馮秀藻，陶炳炎.農(nóng)業(yè)氣象學(xué)原理[M].北京：氣象出版社，1991：127-143.

[16] GAN T Y.Hydroclimatic trends and possible climatic warming in the Canadian Prairies[J].Water resource research，1998，34（11）：3009-3015.

[17] 褚健婷，夏軍，許崇育，等.氣海河流域氣象和水文降水資料對(duì)比分析及時(shí)空變異[J].地理學(xué)報(bào)，2009，64（9）：1083-1092.

[18] 馬柱國(guó)，邵麗娟.中國(guó)北方近百年干濕變化與太平洋年代際振蕩的關(guān)系[J].大氣科學(xué)，2006，30（3）：464-474.

[19] 周丹，張勃，李小亞，等.1961～2010年中國(guó)大陸地面氣候要素變化特征分析[J].長(zhǎng)江流域資源與環(huán)境，2014，23（4）：549-558.

[20] 張璐，周躍.Mann-Kendall檢驗(yàn)及其在河流懸沙濃度時(shí)間序列分析中的運(yùn)用[J].亞熱帶水土保持，2007，19（4）：13-16.