張瑞，汪亞平，潘少明

(1. 淮海工學(xué)院 測繪工程學(xué)院，江蘇 連云港 222005；2. 南京大學(xué) 海岸與海島開發(fā)教育部重點實驗室，江蘇 南京，210093；3. 江蘇省海洋資源開發(fā)研究院，江蘇 連云港 222001)

近50年來長江入海徑流量對太平洋年代際震蕩變化的響應(yīng)

張瑞1,2,3，汪亞平2，潘少明2

(1. 淮海工學(xué)院 測繪工程學(xué)院，江蘇 連云港 222005；2. 南京大學(xué) 海岸與海島開發(fā)教育部重點實驗室，江蘇 南京，210093；3. 江蘇省海洋資源開發(fā)研究院，江蘇 連云港 222001)

太平洋年代際震蕩（Pacific Decadal Oscillation）是類似于ENSO型的具有年代際時間尺度變化的太平洋氣候變率。PDO既對長期的氣候變化趨勢產(chǎn)生擾動，又會對年際變化事件（如ENSO）產(chǎn)生重要影響，可使ENSO事件的頻率和強度改變，以致影響到ENSO和季風(fēng)的關(guān)系。通過對百年來的PDO數(shù)據(jù)進(jìn)行HHT分析，得到8個IMF和1個趨勢項，這表明PDO存在顯著的15 a、30 a、60 a的年代際變化，并且還有準(zhǔn)2 a、6 a的年際波動。同時對近50年來的長江入海月均徑流量序列也進(jìn)行HHT分析，得到7個IMF和1個趨勢項，結(jié)果表明長江入海月均徑流量序列存在有明顯的季節(jié)變化，顯著的準(zhǔn)2 a、4 a、8 a的年際變率，和顯著的16 a、22 a的年代際變率。并且對近50年來的長江入海月均徑流量序列和PDO序列的EMD分量進(jìn)行相關(guān)性分析，分析得知兩者在年代際變化中有較高的相關(guān)性，這表明在長時間尺度的年代際變化中PDO對長江入海徑流有著顯著地影響。

PDO；EMD；徑流量；長江

太平洋年代際震蕩（Pacific Decadal Oscilla--tion）是類似于ENSO型的具有年代際變化的太平洋氣候變率模式，并與ENSO一起對區(qū)域氣候產(chǎn)生顯著的影響[1,2]。PDO指數(shù)通常用北太平洋的 SST的經(jīng)驗正交函數(shù)（EOF）的第一主分量來表示[3,4]。PDO與 ENSO表現(xiàn)形式相似，但是時間尺度卻不同，典型的PDO事件一般持續(xù)時間為20～30年，而ENSO事件持續(xù)的時間尺度為6～18個月；并且在影響的空間區(qū)域上也有所不同，PDO事件主要影響北太平洋和北美地區(qū)，次要影響地區(qū)是熱帶太平洋[1-4]，而ENSO正好相反。PDO可分為冷、暖位相，在PDO暖位相，PDO指數(shù)為正值，熱帶中東太平洋異常暖，北太平洋中部異常冷，而沿北美西岸卻異常暖；反之，則為PDO冷位相。一般認(rèn)為，1890年以來PDO經(jīng)歷了3次顯著的突變，突變點分別是1925年、1947年和1976年，即1890－1925年和1947－1976年為冷位相階段；1926－1946年和 1977年至今為暖位相階段[1,5,6]。PDO對北太平洋和北美地區(qū)的氣候、水文、生態(tài)系統(tǒng)以及漁業(yè)生產(chǎn)產(chǎn)生重要的影響[1,7-10]，并且與ENSO共同對我國的氣候變化產(chǎn)生影響[3,4,11-15]。前人對長江入海徑流的變化特征及其可能的影響因素有不少有益的探討[16-19]，而 PDO對長江入海徑流量的影響研究較少，但是其影響程度值得進(jìn)一步研究。本文利用近50年的長江入海徑流量和近百年的PDO數(shù)據(jù)，應(yīng)用HHT方法來探討長江入海徑流量對PDO的響應(yīng)特征。

1 研究區(qū)域材料和方法

長江是中國第一大河、世界第三大河，全長6 300 km，流域總面積達(dá)1.94×106km2。大通水文站是海洋潮汐上漲所能到達(dá)的河流上界，其下游為感潮河段，它也是長江流域最下游的一個具有長期觀測資料的水文站，控制著 94% 的流域面積[20]。長江流域位于24.5～35.5oN的亞熱帶，氣候濕潤，多年平均降水量1 100 mm a-1左右，約50% 的降水轉(zhuǎn)換為徑流量入海[20]。長江入海徑流量巨大，以大通站為例，多年（1923－2004）平均徑流量為9 156×108m3·a-1。由于受東亞季風(fēng)的影響，長江入海徑流量表現(xiàn)出顯著的季節(jié)性變化特征[20-22]。受氣候變化的影響，長江流域的降水量和徑流量年際變化十分明顯，存在從2－3 a到約35－40 a的多種周期[18]。

2 數(shù)據(jù)資料和方法

2.1 數(shù)據(jù)資料

長江入海徑流量數(shù)據(jù)采用長江大通水文站1950－2004年共計55年的月均流量資料（其中1950－1987年流量來自中華人民共和國水文年鑒[23]）。PDO指數(shù)采用1900－2006年北太平洋20oN以北的月平均SST進(jìn)行經(jīng)驗正交函數(shù)分析（EOF）所得第一模態(tài)的時間系數(shù)。

2.2 分析方法

本文采用希爾伯特－黃變換（Hilbert-Huang Transform，簡稱 HHT）[24]來分析長江入海徑流量與太平洋年代際震蕩的相互關(guān)系，HHT方法的核心部分稱為經(jīng)驗?zāi)B(tài)分解（Empirical Mode Decomposition，簡稱EMD），該方法可以將非線性、非平穩(wěn)過程的數(shù)據(jù)進(jìn)行線性化和平穩(wěn)化的處理，同時保留數(shù)據(jù)本身的特點，具體步驟見參考文獻(xiàn)[24]。數(shù)據(jù)邊界問題是 HHT方法進(jìn)行樣條插值的關(guān)鍵所在，本文采用的端點抑制方法是在端點附近數(shù)據(jù)變化的“平衡位置”附加上兩條平行線段，再進(jìn)行外延而成[25]。由于兩條直線段上的點即為極大值，也是極小值，使得端點處的上、下包絡(luò)線向直線段集中，從而實現(xiàn)抑制。此抑制方法，也解決了Hilbert變換中需要處理的端點效應(yīng)。

3 結(jié) 果

3.1 Hilbert-Huang變換

3.1.1 長江入海徑流量的 Hilbert-Huang變換經(jīng)過EMD分解后長江入海徑流包含有7個內(nèi)在模函數(shù)（圖1a～g）和一個趨勢項Re（圖1 h）。內(nèi)在模函數(shù)經(jīng)過Hilbert變換后，得到了頻率、周期、振幅等統(tǒng)計結(jié)果（表1）。其中第一個內(nèi)在模函數(shù)IMF1是振幅最大，頻率最高，波長最短的一個波動，依次下去的其他內(nèi)在模函數(shù)振幅逐漸變小，頻率逐漸變低 ，波長逐漸變大。第一個模態(tài)到第七個模態(tài)的平均周期有短變長。第一模態(tài)IMF1的平均周期是0.30 a，這表明長江入海徑流存在明顯的年內(nèi)變化；第二模態(tài)IMF2至第四模態(tài)IMF4的平均周期分別是1.01 a 、2.08 a和4.26 a，這表明長江入海徑流具有1～4年的年際變化；第五模態(tài)IMF5至第七模態(tài)IMF7的平均周期分別是8.01 a、15.02 a和20.62 a，表明長江入海徑流存在著長時間尺度的年代際波動。圖 1中（h）是長江入海徑流經(jīng)過經(jīng)驗?zāi)B(tài)分解后得到的趨勢項Re，趨勢項Re顯示出近50年來長江入海徑流呈輕微升高的趨勢。

圖2是各個內(nèi)在模函數(shù)的Hilbert譜，像點以色標(biāo)來表示能量的大小，Hilbert譜右邊的標(biāo)尺表示從下向上能量從弱到強。從圖2可以看出主要的波動能量集中在頻率0.001～0.1范圍內(nèi)，這說明年際變化和年代際變化是長江入海徑流存在的重要變化。除此之外，還有一部分較高的能值集中在頻率0.25～0.45范圍內(nèi)，即表明年內(nèi)的季節(jié)變化也是長江入海徑流具有的重要變化。

圖 1 長江入海徑流的經(jīng)驗?zāi)B(tài)分解.(a)～(g)7個內(nèi)在模函數(shù)(IMF1～7)，(h)趨勢項ReFig. 1 Empirical mode decomposition of the water discharges from Changjiang River. (a)～(g) are seven IMF, and (h) is residue Re

表 1 長江入海徑流量Hilbert變換后內(nèi)在模函數(shù)的統(tǒng)計特征值Tab. 1 Statistic characteristics of IMF of the water discharges from Changjiang River after Hilbert transform

3.1.2 PDO的Hilbert-Huang變換 經(jīng)過EMD分解后可以得出，PDO時間系列包含有8個內(nèi)在模函數(shù)(圖 3(b)～(i))和一個趨勢項 Re(圖 3(j))，其中圖3(a) 是PDO時間系列圖。內(nèi)在模函數(shù)經(jīng)過Hilbert變換后，得到了頻率、周期、振幅等統(tǒng)計結(jié)果(表2)。

圖 2 長江入海徑流量的內(nèi)在模函數(shù)的Hilbert譜Fig. 2 Hilbert spectrum of IMF of the water discharges from Changjiang River

圖 3 PDO的經(jīng)驗?zāi)B(tài)分解.(a) PDO指數(shù)(b)～(i)8個內(nèi)在模函數(shù)(IMF1～8)，(j)趨勢項 ReFig. 3 Empirical mode decomposition of PDO. (a)PDO Index, (b)～(i) are eight IMF, and (j) is residue Re

圖 4 PDO指數(shù)內(nèi)在模函數(shù)的Hilbert譜Fig. 4 Hilbert spectrum of IMF of the PDO Index

圖4是PDO指數(shù)內(nèi)在模函數(shù)的Hilbert譜。圖中的像點以色標(biāo)來表示能量的大小，Hilbert譜右邊處也有集中分布，還有一部分低能值分布在0.25～0.45范圍內(nèi)。這說明10年以上的年代際波動是PDO變化的主要周期。

3.2 長江入海徑流量與PDO相關(guān)性

1950－2004年的長江入海徑流量與PDO的相關(guān)系數(shù)是 0.60，說明長江入海徑流與PDO之間通過海氣相互作用而緊密聯(lián)系。為了深入探討它們之間的相關(guān)性，將長江入海徑流量與 PDO的經(jīng)驗?zāi)B(tài)分解后的內(nèi)在模函數(shù)進(jìn)行相關(guān)性分析。表 3是1950－2004年的長江入海徑流量與同期PDO指數(shù)的各內(nèi)在模函數(shù)的相關(guān)系數(shù)。從表3中可以看出兩者在高頻高能模態(tài)上相關(guān)性較差，只有長江入海徑流量和 PDO的第二個模態(tài)的相關(guān)性顯著，其相關(guān)系數(shù)為0.105 （p＜0.01），而其他高頻模態(tài)的相關(guān)系數(shù)接近為零。在低頻模態(tài)上，兩者的相關(guān)性較高，第四、第五模態(tài)代表兩者的年際波動，其相關(guān)系數(shù)分別是0.358、0.301。PDO的第六、第七、第八模態(tài)和長江入海流量的第六、第七模態(tài)代表了兩者的年代際波動，兩個第六、第七模態(tài)相關(guān)系數(shù)分別是-0.334、-0.688，PDO的第八模態(tài)和長江入海徑流量的第七模態(tài)相關(guān)系數(shù)為0.275；并且PDO的第六模態(tài)和長江入海徑流量的第七模態(tài)相關(guān)系數(shù)為-0.452，PDO的第七模態(tài)和長江入海徑流量的第六模態(tài)相關(guān)系數(shù)為 0.462。將兩者的趨勢項也進(jìn)行相關(guān)性分析發(fā)現(xiàn)，兩者的第七模態(tài)的相關(guān)系數(shù)為-0.688，PDO的低頻模態(tài)和趨勢項與長江入海徑流量的趨勢項之間的相關(guān)系數(shù)分別是-0.325、0.648、0.495，說明PDO和長江入海徑流量波動趨勢上是有一定聯(lián)系的。

4 討 論

長江流域地處亞熱帶氣候區(qū)，氣候濕潤，雨量豐沛，對PDO、ENSO等天氣系統(tǒng)異常響應(yīng)敏感。降水是影響徑流的直接的天氣因素。有研究表明，長江流域的徑流與降水的相關(guān)系數(shù)達(dá)到了 0.835，說明徑流與降水的關(guān)系十分密切[26]。東亞季風(fēng)是影響我國東部降水最直接的因子之一，尤其是東亞夏季風(fēng)是我國東部地區(qū)夏季降水的最主要的水汽來源，東亞和西北太平洋地區(qū)的氣候（降水、溫度等）都存在極明顯的 10 年時間尺度的準(zhǔn)周期振蕩現(xiàn)象，其周期為11～16 a；并且大氣環(huán)流演變也有明顯的10年時間尺度的準(zhǔn)周期振蕩，其周期為12～18 a[27]。PDO與東亞大氣環(huán)流及中國氣候年代際變化關(guān)系密切，PDO暖位相期冬季江淮及長江流域降水偏少，夏季長江流域大部分地區(qū)降水異常偏多，當(dāng)PDO冷位相期時，正好相反[3]。由此看來，我國氣候變率和東亞大氣環(huán)流異常與 PDO的異常波動是有直接關(guān)系的。不僅如此，PDO還對ENSO事件對我國夏季氣候影響具有調(diào)制作用，但不管 PDO處于暖位相還是冷位相，ENSO處于發(fā)展階段還是衰減階段，都會使得長江流域降水異常偏多[3]。同時，PDO也對我國其他區(qū)域如華北、東北、華南等地區(qū)產(chǎn)生重要影響。PDO通過遙相關(guān)影響東亞夏季風(fēng)，進(jìn)而來影響我國華北地區(qū)的降水，造成華北夏季降水異常具有年代際震蕩的特征[14]。長江中下游梅雨的年代際變化與PDO的震蕩也有一定關(guān)系，當(dāng)PDO處于暖位相期時，長江中下游易出現(xiàn)多梅雨的趨勢，反之，長江中下游易出現(xiàn)少梅雨的趨勢[28]。

經(jīng)過HHT分析，PDO有準(zhǔn)15 a、30 a和60 a的年代際波動，而長江流域降水具有準(zhǔn)14 a、22 a、35～40 a的周期[18]，這說明PDO通過海氣相互作用對我國的天氣系統(tǒng)、大氣環(huán)流產(chǎn)生重要影響，使得降水也具有了 PDO的年代際震蕩的特征。而對于長江流域而言，徑流深受降水的影響，因此降水將自身具有的年代際震蕩特征傳遞給了徑流，因此長江入海徑流量也表現(xiàn)出了準(zhǔn)16 a和22 a的年代際波動特征。分析得知，總體上長江入海徑流量的趨勢項呈現(xiàn)出輕微上升的趨勢，而 PDO的趨勢項也表現(xiàn)出上升趨勢，并且兩者的相關(guān)性較高，這可能由于近20 a來的PDO暖位相與東亞大氣環(huán)流的年代際變化同步，由于 PDO的暖位相北太平洋為持續(xù)負(fù)的海溫距平，因此太平洋的負(fù)海溫距平對近距離的鄂海高壓生成可能有很大的貢獻(xiàn)，容易形成阻塞高壓，從而對長江流域的梅雨鋒的穩(wěn)定起較大作用，使得長江中下游的夏季降水和梅雨呈增長的趨勢[29]。因此使得長江入海徑流量在大尺度上的波動特征與PDO是相同的。

人類活動也對長江入海徑流變化產(chǎn)生一定影響，而這種影響可能主要表現(xiàn)在在2個方面，一方面有些人類活動（如毀林開荒）可能會使得徑流在短時間尺度上出現(xiàn)突發(fā)性變化，另一方面有些人類活動還可能影響徑流的趨勢性變化。修建水庫等人為活動可以減小徑流的變異，并且由于水庫的調(diào)蓄功能還可能削弱徑流與降水的相關(guān)性，弱化了氣候變化對長江入海徑流的影響。而近幾十年來的跨流域調(diào)（如南水北調(diào)）、生產(chǎn)生活的取水量增加等人類活動在一定程度上會減少長江入海徑流量，并且還可能使得長江入海徑流呈現(xiàn)下降趨勢，而降水量卻沒有表現(xiàn)出顯著的上升或下降趨勢[18]?？傊祟惢顒訒谝欢ǔ潭壬细蓴_長江入海徑流量對PDO的響應(yīng)，但是這種干擾機制有待進(jìn)一步探討。

表3 1950－2004年的長江入海徑流與PDO指數(shù)的各內(nèi)在模函數(shù)的相關(guān)系數(shù)（*表示p＜0.01）Tab. 3 Coefficient correlation of IMFs and the residues between the discharge and PDO from 1950 to 2004(*denotes p＜0.01)

5 討 論

本文采用HHT方法對1950－2004年的長江月均入海徑流量序列和1900－2006年的PDO序列進(jìn)行分析，研究結(jié)果顯示：

（1）對1950－2004年的長江入海月均徑流量序列進(jìn)行HHT分析得知，長江入海徑流具有7個IMF和1個趨勢項。結(jié)果表明，1950－2004年的長江入海月均徑流量序列存在顯著的準(zhǔn)2 a、4 a、8 a的年際變化，和16 a、22 a的年代際變化，并且還存在較明顯的季節(jié)變化。

（2） 對1900－2006年的PDO序列進(jìn)行HHT分析得知，近百年的PDO序列具有8個IMF和1個趨勢項，結(jié)果表明，近百年的 PDO序列存在顯著的15 a、30 a、60 a的年代際變化和準(zhǔn)2 a、4 a、6 a的年際波動。

（3） 通過分別對兩者的IMF和趨勢項進(jìn)行相關(guān)分析得知，在代表年代際震蕩的低頻模態(tài)上兩者具有良好的相關(guān)性，這表明在長時間尺度的年代際周期變化上 PDO是影響長江入海徑流量的重要因素。

[1] Mantua N J, Hare S R, Zhang Y, et al. A Pacific interdecadal climate oscillation with impacts on salmon production [J]. Bulletin of the American Meteorological Society, 1997, 78(6): 1069-1079.

[2] Zhang Y, Battisi J M. ENSO-like interdecadal variability, 1900-1993[J]. Journal of Climate, 1997, 10: 1004-1020.

[3] 朱益民, 楊修群. 太平洋年代際振蕩與中國氣候變率的聯(lián)系 [J].氣象學(xué)報, 2003, 61(6): 641-654.

[4] Yang H J, Zhang Q. On the decadal and interdecadal variability in the Pacific Ocean [J]. Advances in Atmospheric Sciences, 2003,20(2): 173-184.

[5] Minobe S. A 50-70 year climate oscillation over the North Pacific and North America [J]. Geophysical Research Letters, 1997, 24:683-686.

[6] Hare S R, Mantua N J. Empirical evidence for North Pacific regime shifts in 1977 and 1989 [J]. Process in Oceanography, 2000, 47:103-145.

[7] Bitz C M, Battisti D S. Interannual to decadal variability in climate and glacier mass balance in Washington, Western Canada, and Alaska [J]. Journal of Climate, 1999, 12: 3181-3196.

[8] Benson L, Linsley B, Smoot J, et al. Influence of the pacific decadal oscillation on the climate of the Sierra Nevada, California and Nevada [J]. Quaternary Research, 2003, 59(2): 151-159.

[9] Neal E G, Walter M T, Coffeen C. Linking the pacific decadal oscillation to seasonal stream discharge patterns in Southeast Alaska[J]. Journal of Hydrology, 2002, 263: 188-197.

[10] Rita V A, Mary T K. ENSO-related rainfall anomalies in South America and associated circulation features during warm and cold pacific decadal oscillation regimes [J]. International Journal of Climatology, 2005, 25: 2017-2030.

[11] 于淑秋, 林學(xué)椿. 北太平洋海溫的氣候躍度及其對中國汛期降水的影響 [J]. 熱帶氣象學(xué)報, 1997, 13(3): 265-275.

[12] 王東曉, 謝強, 劉赟, 等. 太平洋年代際海洋變率研究進(jìn)展 [J].熱帶海洋學(xué)報, 2003, 22(1): 76-83.

[13] 秦劍, 田永麗, 任菊章, 等. 太平洋年代際濤動與云南夏季氣溫的年代際變化 [J]. 高原氣象, 2004, 23: 69-76.

[14] 楊修群, 謝倩, 朱益民, 等. 華北降水年代際變化特征及相關(guān)的海氣異常型 [J]. 地球物理學(xué)報, 2005, 48(4): 789-797.

[15] 紀(jì)忠萍, 熊亞麗, 谷德軍, 等. 廣東汛期開始日期的年際和年代際變化研究 [J]. 大氣科學(xué), 2005, 29(2): 292-300.

[16] Zhang Q, Jiang T, Marco G, et al. Precipitation, temperature and discharge analysis from 1951 to 2002 in the Yangtze Catchment,China [J]. Hydrological Sciences Journal, 2005, 50(1): 65-80.

[17] 張強, 陳桂亞, 許崇育, 等. 長江流域水沙周期特征及可能影響原因 [J]. 水科學(xué)進(jìn)展, 2009, 20(1)：80-85.

[18] Yang S L, Liu Z, Dai S B, et al. Temporal variations in water resources in the Yangtze River (Changjiang) over the Industrial Period, based on reconstruction of missing monthly discharges [J].Water Resources Research. 2010. 46: 1-13.

[19] 張瑞, 汪亞平, 潘少明. 長江大通水文站徑流量的時間系列分析[J]. 南京大學(xué)學(xué)報(自然科學(xué)), 2006, 42(4)：423-434.

[20] Yang S L, Zhao Q Y, Belkin I M. Temporal variation in the sediment load of the Yangtze River and the influences of the human activities[J]. Journal of Hydrology. 2002, 263: 56-71.

[21] 張瑞, 汪亞平, 潘少明. 近 50 年來長江入河口區(qū)含沙量和輸沙量的變化趨勢 [J]. 海洋通報, 2008, 27(2):1-9.

[22] Xu K H, Milliman J D. Seasonal variations of sediment discharge from the Yangtze River before and after impoundment of the Three Gorges Dam [J]. Geomorphology, 2009, 104(3-4): 276-283.

[23] 中華人民共和國水利部. 中華人民共和國水文年鑒 [M]. 北京:中國水利電力出版社, 1950-1987: 11.

[24] Norden E H, Shen Z, Long S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis [J]. Proceedings of the Royal Society A 1998, 454(1971): 899-955.

[25] 熊學(xué)軍, 郭炳火, 胡莜敏, 等. EMD方法和Hilbert譜分析法的應(yīng)用與探討 [J]. 黃渤海海洋, 2002, 20(2): 12-21.

[26] Yang S L, Gao A, Hotz H J, et al. Trends in annual discharge from the Yangtze River to the sea (1865–2004) [J]. Hydrological Sciences Journal, 2005, 50: 825–836.

[27] 李崇銀, 廖清海. 東亞和西北太平洋地區(qū)氣候的準(zhǔn) 10年尺度振蕩及其可能機制. 氣候與環(huán)境研究 [J]. 1996, 1(2): 124-133.

[28] 魏鳳英, 宋巧云. 全球海表溫度年代際尺度的空間分布及其對長江中下游梅雨的影響 [J]. 氣象學(xué)報, 2005, 63(4): 477-485.

[29] 王亞非, 高橋清利. 長江中下游降水以及東亞夏季環(huán)流的年代際變化 [J]. 熱帶氣象學(xué)報, 2005, 21(4): 351-358.

Response of discharges from Changjiang River on Pacific Decadal Oscillation in the past 50 years

ZHANG Rui1,2,3, WANG Ya-ping2, PAN Shao-ming2

(1. School of Geodesy and Geomatics Engineering, Huaihai Institute of Technology, Lianyungang 222005, China; 2. MOE Key Laboratory of Coast and Island Development, Nanjing University, Nanjing 210093, China; 3. Jiangsu Marine Resources Development Research Institute, Lianyungang 222001, China)

Pacific Decadal Oscillation (PDO) is a long-term ENSO-like variability of the north Pacific. The index of PDO is denoted by the first principal component of EOF of the north Pacific SST. The PDO not only disturb the long-live trend of climate changes, but also influences significantly annual climate change (e.g. ENSO). The relationship between ENSO and monsoon also is influenced by PDO. Hilbert-Huang Transform (HHT) is used to analyze PDO over one hundred years to obtain 8 IMF and 1 residual. The results show that PDO has the quasi-2 and 6-year variability as well as significant quasi-15 and 30 and 60-year multi-decadal timescales. The discharges from Changjiang River in the last 50 years exit 7 IMF and 1 residual, which indicate that water discharges have seasonal change, quasi-2 and 4 and 8-year variability as well as significant 16 and 22-year multi-decadal variability. It is found that the good correlations between the discharges and PDO exit on the low frequency of long cycle multi-decadal timescales. Moreover, the relationship of both residuals is also significant. It is indicated that PDO influences water discharge from Changjiang River on long-term multi-decadal timescales.

PDO; EMD; Discharge; Changjiang River

TV148

A

1001-6932(2011)05-0572-06

2010-12-29；

2011-03-21

國家重點基礎(chǔ)研究發(fā)展規(guī)劃項目（2002CB412401）；國家自然科學(xué)基金（40106009、40576040）；淮海工學(xué)院科研啟動基金（KQ09041）。

張瑞（1982－），男，山東濟(jì)南人，講師，博士，主要從事陸海相互作用、河口水環(huán)境、近海沉積動力學(xué)等方面的研究。電子郵箱：rzhang_838@163.com。