氣候變化對黃河中上游地區(qū)潛在蒸散影響的估算

王亞俊，李俊，林忠輝，同小娟，邢魯敏

(1.中國科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室，100101，北京;2.中國科學(xué)院大學(xué)，100049，北京;3.北京林業(yè)大學(xué)林學(xué)院，100083，北京;4.禹城市建筑工程公司，251200，山東禹城)

蒸散是陸地水循環(huán)和能量平衡的一個重要過程，其中潛在蒸散表征一個地區(qū)的最大可能蒸散量，在區(qū)域水量平衡、農(nóng)業(yè)需水和氣候變化影響等研究中具有重要意義。黃河作為華北和西北地區(qū)重要供水水源，承擔(dān)著農(nóng)業(yè)灌溉用水、居民和工礦用水任務(wù)。研究表明，包括黃河流域在內(nèi)的中國地區(qū)蒸散和潛在蒸散都呈現(xiàn)減少趨勢［1-3］，對潛在蒸散變化貢獻(xiàn)最大的有最高氣溫［3］、輻射(日照時(shí)間)［4-5］、風(fēng)速［6］和相對濕度［7］等。不同區(qū)域、不同學(xué)者、不同方法得到的研究結(jié)果不盡相同。以往研究多采用相關(guān)分析和多元回歸方法［8-11］分析氣候因子對潛在蒸散的影響，只得到統(tǒng)計(jì)關(guān)系，缺乏機(jī)制基礎(chǔ)。筆者采用偏微分方法［12-13］研究氣候因子變化對潛在蒸散變化的影響，著重考慮變化機(jī)制，計(jì)算結(jié)果明確、直觀、可信。筆者結(jié)合黃河中上游地區(qū)近50年氣候變化的趨勢，采用偏微分方法分析潛在蒸散對氣候因子(實(shí)際水汽壓 ea、日照時(shí)間 SH、風(fēng)速 v、最高氣溫tx、最低氣溫tm)變化的敏感性和氣候因子對潛在蒸散變化的貢獻(xiàn)率，揭示不同時(shí)間尺度上影響黃河中上游地區(qū)潛在蒸散變化的主導(dǎo)氣候要素，為流域水資源管理、區(qū)域可持續(xù)發(fā)展和氣候變化適應(yīng)對策研究提供科學(xué)依據(jù)。

1 研究區(qū)概況

黃河發(fā)源于青海省巴顏喀拉山北麓海拔4 500 m的約古宗列盆地，地勢自西向東呈青海高原、內(nèi)蒙古高原和黃土高原、魯中南低山丘陵區(qū)三級階梯逐級下降。黃河流域(E 96°～119°，N 32°～42°)氣候類型多樣，平均氣溫在-4～14℃之間，年降水量平均為478 mm，降水量從東南到西北遞減，分屬半濕潤、半干旱、干旱地區(qū)。

2 資料與方法

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

本研究中所用氣象數(shù)據(jù)來自中國氣象局國家氣象信息中心，選取黃河中上游地區(qū)62個氣象站點(diǎn)(圖1)，數(shù)據(jù)包括日最高、最低和平均氣溫，風(fēng)速，相對濕度和日照時(shí)間，月值和年值均由日值統(tǒng)計(jì)得出。

2.2 變化趨勢檢驗(yàn)及突變點(diǎn)檢驗(yàn)

用線性擬合的斜率表示氣候要素的變化趨勢，正值表示增加趨勢，負(fù)值表示減小趨勢。變化趨勢的信度檢驗(yàn)采用Mann-Kendall(M-K)趨勢檢測法，該方法適用于水文、氣象等非正態(tài)分布的數(shù)據(jù)，計(jì)算簡便，廣泛應(yīng)用于時(shí)間序列趨勢的非參數(shù)檢驗(yàn)［14］。M-K公式中|Z|≥1.96，即置信水平達(dá)0.05的變化趨勢為顯著(a＜0.05)，其余為變化趨勢不顯著。

氣候突變是氣候變化過程中存在的某種不連續(xù)現(xiàn)象，氣象序列突變點(diǎn)檢驗(yàn)方法眾多，本文選擇應(yīng)用廣泛的M-K突變點(diǎn)檢驗(yàn)方法來檢測氣象要素過去50年時(shí)間序列可能的突變點(diǎn)。

2.3 潛在蒸散

潛在蒸散(E0)是指下墊面足夠濕潤條件(包括了地上植被含水量和土壤含水量均達(dá)到飽和以上的情況)下供應(yīng)的最大蒸發(fā)量。本文采用1998年聯(lián)合國糧農(nóng)組織(FAO)改進(jìn)的Penman-Monteith(P-M)公式，該公式綜合了空氣動力學(xué)的湍流傳輸與能量平衡，考慮了植被的生理特征，在干旱和濕潤條件下準(zhǔn)確性都相對較高，已在世界范圍內(nèi)得到廣泛應(yīng)用［15］。其中參考假想面類似于同一高度、生長旺盛、完全覆蓋地面、水分充足的廣闊綠色植被(高0.12 m，表面阻力為70 s/m，反射率為0.23)［16］。計(jì)算公式為

式中:E0為潛在蒸散量，mm;Rn為凈輻射，MJ/(m2·d);G為土壤熱通量，MJ/(m2·d);t為2 m高處日平均氣溫，℃;v2為2 m高處風(fēng)速，m/s;es為飽和水汽壓，kPa;ea為實(shí)際水汽壓，kPa;Δ為飽和水汽壓曲線斜率，kPa/℃;γ 為干濕表常數(shù)，kPa/℃。Rn、G、Δ、v2可通過氣象臺站觀測資料計(jì)算求得，土壤熱通量G很小，不予考慮。

圖1 黃河中上游地區(qū)氣象站點(diǎn)分布Fig．1 Distribution of meteorological stations in the upper-middle reach of the Yellow River

式(1)中，Rn為太陽輻射 Rs的函數(shù)，即 Rn=f(Rs)。Rs可由下式計(jì)算:

1．2 實(shí)驗(yàn)方法和程序 給全體被試在一安靜房間內(nèi)依次施測區(qū)分愿望任務(wù)、區(qū)分信念任務(wù)、意外內(nèi)容、意外地點(diǎn)任務(wù)。區(qū)分愿望任務(wù)測試了兒童對愿望的理解，總分1分；區(qū)分信念任務(wù)測試兒童是否能區(qū)分自己和他人的信念，總分1分；后兩個任務(wù)則測試了兒童對錯誤信念的理解，總分3分。所有的任務(wù)均用實(shí)物模型演示。

式中:Rs為太陽總輻射，MJ/(m2·d);SH為日照時(shí)間，h;N為最大可能日照時(shí)間，h;Ra為天文輻射，MJ/(m2·d);a和b分別為經(jīng)驗(yàn)系數(shù)，中國地區(qū)一般取a=0.25，b=0.5。本文采用黃河流域內(nèi)外輻射觀測站的數(shù)據(jù)，通過最小二乘法和ArcGIS反距離加權(quán)插值方法，確定整個流域內(nèi)氣象站點(diǎn)的a和b值［3］，采用修正后的 a、b值計(jì)算潛在蒸散。

2.4 濕潤指數(shù)

濕潤指數(shù)(Humidity index，HI)即表征一個地區(qū)干濕程度的指標(biāo)，在地理學(xué)和生態(tài)學(xué)研究中長期應(yīng)用，特別是氣候變化和干旱化等研究中使用的更為廣泛［17］，即

式中:HI為濕潤指數(shù);P為降水量，mm。本文采用王菱等［18］提出的劃分標(biāo)準(zhǔn)來研究黃河流域干濕界限變化，即:HI＜0.03為極端干旱，HI=0.03～0.20為干旱氣候帶，HI=0.20～0.50為半干旱氣候帶，HI=0.05～0.75為半濕潤氣候帶，HI=0.75～1.0為比較濕潤氣候帶，HI＞1.0為濕潤氣候帶。

2.5 敏感性分析

已有很多研究使用偏微分方法檢驗(yàn)?zāi)Ｐ妥宰兞繉σ蜃兞康拿舾行?，本文敏感性分析有助于定量解釋潛在蒸散變化的原因，且已證明其結(jié)果是有效的［19］。其定義為

式中:Sx為潛在蒸散關(guān)于氣侯因子x的敏感系數(shù)，量綱為1。各個氣候變量對潛在蒸散的敏感性計(jì)算方法已經(jīng)得到了廣泛的應(yīng)用［11-13］。使用敏感系數(shù)的優(yōu)點(diǎn)在于量綱為1，方便對不同量綱的影響變量進(jìn)行排序評價(jià)［13］。

2.6 貢獻(xiàn)率評價(jià)

貢獻(xiàn)率評價(jià)能夠定量地確定氣候因子變化在導(dǎo)致潛在蒸散量變化中所占的比例，可通過求偏導(dǎo)數(shù)方法定量確定。假設(shè)各氣候因子相互獨(dú)立，相關(guān)性小，則計(jì)算潛在蒸散的全微分表達(dá)式［12］為

3 結(jié)果與分析

3.1 氣候要素變化特征

黃河中上游地區(qū)相對濕度(RH)總體下降，年變化幅度為-0.048%，但未達(dá)顯著水平;實(shí)際水汽壓(ea)總體上升幅度為0.0005 h Pa/a(a＜0.01)。流域相對濕度春季和秋季變化幅度最大，實(shí)際水汽壓流域春季有微弱的減少趨勢，冬季和夏季增加幅度最大。相對濕度M-K圖顯示突變點(diǎn)發(fā)生在2004年，實(shí)際水汽壓突變點(diǎn)較氣溫提前，發(fā)生在20世紀(jì)80年代中期。

黃河中上游地區(qū)日照時(shí)間(SH)總體下降，變化幅度為-0.009 h/a(a＜0.01)。日照時(shí)間春季略有增加，其他季節(jié)減少，夏季和冬季減少趨勢更明顯，與韓虹等［21］研究結(jié)果相似(表1)。日照時(shí)間在20世紀(jì)70年代后期發(fā)生突變，與中國改革開放時(shí)期相符，這一期間工業(yè)活動大力開展，相應(yīng)的大氣污染排放物增多，氣溶膠含量增加，大氣透明度降低，到達(dá)地表的輻射量減少。

表1 黃河中上游地區(qū)氣候因子年變化率(1961—2010年)Tab．1 Annual change rates of climate factors in the upper-middle reach of the Yellow River(1961—2010)

全球氣溫升高導(dǎo)致大洋環(huán)流、大氣與地表能量交換等發(fā)生變化，進(jìn)而影響各地的風(fēng)速、最大風(fēng)速時(shí)間及颶風(fēng)頻率。大量研究［14，22-23］表明中國大部分地區(qū)風(fēng)速減少，本研究結(jié)果顯示黃河中上游地區(qū)風(fēng)速總體下降，年變化幅度為-0.006 m/s(a＜0.01)，略低于全國-0.012 m/s［24］的年變化率。風(fēng)速在四季都呈顯著減少趨勢，冬季減少程度最大(表1)。風(fēng)速顯著減少可能與冬夏季風(fēng)的不斷減弱有關(guān)［16］，也可能是大氣環(huán)流減弱、城市化等因素引起的［25］。流域年平均風(fēng)速在20世紀(jì)80年代初發(fā)生突變，1980年之前風(fēng)速有微弱的增加趨勢，之后風(fēng)速逐年減少。

黃河中上游地區(qū)降水量(P)總體下降，幅度為-0.119 mm/a(表1)，與徐宗學(xué)等［26］得到的黃河流域的年均變化率-0.129 mm/a相比偏少。秋季降水量顯著減少，變化率達(dá)-0.287 mm/a(a＜0.01)，對年降水量減少貢獻(xiàn)較大，冬季降水量有略微增加趨勢，但不顯著(表1)。流域降水量突變點(diǎn)發(fā)生在20世紀(jì)80年代中期，與實(shí)際水汽壓突變時(shí)間相同，較風(fēng)速和日照時(shí)間的突變時(shí)間有所滯后。

黃河中上游地區(qū)多年平均潛在蒸散量(E0)年變化為0.079 mm/a，變化不顯著，與他人研究結(jié)果［1-2］有出入，可能是未考慮黃河下游的原因。春季和秋季，黃河中上游地區(qū)潛在蒸散呈增加趨勢，夏季和冬季則呈減少趨勢(表1)。潛在蒸散的變化是由氣候因子共同作用的結(jié)果。王幼奇等［27］認(rèn)為自1957年以來，黃河流域潛在蒸散量總體上呈下降—穩(wěn)定—上升的趨勢，結(jié)合M-K檢驗(yàn)結(jié)果得出，潛在蒸散在20世紀(jì)80年代初發(fā)生突變，呈先下降后有微弱上升的趨勢。

3.2 干濕界限變化

圖2 黃河中上游地區(qū)氣候因子年值M-K突變點(diǎn)檢驗(yàn)圖Fig．2 M-K break point test graph of climate factors'annual value in the upper-middle reach of the Yellow River

如圖3所示，黃河流域處于半干旱半濕潤氣候帶，西北小部分地區(qū)處于干旱氣候帶。圖3分別以20世紀(jì)60年代(1961—1970年)為基礎(chǔ)，比較每10年氣候干濕帶的變化。70年代流域中游區(qū)0.5線(即HI=0.5的線，下同)和0.75線向東南方向移動，且60年代北部干旱氣候區(qū)域消失，半干旱氣候區(qū)擴(kuò)大，比較濕潤區(qū)減小，整體有暖干的變化趨勢。而在黃河源區(qū)70年代0.5線消失，0.75線向西北移動，比較濕潤區(qū)有擴(kuò)大(圖3(a));80年代相較60年代的黃河中游地區(qū)0.2、0.5線向東南方移動，0.75線變化不明顯，流域南部出現(xiàn)了1.0線，干旱氣候區(qū)向南有所延伸，半濕潤氣候區(qū)和比較濕潤區(qū)縮小，但南部濕潤氣候區(qū)向北擴(kuò)展，流域北部干旱化，南部有變濕趨勢。黃河源區(qū)0.5、0.75線都向北部移動，出現(xiàn)了1.0濕潤區(qū)，呈現(xiàn)濕潤的變化趨勢(圖3(b))。90年代黃河流域干濕界限變化特點(diǎn)與70年代相似，0.5和0.75線向南移動程度加大，干旱區(qū)南移(圖3(c));21世紀(jì)初10年黃河中游干濕界限的變化也同70年代相似，0.5線向南移動較大，干旱區(qū)范圍擴(kuò)大，0.75線較60年代向南偏移，濕潤區(qū)范圍縮小。黃河源區(qū)0.5線消失，0.75線北移幅度小于南移幅度移動，北部有變濕潤的跡象，南部趨于變干(圖3(d))。

從以上結(jié)果可以看出，近50年流域半干旱區(qū)和半濕潤區(qū)的分界線(0.5線)不斷向東南推進(jìn)，半干旱區(qū)面積不斷擴(kuò)大，氣候趨向干旱化，與王菱等［18］研究結(jié)果基本一致。黃河源區(qū)北部可能是降水變率大于潛在蒸散變率，濕潤指數(shù)增加，氣候持續(xù)變濕，而源區(qū)南部20世紀(jì)70年代、80年代、90年代有變濕趨勢，但2010年與60年代相比變干。M-K突變點(diǎn)顯示在70年代和21世紀(jì)初，且與日照時(shí)間M-K圖走勢一致，與日照時(shí)間相關(guān)性大(R2=0.65，P＜0.01)。從圖3看出，70年代相較60年代的干濕界限的確發(fā)生了較大的變化，干旱化明顯，同楊建平等［28］結(jié)果相似，21世紀(jì)初沒有顯著突變情況。東南季風(fēng)、西南季風(fēng)、高原季風(fēng)、中緯度西風(fēng)環(huán)流、西太平洋副熱帶高壓的年代際變化是過去中國干濕氣候界線波動與氣候干濕變化的根本原因。黃河流域中上游局地而言，濕潤指數(shù)降低，氣候干旱程度增加，主要與降水的減少有關(guān)。

圖3 黃河中上游地區(qū)氣候干濕界限年代際變化Fig．3 Decadal change for the borderline of wet and dry climate zone in the upper-middle reach of the Yellow River

3.3 敏感性分析

3.3.1 全年和各季敏感系數(shù) 潛在蒸散對氣候因子的敏感性隨著時(shí)間和地區(qū)發(fā)生變化，黃河流域中上游四季及全年潛在蒸散對氣候因子的敏感系數(shù)(式4)見表2。年尺度上，各氣候因子相對獨(dú)立的情況下，日照時(shí)間、風(fēng)速、日最高氣溫、日最低氣溫增加10%，潛在蒸散分別增加 1.99%、2.01%、3.93%、0.29%，實(shí)際水汽壓增加10%，潛在蒸散減少6.37%。這表明年潛在蒸散對實(shí)際水汽壓最敏感，其次為日最高氣溫、風(fēng)速、日照時(shí)間，對日最低氣溫度最不敏感，與江瀅等［24］的結(jié)果一致。除少數(shù)站點(diǎn)(青海的貴德站、甘肅的靖遠(yuǎn)站、山西的河曲站，四川的馬爾康站)潛在蒸散對日最高氣溫最敏感外，多數(shù)站點(diǎn)都對實(shí)際水汽壓最敏感。黃河中上游地區(qū)降水變率大，潛在蒸散對實(shí)際水汽壓最敏感，表明水分是限制該地區(qū)農(nóng)業(yè)生產(chǎn)的主要因素［29］。

季節(jié)尺度上，日照時(shí)間和日最高氣溫的敏感系數(shù)都呈現(xiàn)夏季最高，冬季最低，與其本身因子季節(jié)變化特征相似。風(fēng)速敏感系數(shù)是夏季最低，冬季最高，與冬季風(fēng)速變化率最高一致。實(shí)際水汽壓敏感系數(shù)絕對值在秋冬季高，春夏季低。日最低氣溫的敏感系數(shù)夏季最高，冬季為負(fù)值，冬季日最低氣溫的變化方向與潛在蒸散變化方向相反。春季、秋季和冬季實(shí)際水汽壓敏感系數(shù)絕對值仍然是最大值，即春、秋和冬季潛在蒸散對實(shí)際水汽壓最敏感。夏季，敏感系數(shù)最大值是日最高氣溫，這與楊小利等［16］的研究結(jié)果一致。

表2 黃河中上游地區(qū)潛在蒸散對氣候因子變化敏感系數(shù)(1961—2010年)Tab．2 Sensitivity coefficient of climate factors to potential evapotranspiration in the upper-middle reach of the Yellow River(1961—2010)

不同氣候因子對潛在蒸散的敏感程度是隨著時(shí)間發(fā)生變化(表3)。過去50年，潛在蒸散對日照時(shí)間變化的敏感系數(shù)(S(SH))呈顯著下降趨勢(a＜0.01)，對風(fēng)速、實(shí)際水汽壓和最低溫度變化的敏感系數(shù)(分別為S(v)、S(ea)和S(tn))都呈顯著上升趨勢(a＜0.01)，S(tx)有上升趨勢，但變化不顯著。在各氣候因子中，實(shí)際水汽壓(ea)的敏感系數(shù)變化率最大。S(ea)為負(fù)值，其絕對值呈下降趨勢，即潛在蒸散對實(shí)際水汽壓的敏感性減少。

表3 黃河中上游地區(qū)潛在蒸散對氣候因子變化敏感系數(shù)的年變化率(1961—2010年)Tab．3 Annual change rate of climate factors'sensitive coefficient to potential evapotranspiration in upper-middle reach of Yellow River(1961—2010)

3.3.2 敏感系數(shù)與高程關(guān)系 綜合流域內(nèi)氣象站點(diǎn)高程資料，發(fā)現(xiàn)潛在蒸散對敏感系數(shù)與高程有一定關(guān)系。黃河中上游地區(qū)地勢呈3級階梯，高程范圍從100 m到4 000多m，流域內(nèi)高程的變化大，跨越多種氣候區(qū)。潛在蒸散對日照時(shí)間敏感系數(shù)隨著高程增加而增加(P＜0.05)(圖4(a))，高程越高太陽輻射越強(qiáng)，這與潛在蒸散對日照時(shí)間敏感系數(shù)夏季時(shí)最高一致。潛在蒸散對風(fēng)速、最高氣溫和最低氣溫敏感系數(shù)隨著高程增加而減少(P＜0.05)(圖4(a)和(c))，潛在蒸散對氣溫的敏感系數(shù)與高程相關(guān)性最好(圖4(c))。在高海拔地區(qū)，最低氣溫的敏感系數(shù)為較大的負(fù)值，可能是由于氣溫極低時(shí)潛在蒸散對氣溫偏導(dǎo)數(shù)值小，而氣溫絕對值較大，導(dǎo)致敏感系數(shù)數(shù)值較大［30］。潛在蒸散對實(shí)際水汽壓敏感系數(shù)與高程關(guān)系不顯著(圖4(b))。以上分析結(jié)果與劉昌明等［3］的結(jié)果一致。潛在蒸散發(fā)對氣象要素的敏感系數(shù)與高程的線性關(guān)系反映了不同海拔的地形、植被、土壤和大氣之間復(fù)雜的反饋機(jī)制。

圖4 黃河中上游地區(qū)潛在蒸散對氣候因子的敏感系數(shù)與高程的關(guān)系Fig．4 Correlations between the sensitivities and altitudes in the upper-middle reach of the Yellow River

3.4 貢獻(xiàn)率分析

各個氣候因子對潛在蒸散變化的貢獻(xiàn)率可通過潛在蒸散對各氣候因子的偏導(dǎo)數(shù)和各氣候因子的變化斜率估算(式5)得到。結(jié)果表明，區(qū)域內(nèi)各站點(diǎn)的氣候因子對潛在蒸散變化貢獻(xiàn)率之和與潛在蒸散變化率吻合較好，春、夏、秋、冬和全年線性擬合R2分別為0.88、0.96、0.85、0.85 和 0.91，都達(dá)到顯著相關(guān)(P＜0.005)(圖5)。這表明采用式(5)定量計(jì)算該流域氣候因子對潛在蒸散貢獻(xiàn)率的方法是可行的。

圖5 黃河中上游地區(qū)各氣候因子對潛在蒸散變化的貢獻(xiàn)率總和與潛在蒸散變化率比較Fig．5 Comparison between contribution of climate factors to E0change and change rate of E0

年尺度上，風(fēng)速對潛在蒸散變化貢獻(xiàn)率最大，之后依次為日最高氣溫、日氣溫最低、實(shí)際水汽壓，貢獻(xiàn)率最小的是日照時(shí)間(表4)。除了日最高氣溫、日最低氣溫變化對潛在蒸散的影響為正效應(yīng)，其余氣候因子變化對潛在蒸散的影響為負(fù)效應(yīng)(表1)。由于黃河中上游地區(qū)風(fēng)速下降明顯且潛在蒸散對風(fēng)速敏感性較大，風(fēng)速成為過去50年潛在蒸散減少的主導(dǎo)因素，這與尹云鶴等［25］和高歌等［31］的研究結(jié)果一致。雖然潛在蒸散對實(shí)際水汽壓最敏感，但由于實(shí)際水汽壓的年變化率較低，其對潛在蒸散變化的影響較小，使得實(shí)際水汽壓成為對潛在蒸散變化貢獻(xiàn)最小的氣候因子。

表4 黃河中上游流域氣候因子對潛在蒸散變化的貢獻(xiàn)率Tab．4 Contribution of climate factors to E0change in the upper-middle reach of the Yellow River mm/a

季節(jié)尺度上，春季風(fēng)速對潛在蒸散的貢獻(xiàn)率最大。夏季太陽輻射強(qiáng)，日照時(shí)間值最大，且日照時(shí)間夏季變化率較大，敏感系數(shù)較大。日照時(shí)間成為夏季E0變化主導(dǎo)因子，風(fēng)速次之。秋季日最低氣溫對E0變化貢獻(xiàn)最大，風(fēng)速次之。由于日最高氣溫冬季敏感系數(shù)較大，同時(shí)其冬季的變化率最大，日最高氣溫成為冬季E0變化的主導(dǎo)因子，風(fēng)速則次之。在年尺度上，潛在蒸散年際變化貢獻(xiàn)最大的氣候因子為風(fēng)速，之后依次為最高氣溫、最低氣溫、日照時(shí)間和實(shí)際水汽壓。雖然潛在蒸散對實(shí)際水汽壓最敏感，但由于其變化微弱，因此對潛在蒸散的貢獻(xiàn)最小。

4 結(jié)論

1961—2010年黃河中上游地區(qū)氣溫顯著升高，日最低氣溫升高幅度大于日最高氣溫，冬季升溫幅度大于其他季節(jié);日照時(shí)間和風(fēng)速顯著下降;降水量、相對濕度呈下降趨勢，實(shí)際水汽壓和潛在蒸散呈升高趨勢，均未達(dá)顯著水平。流域半干旱區(qū)和半濕潤區(qū)的分界線(HI=0.5線)不斷向東南推進(jìn)，半干旱區(qū)面積不斷擴(kuò)大，黃河中上游地區(qū)氣候呈暖干趨勢變化。

近50年來黃河中上游地區(qū)潛在蒸散對實(shí)際水汽壓最敏感，之后依次為日最高氣溫、風(fēng)速和日照時(shí)間，對日最低氣溫最不敏感。夏季潛在蒸散對日最高氣溫最敏感，其他季節(jié)對實(shí)際水汽壓最敏感。近50年來潛在蒸散對日照時(shí)間敏感性顯著減少，對風(fēng)速、最低氣溫和實(shí)際水汽壓的敏感性顯著增加，對最高氣溫的敏感性有所增加，但未達(dá)顯著水平;潛在蒸散對風(fēng)速、最高氣溫和最低氣溫敏感系數(shù)隨海拔增加而減少，對日照時(shí)間敏感系數(shù)隨海拔增加而增加。

年尺度上，風(fēng)速對潛在蒸散變化的貢獻(xiàn)最大，之后依次為日最高氣溫、日最低氣溫和日照時(shí)間，實(shí)際水汽壓的貢獻(xiàn)最小。對潛在蒸散變化貢獻(xiàn)最大的因子春季為風(fēng)速，夏季為日照時(shí)間，秋季為日最低氣溫，冬季為日最高氣溫。各氣候因子對潛在蒸散變化的貢獻(xiàn)率春夏較高，秋冬較低。定量研究年和季節(jié)尺度上氣候因子變化對潛在蒸散的影響，有助于加深人們對氣候變化影響機(jī)制的認(rèn)識，并為氣候變化適應(yīng)對策研究和農(nóng)業(yè)可持續(xù)發(fā)展提供理論依據(jù)和數(shù)據(jù)支持。

［1］ 謝賢群，王菱.中國北方近50年潛在蒸發(fā)的變化［J］.自然資源學(xué)報(bào)，2007，22(5):683-691

［2］ 馬雪寧，張明軍，王圣杰，等.“蒸發(fā)悖論”在黃河流域的探討［J］.地理學(xué)報(bào)，2012，67(5):645-656

［3］ 劉昌明，張丹.中國地表潛在蒸散發(fā)敏感性的時(shí)空變化特征分析［J］.地理學(xué)報(bào)，2011，66(5):579-588

［4］ Thomas A.Spatial and temporal characteristics of potential evapotranspiration trends over China［J］.International Journal of Climatology，2000，20(4):381-396

［5］ Liu Qiang，Yang Zhifeng.Quantitative estimation of the impact of climate change on actual evapotranspiration in the Yellow River Basin，China［J］.Journal of Hydrology，2010，395:226-234

［6］ Yang Hanbo，Yang Dawen.Climatic factors influencing changing pan evaporation across China from 1961 to 2001［J］.Journal of Hydrology，2012，414/415:184-193

［7］ Li Zhi，Zheng Fenli，Liu Wenzhao.Spatiotemporal characteristic of potential evapotranspiration during 1961-2009 and its projected changes during 2011-2099 on the Loess Plateau of China［J］.Agricultural and Forest Meteorology，2012，154/155:147-155

［8］ Chen Shenbin，Liu Yunfeng，Thomas A.Climatic change on the Tibetan Plateau:Potentialevapotranspiration Trends from 1961—2000 ［J］.Climatic Change，2006，76:291-319

［9］ Zhang Yongqiang，Liu Changming，Tang Yanhong，et al.Trends in pan-evaporation and reference and actual evapotranspiration across the Tibetan Plateau［J］.Journal of Geographical Research，2006，112，D12110，doi:10.1029/2006JD008161

［10］ Xu Chongyu，Gong Lebing，Jiang Tong，et al.Analysis of spatial distribution and temporal trend of potential evapotranspiration and pan evaporation in Changjiang(Yangtze River)catchment［J］.Journal of Hydrology，2006，327:81-93

［11］ Gong Lebing，Xu Chongyu，Chen Deliang，et al.Sensitivity of the Penman-Monteith potential evapotranspiration to key climatic variables in the Changjiang(Yangtze River)basin［J］.Journal of Hydrology，2006，392:620-629

［12］ Zheng Hongxing，Liu Xiaomang，Liu Changming，et al.Assessing contributionsto pan-evaporation trendsin Haihe River Basin，China［J］.Journal of Geophysical Research， 2009， 114， D24105， doi:10.1029/2009JD012203

［13］ Meng Dejuan，Mo Xingguo.Assessing the effect of climate change on mean annual runoff in the Songhua River basin，China ［J］.Hydrological Processes，2012，26:1050-1061

［14］王遵婭，丁一匯，何金海，等.近50年來中國氣候變化特征的再分析［J］.氣象學(xué)報(bào)，2004，62(2):228-236

［15］楊小利，蒲金涌，馬棚里，等.陸面潛在蒸散計(jì)算模型在甘肅省黃河中上游地區(qū)的適用性研究［J］.水土保持通報(bào)，2010，30(2):184-189

［16］Allen R，Pereira L S，Raes D，et al.Crop evapotranspiration:Guidelines for computing crop requirements［M］.Rome:Food and Agricultural Organization of the United Nations，1998:56-57

［17］毛飛，孫涵，楊紅龍，等.干濕氣候區(qū)劃研究進(jìn)展［J］.地理科學(xué)進(jìn)展，2011，30(1):17-26

［18］王菱，謝賢群，李運(yùn)生，等.中國北方地區(qū)40年來濕潤指數(shù)和氣候干濕帶界線的變化［J］.地理研究，2004，23(1):45-54

［19］ Yin Yunhe，Wu Shaohong，Chen Gang，et al.Attribution analysis of potential evapotranspiration changes in China since the 1960s［J］.Theoretical and Applied Climatology，2010，101:19-28

［20］黃沖，張芳敏.全球氣候變化下的中國氣候變化特征分析［C］.杭州:第26屆中國氣象學(xué)會年會氣候變化分會場論文集，2009:1354-1360

［21］韓虹，任國玉，王文，等.黃河中上游地區(qū)太陽輻射時(shí)空演變特征［J］.氣候與環(huán)境研究，2008，13(1):61-66

［22］ McVicar T，Roderick M.Atmospheric science:Wind of change［J］.Nature Geoscience，2012，3(11):747-748

［23］ Xu Ming，Chang Chipei，F(xiàn)u Congbin，et al.Steady decline of east Asian monsoon winds，1969-2000:Evidence from direct ground measurements of wind speed ［J］.Journal of Geophysical Research，2006，111:D24111，doi:10.1029/2006JD007337

［24］江瀅，羅勇，趙宗慈.近50年中國風(fēng)速變化多氣候模式模擬檢驗(yàn)［J］.氣象學(xué)報(bào)，2009，67(6):923-934

［25］尹云鶴，吳紹洪，戴爾阜.1971—2008年我國潛在蒸散時(shí)空演變的歸因［J］.科學(xué)通報(bào)，2010，55(22):2226-2234

［26］徐宗學(xué)，張楠.黃河流域近50年降水變化趨勢分析［J］.地理研究，2006，25(1):27-34

［27］王幼奇，樊軍，邵明安，等.黃河中上游地區(qū)近50年潛在蒸散量變化特征［J］.農(nóng)業(yè)工程學(xué)報(bào)，2008，24(9):6-10

［28］楊建平，丁永健，陳仁升，等.近50年來中國干濕氣候界限的10年波動［J］.地理學(xué)報(bào)，2002，57(6):655-661

［29］劉文泉，王馥棠.黃河中上游地區(qū)地區(qū)農(nóng)業(yè)生產(chǎn)對氣候變化的脆弱性分析［J］.南京氣象學(xué)院學(xué)報(bào)，2002，25(5):620-624

［30］杜加強(qiáng)，舒儉民，劉成程，等.黃河上游潛在蒸散量變化特征及其對氣候變化的響應(yīng)［J］.農(nóng)業(yè)工程學(xué)報(bào)，2012，28(4):92-100

［31］高歌，陳德亮，任國玉，等.1956—2000年中國潛在蒸散量變化趨勢［J］.地理研究，2006，20(3):378-387