邵 駿,錢曉燕,謝 珊,杜 濤,湯瑤瑤,向碧為

(1.長(zhǎng)江水利委員會(huì)水文局,武漢 430010; 2.武漢市水務(wù)科學(xué)研究院,武漢 430014; 3.長(zhǎng)江水利委員會(huì)水文局長(zhǎng)江上游水文水資源勘測(cè)局,重慶 400020; 4.上海勘測(cè)設(shè)計(jì)研究院有限公司,上海 200435)

0 引 言

水文循環(huán)是氣候系統(tǒng)的重要組成部分。全球氣候變化通過(guò)影響氣溫、降水、蒸發(fā)、徑流、土壤水分等水文循環(huán)過(guò)程,對(duì)徑流的時(shí)空分布產(chǎn)生重大改變[1-3]。由于降水是影響河川徑流年際和年代際變化的主要驅(qū)動(dòng)因子,而降水本身又取決于一定的大氣環(huán)流形勢(shì)與局地天氣過(guò)程,因此河川徑流的演變與大氣環(huán)流是聯(lián)系在一起的。ENSO(El Nio Southern Oscillation)是以赤道東太平洋海表溫度異常為特征的周期性現(xiàn)象,是全球氣候系統(tǒng)中最為顯著的年際信號(hào)。ENSO事件對(duì)我國(guó)氣候的影響得到了廣泛關(guān)注和研究,Blender等[4]和Wei等[5]研究表明長(zhǎng)江流域的降水和徑流受東亞夏季季風(fēng)和ENSO事件影響較大;Gong等[6]研究發(fā)現(xiàn)我國(guó)東部冬季、秋季降水量與ENSO事件有顯著的關(guān)系,并在ENSO暖事件(El Nio)年降水量顯著減少;王根緒等[7]研究發(fā)現(xiàn)ENSO事件與黃河源區(qū)降水與徑流年際波動(dòng)變化有很好的相關(guān)性,暖事件對(duì)應(yīng)年份黃河源區(qū)徑流量減少,而冷事件則使其增加;Zhu[8]的研究發(fā)現(xiàn),在ENSO發(fā)展階段中國(guó)北方降水偏少,而在ENSO衰減階段長(zhǎng)江流域降水偏多;彭俊等[9]研究發(fā)現(xiàn)長(zhǎng)江流域徑流量與ENSO事件具有相同的2～8 a周期變化,在El Nio發(fā)生期,徑流量較低,在ENSO冷事件(La Nia)發(fā)生期,徑流量較高;藍(lán)永超等[10]研究發(fā)現(xiàn)El Nio事件發(fā)生當(dāng)年,祁連山東段降水偏少、徑流偏枯,El Nio事件次年祁連山區(qū)東段降水減少、徑流偏估的程度不如El Nio事件當(dāng)年那樣顯著。

長(zhǎng)江源區(qū)近些年來(lái)受全球氣候變化影響巨大,主要表現(xiàn)在氣溫升高、徑流持續(xù)增加并一直延續(xù)至今。長(zhǎng)江源區(qū)位于青藏高原腹地,其河川徑流受青藏高原大尺度氣候影響顯著。研究表明,ENSO事件顯著影響次年青藏高原西南部夏季季節(jié)和月平均降水與溫度[11]。ENSO事件發(fā)生以后,通過(guò)直接影響Walker環(huán)流、Hadley環(huán)流和西北太平洋異常反氣旋(Northwest Pacific Anomalous Anticyclone,NWPAC),從而調(diào)制太平洋副熱帶高壓和西太平洋夏秋季臺(tái)風(fēng),影響孟加拉灣和印度洋水汽向青藏高原的輸送[11],進(jìn)而影響長(zhǎng)江源區(qū)的天氣系統(tǒng)和局地氣候。目前尚未有研究關(guān)注長(zhǎng)江源區(qū)徑流變化對(duì)ENSO事件的響應(yīng)。研究長(zhǎng)江源區(qū)徑流演變與ENSO事件之間的響應(yīng)關(guān)系,可為長(zhǎng)江源區(qū)中長(zhǎng)期徑流預(yù)測(cè)、水資源管理和保護(hù)政策的制定提供重要的參考。

徑流與氣候之間的關(guān)系是非線性的、非平穩(wěn)的,僅通過(guò)分析觀測(cè)到的徑流數(shù)據(jù)的變化規(guī)律,很難將水文變化與氣候變化充分聯(lián)系起來(lái)。同時(shí)由于ENSO事件的發(fā)生和發(fā)展機(jī)理較為復(fù)雜,其對(duì)河川徑流的影響并非存在直接的關(guān)聯(lián)性,因此利用統(tǒng)計(jì)學(xué)方法來(lái)分析長(zhǎng)江源區(qū)徑流對(duì)ENSO事件的響應(yīng)不失為有效方法之一。本文旨在通過(guò)沱沱河水文站、直門達(dá)水文站實(shí)測(cè)徑流數(shù)據(jù),分析長(zhǎng)江源區(qū)徑流演變事實(shí),利用海洋尼諾指數(shù)分析ENSO事件的強(qiáng)度與時(shí)間特征,研究徑流豐枯變化對(duì)ENSO事件的響應(yīng)規(guī)律及聯(lián)系,為揭示氣候變化背景下長(zhǎng)江源區(qū)水循環(huán)機(jī)制提供支撐。

1 研究數(shù)據(jù)及研究方法

1.1 研究數(shù)據(jù)

ENSO是厄爾尼諾-拉尼娜事件以及南方濤動(dòng)的合稱,表征ENSO的指標(biāo)較多,本文選取海洋尼諾指數(shù)(Oceanic Nio Index,ONI)作為研究指標(biāo)。ONI是中東太平洋Nio 3.4區(qū)(5°N—5°S,120°W—170°W)的海表溫度距平的3個(gè)月滑動(dòng)平均值。研究數(shù)據(jù)來(lái)源為美國(guó)國(guó)家海洋和大氣管理局(NOAA)氣候預(yù)測(cè)中心。

根據(jù)NOAA對(duì)ENSO事件的確定標(biāo)準(zhǔn),ONI連續(xù)5個(gè)月≥0.5 ℃(≤0.5 ℃)定義為1次El Nio(La Nia)事件。參考文獻(xiàn)[12]的研究成果,對(duì)ONI指數(shù)累積值ΣONI進(jìn)行標(biāo)準(zhǔn)化處理,將強(qiáng)度劃分為極強(qiáng)、強(qiáng)、中等、弱、極弱5個(gè)等級(jí)(表1)。

表1 ENSO事件的海溫強(qiáng)度指標(biāo)Table 1 Sea surface temperature (SST) intensity indexof ENSO event

本研究水文數(shù)據(jù)采用長(zhǎng)江源區(qū)沱沱河站和直門達(dá)站1960—2020年實(shí)測(cè)月均流量。

1.2 研究方法

本研究主要采用Mann-Kendall非參數(shù)檢驗(yàn)、交叉小波變換[13]等方法,簡(jiǎn)述如下。

1.2.1 Mann-Kendall非參數(shù)檢驗(yàn)

時(shí)間序列數(shù)據(jù)(x1,…,xn)是n個(gè)獨(dú)立的、隨機(jī)變量同分布的樣本,對(duì)于所有的k,j

(1)

其中,

(2)

定義統(tǒng)計(jì)變量:

(3)

其中:

E(sk)=k(k-1)/4 ;

(4)

Var(sk)=k(k-1)(2k+5)/72 。

(5)

式中UFk為標(biāo)準(zhǔn)正態(tài)分布。對(duì)于給定顯著性水平α,若|UFk|>UFα/2,則表明序列存在明顯的趨勢(shì)變化。

1.2.2 交叉小波變換

r(X,Y)=

(6)

通過(guò)小波凝聚譜可以量度時(shí)頻范圍內(nèi)兩者之間局部相關(guān)的密切程度,其表達(dá)式為

式中S表示平滑函數(shù)。本文選擇Morlet小波為小波基函數(shù)。

2 研究結(jié)果與討論

2.1 長(zhǎng)江源區(qū)徑流演變趨勢(shì)

圖1為沱沱河站、直門達(dá)站年徑流量的滑動(dòng)平均曲線和Mann-Kendall曲線。由兩站的滑動(dòng)平均曲線可見(jiàn),沱沱河站在20世紀(jì)70—80年代呈現(xiàn)略微下降的態(tài)勢(shì),在2000年之后呈現(xiàn)出大幅增加的態(tài)勢(shì)。直門達(dá)站在2000年之前基本維持在均值附近小幅震蕩,在2005年之后呈現(xiàn)升高的態(tài)勢(shì)。由兩站的Mann-Kendall曲線可見(jiàn),2002年左右沱沱河站正、逆序列出現(xiàn)交叉,2009年左右直門達(dá)站正、逆序列出現(xiàn)交叉,表明兩站在2002年、2009年均發(fā)生突變?cè)黾拥膽B(tài)勢(shì)。采用Mann-Kendall檢驗(yàn)對(duì)兩站年徑流量進(jìn)行趨勢(shì)檢驗(yàn)(結(jié)果見(jiàn)表2),兩站Mann-Kendall檢驗(yàn)結(jié)果均通過(guò)置信度為95%的顯著性檢驗(yàn),表明兩站均呈現(xiàn)顯著的上升趨勢(shì)。

圖1 沱沱河站、直門達(dá)站年徑流量滑動(dòng)平均曲線與Mann-Kendall曲線Fig.1 Moving-average curves and Mann-Kendall curves of annual runoff at Tuotuohe and Zhimenda hydrological stations

表2 沱沱河站、直門達(dá)站年徑流量Mann-Kendall檢驗(yàn)結(jié)果Table 2 Mann-Kendall test results of annual runoff at Tuotuohe and Zhimenda hydrological stations

2.2 徑流變化對(duì)ENSO事件的響應(yīng)

2.2.1 徑流豐枯變化與ENSO事件的關(guān)系

將1960年1月—2020年12月以來(lái)的ENSO事件強(qiáng)度進(jìn)行劃分,成果如表3所示。同時(shí),統(tǒng)計(jì)沱沱河站和直門達(dá)站相應(yīng)于ENSO事件對(duì)應(yīng)年份的年徑流距平程度,結(jié)果也列于表3。為便于直觀分析,將表征ENSO事件強(qiáng)度的ONI變化過(guò)程線與沱沱河站、直門達(dá)站徑流距平過(guò)程線繪制于圖2。

圖2 1960—2020年ENSO事件及與沱沱河站、直門達(dá)站年徑流量過(guò)程線對(duì)比Fig.2 ENSO events during 1960-2020 and comparison with annual runoff hydrograph of Tuotuohe and Zhimenda hydrological stations

表3 1960—2020年ENSO事件特征值及沱沱河站、直門達(dá)站年徑流距平成果Table 3 Characteristic values of ENSO events and annual runoff anomalies at Tuotuohe and Zhimenda hydrological stations during 1960-2020

從表3可以看出,1960年以來(lái)出現(xiàn)過(guò)多次極強(qiáng)的ENSO暖事件(1982-05—1983-06、1997-05—1998-05、2015-03—2016-05)和極強(qiáng)的ENSO冷事件(1998-07—2001-03)。在極強(qiáng)暖事件發(fā)生當(dāng)年或次年,沱沱河站和直門達(dá)站徑流均在不同程度上表現(xiàn)出偏枯的態(tài)勢(shì);而在極強(qiáng)冷事件的發(fā)展階段(當(dāng)年),沱沱河站呈現(xiàn)出偏豐態(tài)勢(shì)、直門達(dá)站呈現(xiàn)偏枯態(tài)勢(shì),在衰減階段(次年)兩站均呈現(xiàn)偏豐或略偏豐態(tài)勢(shì)。對(duì)于1972-05—1973-03、1986-09—1988-02、1991-06—1992-07出現(xiàn)的強(qiáng)暖事件,兩站徑流基本上均呈現(xiàn)偏枯態(tài)勢(shì);對(duì)于1973-06—1974-07出現(xiàn)的強(qiáng)冷事件,兩站在發(fā)展階段均呈現(xiàn)偏枯態(tài)勢(shì),在衰減階段呈現(xiàn)略偏豐態(tài)勢(shì)。從上述分析可知,當(dāng)ENSO發(fā)生強(qiáng)暖事件和暖事件期間,長(zhǎng)江源區(qū)徑流基本上表現(xiàn)出偏枯的態(tài)勢(shì),而在強(qiáng)冷和冷事件期間,在發(fā)展期呈現(xiàn)偏枯、衰落期呈現(xiàn)偏豐的態(tài)勢(shì)。

同樣,對(duì)于中等、弱、極弱的暖事件和冷事件出現(xiàn)年份和沱沱河站、直門達(dá)站徑流變化規(guī)律進(jìn)行分析,可以得出以下結(jié)論:對(duì)于中等暖事件長(zhǎng)江源區(qū)兩站基本處于偏豐態(tài)勢(shì),對(duì)于弱和極弱暖事件兩站基本處于偏枯態(tài)勢(shì)。而對(duì)于中等、弱和極弱冷事件而言,根據(jù)上節(jié)結(jié)論,以2000年為分界點(diǎn):1960—2000年對(duì)于中等冷事件,在發(fā)展階段兩站徑流處于偏枯、衰落階段處于偏豐態(tài)勢(shì);對(duì)于弱和極弱冷事件,兩站徑流基本處于偏枯態(tài)勢(shì)。2000—2020年受長(zhǎng)江源區(qū)氣候變暖等因素影響,無(wú)論對(duì)于中等、弱還是極弱冷事件,基本上都呈現(xiàn)出偏豐態(tài)勢(shì)。

2.2.2 徑流豐枯變化對(duì)ENSO事件的響應(yīng)

根據(jù)ENSO事件發(fā)生年份,以2000年為節(jié)點(diǎn),分兩段(1960—2000年、2000—2020年)統(tǒng)計(jì)長(zhǎng)江源區(qū)沱沱河站、直門達(dá)站徑流的距平程度及豐枯概率。其中ENSO事件發(fā)生年份分為發(fā)展階段(當(dāng)年)和衰減階段(次年),豐枯程度按照距平正值(+)、距平負(fù)值(-)分為偏豐和偏枯,成果列于表4。

表4 ENSO事件發(fā)生年份沱沱河站、直門達(dá)站徑流豐枯概率統(tǒng)計(jì)Table 4 Probability of dry or abundant runoff at Tuotuo-he Station and Zhimenda Station in ENSO event years

從表4成果可以得出以下結(jié)論:

(1)1960—2000年期間,暖事件發(fā)生年份,無(wú)論是當(dāng)年還是次年,長(zhǎng)江源區(qū)沱沱河站、直門達(dá)站出現(xiàn)徑流偏枯的概率較高。對(duì)于暖事件發(fā)生當(dāng)年,兩站徑流偏枯的概率均為76.9%;次年偏枯的概率分別為84.6%、76.9%。

(2)1960—2000年期間,冷事件發(fā)生年份,兩站均呈現(xiàn)出在當(dāng)年偏枯、次年豐枯概率基本相當(dāng)?shù)膽B(tài)勢(shì)。兩站當(dāng)年偏枯的概率均為77.8%;次年豐枯概率分別為4∶5和5∶4,出現(xiàn)概率基本相當(dāng)。

(3)2000年以后,長(zhǎng)江源區(qū)徑流偏豐態(tài)勢(shì)尤為顯著。對(duì)于暖事件發(fā)生年份,沱沱河站表現(xiàn)為當(dāng)年偏豐(66.7%概率)、次年豐枯相當(dāng)(各50%概率),直門達(dá)站在當(dāng)年徑流仍偏枯(66.7%概率)、次年豐枯相當(dāng)(各50%概率)。對(duì)于冷事件發(fā)生年份,無(wú)論是當(dāng)年還是次年,兩站均表現(xiàn)為偏豐概率大于偏枯概率。

朱鈺等[14]研究發(fā)現(xiàn),對(duì)于地處青藏高原的三江源區(qū),降水較少的冬半年多伴隨著暖事件(El Nio)的出現(xiàn);而降水峰值多伴隨著冷事件(La Nia)的出現(xiàn)。薛殷宗[15]在研究青海湖水位變化規(guī)律時(shí)發(fā)現(xiàn),當(dāng)暖事件(El Nio)增強(qiáng)時(shí),青海湖地區(qū)降水量明顯減少;冷事件(La Nia)增強(qiáng)時(shí),降水增加趨勢(shì)明顯。上述研究成果與本研究結(jié)論基本吻合。

2.3 徑流與ENSO事件的多尺度相關(guān)特征及響應(yīng)機(jī)制

采用交叉小波變換進(jìn)一步分析長(zhǎng)江源區(qū)徑流與ENSO事件之間的多尺度相關(guān)特征,結(jié)果見(jiàn)圖3。交叉小波變換重點(diǎn)突出時(shí)頻域中高能量區(qū)的相互關(guān)系,小波相干則重點(diǎn)突出時(shí)頻域中低能量區(qū)的相互關(guān)系?，F(xiàn)有研究表明,ENSO事件存在2～8 a時(shí)間尺度上的年際周期變化[16],這也直接體現(xiàn)在ENSO事件與沱沱河站、直門達(dá)站的交叉小波功率譜上。由圖3可見(jiàn),徑流與ENSO事件兩者之間在整個(gè)1960—2020年時(shí)間軸上具有8～16個(gè)月時(shí)間尺度上的共振周期,其中1960—1970年時(shí)間軸上,兩者基本呈現(xiàn)同頻同位相變化態(tài)勢(shì);在1970—2000年共振周期位相關(guān)系復(fù)雜,既有同位相、也有反位相變化態(tài)勢(shì);2000年以后,兩者之間呈現(xiàn)近似的同位相變化態(tài)勢(shì)。綜上分析可知,在1970年以前和2000年以后,長(zhǎng)江源區(qū)徑流與ENSO事件之間的時(shí)頻結(jié)構(gòu)具有一定程度的正相關(guān)性。

圖3 沱沱河站和直門達(dá)站月徑流與ONI之間的交叉小波功率譜與小波凝聚譜Fig.3 Cross wavelet transform and wavelet coherence between monthly runoff and Oceanic Nino Index at Tuotuohe station and Zhimenda station

ENSO事件作為大尺度海洋與大氣交互作用的事件,可通過(guò)海-氣相互耦合直接影響熱帶大氣環(huán)流,并通過(guò)大氣遙相關(guān)影響青藏高原季風(fēng)產(chǎn)生影響,從而進(jìn)一步影響長(zhǎng)江源區(qū)降水和徑流的時(shí)空分布。ENSO事件發(fā)生時(shí),赤道中、東太平洋增暖,赤道太平洋西部與東部之間的溫差減小,減弱甚至破壞緯向Walker環(huán)流,導(dǎo)致赤道太平洋對(duì)流活躍區(qū)東移到中太平洋,西太平洋和印度洋因海溫下降、大氣對(duì)流活動(dòng)減弱,降水減少[17-18]。Walker環(huán)流和Hadley環(huán)流強(qiáng)度變化之間是相互制約的[19]。當(dāng)赤道東太平洋增暖時(shí),緯向Walker環(huán)流減弱,而經(jīng)向Hadley環(huán)流加強(qiáng),從而增加了向極地方向的熱量輸送;Hadley環(huán)流加強(qiáng)后,定常渦旋向極經(jīng)向動(dòng)量輸送加強(qiáng),中緯度西風(fēng)加強(qiáng),進(jìn)而副高減弱,并向東向南收縮,東亞夏季風(fēng)減弱[20]。中緯度地區(qū)西風(fēng)帶位置南壓,西太平洋夏秋季臺(tái)風(fēng)減少。王紹武和龔道溢在研究ENSO事件影響機(jī)制時(shí)[21-23],提出ENSO事件的影響模式:赤道中、東太平洋海溫低(高)→南方濤動(dòng)強(qiáng)(弱)→Walker環(huán)流強(qiáng)(弱)→Hadley環(huán)流西強(qiáng)東弱(東強(qiáng)西弱)→太平洋副熱帶高壓西強(qiáng)東弱(東強(qiáng)西弱)→西太平洋夏秋季臺(tái)風(fēng)不活躍(活躍)。

青藏高原河川徑流主要受夏季的降水影響。青藏高原的水汽來(lái)源有3條路徑:一是阿拉伯海和孟加拉灣,二是中國(guó)南海,三是中緯度西風(fēng)帶[24]。太平洋海溫異常遠(yuǎn)距離影響著南亞地區(qū)的降水以及青藏高原的水汽供給,尤其是印度洋偶極子(Indian Ocean Dipole,IOD)和ENSO事件[25]。正事件(即印度洋正偶極子、El Nio以及二者共同作用)會(huì)使得源于阿拉伯海、孟加拉灣、南海和東南季風(fēng)輸送到青藏高原的水汽輸送減弱,導(dǎo)致青藏高原降水量減少,而負(fù)事件(印度洋負(fù)偶極子、La Nia以及二者共同作用)會(huì)使得水汽輸送加強(qiáng)[26]。薛殷宗[15]研究發(fā)現(xiàn),ENSO事件發(fā)生后引起緯向和經(jīng)向大氣環(huán)流異常,改變青藏高原氣壓場(chǎng)和風(fēng)場(chǎng)的位置以及強(qiáng)度,使得海洋向高原熱量和水汽輸送產(chǎn)生變化,從而影響青藏高原降水。梁涵洲等[11]研究表明,ENSO事件主要有3條途徑顯著影響青藏高原夏季降水和溫度:第1條途徑直接通過(guò)改變Walker環(huán)流、Hadley環(huán)流和NWPAC異常響應(yīng),影響南亞和印度洋向青藏高原的水汽輸送;第2條途徑直接通過(guò)激發(fā)東亞-太平洋(EAP)遙相關(guān)型波列,使得青藏高原東北部夏季受EAP北部反氣旋影響,不利于降水發(fā)生;第3條途徑在La Nia發(fā)展年夏季印度降水增加,造成中西亞地區(qū)西北側(cè)對(duì)流層中高層異常反氣旋環(huán)流,通過(guò)Rossby頻散沿著夏季急流激發(fā)環(huán)球遙相關(guān)(Circumglobal Teleconnection,CGT),進(jìn)而影響整個(gè)高原降水與溫度。

3 結(jié)論與展望

針對(duì)近60 a來(lái)長(zhǎng)江源區(qū)徑流變化事實(shí),采用沱沱河站、直門達(dá)站1960—2020年實(shí)測(cè)徑流系列研究長(zhǎng)江源區(qū)年徑流的演變趨勢(shì),利用ONI指數(shù)分析近60 a來(lái)ENSO事件的強(qiáng)度與時(shí)間特征,研究?jī)烧緩搅髫S枯變化及其對(duì)ENSO事件的響應(yīng),探討ENSO事件對(duì)徑流的可能影響機(jī)制。通過(guò)上述研究得到以下結(jié)論:

(1)長(zhǎng)江源區(qū)沱沱河站與直門達(dá)站的徑流變化趨勢(shì)基本一致,即總體上均呈現(xiàn)震蕩上升的趨勢(shì),尤其是在2000年以后呈現(xiàn)出大幅度增加的趨勢(shì)。Mann-Kendall檢驗(yàn)結(jié)果表明,兩站在進(jìn)入2000年以后均發(fā)生突變,并呈現(xiàn)出顯著的上升趨勢(shì)。

(2)1960—2000年期間,暖事件發(fā)生年份,無(wú)論是當(dāng)年還是次年,長(zhǎng)江源區(qū)沱沱河站、直門達(dá)站出現(xiàn)徑流偏枯的概率較高。冷事件發(fā)生年份,兩站均呈現(xiàn)出在當(dāng)年偏枯、次年豐枯概率基本相當(dāng)?shù)膽B(tài)勢(shì)。2000年以后,長(zhǎng)江源區(qū)徑流偏豐態(tài)勢(shì)尤為顯著。對(duì)于暖事件發(fā)生年份,豐枯基本相當(dāng)。對(duì)于冷事件發(fā)生年份,無(wú)論是當(dāng)年還是次年,兩站均表現(xiàn)為偏豐概率大于偏枯概率。

(3)徑流與ENSO事件兩者之間在整個(gè)1960—2020年時(shí)間軸上具有8～16個(gè)月時(shí)間尺度上的共振周期,其中1960—1970年時(shí)間軸上,兩者基本呈現(xiàn)同頻同位相變化態(tài)勢(shì);在1970—2000年共振周期位相關(guān)系復(fù)雜,既有同位相、也有反位相變化態(tài)勢(shì);2000年以后,兩者之間呈現(xiàn)近似的同位相變化態(tài)勢(shì)。

(4)青藏高原的水汽來(lái)源有3條路徑,一是阿拉伯海和孟加拉灣,二是中國(guó)南海,三是中緯度西風(fēng)帶。ENSO事件發(fā)生后引起緯向和經(jīng)向大氣環(huán)流異常,改變青藏高原氣壓場(chǎng)和風(fēng)場(chǎng)的位置以及強(qiáng)度,使得海洋向高原熱量和水汽輸送產(chǎn)生變化,從而影響青藏高原降水。

ENSO事件與徑流變化之間的作用機(jī)制是十分復(fù)雜的,本文初步分析了長(zhǎng)江源區(qū)徑流變化及其與ENSO事件之間的響應(yīng)關(guān)系,并通過(guò)現(xiàn)有研究成果,初步梳理了ENSO事件對(duì)青藏高原氣候的可能影響機(jī)制。然而,在全球氣候變暖的背景下,氣溫升高對(duì)青藏高原降水和源區(qū)徑流的影響遠(yuǎn)大于ENSO事件對(duì)大氣的遙相關(guān)影響。一方面氣溫升高導(dǎo)致冰川和融雪補(bǔ)給增加,直接影響源區(qū)徑流;另一方面氣溫升高導(dǎo)致青藏高原熱低壓和夏季風(fēng)的加強(qiáng),強(qiáng)化了阿拉伯海和孟加拉灣水汽向青藏高原輸送,間接影響源區(qū)徑流[27]。本文研究結(jié)果也顯示,2000年以后長(zhǎng)江源區(qū)徑流偏豐態(tài)勢(shì)尤為顯著,其豐枯變化機(jī)制并不遵循2000年以前的變化規(guī)律,也從側(cè)面印證長(zhǎng)江源區(qū)作為青藏高原腹地對(duì)全球氣候變暖的響應(yīng)最為敏感。今后還將進(jìn)一步深入研究全球氣候變化對(duì)青藏高原水循環(huán)和徑流形成機(jī)制。