太白山保護(hù)區(qū)氣溫變化及其空間差異性分析*

劉榮娟，白紅英，馬新萍，馮林林，秦 進(jìn)，程德強(qiáng)

（西北大學(xué)城市與環(huán)境學(xué)院，西安 710127）

?

太白山保護(hù)區(qū)氣溫變化及其空間差異性分析*

劉榮娟，白紅英**，馬新萍，馮林林，秦 進(jìn)，程德強(qiáng)

（西北大學(xué)城市與環(huán)境學(xué)院，西安 710127）

摘要：根據(jù)太白山保護(hù)區(qū)周邊15個氣象站1959-2013年以及保護(hù)區(qū)內(nèi)11個高山氣象站2012-2014年的數(shù)據(jù)，采用Kriging插值法、一元線性回歸分析法、Mann-Kendall突變檢驗法等，結(jié)合ArcGIS分析太白山保護(hù)區(qū)的氣溫變化趨勢和突變特征，并對比保護(hù)區(qū)南北坡和高海拔區(qū)域（海拔≥2800m）與保護(hù)區(qū)整體的氣溫變化差異。結(jié)果表明：（1）1959-2013年太白山保護(hù)區(qū)的氣溫呈波動上升趨勢，線性傾向率為0.329℃·10a?1（P＜0.05）。氣溫變化空間差異明顯，東南部和西部地區(qū)氣溫上升趨勢顯著（P＜0.05或P＜0.01），最高達(dá)0.383℃·10a?1；四季氣溫均呈上升趨勢，冬季升溫最快，春季最慢。（2）M-K檢驗表明，1998年前后保護(hù)區(qū)整體氣溫發(fā)生增溫突變。突變后氣溫上升加快，上升趨勢達(dá)顯著水平的區(qū)域為86.76%。（3）空間差異分析顯示，55a來太白山保護(hù)區(qū)北坡年平均氣溫和升溫速率明顯低于南坡，北坡和南坡年平均氣溫線性傾向率分別為0.247和0.291℃·10a?1。高海拔區(qū)域的年平均氣溫線性傾向率為0.244℃·10a?1，低于保護(hù)區(qū)整體，其氣溫突變發(fā)生在2001年，突變前東南部升溫最快，突變后西部升溫速率降低，中部地區(qū)升溫明顯加快。

關(guān)鍵詞：太白山保護(hù)區(qū)；氣溫變化；空間差異；Kriging插值；氣溫變化傾向率

劉榮娟，白紅英，馬新萍，等.太白山保護(hù)區(qū)氣溫變化及其空間差異性分析[J].中國農(nóng)業(yè)氣象，2016，37（1）:1-10

氣候變化將嚴(yán)重影響全球各地的生產(chǎn)生活[1]。IPCC第五次評估報告指出[2]，全球氣候系統(tǒng)變暖的事實毋庸置疑，1951年以來全球地表溫度的平均上升速率為0.12（0.08～0.14）℃·10a?1。第二、第三次氣候變化國家評估報告[3-4]指出，1951-2009年中國陸地表面平均氣溫上升了1.38℃，明顯高于同期全球或北半球其它區(qū)域的平均增溫速率，尤其近15a處于近百年來升溫最快的階段。有關(guān)氣候變化的研究逐漸受到重視。

研究區(qū)域氣候在全球變化背景下的變化特征具有重要意義，而利用GIS技術(shù)生成高分辨率的氣候數(shù)據(jù)，分析氣候資源的空間分布及趨勢模擬，能達(dá)到較好的效果[5-6]。段春鋒等[7-10]對中國整體及區(qū)域氣溫變化的研究表明，各區(qū)域氣溫變化并不一致。施雅風(fēng)等[11]認(rèn)為，西北氣候由暖干型向暖濕轉(zhuǎn)型；吳昊旻等[12]對長三角地區(qū)的季節(jié)氣候變化特征進(jìn)行了研究，認(rèn)為該區(qū)域春夏季時間延長，秋冬季縮短，氣溫呈上升趨勢；楊再禹等[13]研究發(fā)現(xiàn)，黔東南地區(qū)1961-2009年氣候向暖干化發(fā)展。太白山位于秦嶺山脊線處，高差大，坡度陡，地貌特征復(fù)雜，是對全球變化響應(yīng)敏感的地區(qū)之一。其所處的秦嶺地帶是中國南北重要的地理分界線，引起眾多關(guān)注[14-16]。對秦嶺氣候分界線[17]、南北氣候變化差異[18]等的研究表明，秦嶺南北坡的氣溫具有同步增加的趨勢[19]。但目前關(guān)于太白山區(qū)氣候的研究報道鮮見，尤以針對太白山復(fù)雜的地形內(nèi)氣候變化差異的相關(guān)研究缺乏。本文以太白山周邊氣象站點氣象數(shù)據(jù)為基礎(chǔ)，基于Kriging插值法和數(shù)理統(tǒng)計方法，對全球氣溫明顯升高的近半個世紀(jì)以來太白山保護(hù)區(qū)的氣候變化特征進(jìn)行研究，將其與秦嶺山地和西北地區(qū)進(jìn)行比較，并探討南北坡之間、高海拔區(qū)域與保護(hù)區(qū)整體之間的氣溫變化特征的異同。

1 資料與方法

1.1 研究區(qū)概況

太白山國家級自然保護(hù)區(qū)地處陜西秦嶺中段，總面積563.25km2，位于33°49′-34°08′N，107°22′-107°52′E，地跨陜西眉縣、太白和周至3個縣[20]。海拔落差大，最高處3771.2m，位于主峰拔仙臺，最低處1028m。屬內(nèi)陸季風(fēng)氣候，年平均降水量500～1100mm，年平均氣溫5.9～7.5℃，南北坡氣溫變化有所差異。

1.2 數(shù)據(jù)來源

25m分辨率的DEM高程數(shù)據(jù)來自國家測繪局。太白山保護(hù)區(qū)周邊15個基本氣象站點1959-2013年以及保護(hù)區(qū)內(nèi)11個高山氣象站點2012-2014年氣象數(shù)據(jù)來自陜西省氣象局，站點分布見圖1。基本氣象站長期氣象資料主要用于氣溫變化特征分析，2012-2014年高山氣象站氣象資料主要用于山地氣溫垂直遞減率的計算和空間插值結(jié)果的驗證。

圖1 研究區(qū)及其周邊各氣象站點位置(a)和太白山保護(hù)區(qū)地形圖(b)Fig.1 Meteorological stations(a) and topography(b) in the protected area of Taibai Mountain注：圖b中外部實線包圍的是太白山保護(hù)區(qū)即I區(qū)，內(nèi)部填充部分為高海拔區(qū)域即II區(qū)。下同Note: In Fig.b， the range in outside solid line is I area which is Taibai Mountain protected area.The inside filling range is II area which is high altitude area （altitude ≥2800m）.The same as below

1.3 研究方法

1.3.1 一元線性回歸分析法

該方法可以對氣溫柵格圖內(nèi)每個像元的變化趨勢進(jìn)行模擬，采用最小二乘法計算其傾向值，以估計時間序列與氣溫序列的回歸系數(shù)，計算式為

式中，θslope為趨勢斜率（℃·10a?1），體現(xiàn)像元在一定時間內(nèi)的變化趨勢，n為年份累加，Ti為第i年的年平均氣溫（℃），θslope＞0說明氣溫升高；θslope＜0則氣溫降低。

1.3.2 Mann-Kendall檢驗法

均值突變是氣候突變的一種，能一定程度上反應(yīng)基本氣候狀況的變化，而M-K檢驗則可有效檢測氣溫均值突變。其步驟為先將原始時間序列定義為X1，X2，…，Xn，ni表示時間序列中第i個樣本Xi大于Xj的累計數(shù)，dk為統(tǒng)計量，即

式中，dk的均值E（dk）和方差Var（dk）分別為

將dk標(biāo)準(zhǔn)化后得

U（dk）所得曲線UF若超過信度線（u0.05＝±1.96，u0.01＝±2.56），表示變化趨勢顯著，UF曲線與反序列計算得出的UB曲線在置信區(qū)間內(nèi)的交點即為突變的開始點[21-22]。

1.3.3 Kriging插值法及其驗證

利用Kriging插值法[23]對氣溫要素進(jìn)行空間插值，變異函數(shù)均選擇球面函數(shù)。為保證山地氣溫插值結(jié)果的準(zhǔn)確性，插值過程中需考慮高程因素的影響[24]。首先確定氣溫垂直遞減率，對26個氣象站數(shù)據(jù)進(jìn)行選擇整理，結(jié)合海拔高度分別計算南北坡年、季、月的氣溫垂直遞減率[25]，保護(hù)區(qū)整體的氣溫垂直遞減率為南北坡平均值，而在保護(hù)區(qū)南北坡分別進(jìn)行插值時則使用各自的氣溫垂直遞減率?？紤]到高海拔區(qū)域與保護(hù)區(qū)整體的差異性，上述計算的遞減率并不適用。太白山11個高山站點平均海拔2272m，更接近保護(hù)區(qū)高海拔地區(qū)的氣溫特征，利用2012-2014年氣象數(shù)據(jù)計算其氣溫遞減率，再次對高海拔地區(qū)的氣溫進(jìn)行插值。表1為太白山高海拔區(qū)域與保護(hù)區(qū)整體氣溫垂直遞減率，1月和7月均為兩個區(qū)域氣溫遞減率最低和最高的月份，由表可見，高海拔區(qū)域各個月份、季節(jié)和年平均氣溫的遞減率均高于保護(hù)區(qū)整體，氣溫隨海拔升高而降低的現(xiàn)象在高海拔地區(qū)更為明顯，且時間上以冬季差異最大，空間上表現(xiàn)為北坡差異更大。

利用表1中的垂直遞減率將各站點實測氣溫值訂正至海平面氣溫，在ArcGIS10.1中進(jìn)行插值，將插值后的柵格數(shù)據(jù)與區(qū)域DEM數(shù)據(jù)相疊加，利用柵格計算工具，最終得到具有地形特征的太白山保護(hù)區(qū)整體、保護(hù)區(qū)南北坡和高海拔區(qū)域的氣溫柵格數(shù)據(jù)，即

式中，T0為插值后生成的海平面的溫度（℃），A為研究區(qū)氣溫垂直遞減率（℃·m?1），E為海拔（m） ，T為最終生成的研究區(qū)氣溫（℃）。

插值后輸出各站點的交叉驗證結(jié)果，分析插值誤差大小，交叉驗證結(jié)果顯示平均標(biāo)準(zhǔn)差為0.877。將太白山保護(hù)區(qū)內(nèi)4個站點2012、2013年的實測溫度數(shù)據(jù)與插值結(jié)果進(jìn)行對比，其結(jié)果顯示（表2），插值誤差絕大部分在±1℃以內(nèi)，表明該插值結(jié)果具有較高的可信性，可以用來分析太白山保護(hù)區(qū)的氣溫變化特征。

表1 太白山保護(hù)區(qū)（I）和高海拔區(qū)域（II）內(nèi)南北坡氣溫的垂直遞減率（℃·m?1）Table 1 Temperature lapse rate of southern and northern slopes in I area and II area（℃·m?1）

表2 太白山區(qū)域站點2012、2013年不同時段平均氣溫插值結(jié)果與實測數(shù)據(jù)的對比(℃)Table 2 Contrast between average temperature interpolation results with measured data in regional sites of Taibai Mountain during different periods in 2012 and 2013(℃)

2 結(jié)果與分析

2.1 太白山保護(hù)區(qū)氣溫變化趨勢分析

2.1.1 年平均氣溫

根據(jù)55a插值后的柵格數(shù)據(jù)，太白山保護(hù)區(qū)年平均氣溫為3.540℃，其中高海拔區(qū)域（海拔≥2800m）為0.225℃，各區(qū)域年平均氣溫的空間分布和變化趨勢見圖2。由圖2b可見，55a來兩個研究區(qū)域內(nèi)年平均氣溫均表現(xiàn)為波動上升趨勢，整個保護(hù)區(qū)年平均氣溫的線性傾向率為0.329℃·10a?1（P＜0.05），高海拔區(qū)域為0.244℃·10a?1（P＜0.05）。由圖3可見，各柵格的氣溫變化率存在明顯的空間差異，太白山保護(hù)區(qū)內(nèi)大部分地區(qū)線性變化較快，年平均氣溫呈顯著（P＜0.05）上升趨勢，東南部和西部地區(qū)呈極顯著（P＜0.01）上升趨勢，最高達(dá)0.383℃·10a?1，其它地區(qū)的線性變化趨勢不顯著。

2.1.2 季平均氣溫

由圖4可見，太白山保護(hù)區(qū)四季平均氣溫變化率存在不同的空間差異。除春季有小范圍不顯著下降外，其余均呈上升趨勢。其中冬季平均氣溫上升最快，其線性傾向率為0.468℃·10a?1（P＜0.05），最高達(dá)0.644℃·10a?1，春季氣溫上升最慢。秋、冬季氣溫變化率的空間分布相似，東南部線性變化較快，氣溫呈顯著（P＜0.05）或極顯著（P＜0.01）上升趨勢；春夏兩季線性變化較緩，其中北部變化最快，上升趨勢達(dá)到顯著或極顯著水平（P＜0.05或P＜0.01）。西部地區(qū)四季升溫均不顯著。

圖2 1959-2013年太白山保護(hù)區(qū)（I區(qū)和II區(qū)）年平均氣溫的空間分布(a)及其年際變化(b)Fig.2 Spatial distribution(a) and variation (b) of annual mean temperature in I area and II area during 1959-2013

圖3 1959-2013年太白山保護(hù)區(qū)年平均氣溫變化率的空間分布(a)及其F檢驗結(jié)果(b)Fig.3 Spatial distribution(a) of annual mean temperature change rate and its F-test result(b) in protected area during 1959-2013

圖4 1959-2013年太白山保護(hù)區(qū)季平均氣溫變化趨勢的空間分布及F檢驗結(jié)果Fig.4 Spatial distribution and its F-test result of seasonal mean temperature change rate in protected area during 1959-2013

2.2 氣溫突變特征分析

由圖5可見，研究期內(nèi)（1959-2013年），太白山保護(hù)區(qū)年平均氣溫在1984-1995年有所下降，其余時段內(nèi)均呈上升趨勢；UF（I）與UB（I）在u0.05＝±1.96（P＜0.05）的兩條臨界線內(nèi)相交于1998年，說明太白山保護(hù)區(qū)年平均氣溫于1998年前后發(fā)生增溫突變，并在2002年超過1.96的臨界線（P＜0.05），達(dá)到顯著增溫水平；高海拔地區(qū)年平均氣溫在1986年之前呈波動下降或上升趨勢，UF（II）與UB（II）在u0.05＝±1.96（ P＜0.05）的兩條臨界線內(nèi)相交于2001年，說明太白山高海拔區(qū)域年平均氣溫于2001年前后發(fā)生增溫突變，并在2003年超過1.96的臨界線（ P＜0.05），達(dá)到顯著增溫水平。

根據(jù)上述結(jié)果，將太白山保護(hù)區(qū)氣溫分為1959-1998年和1999-2013年2個時段進(jìn)行分析。由圖6可見，突變前后氣溫均呈上升趨勢，但突變前的氣溫變化率明顯小于突變后，且其空間分布差異較大。1959-1998年除西部區(qū)域升溫顯著外，其它區(qū)域平均氣溫變化均不顯著；1999-2013年太白山保護(hù)區(qū)平均氣溫線性傾向率為0.348℃·10a?1（P＜0.05），以中部地區(qū)最高，達(dá)到0.385℃·10a?1（P＜0.05）。經(jīng)統(tǒng)計突變年前后氣溫變化顯著性水平的柵格數(shù)，發(fā)現(xiàn)全保護(hù)區(qū)年平均氣溫上升顯著的區(qū)域所占比例由突變前的43.59%上升至突變后的86.76%。

圖5 1959-2013年太白山保護(hù)區(qū)（I）與高海拔區(qū)域（∏）內(nèi)年平均氣溫M-K突變檢驗Fig.5 M-K test of annual temperature in protected area(I)and high altitude area(II)during 1959-2013

圖6 1959-2013年氣溫突變前(a)、后(b)太白山保護(hù)區(qū)（I區(qū)）年平均氣溫變化率的空間分布及其F檢驗結(jié)果Fig.6 Spatial distribution of annual mean temperature change rate and its F-test result in protected area (I area) before(a) and after(b) mutation year

對2012-2014年11個高山站氣象數(shù)據(jù)進(jìn)行插值并計算變化率后可知，高海拔區(qū)域（海拔≥2800m）突變前（1959-2001年）年平均氣溫為-0.108℃，突變后（2002-2013年）年平均氣溫為0.959℃。圖7顯示了突變前后該區(qū)域氣溫變化率的空間分布。由圖可見，突變前（圖7a）氣溫變化趨勢不顯著，以東南部上升速度最快；突變后（圖7b）氣溫變化線性傾向率達(dá)到0.279℃·10a?1（P＜0.05），其中西部升溫速率降低，中部地區(qū)升溫明顯加快，最高達(dá)到0.470℃·10a?1（P＜0.05）。

2.3 太白山保護(hù)區(qū)南北坡氣溫變化差異分析

太白山地跨秦嶺分水嶺南北兩側(cè)，南坡平均海拔2773m，北坡平均海拔2304m，氣候特征有所差異。對Kriging插值結(jié)果計算表明，北坡年平均氣溫（5.179℃）高于南坡（2.427℃），但突變前后北坡增溫速率（0.247℃·10a?1）均明顯低于南坡（0.291℃·10a?1）。由圖8可知，南北坡氣溫線性傾向率在突變后均有所增加，且南北坡增溫速率差異加大；空間分布上北坡變化率分布具有明顯的海拔變化特征，由北向南隨著海拔升高，升溫速率逐漸增加，而南坡主要表現(xiàn)為東西向的差異。

圖7 1959-2013年突變前后太白山高海拔區(qū)域（II區(qū)）年平均氣溫變化率的空間分布Fig.7 Spatial distribution of annual mean temperature change rate in high altitude area (II area) before(a) and after(b) mutation year

圖8 1959-2013年太白山保護(hù)區(qū)突變前(a)后(b)南北坡年平均氣溫變化率的空間分布Fig.8 Spatial distribution of annual mean temperature change rate of northern and southern slopes in protected area of Taiabai Mountain before(a) and after(b) mutation during 1959-2013

3 討論與結(jié)論

（1）1959-2013年，太白山保護(hù)區(qū)年平均氣溫呈波動上升趨勢，線性斜率為0.329℃·10a?1（P＜0.05），可以分為3個不同階段：1959-1979年的波動增溫期，1979-1998年的平緩升溫期和1999-2013年的迅速升溫期。四季平均氣溫均呈上升趨勢，冬季升溫最快，春季最慢。太白山保護(hù)區(qū)無論年或季節(jié)的平均升溫速率均高于中國平均速率（0.04℃·10a?1），且冬季氣溫升高對保護(hù)區(qū)的氣溫變化影響最大，這與全球變暖的特征頗為相似。同時段內(nèi)，秦嶺關(guān)中、陜南地帶[26]氣溫傾向率（0.21℃·10a?1、0.12℃·10a?1）均稍低于此；而干旱區(qū)[27]（0.39℃·10a?1）的氣溫傾向率則略高于此。王勁松等[28]的研究結(jié)果表明，西北區(qū)1960-2005年氣溫變化幅度為0.37℃·10a?1，1985年之后是該地區(qū)增溫的主要階段，并以冬季增溫率最高。馬新萍[29]研究發(fā)現(xiàn)，1959-2013年秦嶺山地氣溫變化率為0.174℃·10a?1，太白山在秦嶺山地內(nèi)屬升溫顯著地區(qū)，與本文研究結(jié)論相符合。從全球變化的大環(huán)境來看，太白山保護(hù)區(qū)將繼續(xù)保持升溫趨勢。

（2）M-K檢驗結(jié)果顯示，1998年前后保護(hù)區(qū)氣溫發(fā)生增溫突變，而高海拔地區(qū)在2001年發(fā)生突變。突變前年氣溫平均變化率明顯小于突變后。研究時段內(nèi)，西部地區(qū)由升溫最快變?yōu)樯郎刈盥膮^(qū)域之一，東部地區(qū)和中部高海拔區(qū)域升溫速率呈加大趨勢；且氣溫上升趨勢達(dá)到顯著水平的區(qū)域比例大幅增加。

關(guān)于區(qū)域氣溫的突變特征前人做了許多研究，但結(jié)論并不一致。對秦嶺地區(qū)氣溫突變的研究通常分南北坡分別進(jìn)行，研究表明其突變基本均發(fā)生在20世紀(jì)90年代，這與全球氣候變化的大形勢相吻合。宋佃星等[18]認(rèn)為秦嶺北部在1989 年發(fā)生突變，升溫幅度為0.34℃·10a?1；秦嶺南部在1992 年發(fā)生突變，升溫幅度為0.14℃·10a?1，兩地同樣夏冬兩季升溫幅度大。李雙雙等[26]和白晶等[30]都將秦嶺分為關(guān)中和陜南進(jìn)行分析，前者認(rèn)為1959-2009年秦嶺南北氣候變化具有同步性，增溫趨勢明顯，關(guān)中地區(qū)氣溫突變（1995 年）早于陜南（1998 年）；后者同樣認(rèn)為關(guān)中地區(qū)氣溫突變時期（1990 年）早于陜南（1998年），兩者時間上略有差異，但均與本文結(jié)論類似。20世紀(jì)80年代北半球發(fā)生變暖的突變[31]，從90年代初期，中國東北、華北、青藏高原及西北地區(qū)先后進(jìn)入增暖期[32]，而秦嶺氣溫突變發(fā)生時間滯后，太白山保護(hù)區(qū)及其高海拔區(qū)域突變滯后現(xiàn)象則尤其明顯。

（3）太白山保護(hù)區(qū)內(nèi)氣溫變化率存在明顯的空間差異。北坡氣溫高于南坡，但其升溫速率較低，且北坡基本表現(xiàn)為由北向南升溫速率逐漸增加，而南坡主要表現(xiàn)為東西向的差異。太白山區(qū)2800m以上的高海拔區(qū)域主要分布有太白紅杉、灌叢、高山草甸等植被類型[33]，該區(qū)域幾乎沒有人類活動的影響，其氣溫變化特征與低海拔區(qū)域不同，更能體現(xiàn)出排除人為干擾和城市熱島效應(yīng)的區(qū)域氣候變化特征。高海拔區(qū)域氣溫遞減率高于保護(hù)區(qū)整體，但二者波動趨勢一致。高海拔區(qū)域南北坡均會受到冬季風(fēng)強(qiáng)烈影響，冬季遞減率與保護(hù)區(qū)整體差異尤為明顯，又因南坡氣溫普遍較低，所以北坡氣溫垂直遞減率普遍高于南坡，且與保護(hù)區(qū)整體的差異更大。高海拔區(qū)域突變前后氣溫變化率均低于太白山保護(hù)區(qū)整體，且空間分布有較大差異，西部升溫速率降低，中部地區(qū)升溫顯著。

本文分析了不同海拔不同坡向的山地氣溫變化特征以及氣溫垂直遞減率的差異，但氣溫變化的影響因素較為復(fù)雜，太白山保護(hù)區(qū)尤其高海拔地區(qū)的氣溫變化是否與大氣環(huán)流等有關(guān)還需要進(jìn)一步探討。

參考文獻(xiàn)References

[1]Nielsen J ?，D’haen A L.Asking about climate change: reflections on methodology in qualitative climate change research published in global environmental change since 2000[J].Global Environmental Change，2014，（24）:402-409.

[2]IPCC.Working group I contribution to the IPCC fifth assessment report climate change 2013:the physical science basis summary for policymakers[R/OL].Http:// www.climatechange2013.org/images/uploads/WGIA5 SPM_Approved27Sep 2013.pdf，2015-1-13.

[3]《第二次氣候變化國家評估報告》編寫委員會.第二次氣候變化國家評估報告[M].北京: 科學(xué)出版社，2011:80-105. Editorial Board on Second Chinese National Assessment Report on Climate Change.The second national assessment report on climate change[M].Beijing:Science Press，2011:80-105.（in Chinese）

[4]侯露露，王海林.第三次氣候變化國家評估報告 [EB/OL]. http://www.cma.gov.cn/2011xwzx/2011xqxxw/ 2011xqxyw/ 201412/t20141207_269047.html，2015-01-13. Hou L L，Wang H L.China published national climate change third assessment repor[EB/OL].http://www.cma.gov.cn/ 2011 xwzx/2011xqxxw/2011xqxyw/201412/t20141207_269047.html，2015-01-13.（in Chinese ）

[5]王瑩，蘇永秀，李政.1961-2010年廣西氣溫變化對全球變暖的響應(yīng)[J].自然資源學(xué)報，2013，28（10）:1707-1717. Wang Y，Su Y X，Li Z.Responses of temperature change in Guangxi to global warming during 1961-2010[J].Journal of Natural Resource，2013，28（10）:1707-1717.（ in Chinese ）

[6]蔡福，張玉書，陳鵬獅，等.近50a遼寧熱量資源時空演變特征分析[J].自然資源學(xué)報，2009，24（9）:1635-1646. Cai F，Zhang Y S，Chen P S，et al.Spatio-temporal characteristics of Liaoning thermal resources in recent 50 years[J].Journal of Natural Resources，2009，24（9）:1635-1646. （ in Chinese ）

[7]段春鋒，繆啟龍，曹雯，等.以高山站為背景研究城市化對氣溫變化趨勢的影響[J].大氣科學(xué)，2012，36（4）:811-822. Duan C F，Miao Q L，Cao W，et al.Effect of urbanization on variation trends of air temperatures based on mountain stations[J].Chinese Journal of Atmospheric Sciences，2012，36（4）:811-822.（in Chinese）

[8]Zhang Y J，Kang S C，Qin D H，et al.Seasonal air temperature variations retrieved from a Geladaindong ice core，Tibetan Plateau[J].Journal of Geographical Sciences，2007，17（4）:431-441.

[9]李林，陳曉光，王振宇，等.青藏高原區(qū)域氣候變化及其差異性研究[J].氣候變化研究進(jìn)展，2010，6（3）:181-186. Li L，Chen X G，Wang Z Y，et al.Climate change and its regional differences over the Tibetan Plateau [J].Advantage in Climate Change Research，2010，6（3）:181-186.（in Chinese）

[10]郭渠，孫衛(wèi)國，程炳巖，等.我國西北地區(qū)氣候變化與北極濤動的交叉小波分析[J].南京氣象學(xué)院學(xué)報，2008，31（6）:811-818. Guo Q，Sun W G，Cheng B Y，et al.Cross-Wavelet analysis of the arctic oscillation and climate change in Northwest China[J].Journal of Nanjing Institute of Meteorology，2008，31（6）:811-818.（in Chinese）

[11]施雅風(fēng)，沈永平，李棟梁，等.中國西北氣候由暖干向暖濕轉(zhuǎn)型的特征和趨勢探討[J].第四紀(jì)研究，2003，23（2）:152-164. Shi Y F，Shen Y P，Li D L，et al.Discussion on the present climate change from warm-dry to warm wet in Northwest China[J].Quaternary Sciences，2003，23（2）: 152-164.（in Chinese）

[12]吳昊旻，黃安寧，黃旋旋.近50年長三角地區(qū)季節(jié)的氣候變化特征[J].中國農(nóng)業(yè)氣象，2012，32（3）:317-324. Wu H M，Huang A N，Huang X X.Variation characteristics ofseasons in the Yangtze River delta during the past 50 years[J].Chinese Journal of Agrometeorology，2012，33（3）: 317-324.（in Chinese）

[13]楊再禹，龍先菊.黔東南地區(qū)近49年的氣候變化特征[J].中國農(nóng)業(yè)氣象，2011，32（3）:338-345. Yang Z Y，Long X J.Analysis on characteristics of climate change in last 49 years in Southeastern Guizhou [J].Chinese Journal of Agrometeorology，2011，32（3）:338-345.（in Chinese）

[14]周旗，卞娟娟，鄭景云.秦嶺南北1951-2009年的氣溫與熱量資源變化[J].地理學(xué)報，2011，66（9）:1211-1218. Zhou Q，Bian J J，Zheng J Y.Variation of air temperature and thermal resources in the northern and southern regions of the Qinling Mountains from 1951 to 2009[J].Acta Geographica Sinica，2011，66（9）:1211-1218.（in Chinese）

[15]延軍平，鄭宇.秦嶺南北地區(qū)環(huán)境變化響應(yīng)比較研究[J].地理研究，2001，20（5）:576-582. Yan J P，Zheng Y.A comparative study on environmental change response over the northern and the southern regions of the Qinling Mountains[J].Geographical Research，2001，20（5）:576-582.（in Chinese）

[16]萬佳，延軍平.秦嶺南北生態(tài)地理界線對氣候變化的響應(yīng)[J].干旱區(qū)資源與環(huán)境，2013，27（9）:104-108. Wan J，Yan J P.Interdecadal changes of eco-geographical boundaries in Qinling Mountains[J].Journal of Arid Land Resources and Environment，2013，27（9）:104-108.（in Chinese）

[17]白紅英，馬新萍，高翔，等.基于DEM的秦嶺山地1月氣溫及0℃等溫線變化[J].地理學(xué)報，2012，67（11）:1443-1450. Bai H Y，Ma X P，Gao X，et al.Variations in January temperature and 0℃isothermal curve in Qinling Mountains based on DEM[J].Acta Geographica Sinica，2012，67（11）: 1443-1450.（in Chinese）

[18]宋佃星，延軍平，馬莉.近 50 年來秦嶺南北氣候分異研究[J].干旱區(qū)研究，2011，28（3）: 492-499. Song D X，Yan J P， Ma L.Study on climatic differentiation in the south and north Qinling Mountains in recent 50 years[J].Arid Zone Research，2011，28（3）:492-499.（in Chinese）

[19]孫華，白紅英，張清雨，等.秦嶺南北地區(qū)植被覆蓋對區(qū)域環(huán)境變化的響應(yīng)[J].環(huán)境科學(xué)學(xué)報，2009，29（12）:2635-2641. Sun H，Bai H Y.Zhang Q Y，et al.Responses of vegetation coverage to regional environmental change in the northern and southern Qinling Mountains[J].Acta Scientiae Circumstantiae，2009，29（12）:2635-2641.（in Chinese）

[20]任毅，劉明時，田聯(lián)會，等.太白山自然保護(hù)區(qū)生物多樣性研究與管理[M].北京:中國林業(yè)出版社，2006:1. Ren Y，Liu M S，Tian L H，et al.Biodiversity conservation and management of Taibai nature reserve[M].Beijing:China Forestry Press，2006:1.（in Chinese）

[21]曹潔萍，遲道才，武立強(qiáng)，等.Mann-Kendall檢驗方法在降水趨勢分析中的應(yīng)用研究[J].農(nóng)業(yè)科技與裝備，2008，（5）:35-37，40. Cao J P，Chi D C，Wu L Q，et al.Mann-Kendall examination and applicayion in the analysis of precipitation trend [J].Agriculture Science Technology and Equipment，2008，（5）:35-37，40.（in Chinese）

[22]符淙斌，王強(qiáng).氣候突變的定義和檢測方法[J].大氣科學(xué)，1992，16（4）:482-493. Fu C B，Wang Q.The definition and detection of the abrupt climatic change[J].Scientia Atmospheric Sinica，1992，16（4）: 482-493.（in Chinese）

[23]彭彬，周艷蓮，高蘋，等.氣溫插值中不同空間插值方法的適用性分析:以江蘇省為例[J].地球信息科學(xué)學(xué)報，2011，13（4）: 539-548. Peng B，Zhou Y L，Gao P，et al.Analysis of the use of different spatial interpolation methods in the temperature interpolation taking Jiangsu Province as an example[J]. Journal of Earth Information Science，2011，13（4）:539-548. （in Chinese）

[24]莫申國，張百平.基于DEM的秦嶺溫度場模擬[J].山地學(xué)報，2007，25（4）:406-411. Mo S G，Zhang B P.Simulation of temperature fields based on DEM in Qinling Mountains[J].Journal of Mountain Science， 2007，25（4）:406-411.（in Chinese）

[25]何紅.基于RS和DEM的太白山高山林線對氣溫變化響應(yīng)的研究[D].西安:西北大學(xué)，2014. He H.Research of timberline response to temperature changes in Taibai mountain and based DEM and RS[D].Xi'an: Northwest University，2014.（in Chinese）

[26]李雙雙，延軍平，萬佳.全球氣候變化下秦嶺南北氣溫變化特征[J].地理科學(xué)，2012，32（7）:853-858. Li S S，Yan J P，Wan J.The characteristics of temperature change in Qinling mountains[J].Scientia Geographica Sinica，2012， 32（7）:853-858.（in Chinese）

[27]張雪芹，孫楊，毛煒嶧，等.中國干旱區(qū)氣溫變化對全球變暖的區(qū)域響應(yīng)[J].干旱區(qū)研究，2010，27（4）:592-599. Zhang X Q，Sun Y，Mao W Y，et al.Regional response of temperature change in the arid regions of China to global warming[J].Arid Zone Reserch，2010，27（4）:592-599.（in Chinese）

[28]王勁松，費曉玲，魏鋒.中國西北近50a來氣溫變化特征的進(jìn)一步研究[J].中國沙漠，2008，28（4）:724-732. Wang J S，F(xiàn)ei X L，Wei F，et al.Further study of temperature change in Northwest China in recent 50 years[J].Journal of Desert Research，2008，28（4）:724-732.（in Chinese）

[29]馬新萍.秦嶺林線及其對氣候變化的響應(yīng)[D].西安:西北大學(xué)，2015. Me X P.Qinling mountains timberline and its response to climate change[D].Xi'an:Northwest University，2015.（in Chinese）

[30]白晶，延軍平，蘇坤慧.1958-2007年秦嶺南北氣候變化的差異性分析[J].陜西師范大學(xué)學(xué)報（自然科學(xué)版），2010，38（6）: 98-105. Bai J，Yan J P，Su K H.Differential analysis of abrupt climate change between southern and northern Qinling Mountains in the past 50 years[J].Journal of Shaanxi Normal University （Natural Science Edition），2010，38（6）:98-105.（in Chinese）

[31]丁一匯，張莉.青藏高原與中國其他地區(qū)氣候突變時間的比較[J].大氣科學(xué)，2008，32（4）:794-805. Ding Y H，Zhang L.Inter comparison of the time for climate abrupt change between the Tibetan Plateau and other regions in China[J].Chinese Journal of Atmospheric Sciences，2008， 32（4）:794-805.（in Chinese）

[32]宋燕，季勁鈞.氣候變暖的顯著性檢驗以及溫度場和降水場的時空分布特征[J].氣候與環(huán)境研究，2006，10（2）:157-165. Song Y，Ji J J.The remarkable test of abrupt climatic warming and spatiotemporal distribution features of temperature and precipitation fields[J].Climatic and Environmental Research， 2006，10（2）:157-165.（in Chinese）

[33]陜西省地方志編纂委員會.陜西省植被志[M].西安:西安地圖出版社，2011:488-489. Shaanxi Province Local History Compilation Committee. Vegetation in Shaanxi province[M].Xi'an:Xi'an Map Press，2011:488-489.（in Chinese）

Temperature Variation and Spatial Difference in the Protected Area of Taibai Mountain

LIU Rong-juan, BAI Hong-ying, MA Xin-ping, FENG Lin-lin, QIN Jin, CHENG De-qiang
（Department of Urban and Environment， Northwestern University， Xi'an， Shaanxi 710127， China）

Abstract：Based on the data of 15 meteorological stations from the period of 1959 to 2013 around Taibai Mountain and the data of 11 alpine meteorological stations from the period of 2012 to 2014 inside Taibai Mountain protected area， authors used the Kriging， linear regression and Mann-Kendall mutation test method combined with ArcGIS to analyze the variation trend and mutation characteristics of temperature in the protected area， and contrasted the differences in climate change between southern and northern aspect of the conservation area, high altitude area （altitude≥2800m） and overall the protected area of Taibai Mountain. The results showed: （1）there was an upward trend of annual mean temperature in the protected area of Taibai Mountain from1959 to 2013， and linear trend rate was 0.329℃·10y?1. Spatial differences of temperature variation was very significant， southeastern and western regions showed a significant warming trend （P＜0.05， P＜0.01）， up to 0.383℃·10y?1. Every seasonal mean temperature was increased; winter warming rate was the fastest while spring was the slowest. （2）According to result of M-K test， a warm abrupt change of the annual mean temperature in overall protected area was occurred in 1998. After 1998， temperature variation rate was increasing rapidly. The proportion of the Area’s significant warming trend reached 86.76%. （3）Analysis of spatial difference indicated: annual mean temperature and change rate of northern aspect was lower than southern aspect in the protected area of Taibai Mountain. In latest 55 years， their linear trend rate were 0.247℃·10y?1and 0.291℃·10y?1. Temperature linear trend rate of high altitude area（0.244℃·10y?1） was lower than overall protected area. A warm mutation of the annual mean temperature was appeared in 2001.book=2,ebook=5Temperature of southeastern region raised fast before 2001. Warming rate of western region was declined while central region’s was significantly accelerated after 2001.

Key words：Taibai Mountain; Temperature variation; Spatial difference; Kriging; Temperature change trend rate

作者簡介：劉榮娟（1990-），女，碩士生，主要從事植被變化研究。E-mail：lrj61887@163.com

基金項目：林業(yè)公益性行業(yè)科研專項“秦嶺天然林對氣候變化的時空響應(yīng)及管理對策”（201304309）

* 收稿日期：2015-06-03**通訊作者。E-mail：hongyingbai@163.com

doi:10.3969/j.issn.1000-6362.2016.01.0011