康　磊, 劉世榮,*, 劉憲釗

1 中國(guó)林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所，北京　100091

2 中國(guó)林業(yè)科學(xué)研究院資源信息研究所，北京　100091

岷江上游水文氣象因子多尺度周期性分析

康磊1, 劉世榮1,*, 劉憲釗2

1 中國(guó)林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所，北京100091

2 中國(guó)林業(yè)科學(xué)研究院資源信息研究所，北京100091

摘要:采用Dmey小波變換法對(duì)岷江上游雜谷腦流域1959年至2006年月徑流量、月平均氣溫和月降水量不同時(shí)間尺度下的變化周期進(jìn)行分析，探討三者在長(zhǎng)時(shí)間序列周期性變化中的相互響應(yīng)，并根據(jù)主周期預(yù)測(cè)未來(lái)氣溫、降水和徑流的變化趨勢(shì)，結(jié)果表明：氣溫、徑流和降水存在多尺度周期性變化，在不同的尺度周期中，表現(xiàn)出不同的冷暖、豐枯和干濕的振蕩規(guī)律，總體表現(xiàn)為由小尺度無(wú)明顯規(guī)律的劇烈振蕩向大尺度有明顯規(guī)律的振蕩變化。3個(gè)要素同以8—12個(gè)月的小尺度為周期劇烈振蕩，在較大時(shí)間尺度上，具有明顯規(guī)律振蕩變化的周期分別為氣溫500個(gè)月、徑流150個(gè)月和降水120個(gè)月。受森林砍伐的影響，研究區(qū)域在1962—1988年期間徑流對(duì)降水變化的周期性響應(yīng)遲鈍，而在1988—2006年期間基本同步。根據(jù)大尺度周期性波動(dòng)趨勢(shì)預(yù)測(cè)，未來(lái)十幾年研究區(qū)域處于偏暖的年代際背景下，未來(lái)6—7a為多雨期，但徑流量偏少。

關(guān)鍵詞：多尺度；小波分析；長(zhǎng)時(shí)間序列；岷江上游

氣候變化將改變?nèi)蛩难h(huán)的現(xiàn)狀，引起水資源在時(shí)空上的重新分配，并對(duì)降水、蒸散、徑流等造成直接影響[1- 6]。根據(jù)現(xiàn)有研究成果，氣候系統(tǒng)的復(fù)雜與多變不僅已由眾多研究所證實(shí)，其時(shí)間和空間變化上的區(qū)域性、關(guān)聯(lián)性、周期性和不確定性，以及氣候變化對(duì)生態(tài)、農(nóng)業(yè)生產(chǎn)、社會(huì)經(jīng)濟(jì)和可持續(xù)發(fā)展的影響，成為近年來(lái)國(guó)內(nèi)外研究的熱點(diǎn)問(wèn)題[7- 11]。

水文氣象變量的時(shí)間序列同樣具有尺度特征，且大時(shí)間尺度包含著小時(shí)間尺度，不同時(shí)間尺度隱含著不同的氣候變化規(guī)律和水資源變化趨勢(shì)。小波分析中的多分辨率分析方法可以在不同的時(shí)間尺度上分析水文氣象時(shí)間序列的周期性和演變趨勢(shì)，并可借助小波變換的時(shí)頻局部化優(yōu)勢(shì)準(zhǔn)確地找到時(shí)間序列的大小時(shí)間尺度和突變點(diǎn)所在的位置[12]。目前小波分析已廣泛用于氣溫變化分析、降水變化分析、多尺度分析等，并取得了一些成果。邵駿等對(duì)岷江上游年徑流變化進(jìn)行了多時(shí)間尺度分析，將小波分析理論應(yīng)用于水文序列的多尺度分析中，表明岷江上游的地表水資源具有多時(shí)間尺度變化，具有6—7a、13a、16a左右的周期變化，利用小波分析可以準(zhǔn)確分析它的變化階段和序列特征[13]。王文圣等利用Marr小波分析作為母函數(shù)分析了長(zhǎng)江宜昌站98a的年平均流量資料，并根據(jù)小波變換系數(shù)分析了在4a、16a和32a 3個(gè)時(shí)間尺度上周期變化以及豐水期和枯水期交替變化對(duì)應(yīng)的突變點(diǎn)[14]。樊高峰等應(yīng)用小波分析對(duì)杭州市1951年以來(lái)的夏季氣溫進(jìn)行分析，發(fā)現(xiàn)杭州夏季平均氣溫發(fā)生過(guò)4次轉(zhuǎn)折，每次年代際轉(zhuǎn)換中都伴隨著劇烈振蕩[15]。

由于農(nóng)業(yè)不斷發(fā)展，人口爆炸式增長(zhǎng)和不合理的資源開發(fā)利用，我國(guó)的水資源狀況日益惡化，并帶來(lái)一系列的生態(tài)和環(huán)境問(wèn)題。在中國(guó)許多地區(qū)，水資源已成為制約地區(qū)可持續(xù)發(fā)展的瓶頸，而水資源的合理利用和配置成為經(jīng)濟(jì)可持續(xù)發(fā)展的前提之一[16- 18]。岷江上游地處四川盆地向青藏高原的過(guò)渡區(qū)，是長(zhǎng)江上游重要的水源涵養(yǎng)生態(tài)功能區(qū)，也是我國(guó)一個(gè)重要的大尺度、復(fù)合型生態(tài)過(guò)渡帶和生態(tài)系統(tǒng)脆弱區(qū)，其自然環(huán)境的復(fù)雜性、生態(tài)系統(tǒng)的脆弱性、經(jīng)濟(jì)發(fā)展的邊緣性和社會(huì)文化的過(guò)渡性，在我國(guó)都具有典型的代表性。隨著岷江流域上游人口的增加和20 世紀(jì)中后葉的大規(guī)模天然林采伐, 該地區(qū)森林覆蓋率顯著下降, 森林植被嚴(yán)重退化。而氣候變化和土地利用格局的變化引起流域水文狀況明顯變化, 進(jìn)而影響岷江下游的水資源管理和農(nóng)業(yè)可持續(xù)發(fā)展[19- 20]。

本文旨在利用小波分析工具通過(guò)對(duì)1959—2006 年氣溫、降水和徑流時(shí)間序列變化規(guī)律進(jìn)行深入分析，揭示岷江上游流域水文氣象長(zhǎng)期變化規(guī)律和周期性特征，預(yù)測(cè)流域氣候變化趨勢(shì)，為全面認(rèn)識(shí)岷江上游氣候變化特征以及各氣象因子變化的響應(yīng)特征、科學(xué)制定水資源管理和天然林保護(hù)策略提供技術(shù)參考。

圖1　研究區(qū)的地理位置Fig.1　Geographical location of the study area

1研究區(qū)域概況

雜古腦流域是岷江上游地區(qū)一個(gè)重要流域，岷江就水量而言，是長(zhǎng)江上游最大的支流。岷江發(fā)源于岷山南麓的貢杠嶺和郎架山，由北向南流經(jīng)四川省阿壩州、都江堰市、成都市、樂(lè)山市、于宜賓市與金沙江匯合，流入長(zhǎng)江。雜古腦河是岷江上游的主要的支流之一，位于四川省理縣境內(nèi)(102°6′—103°6′E, 31°2′—31 °9′N),總長(zhǎng)約168 km，流域面積約4632 km2。雜古腦流域以雜古腦水文站(103°10′E, 31°26′N)為界又可分為兩個(gè)部分：雜古腦上游流域和雜古腦下游流域,流域出口分別為雜古腦水文站(103°10′E, 31°26′N)和桑坪水文站(103°35′E, 31°29′N)。鑒于雜古腦下游植被數(shù)據(jù)缺乏，本研究不涉及雜古腦下游流域，研究區(qū)域僅為雜古腦上游流域(圖1)。雜古腦上游主河道長(zhǎng)度為113 km,河寬介于4—36 m，平均河道比降為23 ‰，流域面積為2528 km2。在雜古腦上游流域設(shè)有雜古腦水文站、米亞羅雨量站、雜古腦雨量站和理縣氣象站等觀測(cè)臺(tái)站，雜古腦流域上游森林采伐和更新經(jīng)營(yíng)歷史在整個(gè)雜古腦流域具有非常典型的代表性，森林采伐和更新資料比較齊全。

雜古腦上游流域氣候受高原地形影響，為冬寒夏涼的高山氣候， 1月份月平均氣溫-8 ℃，7月份月平均氣溫12.6 ℃，≥10 ℃的年積溫為1200—1400 ℃。全年降水量在627.5—1478.0 mm，平均降水量為1067.6 mm。降水主要發(fā)生在5—7月。

2研究方法

氣象數(shù)據(jù)取自流域內(nèi)雜古腦水文站、米亞羅雨量站、雜古腦雨量站和理縣氣象站，雜谷腦流域徑流數(shù)據(jù)的收集包括雜谷腦水文站(31.26′E，103.10′N)1959年—2006年的日徑流量、日最大徑流量、日最小徑流量等實(shí)測(cè)徑流數(shù)據(jù)。本文利用小波分析，對(duì)流域內(nèi)月徑流量、月平均氣溫和月降水量的長(zhǎng)時(shí)間序列的周期性特征進(jìn)行分析。

小波是具有震蕩特性、能迅速衰減到零的一類函數(shù)，即∫RΨ(t)dt=0，由Ψ(t) 的伸縮和平移構(gòu)成一簇函數(shù)系：

式中,Ψa,b(t) 為子小波，a為尺度因子或頻率因子，b為時(shí)間因子或平移因子。對(duì)于能量有限信號(hào)f(t),其連續(xù)小波變換定義為：

任何一個(gè)時(shí)間序列或信號(hào)都由長(zhǎng)期趨勢(shì)、各種周期性成分以及隨機(jī)成分組成，水文氣象要素時(shí)間序列也不例外。通過(guò)小波變換系數(shù)和小波方差可以了解時(shí)間序列存在的周期性成分，小波方差可以由下式計(jì)算：

Var(a)=∫R|Wf(a,b)|2db

式中，Var(a)是指在尺度a下的小波方差，Wf(a,b)為小波變換系數(shù)。

小波方差隨尺度a變化過(guò)程稱小波方差圖。它反映了波動(dòng)的能量隨尺度的分布，根據(jù)峰值對(duì)應(yīng)的時(shí)間尺度可以確定一個(gè)時(shí)間序列中存在的主要周期成分。小波方差可以確定時(shí)間序列的周期性成分組成情況，確定整個(gè)時(shí)間序列中存在的各種主要周期。任何信號(hào)或者時(shí)間系列都可以用理想的低通濾波器和高通濾波器分解為高頻部分d1和低頻部分a1，其中低頻部分a1可以用低通濾波器和高通濾波器進(jìn)一步分解為低頻部分a2和高頻部分d2,同樣a2可以分解為低頻部分a3和高頻部分d3，如此重復(fù)進(jìn)行，可以獲得任意尺度上時(shí)間序列的概貌(低頻部分)和細(xì)節(jié)成分(高頻部分)。對(duì)于大多數(shù)信號(hào)來(lái)說(shuō)，低頻部分往往是最重要的，給出了信號(hào)的特征。

3結(jié)果分析

3.1流域月徑流量長(zhǎng)時(shí)間序列多尺度周期性變化特征

為分析雜古腦流域徑流動(dòng)態(tài)的周期性，對(duì)雜古腦流域標(biāo)準(zhǔn)化徑流時(shí)間序列進(jìn)行Dmey小波連續(xù)變換，得到各個(gè)尺度上的小波變換系數(shù)，在此基礎(chǔ)上計(jì)算小波方差(圖2)。由圖2可以看出雜古腦流域徑流，方差在8—12個(gè)月時(shí)達(dá)到最大，即雜古腦流域徑流主要周期為8—12個(gè)月。雜谷腦流域徑流時(shí)間序列在10個(gè)月尺度上的方差略高于12個(gè)月。在更大時(shí)間尺度上，雜谷腦流域徑流時(shí)間序列主要周期分別是40、65和150個(gè)月。這意味著用年徑流量來(lái)分析森林植被變化對(duì)雜古腦流域的徑流影響時(shí)所采用的時(shí)間尺度并不是最佳尺度。

圖2　雜谷腦徑流時(shí)間序列小波方差Fig.2　Wavelet transform variation of zagunao monthly runoff time serial

由圖3可根據(jù)小波變化系數(shù)確定不同時(shí)期的主導(dǎo)周期成分。圖中顏色越亮，對(duì)應(yīng)的時(shí)間尺度在此時(shí)間序列上越明顯。由圖3可以看出，8—12個(gè)月的時(shí)間尺度在整條時(shí)間序列上都很明顯，即8—12個(gè)月是整個(gè)研究時(shí)間序列的主導(dǎo)周期。更大時(shí)間尺度主導(dǎo)周期僅在部分時(shí)間序列中具有明顯的周期性變化，40個(gè)月和65個(gè)月時(shí)間尺度具有相似的周期性信息，在M100—M450(對(duì)應(yīng)1967—1996年)和M500—M570(對(duì)應(yīng)2000—2006年)區(qū)間上具有顯著的周期性。150個(gè)月尺度在M200—M480(對(duì)應(yīng)1975—1998年)區(qū)間上具有顯著周期性。

圖3　雜谷腦徑流時(shí)間序列徑流小波變換二維系數(shù)Fig.3　Wavelet transform Coefficients of Zagunao normalized monthly runoff time serial

利用Dmey小波變換函數(shù)在水平5(a5)、水平6(a6)和水平7(a7)(相應(yīng)的尺度為25、26和27)上將徑流時(shí)間序列進(jìn)行低通濾波，獲得雜古腦流域地表徑流在3個(gè)尺度上的低頻概貌成分(圖4)。從3個(gè)尺度上都可看出，雜谷腦流域徑流在1959年至1970年之間波動(dòng)幅度較小，趨于平緩，從1971年開始徑流波動(dòng)變大，且在1986年達(dá)到了最低值，即枯水極值。

圖4　雜古腦徑流時(shí)間序列多尺度分析及在尺度40個(gè)月、65個(gè)月和150個(gè)月上的小波變換系數(shù)Fig.4　Multi-resolution analysis and WT(40,b),WT(65,b) and WT(150,b) of Zagunao monthly runoff time serial

小波方差分析表明雜古腦流域徑流存在的主要周期(年尺度以上的)是150個(gè)月、65個(gè)月和40個(gè)月。為了更準(zhǔn)確地揭示雜古腦流域徑流動(dòng)態(tài)規(guī)律，利用Dmey小波對(duì)標(biāo)準(zhǔn)化時(shí)間序列進(jìn)行連續(xù)小波變換，得到雜谷腦流域3個(gè)主要時(shí)間尺度上(a=40，65和150)小波變換系數(shù)圖(圖4)。不同時(shí)間尺度下的小波變換系數(shù)可以反映不同時(shí)間尺度下徑流變化特征，正的小波變換系數(shù)對(duì)應(yīng)于豐水期，負(fù)的小波變換系數(shù)對(duì)應(yīng)于枯水期，小波變換系數(shù)絕對(duì)值越大，表明該時(shí)間尺度變化越明顯。在1959—1962年之間，主導(dǎo)周期為40個(gè)月，在1963—1980年主導(dǎo)周期為65個(gè)月，在1981—1998年主導(dǎo)周期為150個(gè)月。40個(gè)月與65個(gè)月具有相似的周期性波動(dòng)，隨著時(shí)間尺度的縮小，所包含的水文信息也逐漸增多。在150個(gè)月尺度上，雜谷腦流域徑流在1965—2006年期間共經(jīng)歷了7次豐枯交替，分別為1959 —1962年枯水期、1963—1968年豐水期、1969—1975年枯水期、1976—1982年豐水期、1983—1991年枯水期、1992—1998年豐水期和1999—2006年枯水期。在65個(gè)月尺度上則經(jīng)歷了15次豐枯交替。

3.2流域月降水量長(zhǎng)時(shí)間序列多尺度周期性變化特征

由圖5和圖6可知，雜谷腦流域降水量與徑流量相同，8—12個(gè)月依然是最主要周期。在更大尺度上的主要周期分別為20個(gè)月、40個(gè)月和120個(gè)月。20個(gè)月尺度的周期性主要發(fā)生在M110—M180(對(duì)應(yīng)1968—1973年)、M400—M440(對(duì)應(yīng)1992—1995年)和M470—M520(對(duì)應(yīng)1998—2002年)。40個(gè)月尺度的周期性主要發(fā)生在1973—1979年、1986—1992年和2000—2006年，120個(gè)月尺度的周期性主要發(fā)生在1989—2006年。

圖5　雜谷腦降水量時(shí)間序列小波方差Fig.5　Wavelet transform variation of zagunao monthly precipitation time serial

在120個(gè)月時(shí)間尺度上(圖7)，降水量整個(gè)時(shí)間序列上呈現(xiàn)少雨期和多雨期的交替，少雨期為：1964—1971年、1977—1989年和1996—2002年。多雨期為：1959—1963年、1972—1976年、1990—1995年和2003—2006年。

圖7　雜古腦降水量時(shí)間序列多尺度分析及在尺度20個(gè)月、40個(gè)月和120個(gè)月上的小波變換系數(shù)Fig.7　Multi-resolution analysis and WT(20,b),WT(40,b) and WT(120,b) of Zagunao monthly precipitation time serial

3.3流域月平均氣溫長(zhǎng)時(shí)間序列多尺度周期性變化特征

由氣溫小波方差分析圖可知，在整個(gè)時(shí)間序列上，小波方差在10—12個(gè)月是達(dá)到最大，10個(gè)月的方差較12個(gè)月略大，由于雜谷腦流域?qū)儆诟吆貐^(qū)，春季較短，冬季和春季氣溫特點(diǎn)相似，導(dǎo)致10個(gè)月的周期性較12個(gè)月明顯。在較大尺度上，月平均氣溫的小波方差隨著時(shí)間尺度的增加逐漸變大，在500個(gè)月左右達(dá)到最大后開始變小，即500個(gè)月是較大時(shí)間尺度的主要周期。

由氣溫在60個(gè)月、120個(gè)月和500個(gè)月尺度上的小波變換曲線(圖9)可知，60個(gè)月和120個(gè)月尺度的周期性波動(dòng)情況大致相同，1990年之前氣溫波動(dòng)不大，較為平緩，1990年之后，氣溫波動(dòng)變大。小尺度的變化過(guò)程反映了大尺度背景下氣溫的詳細(xì)變化，而大尺度反映了氣溫變化的年代際背景。在500個(gè)月尺度上，小波變換曲線波動(dòng)較60個(gè)月和120個(gè)月明顯。從500個(gè)月尺度小波變換系數(shù)曲線可以看出，整個(gè)研究時(shí)間序列有兩個(gè)氣溫突變點(diǎn)，即1966年和1995年。1950—1966年為偏暖期，1966—1995年為偏冷期，1995—2011年為偏暖期。

圖9　雜古腦氣溫時(shí)間序列在尺度60個(gè)月、120個(gè)月和500個(gè)月上的小波變換系數(shù)Fig.9　WT(60,b),WT(120,b) and WT(500,b) of Zagunao monthly precipitation time serial

圖10　雜古腦氣溫、徑流及降水量時(shí)間序列在長(zhǎng)時(shí)間尺度的小波變換系數(shù)Fig.10　Long-time scale wavelet transform coefficients of Zagunao monthly temperature, runoff and precipitation time serial

4結(jié)論與討論

(1)本研究在總結(jié)已有研究基礎(chǔ)上，通過(guò)完善研究資料，從雜谷腦流域水文氣象因子整體考慮，系統(tǒng)分析徑流、氣溫、降水的長(zhǎng)時(shí)間序列的周期性變化及三者之間的相互響應(yīng)關(guān)系。本文在周期性分析中，使用的均為標(biāo)準(zhǔn)化后的時(shí)間序列，主要為避免在求算小波變化系數(shù)時(shí)出現(xiàn)大數(shù)吃小數(shù)及除小數(shù)的問(wèn)題[21]。

(2)雜谷腦流域月平均氣溫、月徑流量和月降水量的周期性變化過(guò)程中，存在著多重時(shí)間尺度上得復(fù)雜嵌套結(jié)構(gòu)，在不同的尺度周期中，表現(xiàn)出不同的冷暖、豐枯和干濕的振蕩規(guī)律，總體表現(xiàn)為由小尺度無(wú)明顯規(guī)律的劇烈振蕩向大尺度有明顯規(guī)律的振蕩變化。在較小時(shí)間尺度上，氣溫、徑流和降水均以8—12個(gè)月為主要周期高頻振蕩，12個(gè)月的周期性沒有特別突出。究其原因可能與雜谷腦流域特殊的氣溫和降水特點(diǎn)有關(guān)。雜谷腦流域具有平均海拔高、年平均氣溫低和積雪時(shí)間長(zhǎng)的特點(diǎn)，6月底冰雪面積比例為3.53%，因此，冰雪融化是該地區(qū)流域水文研究必須考慮的水文過(guò)程。氣溫對(duì)該流域徑流的直接影響主要反映在氣溫對(duì)融雪徑流的影響上，氣溫增加融雪徑流，可能是研究結(jié)果與其他研究不一致的主要原因。由于雜谷腦流域?qū)儆诟吆貐^(qū)，春季比較短，冬季和春季氣溫和降水特點(diǎn)比較相似，其中11—12月與1—4月都屬于非生長(zhǎng)季節(jié)，降水形式以降雪為主，這無(wú)疑將對(duì)徑流動(dòng)態(tài)的周期變化產(chǎn)生影響，從而形成與低海拔地區(qū)徑流動(dòng)態(tài)不同的周期變化規(guī)律。本文主要側(cè)重于大尺度下流域徑流的動(dòng)態(tài)變化機(jī)制，從而為雜谷腦流域乃至整個(gè)長(zhǎng)江上游的水資源管理和水土流失治理提供理論參考。

(3)同一氣象或水文要素隨主周期尺度的變化均呈現(xiàn)出長(zhǎng)周期主導(dǎo)下的長(zhǎng)短周期振蕩疊加，振蕩幅度隨主周期尺度增加逐漸減弱，這種多重尺度和幅度的周期振蕩特征與氣候過(guò)程的非線性、非平穩(wěn)性特征相一致。由大尺度下月徑流量、月降水量和月均氣溫的小波變換系數(shù)曲線(圖10)可知，1988—2006年，研究區(qū)的徑流量和降水量的周期性變化基本同步，而在1962—1988年徑流量對(duì)降水量變化的周期性響應(yīng)變得遲鈍，徑流波動(dòng)變化相對(duì)降水具有2—3a左右的滯后，究其原因，除了融雪徑流的影響外，1962—1978年，研究區(qū)森林被大規(guī)模砍伐，植被遭到嚴(yán)重破壞，造成土壤干旱，土壤水減少，降水產(chǎn)流量降低。另外，徑流量的低頻振蕩所反映的是總徑流中來(lái)自土壤中水徑流和存在于裂隙中的地下水徑流的周期變化，高頻振蕩反映的是總徑流中來(lái)自于地面徑流的周期性變化規(guī)律。研究區(qū)徑流量的年內(nèi)變化和年際高頻振蕩主要是降水量變化所引起的，而徑流量時(shí)間序列150個(gè)月周期變化是由于土壤和裂隙對(duì)降水量時(shí)間序列周期增長(zhǎng)和放大作用形成的。

(4)本文研究的雜谷腦流域3個(gè)水文氣象因子是其大時(shí)間尺度上的動(dòng)態(tài)機(jī)制，大尺度上生態(tài)過(guò)程和水文過(guò)程都有更好的可預(yù)測(cè)性和規(guī)律性，小尺度上許多局部異質(zhì)性和非線性在大尺度上將被平滑掉。依據(jù)主周期變化趨勢(shì)，在2011—2030 年，氣溫仍處于偏暖期的年代際背景下，但在此期間，氣溫會(huì)以2015 年和2021 年為突變點(diǎn)發(fā)生冷暖的周期性波動(dòng)。2011—2015 年為偏暖期，2016—2021為偏冷期，2022—2027年為偏暖期。研究區(qū)徑流在150個(gè)月周期波動(dòng)下，2014—2020年為枯水期。降水在120個(gè)月周期波動(dòng)下，從波動(dòng)曲線的未來(lái)走勢(shì)可預(yù)測(cè)，2014—2020 年為多雨期。

參考文獻(xiàn)(References):

[1] 石培禮, 李文華. 森林植被變化對(duì)水文過(guò)程和徑流的影響效應(yīng). 自然資源學(xué)報(bào), 2001, 16(5): 481- 487.

[2]姜彤, 蘇布達(dá), 王艷君, 張強(qiáng), 秦年秀, 施雅風(fēng). 四十年來(lái)長(zhǎng)江流域氣溫、降水與徑流變化趨勢(shì). 氣候變化研究進(jìn)展, 2005, 1(2): 65- 68.

[3]王順久. 長(zhǎng)江上游川江段氣溫、降水及徑流變化趨勢(shì)分析. 資源科學(xué), 2009, 31(7): 1142- 1149.

[4]李春強(qiáng), 杜毅光, 李保國(guó). 1965—2005年河北省降水量變化的小波分析. 地理科學(xué)進(jìn)展, 2010, 29(11): 1340- 1344.

[5]李瑞平, 史海濱, 李彰俊, 劉宏云. 連續(xù)小波變換在氣溫和降水變化分析中的應(yīng)用. 灌溉排水學(xué)報(bào), 2008, 27(1): 86- 89.

[6]賀偉, 布仁倉(cāng), 熊在平, 胡遠(yuǎn)滿. 1961—2005年?yáng)|北地區(qū)氣溫和降水變化趨勢(shì). 生態(tài)學(xué)報(bào), 2013, 33(2): 519- 531.

[7]張濟(jì)世, 康爾泗, 藍(lán)永超, 陳仁升, 姚盡忠, 蒲瑞豐, 陳滿祥. 50a來(lái)洮河流域降水徑流變化趨勢(shì)分析. 冰川凍土, 2003, 25(1): 77- 82.

[8]張遠(yuǎn)東. 川西亞高山森林恢復(fù)格局、林地水文效應(yīng)及經(jīng)營(yíng)規(guī)劃[D]. 北京: 中國(guó)林業(yè)科學(xué)研究院博士后出站報(bào)告, 2004.

[9]曾濤, 郝振純, 王加虎. 氣候變化對(duì)徑流影響的模擬. 冰川凍土, 2004, 26(3): 324- 332.

[10]朱道光, 蔡體久, 姚月峰, 琚存勇. 小興安嶺森林采伐對(duì)河川徑流的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2005, 16(12): 2259- 2262.

[11]劉俊萍, 田峰巍, 黃強(qiáng), 佟春生. 基于小波分析的黃河河川徑流變化規(guī)律研究. 自然科學(xué)進(jìn)展, 2003, 13(4): 383- 387.

[12]李崇巍. 岷江上游生態(tài)系統(tǒng)空間格局與生態(tài)水文功能模擬[D]. 北京: 北京師范大學(xué), 2005.

[13]邵駿, 袁鵬, 李秀峰. 岷江上游年徑流變化的多時(shí)間尺度分析. 東北水利水電, 2007, 25(6): 52- 54.

[14]王文圣, 向紅蓮, 黃偉軍, 丁晶. 基于連續(xù)小波變換的徑流分維研究. 水利學(xué)報(bào), 2005, 36(5): 598- 601.

[15]樊高峰, 苗長(zhǎng)明. 用小波分析方法診斷杭州近50a夏季氣溫變化. 氣象科學(xué), 2008, 28(4): 431- 434.

[16]馮亞文, 任國(guó)玉, 劉志雨, 吳吉東, 張雷. 長(zhǎng)江上游降水變化及其對(duì)徑流的影響. 資源科學(xué), 2013, 35(6): 1268- 1276.

[17]秦大河. 中國(guó)氣候與環(huán)境演變(上卷). 北京: 科學(xué)出版社, 2005.

[18]Guo J, Chen H, Xu C Y, Guo S L, Guo J L. Prediction of variability of precipitation in the Yangtze River Basin under the climate change conditions based on automated statistical downscaling. Stochastic Environmental Research and Risk Assessment, 2012, 26(2): 157- 176.

[19]夏軍, 談戈. 全球變化與水文科學(xué)新的進(jìn)展與挑戰(zhàn). 資源科學(xué), 2002, 24(3): 1- 7.

[20]李林, 王振宇, 秦寧生, 馬玉才. 長(zhǎng)江上游徑流量變化及其與影響因子關(guān)系分析. 自然資源學(xué)報(bào), 2004, 19(6): 694- 700.

[21]徐萃薇, 孫繩武. 計(jì)算方法引論(第二版). 北京: 高等教育出版社, 2002: 11- 15.

Multiresolution and periodicity analysis of hydrological and meteorological factors in upper reaches of Minjiang River

KANG Lei1, LIU Shirong1,*, LIU Xianzhao2

1ResearchInstituteofForestEcology,EnvironmentandProtection,ChineseAcademyofForestry,Beijing100091,China2ResearchInstituteofForestResourceInformationTechniques,ChineseAcademyofforestry,Beijing100091,China

Abstract：In this study, the Dmey transformation approach was used to analyze monthly runoff, mean temperature, and precipitation from 1959 to 2006 in the Zagunao basin in the upper reaches of the Minjiang River, to explore the mutual responses to periodic changes in these three variables and to attempt to predict the trends of future temperature, precipitation, and runoff. This study is focused on dynamic-change mechanisms of runoff on a large scale, and should provide helpful insights into water management and control of soil water losses in the Zagunao watershed. The results showed that temperature, runoff, and precipitation changed periodically in several scales, and the oscillation raw data on cold and warm temperatures, plentiful and scarce runoff, and wet and dry periods differed. Overall, the presentation changed from nonobvious regular violent oscillation on a small scale to obvious regular oscillation on a large scale. All three elements showed the small-scale periodic violent oscillation of 8 to 12 months. On a larger time scale, the periods of obvious regular oscillation of temperature consisted of 500 months, and the runoff 150 months, the precipitation 120 months. The Zagunao watershed is characterized by high average elevation, low annual average temperature, and a prolonged snow cover. The percentage of the snow area is 3.53% at the end of June. Therefore, snow melting is a hydrological process that must be taken into account in this area. Temperature is the main factor that influences snowmelt runoff. The latter might be the main reason for the differences between this study and others. Because the Zagunao watershed belongs to an alpine region that has a short spring season, the characteristics of temperature and precipitation in the spring are similar to those of winter. November to December and January to April are nongrowth seasons. The precipitation types are dominated by snow, which should undoubtedly have an impact on the dynamic and periodic changes of runoff, thereby resulting in periodic variation, which is different from that in a low-elevation area. From 1962 to 1988, the periodic response of runoff to precipitation changes was insensitive, but it changed and became synchronized with precipitation from 1988 to 2006. Deforestation and vegetation that were damaged in the study area between 1962 and 1978, as well as the low-frequency oscillation of runoff, might have resulted in soil drought and a reduction in soil water. The high-frequency oscillation of runoff was the reflection of the periodic changes in surface runoff. The intra-annual variation and interannual high-frequency oscillation of runoff in this area were caused by variation of precipitation. The focus of this research was the dynamic mechanisms of these three hydrological and meteorological factors on a large scale. Ecological and hydrological processes on a large scale showed better predictability and regularity. Some local heterogeneity and nonlinearity on a small scale are expected to be smoothed out on a large scale. According to the trend of large-scale periodic fluctuation, the study area should be in a warmer background era in the coming decades, and the precipitation is expected to be plentiful in the next 6—7 years; in contrast, the runoff is expected to be scanty.

Key Words:multiresolution; wavelet analysis; lang-time series; upper reaches of Minjiang River

基金項(xiàng)目:林業(yè)公益行業(yè)科研專項(xiàng)(201404201); 國(guó)家自然科學(xué)基金項(xiàng)目(31400538); 國(guó)家自然科學(xué)基金項(xiàng)目(31290223)

收稿日期:2014- 07- 04; 網(wǎng)絡(luò)出版日期:2015- 07- 22

DOI:10.5846/stxb201407041378

*通訊作者Corresponding author.E-mail: liusr@caf.ac.cn

康磊, 劉世榮, 劉憲釗.岷江上游水文氣象因子多尺度周期性分析.生態(tài)學(xué)報(bào),2016,36(5):1253- 1262.

Kang L, Liu S R, Liu X Z.Multiresolution and periodicity analysis of hydrological and meteorological factors in upper reaches of Minjiang River.Acta Ecologica Sinica,2016,36(5):1253- 1262.