1981～2016年渭河流域上中下游極端氣溫差異研究

張海寧，任源鑫，張新弟，周 旗*

(1.寶雞文理學院 地理與環(huán)境學院，陜西 寶雞 721013；2.陜西省寶雞市馮家山水庫管理局，陜西 寶雞 721300)

政府間氣候變化組織(IPCC)第五次評估報告指出：1880～2012年，全球氣候的平均地表溫度升高了約0.85 ℃；1951～2012年，全球的平均地表溫度升溫速率達到0.12 ℃/10年，約為1880年以來升溫速率的兩倍[1，2]。全球變暖的趨勢愈加嚴重，極端氣候事件的發(fā)生頻率明顯提高，因此極端氣候事件的發(fā)生發(fā)展及其變化規(guī)律備受科學家們的關注。

國內外學者已采用多種方法對極端氣溫事件進行了綜合評價分析。線性回歸算法、Mann-Kendall非參數(shù)檢驗、相關分析法[3，4]等傳統(tǒng)方法在極端事件研究中使用較多。在全球范圍內，Song等[5]對2010年前近30年的世界極端氣溫事件進行分析，發(fā)現(xiàn)全球極端冷暖事件發(fā)生頻率均呈增加趨勢。在國內，眾多學者的研究成果反映出近幾十年國內各流域、地區(qū)極端氣溫事件的發(fā)生呈增加趨勢?；春恿饔騕6,7]、珠江流域[8]、長江流域[9]、中國沿海地區(qū)[10]以及中國西北西南地區(qū)[11-13]冷指數(shù)均呈下降趨勢，暖指數(shù)呈上升趨勢，突變發(fā)生在20世紀80年代至21世紀初期，整體氣溫有所升高，極端氣溫事件發(fā)生頻率有所增加。自1980年，國內相關學者已經開始了渭河流域極端氣溫事件的相關研究。蔡新玲等[14]的研究表明，近50 a陜西省年極端氣溫事件在空間分布上具有較好的一致性。馮星等[15]研究表明，渭河流域近55 a氣溫總體呈小幅增加趨勢，20世紀80年代中期到90年代初期氣溫增長幅度較大。王帥等[16]研究表明，渭河流域各年、季的平均氣溫均表現(xiàn)為顯著上升趨勢，1993年為氣溫突變年。董晴晴等[17]研究發(fā)現(xiàn)，渭河流域關中地區(qū)的平均氣溫呈顯著上升趨勢，冬季對流域平均氣溫的變暖貢獻最大；歷年及各季平均氣溫的突變時間均為20世紀90年代中期。

總體來說，諸多學者偏重于渭河流域極端氣溫事件整體趨勢格局的相關研究，少有學者對渭河流域上、中、下游極端氣溫事件進行差異評價。

本文以渭河流域為研究區(qū)域，定義林家村站之上為渭河上游流域，林家村站至咸陽站之間為渭河中游流域，咸陽站之下為渭河下游流域。借助氣候變化監(jiān)測和指數(shù)專家組 ETCDI(Expect Team on Climate Change Detection and Indices)定義的核心指數(shù)，結合渭河流域的氣候特點，選取可以表征渭河流域極端氣溫事件的12個指數(shù)；基于R ClimDex軟件[18]計算了1981～2016年渭河流域66個氣象站點的氣溫數(shù)據并建立了極端氣溫指數(shù)時間序列。首先運用線性回歸法、Mann-Kendall突變檢驗法對氣溫指數(shù)時間序列進行處理，同時結合累計距平法分析渭河流域上、中、下游36 a氣溫趨勢差異和突變情況。其次結合任國玉等[19]對極端氣溫指數(shù)的分類，對12個極端氣溫指數(shù)進行顯著性檢驗，篩選出能夠表征“相對閾值極端高低溫(日數(shù))、持續(xù)性高低溫事件、年較差”的極端氣溫指數(shù)，對其進行相關分析和R/S分析，研究渭河上、中、下游流域未來的氣溫變化趨勢，為研究區(qū)未來極端氣溫事件的研究提供科學預測和理論參考。

1 研究區(qū)概況

渭河是黃河的第一大支流，其干流全長約818 km，流域覆蓋總面積約13.5萬km2。發(fā)源地位于甘肅省定西市渭源縣，流經甘肅省天水市、陜西省關中平原(主要包括寶雞、咸陽、西安、渭南)等地區(qū)，并在渭南市的潼關縣匯入黃河。渭河流域位于大陸性季風氣候區(qū)，并處于干旱、濕潤地區(qū)的過渡交錯地帶；主要受西太平洋副熱帶高壓的影響，年降水較少，同時蒸發(fā)量較大。渭河流域呈現(xiàn)西高東低的地勢：其中渭河流域的西部為黃土丘陵溝壑區(qū)，東部為關中平原區(qū)，南部為秦嶺山脈，北部為黃土高原地區(qū)。渭河流域春季氣溫不穩(wěn)定，夏季炎熱多雨，秋季涼爽濕潤，冬季寒冷且干燥。流域內多年平均氣溫在7.8～13.5 ℃。極端最高氣溫為42.8 ℃，極端最低氣溫約為-28.1 ℃。

2 研究內容及方法

2.1 數(shù)據來源與處理

為研究渭河流域極端氣溫事件的變化規(guī)律，本文選取渭河流域數(shù)據較為完整、空間分布較為均勻的66個氣象站點(圖1)的逐日氣溫資料。同時對資料進行質量控制，對短期缺測數(shù)據進行插補。其中采用的逐日氣溫資料主要來源于陜西、甘肅兩省氣候中心，部分數(shù)據通過國家氣候數(shù)據共享網(https://data.cma.cn)進行補充，氣溫資料時段取在1981年1月1日～2016年12月31日。

圖1 渭河流域及主要站點示意圖

2.2 研究方法

本文選取 ETCDI定義的12個極端氣溫指數(shù)來分析渭河流域極端氣溫事件，所選指標定義及表示方式見表1；同時采用線性趨勢法、Mann-Kendall突變檢驗和R/S分析等方法對數(shù)據進行處理分析。其中Mann-Kendall突變檢驗[7,19]是分析氣象數(shù)據時間序列趨勢的常用方法，能夠很好地體現(xiàn)氣象數(shù)據的變化趨勢和突變點。R/S分析法又稱為重標極差分析法(Rescaled Range Analysis)，其中Hurst指數(shù)能很好地揭示時間序列中的趨勢成分，證明未來變化的總體趨勢與過去的關系。

表1 極端氣候事件類型及其定義

3 結果與分析

3.1 1981～2016年渭河流域極端氣溫指數(shù)變化趨勢

1981～2016年渭河流域上中下游地區(qū)極端氣溫指數(shù)的變化趨勢見圖2、表2。冷持續(xù)指數(shù)(CSDI)在渭河流域上、下游呈下降趨勢，下降速率分別為0.079 d/10 a、0.12 d/10 a，中游以0.103 d/10 a的速率上升，全流域以0.063 d/10 a的速率下降，但變化幅度均較小，歷年變化較平穩(wěn)；結合圖2a、表2，2008年冷持續(xù)指數(shù)在流域內的變化均出現(xiàn)了極大值。暖持續(xù)指數(shù)(WSDI)在渭河流域內均呈波動上升趨勢，上升速率分別為2.579 d/10 a、2.116 d/10 a、2.052 d/10 a及2.133 d/10 a，上升速率從上游至下游逐漸減小，21世紀初期暖持續(xù)指數(shù)上升速率較快；結合圖2b、表2，1998年后該指數(shù)上升較顯著，表明36 a來渭河流域最高氣溫大于90%分位值的日數(shù)逐漸增加，渭河流域內呈變暖趨勢，且變暖速率從上游至下游逐漸減慢。

生物生長季(GSL)在渭河流域內均呈上升趨勢，上升速率分別為8.109 d/10 a、7.401 d/10 a、7.728 d/10 a和7.721 d/10 a，結合圖2c、表2，21世紀初期變化速率明顯加快，2005年后該指數(shù)上升較為顯著。年氣溫日較差(DTR)在上游以0.162 ℃/10 a的速率下降，中、下游及全流域分別以0.178 ℃/10 a、0.199 ℃/10 a、0.207 ℃/10 a的速率上升，下游的變化速率大于上游；表明渭河流域年氣溫日較差從上游至下游逐漸變大，下游的年氣溫變化波動更大。

熱夜日數(shù)(TR20)在渭河流域內均呈上升趨勢，中游上升速率最大，為5.784 d/10 a，其次是下游和上游，全流域與各區(qū)域的變化趨勢一致。霜凍日數(shù)(FD0)在渭河流域內均呈下降趨勢，中游地區(qū)下降速率最大，為6.121 d/10 a，其次是下游和上游，全流域霜凍日數(shù)的變化趨勢與各區(qū)域變化趨勢相符。冰凍日數(shù)(ID0)在渭河流域內均呈下降趨勢，上游地區(qū)下降速率最大，表明渭河流域日最高氣溫小于0 ℃的日數(shù)在36 a內逐漸減少，即每年的日最高氣溫整體上在升高，整個流域呈現(xiàn)出變暖的趨勢。

夏日日數(shù)(SU25)在上、中、下游均呈上升趨勢，上升速率從上游至下游逐漸減小，上游的上升速率最快，約每10年增加8 d，表明渭河流域各地區(qū)每年日最高氣溫大于25 ℃的天數(shù)呈現(xiàn)逐年增多的趨勢，上游地區(qū)變暖趨勢更加明顯。

冷晝日數(shù)(TX10p)在渭河流域內均呈下降趨勢；暖晝日數(shù)(TX90p)在渭河流域內均呈上升趨勢，結合圖2i、圖2j及表2，冷晝日數(shù)和暖晝日數(shù)的變化速率均為從上游至下游逐漸減小，表明上游地區(qū)極端高溫天氣情況變化更為顯著。

圖2 1981～2016年渭河流域極端氣溫指數(shù)變化趨勢

冷夜日數(shù)(TN10p)在渭河流域上、中、下游均呈下降趨勢，暖夜日數(shù)(TN90p)在上、中、下游均呈上升趨勢。結合表2得知，兩個指數(shù)中游地區(qū)的變化速率最大，其次是上游地區(qū)；全流域內冷夜日數(shù)以1.629 d/10 a的速率呈下降趨勢，暖夜日數(shù)以3.363 d/10 a的速率呈上升趨勢，暖夜日數(shù)的變化速率約是冷夜日數(shù)的2倍。

表2 渭河流域極端氣溫指數(shù)線性趨勢

渭河流域極端冷指標的整體變化呈下降趨勢，而極端暖指標的整體變化呈上升趨勢；冷指標的變化速率為中游最快，暖指標的變化速率為上游最快，其中暖持續(xù)指數(shù)上升速率從上游至下游逐漸減小，上游地區(qū)連續(xù)高溫的日數(shù)逐漸增多；說明渭河流域整體上表現(xiàn)為增暖的趨勢，與謝志祥等[6,9]的研究結果相吻合。

3.2 1981～2016年渭河流域極端氣溫指數(shù)突變分析

突變分析的目的在于找到突變點所對應的時間，因此首先采用M-K檢驗對1981～2016年渭河流域極端氣溫指數(shù)進行突變檢驗，同時結合累計距平法分析其變化趨勢，確定其突變點及突變年份，突變檢驗結果見表3。

表3 渭河流域極端氣溫指數(shù)的突變年份

1981～2016年，冷持續(xù)指數(shù)在1997年左右發(fā)生了由多到少的突變，在2008年左右發(fā)生了由少到多的突變，表明在1997年連續(xù)低溫日數(shù)短時間快速增多。霜凍日數(shù)在1999年左右發(fā)生了由少到多的突變，表明在20世紀末日最低溫度低于0 ℃的天氣日數(shù)突然增多；冰凍日數(shù)在1993年左右發(fā)生了由多到少的突變，表明日最高溫度小于0 ℃的日數(shù)突然減少。冷夜日數(shù)在1997年左右發(fā)生了由多到少的突變，冷晝日數(shù)在1997年左右發(fā)生了由少到多的突變。綜上所述，20世紀90年代末期為渭河流域冷指標的突變時期。

暖持續(xù)指數(shù)在1997年左右發(fā)生了由多到少的突變，熱夜日數(shù)在1993年左右發(fā)生了由少到多的突變，表明此階段日最低氣溫大于20 ℃的日數(shù)突然增加。暖晝日數(shù)在1994年左右發(fā)生了由多到少的突變，暖夜日數(shù)在1997年左右發(fā)生了由少到多的突變。夏日日數(shù)在1993年左右發(fā)生了由少到多的突變，表明此階段日最高氣溫大于25 ℃日數(shù)的上升速率顯著加快，高溫日數(shù)增多。綜上所述，1993～1997年為渭河流域暖指標的突變時期。

綜合得知，1981～2016年渭河流域極端氣溫指數(shù)的突變時期為20世紀90年代，尤其在1993年和1997年突變情況最為明顯；暖指標的突變時期略早于冷指標。

3.3 渭河流域全流域與上、中、下游氣溫的相關分析情況

對12個指數(shù)進行顯著性檢驗，同時結合任國玉等[19]對極端氣溫指數(shù)的分類，篩選出顯著性高且能夠表征“相對閾值極端高低溫(日數(shù))、持續(xù)性高低溫事件、年較差”的6個極端氣溫指數(shù)。其中相對閾值極端高低溫(日數(shù))選取冷晝日數(shù)(TX10p)和暖夜日數(shù)(TN90p)；持續(xù)性高低溫事件選取冷持續(xù)指數(shù)(CSDI)、暖持續(xù)指數(shù)(WSDI)和霜凍日數(shù)(FD0)；年較差選取氣溫日較差(DTR)。

對上、中、下游的6個指數(shù)與渭河全流域進行相關分析，結果見表4。結果表明渭河流域上中下游地區(qū)的相關分析結果均接近1，說明上中下游對渭河流域極端氣溫指數(shù)變化的影響均較大；尤其是下游地區(qū)冷晝日數(shù)(TX10p)、暖夜日數(shù)(TN90p)、氣溫日較差(DTR)的相關分析結果更大，說明下游地區(qū)與渭河流域的極端氣溫指數(shù)相關性更高，對渭河流域極端氣溫指數(shù)影響最大的地區(qū)是下游地區(qū)。

表4 渭河流域極端氣溫指數(shù)相關分析結果

3.4 1981～2016年渭河流域極端氣溫指數(shù)R/S分析

為研究渭河流域上中下游極端氣溫指數(shù)的未來變化趨勢，對顯著性高的6個極端氣溫指數(shù)進行R/S分析，計算極端氣溫指數(shù)的Hurst指數(shù)。Hurst指數(shù)的大小反映了變化趨勢的持續(xù)性(反持續(xù)性)的強弱[20]。當H=0.5時，表示氣候序列是獨立的隨機序列，現(xiàn)在不會影響未來。當0.5

由表5可知，渭河流域上中下游極端氣溫指數(shù)具有Hurst現(xiàn)象，且Hurst指數(shù)均大于0.5，表明未來變化趨勢與1981～2006年間的變化趨勢一致，且持續(xù)性強度較強。過去36 a內，冷指標為下降趨勢且Hurst指數(shù)均大于0.5，表明未來冷指標將繼續(xù)呈下降趨勢；暖指標、氣溫日較差和生物生長季為下降趨勢，且Hurst指數(shù)均大于0.5，表明未來暖指標將繼續(xù)呈上升趨勢。

表5 渭河流域極端氣溫指數(shù)的Hurst指數(shù)

整體來看，極端氣溫指數(shù)的未來總體趨勢具有強持續(xù)性，未來渭河流域變暖趨勢的可能性較大。上游地區(qū)極端氣溫指數(shù)的持續(xù)性程度較中、下游更強，變暖趨勢更明顯。

4 結論

通過對12個極端氣溫指數(shù)的研究，得出關于渭河流域極端氣溫事件的主要結論如下：1981～2016年間，渭河流域暖指標(暖持續(xù)指數(shù)，熱夜日數(shù)，夏日日數(shù)，暖晝日數(shù)，暖夜日數(shù))呈波動上升趨勢變化，變化速率從上游到下游逐漸減小，上游地區(qū)連續(xù)高溫和極端高溫出現(xiàn)的日數(shù)逐漸增多；冷指標(冷持續(xù)指數(shù)，霜凍指數(shù)，冰凍日數(shù)，冷晝日數(shù)，冷夜日數(shù))呈顯著下降趨勢；說明過去36 a渭河流域呈變暖趨勢，上游變暖趨勢最明顯，上游極端氣溫事件的發(fā)生頻率大于下游。

渭河流域極端氣溫指數(shù)的突變時期集中在20世紀90年代，尤其在1993年和1997年突變情況最為明顯；暖指標的突變時間略早于冷指標。對渭河流域極端氣溫變化影響最大的區(qū)域為渭河下游區(qū)域。冷指標在未來將繼續(xù)呈下降趨勢，暖指標在未來將繼續(xù)呈上升趨勢；上游地區(qū)極端氣溫指數(shù)的持續(xù)性程度較中、下游更強，未來變暖的趨勢更明顯。