卓瑩瑩 趙慧霞 魏敏 隋瀟 劉厚鳳

摘 要：為分析氣候變化背景下黃河流域蒸發(fā)量變化特征及其影響因素，采用線性回歸、Mann-Kendall突變檢驗(yàn)分析和經(jīng)驗(yàn)正交分解等方法對(duì)1961—2019年黃河流域95個(gè)氣象觀測(cè)站蒸發(fā)皿蒸發(fā)量變化進(jìn)行分析。結(jié)果表明：近59 a黃河流域蒸發(fā)量呈顯著下降趨勢(shì)，全流域平均氣候傾向率為-60.4 mm/10 a;2004年以后，蒸發(fā)量下降趨勢(shì)加?。夂騼A向率為-219.4 mm/10 a）;夏季下降趨勢(shì)相對(duì)明顯，春、夏、秋、冬氣候傾向率分別為-12.7、-24.5、-10.4、-1.9 mm/10 a;不同河段蒸發(fā)量下降趨勢(shì)不同，上、中、下游的氣候傾向率分別為-59.3、-49.2、-74.9 mm/10 a。流域內(nèi)蒸發(fā)量氣候變化趨勢(shì)存在兩種空間分布模態(tài)：一致性分布模態(tài)和中部與東部、西部反位相的分布模態(tài)。通過(guò)偏相關(guān)、灰色關(guān)聯(lián)度和多元回歸分析發(fā)現(xiàn)，影響流域年蒸發(fā)皿蒸發(fā)量的主要因素有年均氣溫日較差、年日照時(shí)數(shù)、年降水量和相對(duì)濕度。

關(guān)鍵詞：蒸發(fā)皿蒸發(fā)量;時(shí)空變化特征;氣候因子;Mann-Kendall;正交經(jīng)驗(yàn)分解;黃河流域

中圖分類號(hào)：P467;TV882.1 文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2021.07.006

引用格式：卓瑩瑩，趙慧霞，魏敏，等.近59 a黃河流域蒸發(fā)量變化規(guī)律及影響因素[J].人民黃河，2021，43（7）：28-34，77.

Abstract： In order to analyze the variation of pan evaporation and the influencing factors of the Yellow River basin in the context of climate change， this paper used linear regression， Mann-Kendall and EOF to analyze the pan evaporations of 95 meteorological stations in the Yellow River basin from 1961 to 2019. The results show that the pan evaporation has shown a significant downward trend in the past 59 years. The climate inclination rate of the whole basin is -60.4 mm/10 a. However， the decline in evaporation has been accelerated， with an evaporation value of -219.4 mm / 10 a since 2004. For seasonal variation， the climate inclination rates in spring， summer， autumn and winter are -12.7， -24.5， -10.4 and -1.9 mm/10 a respectively. The climate inclination rate in the upper reaches （-59.3 mm/10 a） is the highest， followed by the middle reaches （-49.2 mm/10 a） and the lower reaches （-74.9 mm/10 a）. Two spatial distribution modes in the basin are identified， including consistent distribution and the anti-phase distribution of central and eastern and western. Based on partial correlation， grey correlation and multiple regression analysis， the main factors influencing the pan evaporation in the Yellow River basin are the annual mean daily temperature range， annual sunshine duration， annual precipitation and relative humidity.

Key words： pan evaporation; spatial-temporal distribution; climatic factors; Mann-Kendall; EOF; Yellow River basin

蒸發(fā)是水循環(huán)的重要環(huán)節(jié)之一，對(duì)水資源利用有著重要的影響。黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展是我國(guó)重大國(guó)家戰(zhàn)略，研究黃河流域蒸發(fā)量是黃河流域水資源和生態(tài)環(huán)境研究的基礎(chǔ)。目前已有學(xué)者對(duì)黃河流域的蒸發(fā)量進(jìn)行了大量研究[1-7]，其中：柳春等[4]對(duì)1961—2010年黃河流域蒸發(fā)皿蒸發(fā)量研究發(fā)現(xiàn)，該時(shí)段黃河流域蒸發(fā)量呈顯著下降趨勢(shì)，且四季中夏季的下降趨勢(shì)最顯著;童瑞等[6]對(duì)1961—2012年黃河流域蒸散發(fā)進(jìn)行研究，認(rèn)為黃河流域潛在蒸散發(fā)量亦呈減少趨勢(shì)，在不同河段的增減趨勢(shì)不同;邱新法等[5]研究認(rèn)為1960—2000年黃河流域氣溫雖呈上升趨勢(shì)，但蒸發(fā)皿蒸發(fā)量呈減少趨勢(shì)。

國(guó)內(nèi)外學(xué)者對(duì)蒸發(fā)量的影響因素也進(jìn)行了一系列的研究，認(rèn)為風(fēng)速、相對(duì)濕度、氣溫日較差、日照時(shí)數(shù)等是蒸發(fā)量變化的主要影響因素[8-16]。蒸發(fā)量對(duì)氣溫的響應(yīng)表現(xiàn)出負(fù)相關(guān)，在全球氣候變暖的背景下，我國(guó)的蒸發(fā)量出現(xiàn)減少的趨勢(shì)[17]，即存在“蒸發(fā)悖論”[18]的現(xiàn)象。馬雪寧等[19]的研究結(jié)果表明黃河流域總體上存在“蒸發(fā)悖論”，風(fēng)速的明顯減小是導(dǎo)致黃河流域潛在蒸發(fā)量減少的主要原因;柳春等[4]的研究也表明風(fēng)速減小是導(dǎo)致黃河流域蒸發(fā)量減少的主要原因;劉波等[20]研究認(rèn)為氣溫日較差和風(fēng)速是影響1960—2004年我國(guó)北方蒸發(fā)皿蒸發(fā)量變化的主要因素;Shen等[21]研究認(rèn)為風(fēng)速、氣溫日較差和降水量是影響1955—2001年我國(guó)干旱區(qū)蒸發(fā)皿蒸發(fā)量變化的主要因素;Fu等[22]研究認(rèn)為近50 a來(lái)太陽(yáng)輻射、蒸氣壓差和風(fēng)速是影響蒸發(fā)皿蒸發(fā)量變化的三大潛在因素。

已有對(duì)黃河流域蒸發(fā)皿蒸發(fā)量時(shí)空變化特征及其影響因素的研究時(shí)段多集中于2012年之前，對(duì)于更長(zhǎng)時(shí)間尺度的研究尚未見(jiàn)報(bào)道。對(duì)黃河流域長(zhǎng)時(shí)間序列（1961—2019年）蒸發(fā)皿蒸發(fā)量開(kāi)展研究，可闡明全球氣候變暖情境下黃河流域蒸發(fā)量的最新變化趨勢(shì)及影響因素，對(duì)黃河流域水資源平衡、水資源高效利用、生態(tài)環(huán)境保護(hù)及經(jīng)濟(jì)高質(zhì)量發(fā)展具有重要參考意義。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

黃河流域位于我國(guó)中緯度地帶（東經(jīng)96°—119°、北緯32°—42°），流域總面積為79.5萬(wàn)km2。黃河流域地貌特征復(fù)雜多樣，東西橫跨4個(gè)地貌單元，分別是青藏高原、內(nèi)蒙古高原、黃土高原和黃淮海平原。黃河流域氣候以溫帶大陸性氣候?yàn)橹?，具有氣溫年較差和日較差大、降水時(shí)間集中且年際變化大、蒸發(fā)量較大等特征。

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

氣象數(shù)據(jù)來(lái)自中國(guó)氣象局國(guó)家氣象信息中心，主要為1961—2019年黃河流域95個(gè)國(guó)家氣象觀測(cè)站（見(jiàn)圖1）的逐月蒸發(fā)皿蒸發(fā)量、平均氣溫、平均氣溫日較差、平均相對(duì)濕度、平均2 min風(fēng)速、日照時(shí)數(shù)和降水量數(shù)據(jù)。其中2001—2010年有76個(gè)站點(diǎn)、2011—2019年有91個(gè)站點(diǎn)非結(jié)冰期只有大型蒸發(fā)量數(shù)據(jù)，無(wú)小型蒸發(fā)量數(shù)據(jù)，為了確保蒸發(fā)量數(shù)據(jù)的連續(xù)性和可比性，采用比值法[23]將該時(shí)段的大型蒸發(fā)量數(shù)據(jù)折算為小型蒸發(fā)量數(shù)據(jù)。地圖底圖數(shù)據(jù)來(lái)源于標(biāo)準(zhǔn)地圖服務(wù)系統(tǒng)（http：//bzdt.ch.mnr.gov.cn/）。

1.3 研究方法

采用線性回歸法和Mann-Kendall突變檢驗(yàn)法分析黃河流域蒸發(fā)皿蒸發(fā)量隨時(shí)間的變化趨勢(shì)，運(yùn)用經(jīng)驗(yàn)正交函數(shù)分解法分析黃河流域蒸發(fā)皿蒸發(fā)量的空間變化特征，利用SPSS23.0對(duì)蒸發(fā)皿蒸發(fā)量的影響因子進(jìn)行相關(guān)性分析。

1.3.1 線性回歸法

用yi表示樣本量為n的某一氣象要素變量，用xi表示yi對(duì)應(yīng)的時(shí)間段，建立yi和xi之間的一元線性回歸方程：

式中：i為回歸計(jì)算值;a、b分別為回歸常數(shù)、回歸系數(shù)，a和b可用最小二乘法估算。

1.3.2 Mann-Kendall突變檢驗(yàn)法

Mann-Kendall（M-K）突變檢驗(yàn)法是一種氣象學(xué)上用于判斷變化趨勢(shì)及是否發(fā)生突變的非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，該方法的優(yōu)點(diǎn)是樣本數(shù)據(jù)不需要遵循一定的分布，也不受少數(shù)異常值的影響[24]。

1.3.3 經(jīng)驗(yàn)正交函數(shù)分解法

經(jīng)驗(yàn)正交函數(shù)分解法（EOF）也被稱為主成分分析法，是氣象診斷中經(jīng)常使用的方法之一。該方法的原理是將原變量場(chǎng)分解為幾個(gè)互不相關(guān)的典型模態(tài)來(lái)代替原變量場(chǎng)，在典型模態(tài)中盡可能多地保留原變量場(chǎng)的相關(guān)信息。該方法具有易于使用、不受站點(diǎn)分布限制、原始信息保留較完整和快速展開(kāi)收斂等特點(diǎn)。

2 黃河流域蒸發(fā)皿蒸發(fā)量隨時(shí)間變化規(guī)律

2.1 年際變化特征

統(tǒng)計(jì)黃河流域95個(gè)氣象站點(diǎn)1961—2019年月蒸發(fā)皿蒸發(fā)量數(shù)據(jù)發(fā)現(xiàn)，流域年均蒸發(fā)量為1 666.6 mm，年際波動(dòng)大（年變異系數(shù)為0.08），最大值出現(xiàn)在1972年，為1 955.5 mm，最小值出現(xiàn)在2019年，為1 267.4 mm，二者相差688.1 mm。通過(guò)對(duì)黃河流域上、中、下游數(shù)據(jù)分析發(fā)現(xiàn)，年均蒸發(fā)皿蒸發(fā)量最小的是上游（1 638.3 mm），中游居中，為1 660.4 mm，最大的是下游（1 730.6 mm）。由1961—2019年黃河流域年蒸發(fā)皿蒸發(fā)量變化趨勢(shì)（見(jiàn)圖2）分析得，1961—2019年蒸發(fā)皿蒸發(fā)量總體呈下降趨勢(shì)，平均氣候傾向率為-60.4 mm/10 a（通過(guò)α=0.05顯著性水平檢驗(yàn)）。2004年以后，蒸發(fā)量下降趨勢(shì)加劇，氣候傾向率達(dá)-219.4 mm/10 a（通過(guò)α=0.05顯著性水平檢驗(yàn)）。1961—2019年上、中、下游蒸發(fā)皿蒸發(fā)量氣候傾向率分別為-59.3、-49.2、-74.9 mm/10 a，均通過(guò)α=0.05顯著性水平檢驗(yàn)。這與Xu等[25]、Ji等[26]研究結(jié)論一致。

圖2 1961—2019年黃河流域年蒸發(fā)皿蒸發(fā)量變化趨勢(shì)

對(duì)黃河流域年蒸發(fā)皿蒸發(fā)量進(jìn)行Mann-Kendall突變檢驗(yàn)（見(jiàn)圖3）可知：1975—2019年（除部分年份外），年蒸發(fā)皿蒸發(fā)量整體呈減少趨勢(shì)，且在1982年發(fā)生減少性突變。UF統(tǒng)計(jì)量在2009—2019年超過(guò)臨界值U0.05=±1.96，說(shuō)明在2009—2019年年蒸發(fā)皿蒸發(fā)量減少趨勢(shì)顯著。

2.2 季節(jié)變化特征

1961—2019年，黃河流域蒸發(fā)皿蒸發(fā)量存在明顯的季節(jié)差異（見(jiàn)圖4）。夏季蒸發(fā)量最大，平均值為661.1 mm，冬季最小，平均值為156.2 mm，夏季蒸發(fā)量為冬季的4.2倍。流域內(nèi)蒸發(fā)主要集中在春、夏兩個(gè)季節(jié)，占全年蒸發(fā)量的73.4%。四季蒸發(fā)皿蒸發(fā)量均呈降低的趨勢(shì)，春、夏、秋、冬氣候傾向率分別為-12.7、-24.5、-10.4、-1.9 mm/10 a，其中春、夏、秋三季降低趨勢(shì)顯著，冬季蒸發(fā)量減少不明顯。夏季蒸發(fā)皿蒸發(fā)量的變化幅度遠(yuǎn)大于其他3個(gè)季節(jié)，對(duì)全年蒸發(fā)皿蒸發(fā)量變化趨勢(shì)起了主導(dǎo)作用。

黃河流域各河段各季節(jié)蒸發(fā)皿蒸發(fā)量均呈顯著減少趨勢(shì)，其氣候傾向率見(jiàn)表1，四季蒸發(fā)量減少幅度最大的河段均為下游，各河段不同季節(jié)蒸發(fā)量減少幅度最大的是夏季。

對(duì)黃河流域四季蒸發(fā)皿蒸發(fā)量進(jìn)行M-K突變檢驗(yàn)（見(jiàn)圖5）可知：春季蒸發(fā)皿蒸發(fā)量在1980年前后發(fā)生減少性突變，1987—2019年減少趨勢(shì)顯著。夏季蒸發(fā)皿蒸發(fā)量在1976—2019年呈減少趨勢(shì)，在1976年發(fā)生減少性突變。秋季蒸發(fā)皿蒸發(fā)量在2012年前后發(fā)生減少性突變，2015—2019年UF統(tǒng)計(jì)量超過(guò)U0.05=±1.96，說(shuō)明減少趨勢(shì)顯著。冬季蒸發(fā)皿蒸發(fā)量在1964年發(fā)生減少性突變。

3 黃河流域蒸發(fā)皿蒸發(fā)量空間變化規(guī)律

采用EOF法對(duì)1961—2019年黃河流域蒸發(fā)皿蒸發(fā)量進(jìn)行時(shí)空分解，并進(jìn)行顯著性檢驗(yàn)，得到黃河流域年蒸發(fā)場(chǎng)主要時(shí)空分布形式（見(jiàn)圖6和圖7）。特征向量和特征值誤差范圍分析結(jié)果表明，蒸發(fā)場(chǎng)前6個(gè)特征向量累計(jì)貢獻(xiàn)率達(dá)82.6%，且通過(guò)North顯著性檢驗(yàn)，因此前6個(gè)特征向量能較好地表達(dá)1961—2019年黃河流域蒸發(fā)皿蒸發(fā)量的空間分布特征。本文選用特征向量貢獻(xiàn)率高的第一特征向量和第二特征向量進(jìn)行分析。

第一特征向量貢獻(xiàn)率為41.1%，遠(yuǎn)高于其他特征向量的貢獻(xiàn)率，為黃河流域蒸發(fā)皿蒸發(fā)量的主要空間分布形式（見(jiàn)圖6（a））。該特征向量的各分量均為正值，這表明黃河流域各區(qū)域變化趨勢(shì)基本一致，即黃河流域均為蒸發(fā)量偏大或蒸發(fā)量偏小的分布。墾利、泰山和三門峽為數(shù)值較大處，即變化量大的地區(qū)。

第二特征向量的貢獻(xiàn)率為1.7%，也是1961—2019年黃河流域蒸發(fā)皿蒸發(fā)量的主要空間分布形式（見(jiàn)圖6（b））。該特征向量的各分量正負(fù)值均有，說(shuō)明黃河流域有兩種分布模式。流域的東部和西部為負(fù)值地區(qū)，流域的中部為正值地區(qū)。西寧和墾利為負(fù)值地區(qū)的高值區(qū)，吳忠和河曲為正值地區(qū)的高值區(qū)，這說(shuō)明上述地區(qū)為蒸發(fā)量波動(dòng)變化大的地區(qū)。

根據(jù)EOF法分解的黃河流域蒸發(fā)量空間形式可知，黃河流域蒸發(fā)場(chǎng)存在4種形式：第一特征向量決定全年全流域蒸發(fā)量偏多或偏少的分布特點(diǎn);第二特征向量決定全年流域東部和西部蒸發(fā)量偏多、中部蒸發(fā)量偏少或者東部和西部蒸發(fā)量偏少、中部蒸發(fā)量偏多的分布特點(diǎn)。

特征向量所對(duì)應(yīng)的時(shí)間系數(shù)代表了這一地區(qū)由特征向量所表征的分布形式隨時(shí)間變化特征。時(shí)間系數(shù)為正值代表與空間模態(tài)同方向，負(fù)值則相反。時(shí)間系數(shù)的絕對(duì)值越大，表示這種分布形式越顯著。第一特征向量時(shí)間系數(shù)線性趨勢(shì)預(yù)測(cè)呈下降趨勢(shì)（見(jiàn)圖7（a）），說(shuō)明黃河流域59 a來(lái)蒸發(fā)量有減少趨勢(shì)。在1961—2019年，時(shí)間系數(shù)為正值的有31 a，這31 a全流域蒸發(fā)量偏少，剩下的28 a蒸發(fā)量偏多。第二特征向量時(shí)間系數(shù)線性趨勢(shì)預(yù)測(cè)呈上升趨勢(shì)（見(jiàn)圖7（b）），說(shuō)明59 a來(lái)流域中部蒸發(fā)量呈增加趨勢(shì)，東部和西部蒸發(fā)量呈減少趨勢(shì)。

4 影響蒸發(fā)皿蒸發(fā)量的氣候因子

影響蒸發(fā)量的氣候因子有很多，按影響因子的性質(zhì)可分為熱力因子、水分因子和動(dòng)力因子三類。選用以氣溫、氣溫日較差和日照時(shí)數(shù)為代表的熱力因子，以相對(duì)濕度和降水量為代表的水分因子，以風(fēng)速為代表的動(dòng)力因子，對(duì)它們與蒸發(fā)皿蒸發(fā)量的關(guān)系進(jìn)行偏相關(guān)分析。

4.1 氣候因子隨時(shí)間變化規(guī)律

蒸發(fā)皿蒸發(fā)量與影響因子之間的量綱不同，為進(jìn)一步準(zhǔn)確分析其相關(guān)性，對(duì)所有數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，結(jié)果見(jiàn)圖8。黃河流域蒸發(fā)皿蒸發(fā)量與氣溫日較差和日照時(shí)數(shù)的變化趨勢(shì)相似，與氣溫的變化趨勢(shì)相反，與相對(duì)濕度、降水量波動(dòng)趨勢(shì)相反。

各氣候因子在四季的變化趨勢(shì)不一（見(jiàn)表2）。近59 a來(lái)，黃河流域夏季氣溫上升趨勢(shì)最顯著，冬季呈顯著降低趨勢(shì);熱力因子中的氣溫日較差、日照時(shí)數(shù)除春季外其他季節(jié)均呈減小趨勢(shì);風(fēng)速和相對(duì)濕度各季節(jié)均呈減小趨勢(shì);降水量在春、秋季呈不顯著減少趨勢(shì)，在夏、冬季呈不顯著增加趨勢(shì)。

通過(guò)相關(guān)系數(shù)矩陣（見(jiàn)表3）分析發(fā)現(xiàn)，研究區(qū)域各氣候因子間也存在一定相關(guān)性，相關(guān)性最顯著的是日照時(shí)數(shù)與氣溫日較差，說(shuō)明日照時(shí)數(shù)通過(guò)影響太陽(yáng)輻射量對(duì)氣溫日較差產(chǎn)生影響。降水量與氣溫日較差、日照時(shí)數(shù)均為顯著負(fù)相關(guān)關(guān)系，氣溫日較差和日照時(shí)數(shù)與相對(duì)濕度也均為顯著負(fù)相關(guān)關(guān)系，氣候因子間的相互關(guān)系都將影響蒸發(fā)皿蒸發(fā)量的變化。

4.2 蒸發(fā)皿蒸發(fā)量與氣候因子的相關(guān)性

黃河流域蒸發(fā)皿蒸發(fā)量與氣候影響因子分析（見(jiàn)表4）表明流域蒸發(fā)量與氣溫日較差相關(guān)性最顯著。從全流域看，春、秋、冬季蒸發(fā)皿蒸發(fā)量的最大影響因子均為氣溫日較差，相關(guān)系數(shù)分別為0.797、0.911、0.844，除此之外秋季蒸發(fā)皿蒸發(fā)量受到日照時(shí)數(shù)的影響較大，相關(guān)系數(shù)為0.865。春季和冬季蒸發(fā)皿蒸發(fā)量受相對(duì)濕度的影響較大，相關(guān)系數(shù)分別為-0.761和-0.791。全年和夏季蒸發(fā)皿蒸發(fā)量受降水量因素影響也較大，相關(guān)系數(shù)分別為-0.701和-0.636。

黃河流域東西跨度大，不同河段之間的大氣環(huán)流具有一定的差異，導(dǎo)致氣溫、氣溫日較差、風(fēng)速、相對(duì)濕度等因子的變化規(guī)律在各河段也表現(xiàn)出不同的特征（見(jiàn)表4）。采用偏相關(guān)分析法分析不同河段蒸發(fā)皿蒸發(fā)量的主要影響因素，結(jié)果表明各河段全年及四季的蒸發(fā)量與氣溫日較差相關(guān)性最顯著，但不同河段、不同季節(jié)其與各影響因子之間的相關(guān)系數(shù)不同。夏季，上、中、下游蒸發(fā)皿蒸發(fā)量的主要影響因子為水分因子，濕度和降水量較大導(dǎo)致蒸發(fā)量減少，黃河流域大部分處于季風(fēng)氣候區(qū)，夏季是黃河流域降水量和云量最多的季節(jié)，云量多少通過(guò)影響日照時(shí)數(shù)影響太陽(yáng)輻射的強(qiáng)弱，進(jìn)而影響蒸發(fā)量。秋、冬季蒸發(fā)量的主要影響因子為氣溫日較差（熱力因子），春季蒸發(fā)量的主要影響因子上游為氣溫日較差（熱力因子），中、下游為相對(duì)濕度。

在目前針對(duì)黃河流域蒸發(fā)皿蒸發(fā)量影響因素的研究中，不同學(xué)者的結(jié)論有所不同，原因是已有研究在影響因子的選擇上只有柳春等[4]將氣溫日較差納入分析，童瑞等[6]、馬雪寧等[19]、Xu等[25]和Ji等[26]未將該因子納入分析。同時(shí)，研究的時(shí)間序列不同也可能會(huì)導(dǎo)致分析結(jié)果的不同。

進(jìn)一步對(duì)年蒸發(fā)皿蒸發(fā)量相關(guān)因子進(jìn)行灰色關(guān)聯(lián)度和多元回歸分析，結(jié)果表明，其重復(fù)的因子為年均氣溫日較差、年日照時(shí)數(shù)和年降水量，與前述偏相關(guān)分析結(jié)果一致，這說(shuō)明年均氣溫日較差、年日照時(shí)數(shù)和年降水量是黃河流域年蒸發(fā)量變化的主要影響因子。影響蒸發(fā)量的因素有許多，不同的影響因子不僅受到不同作用機(jī)制的影響，同時(shí)影響因子間也存在相互影響，因此在該方面還有待進(jìn)行進(jìn)一步研究。

5 結(jié) 論

（1）1961—2019年黃河流域年蒸發(fā)皿蒸發(fā)量呈顯著下降趨勢(shì)，2004年以來(lái)下降趨勢(shì)加劇，春、夏、秋、冬氣候傾向率分別為-12.7、-24.5、-10.4、-1.9 mm/10 a。上、中、下游的氣候傾向率分別為-59.3、-49.2、-74.9 mm/10 a。

（2）根據(jù)EOF法分析，前兩個(gè)特征向量的累計(jì)貢獻(xiàn)率達(dá)42.8%，可以較好地反映流域內(nèi)蒸發(fā)量變化特征的空間分布模態(tài)。第一特征向量場(chǎng)呈一致性分布，第二特征向量場(chǎng)表現(xiàn)為中部與東部、西部反位相的分布格局。

（3）對(duì)黃河流域年蒸發(fā)皿蒸發(fā)量影響較大的因子有氣溫日較差、日照時(shí)數(shù)、相對(duì)濕度和年降水量，其中氣溫日較差是最重要的影響因子。不同的時(shí)間尺度以及不同的河段其主要影響因子有所不同，全流域、中游、下游的全年和四季影響蒸發(fā)皿蒸發(fā)量變化的主要因子均為氣溫日較差、日照時(shí)數(shù)和相對(duì)濕度;上游地區(qū)影響蒸發(fā)皿蒸發(fā)量變化的主要因子為氣溫日較差、日照時(shí)數(shù)和降水量。

參考文獻(xiàn)：

[1] 丁磊，劉廷璽，張海霞，等.黃河內(nèi)蒙古段干流水面蒸發(fā)量變化特征[J].人民黃河，2012，34（6）：38-40.

[2] 李高偉，韓美，張東啟.1961—2013年黃河三角洲氣候變化趨勢(shì)研究[J].人民黃河，2017，39（1）：30-37.

[3] 劉猛，夏自強(qiáng)，李俊芬，等.全球氣候變暖背景下黃河源區(qū)蒸發(fā)變化研究[J].人民黃河，2008，30（9）：30-31，108.

[4] 柳春，王守榮，梁有葉，等.1961—2010年黃河流域蒸發(fā)皿蒸發(fā)量變化及影響因子分析[J].氣候變化研究進(jìn)展，2013，9（5）：327-334.

[5] 邱新法，劉昌明，曾燕.黃河流域近40年蒸發(fā)皿蒸發(fā)量的氣候變化特征[J].自然資源學(xué)報(bào)，2003，18（4）：437-442.

[6] 童瑞，楊肖麗，任立良，等.黃河流域1961—2012年蒸散發(fā)時(shí)空變化特征及影響因素分析[J].水資源保護(hù)，2015，31（3）：16-21.

[7] 于蘭蘭，劉健.近50 a黃河三角洲潛在蒸發(fā)變化特征分析[J].人民黃河，2012，34（12）：46-48.

[8] HE Y， LIN K， CHEN X， et al. Classification-Based Spatiotemporal Variations of Pan Evaporation Across the Guangdong Province， South China [J]. Water Resources Management， 2015， 29（3）： 901-912.

[9] LIU Y J， CHEN J， PAN T. Analysis of Changes in Reference Evapotranspiration， Pan Evaporation， and Actual Evapotranspiration and Their Influencing Factors in the North China Plain During 1998-2005 [J]. Earth Space Sci， 2019， 6（8）： 1366-1377.

[10] 曾燕，邱新法，劉昌明，等.1960—2000年中國(guó)蒸發(fā)皿蒸發(fā)量的氣候變化特征[J].水科學(xué)進(jìn)展，2007，18（3）：311-318.

[11] 劉敏，沈彥俊，曾燕，等.近50年中國(guó)蒸發(fā)皿蒸發(fā)量變化趨勢(shì)及原因[J].地理學(xué)報(bào)，2009，64（3）：259-269.

[12] 祁添垚，張強(qiáng)，王月，等.1960—2005年中國(guó)蒸發(fā)皿蒸發(fā)量變化趨勢(shì)及其影響因素分析[J].地理科學(xué)，2015，35（12）：1599-1606.

[13] 申雙和，盛瓊.45年來(lái)中國(guó)蒸發(fā)皿蒸發(fā)量的變化特征及其成因[J].氣象學(xué)報(bào)，2008，66（3）：452-460.

[14] 陳旭，戴安全，鞠彬.哈巴河地區(qū)潛在蒸發(fā)變化特征及影響因素分析[J].人民黃河，2016，38（11）：54-57.

[15] 廖曉芳，錢勝，彭彥銘，等.蒸發(fā)皿蒸發(fā)和潛在蒸散發(fā)對(duì)氣候變化的響應(yīng)[J].人民黃河，2010，32（11）：42-44，152.

[16] 嚴(yán)登明，翁白莎，王俊鋒，等.皖北地區(qū)“蒸發(fā)悖論”探討[J].人民黃河，2017，39（4）：21-24.

[17] 左洪超，李棟梁，胡隱樵，等.近40 a中國(guó)氣候變化趨勢(shì)及其同蒸發(fā)皿觀測(cè)的蒸發(fā)量變化的關(guān)系[J].科學(xué)通報(bào)，2005，50（11）：1125-1130.

[18] RODERICK M L， FARQUHAR G D. The Cause of Decreased Pan Evaporation over the Past 50 Years [J]. Science， 2002， 298（5597）： 1410-1411.

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

[20] 劉波，馬柱國(guó)，丁裕國(guó).中國(guó)北方近45年蒸發(fā)變化的特征及與環(huán)境的關(guān)系[J].高原氣象，2006，25（5）：840-848.

[21] SHEN Y J， LIU C M， LIU M， et al. Change in Pan Evaporation over the Past 50 Years in the Arid Region of China [J]. Hydrol Process， 2010， 24（2）： 225-231.

[22] FU G B， CHARLES S P， YU J J. A Critical Overview of Pan Evaporation Trends over the Last 50 Years [J]. Climate Change， 2009，97（1-2）：193-214.

[23] 褚榮浩，申雙和，李萌，等.小型與E-601型蒸發(fā)皿蒸發(fā)量對(duì)比分析及其折算系數(shù)：以江蘇省為例 [J].氣象科學(xué)，2018，38（2）：247-257.

[24] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京：氣象出版社，2007：69-72.

[25] XU S Q， YU Z B， YANG C G， et al. Trends in Evapotranspiration and Their Responses to Climate Change and Vegetation Greening over the Upper Reaches of the Yellow River Basin[J]. Agric For Meteorol， 2018， 263：118-129.

[26] JI X J， WANG J J， GU W L， et al. Trends in Annual and Seasonal Pan Evaporation in the Lower Yellow River Basin from 1961 to 2010 [J]. Advances in Climate Change Research，2012，3（4）：195-204.

【責(zé)任編輯 張 帥】