張　兵，王中良（天津師范大學(xué) 天津市水資源與水環(huán)境重點(diǎn)實(shí)驗(yàn)室，天津 300387）

天津地區(qū)降水和氣溫的變化趨勢(shì)及多尺度交叉小波分析

張兵，王中良
（天津師范大學(xué) 天津市水資源與水環(huán)境重點(diǎn)實(shí)驗(yàn)室，天津 300387）

摘要：由于降水量和氣溫的變化對(duì)水循環(huán)研究與水資源規(guī)劃和利用具有重要意義，運(yùn)用線性傾向估計(jì)、Mann-Kendall檢驗(yàn)和交叉小波的方法分析了1954年—2013年天津地區(qū)主要?dú)庀笳军c(diǎn)的年降水量和年均氣溫.結(jié)果表明：天津站點(diǎn)的降水傾向率為-16.92～-14.11mm/（10a），年平均氣溫的線性傾向率為0.181～0.309℃/（10a）；天津地區(qū)氣溫在1954年—1982年波動(dòng)變化，1982年發(fā)生突變呈上升趨勢(shì)，且在1991年后上升顯著；天津地區(qū)降水量和氣溫的周期性具有相似性，降水量和氣溫在1960年—1970年呈負(fù)相關(guān)關(guān)系，其共振周期高能量區(qū)主要分布在1965年前后的2a；天津地區(qū)的降水量和氣溫在多時(shí)間尺度上主要呈非線性關(guān)系，降水量減少而氣溫升高，總體呈現(xiàn)干旱化趨勢(shì).

關(guān)鍵詞：年降水量；平均氣溫；突變檢驗(yàn)；交叉小波；天津地區(qū)

天津位于渤海西岸，地處海河流域的九河下梢，具有重要的戰(zhàn)略地位.然而，天津卻是一座資源型和水質(zhì)型缺水城市［1］，2013年天津水資源總量為14.64× 108m3，地表水供水總量為16.23×108m3，人均水資源量?jī)H101 m3［2］，遠(yuǎn)低于世界人均占有量1 000 m3的缺水警戒線，屬于重度缺水地區(qū).在氣候變化和人類活動(dòng)的大背景下，流域的降水量［3］、極端降水事件［4］和潛在蒸散發(fā)［5］等水循環(huán)要素的變化對(duì)水循環(huán)過(guò)程和水資源管理產(chǎn)生影響［6］.降水和氣溫是氣候的重要因素，研究降水和氣溫的變化規(guī)律是研究水循環(huán)過(guò)程變化的基礎(chǔ)，對(duì)水資源的合理規(guī)劃和利用具有重要意義［7-9］.

目前，已有學(xué)者運(yùn)用線性回歸、突變分析和小波分析的方法研究了天津所在區(qū)域降水量和氣溫的變化規(guī)律.袁再健等［7］研究表明海河流域年降水量呈下降趨勢(shì)，氣溫呈上升趨勢(shì)；降水年際變化較為明顯，而氣溫年際變化并不明顯.牛存穩(wěn)等［10］利用小波分析表明，降水量在20世紀(jì)50年代中期到20世紀(jì)70年代中期存在2～3 a的年際變化，15 a左右的年際變化發(fā)生在70年代中期和90年代中后期.王永財(cái)?shù)妊芯勘砻骱：恿饔虻哪昶骄鶜鉁卦?988年發(fā)生突變，降水量在1997年發(fā)生突變［11］.張健等［12］認(rèn)為京津冀降水量呈下降趨勢(shì)，氣候傾向率為15 mm/（10 a），夏季降水量在1996年初減少是一突變現(xiàn)象，其他季節(jié)降水量無(wú)突變現(xiàn)象.這些研究確定了降水和氣溫的突變趨勢(shì)，探討了降水的周期振蕩，但在時(shí)間域和頻率域中，有關(guān)降水和氣溫間多時(shí)間尺度的相關(guān)關(guān)系研究較少.

交叉小波變換可以發(fā)揮小波變換在時(shí)間域和頻率域中均可以表征氣候信號(hào)局部化特征的性能，對(duì)要素間相關(guān)關(guān)系的分析更具優(yōu)越性［13-14］，已應(yīng)用于水文要素和區(qū)域氣象要素多時(shí)間尺度特征的研究［15］.本研究以天津地區(qū)為研究對(duì)象，分析降水和氣溫長(zhǎng)時(shí)間序列的變化特征，運(yùn)用連續(xù)小波和交叉小波變換分析氣溫和降水間的關(guān)系，探討降水和氣溫的多時(shí)間尺度相關(guān)關(guān)系及其所包含的周期性特征，為深入研究天津地區(qū)的氣候變化和水文循環(huán)過(guò)程以及水資源管理和生態(tài)環(huán)境改善提供理論基礎(chǔ).

1　研究區(qū)概況

天津市東臨渤海，北依燕山，南北長(zhǎng)189 km，東西寬117 km，城市面積為11 917.3 km2.天津地區(qū)屬暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候，年平均氣溫為12℃左右，年降水量約為500～700 mm，年蒸發(fā)量為1 029 mm，年日照時(shí)數(shù)為2 610～3 090 h［16］.天津地貌總體為西北高、東南低，海拔由北向南逐漸下降；有山地、丘陵和平原3種地形，平原面積約占93%，如圖1所示.

圖1　天津地區(qū)氣象觀測(cè)點(diǎn)分布情況Fig.1　Distribution of meteorological observation sites in Tianjin area

天津地區(qū)土壤分布由西北向東南依次為潮土、鹽化潮土、沼澤土、鹽化濕潮土和濱海鹽土，海河的五大支流（北運(yùn)河、永定河、大清河、子牙河和南運(yùn)河）在此處匯合入海.天津境內(nèi)河網(wǎng)密布，水庫(kù)湖泊眾多，濕地類型豐富，有濱海濕地、河流濕地和城市濕地.近年來(lái)，受人類活動(dòng)和自然因素變化影響，天津市濕地面積和景觀格局發(fā)生巨大變化［16-17］.在自然因素中，年平均氣溫和降水量對(duì)濕地的影響顯著［16］.

2　研究方法

2.1資料來(lái)源

在中國(guó)地面氣象觀測(cè)站網(wǎng)中，有4個(gè)臺(tái)站位于天津地區(qū)，分別是天津、塘沽、渤海A平臺(tái)和寶坻，如圖1所示，其中天津和塘沽臺(tái)站的觀測(cè)時(shí)間序列最長(zhǎng)，而渤海A平臺(tái)和寶坻的觀測(cè)時(shí)間較短.各觀測(cè)點(diǎn)的相關(guān)數(shù)據(jù)如表1所示，其中年降水量和平均氣溫?cái)?shù)據(jù)來(lái)源于中國(guó)氣象數(shù)據(jù)網(wǎng)，數(shù)據(jù)經(jīng)過(guò)初步質(zhì)量檢測(cè).

表1　天津地區(qū)氣象觀測(cè)站點(diǎn)基本情況Tab.1　Basic case meteorological observation sites in Tianjin area

2.2趨勢(shì)分析和突變檢驗(yàn)

降水量和氣溫的趨勢(shì)分析采用線性傾向估計(jì)方法.用yi表示樣本量為n的某一變量，用ti表示yi所對(duì)應(yīng)的時(shí)間，建立一元線性回歸方程

yi=ati+b（i=1，2，…，n）（1）

式（1）中：回歸系數(shù)a表示變量y的趨勢(shì)傾向，a＞0說(shuō)明y隨時(shí)間t的增加呈上升趨勢(shì)，a＜0則表示y隨時(shí)間呈下降趨勢(shì)，通常稱a為傾向值；b為回歸常數(shù).

降水量和氣溫時(shí)間序列的突變點(diǎn)檢驗(yàn)采用Mann-Kendall（M-K）檢驗(yàn)法.Mann-Kendall檢驗(yàn)法是目前應(yīng)用較普遍的一種非參數(shù)檢驗(yàn)方法，對(duì)于時(shí)間序列x（n），先構(gòu)建秩序列Sk，Sk是i時(shí)刻數(shù)值大于j時(shí)刻數(shù)值個(gè)數(shù)的累計(jì)數(shù)，再計(jì)算序列的均值E（Sk）、方差Var（Sk）和統(tǒng)計(jì)量UFk，其中UF1＝0［18-19］.

按時(shí)間序列x逆序重復(fù)上述過(guò)程，使UBk=-UFk（k=n，n-1，…，1），UB1＝0.繪制UFk和UBk曲線圖，若UFk或UFk值大于0，表明序列呈上升趨勢(shì)，小于0則呈下降趨勢(shì).當(dāng)兩者曲線超過(guò)臨界直線時(shí)，表明上升或下降趨勢(shì)顯著；超過(guò)臨界線的范圍確定為出現(xiàn)突變的時(shí)間區(qū)域.如UFk和UBk2條曲線在臨界線之間出現(xiàn)交叉點(diǎn)，則交點(diǎn)對(duì)應(yīng)的時(shí)刻是突變開(kāi)始的時(shí)間［19］.

2.3小波變換

目前，小波分析可用于多時(shí)間尺度的變化特征研究中，連續(xù)小波分析（continuouswavelettransform，CWT）和交叉小波變換（cross wavelet transform，XWT）的分析方法和程序主要參考已有研究成果［14-15，20］.連續(xù)小波分析（CWT）采用復(fù)Morlet小波：

式（5）中：ω0為頻次；η為時(shí)間.當(dāng)ω0=6時(shí)，認(rèn)為小波的尺度參數(shù)幾乎等于傅里葉周期［14］.時(shí)間序列（xn，n= 1，…，N）的連續(xù)小波變換的卷積和小波標(biāo)準(zhǔn)化可定義為

小波影響錐（cone of influence，COI）表示小波譜區(qū)域以及相應(yīng)的邊緣效應(yīng)，小波譜值在COI的邊緣會(huì)下降e-2［14］.背景功率譜采用紅噪聲檢驗(yàn)，紅噪聲檢驗(yàn)過(guò)程采用一階自回歸方程.背景紅噪聲功率譜為式（7）中：α為紅噪聲功率譜中自回歸方程的相關(guān)系數(shù)；k為傅里葉頻率系數(shù).

交叉小波分析是結(jié)合交叉譜分析與小波變換2種方法，在時(shí)頻域中分析2個(gè)信號(hào)相關(guān)性的分析方法［21］.2個(gè)時(shí)間序列xn和yn的交叉小波譜（XWT）定義為WXY＝WXWY*，其中*為復(fù)共扼，對(duì)應(yīng)交叉小波譜為|WXY|.時(shí)間序列xn和yn的背景功率譜和定義為

式（8）中：Z（υp）為概率p的置信水平，顯著性檢驗(yàn)的標(biāo)準(zhǔn)譜選擇2個(gè)χ2分布積的平方根分布，復(fù)Morlet小波（自由度ν=2）的95%置信度下的置信水平Z（295%）= 3.999［14］.

小波相干（wavelet coherence，WTC）是反映2個(gè)小波變換在時(shí)頻域相干程度的量，定義為

式（9）中：S為平滑算子.小波凝聚譜反映2個(gè)小波變換在時(shí)頻域中的相干程度.交叉小波相位角反映兩序列在不同時(shí)域的滯后性特征，根據(jù)相位角正負(fù)向可以分析時(shí)頻域內(nèi)兩序列間的相關(guān)性［14-15，21］.

3　結(jié)果與分析

3.1降水的變化趨勢(shì)

基于天津地區(qū)氣象站點(diǎn)年降水量的觀測(cè)數(shù)據(jù)，得到降水量的變化趨勢(shì)，結(jié)果如圖2所示.

圖2　天津地區(qū)氣象站點(diǎn)年降水量線性趨勢(shì)Fig.2　Linear trend of annual precipitation at the meteorological observation sites in Tianjin area

由圖2（a）可知，天津站點(diǎn)近60年的平均年降水量為542.9 mm，標(biāo)準(zhǔn)差為148.3 mm，年降水量最多為975 mm（1977年），降水量最少為268.9 mm（1968年）.塘沽站點(diǎn)（圖2（b））的平均年降水量為577.9 mm，標(biāo)準(zhǔn)差為180.9mm，年降水量最多為1080.7mm（1964年），降水量最少為277.3 mm（1968年）.短時(shí)間序列的渤海A平臺(tái)和寶坻站點(diǎn)的平均年降水量分別為364.2 mm和591.2 mm.

通過(guò)線性傾向估計(jì)分析可知，天津地區(qū)年降水量呈減少趨勢(shì)，但均未達(dá)到顯著水平.天津站點(diǎn)的降水傾向率為-16.92 m/（10 a）（n=60，r=-0.20），塘沽為-14.11mm/（10a）（n=60，r=-0.14），即降水量每10年減少14.11～16.92 mm.這與塘沽站降水量減少率為每10年15.86 mm的已有研究成果［22］基本相同.

根據(jù)天津地區(qū)年降水量，運(yùn)用Mann-Kendall方法，繪制正序時(shí)間序列（UF）和逆序時(shí)間序列（UB）曲線，并進(jìn)行突變分析，結(jié)果如圖3所示.

圖3　天津地區(qū)年降水量Mann-Kendall檢驗(yàn)Fig.3　Mann-KendalltestofannualprecipitationinTianjinarea

由圖3可知，天津地區(qū)降水量波動(dòng)變化，有逐年減少的趨勢(shì).在天津站點(diǎn)年降水量的統(tǒng)計(jì)曲線中UF＜0，表明降水量呈減少趨勢(shì).在1957年—1960年和1963年，天津站點(diǎn)降水減少的趨勢(shì)達(dá)到0.05的顯著水平，而其他年份的減少趨勢(shì)均不顯著.20世紀(jì)60年代至80年代，塘沽站點(diǎn)的降水量波動(dòng)變化，多呈現(xiàn)增加趨勢(shì)，如1954年—1955年和1964年—1980年，降水量的UF＞0，表明降水有增加趨勢(shì).然而，在1990年后，降水量的UF＜0，表明降水量一直呈減少趨勢(shì)，但趨勢(shì)性均沒(méi)有達(dá)到0.05的顯著水平.

分析UF和UB統(tǒng)計(jì)曲線在臨界線間的交點(diǎn)可知，天津地區(qū)降水量波動(dòng)明顯，出現(xiàn)多次突變.天津站點(diǎn)降水量在1965年—1970年的波動(dòng)較大，在1980年和1990年左右也出現(xiàn)多次變化.1990年后，天津站點(diǎn)降水量呈減少趨勢(shì)，沒(méi)有發(fā)生突變現(xiàn)象.塘沽站點(diǎn)降水量在1960年發(fā)生突變，隨后呈現(xiàn)先減少后增加的變化趨勢(shì)，并在1978年、1981年和2010年分別出現(xiàn)突變現(xiàn)象，降水量波動(dòng)明顯.

3.2氣溫的變化趨勢(shì)

與降水量減少的趨勢(shì)相反，天津地區(qū)年平均氣溫呈現(xiàn)明顯增加的趨勢(shì)，各站點(diǎn)的統(tǒng)計(jì)結(jié)果如圖4所示.天津站點(diǎn)的平均氣溫為12.6℃，標(biāo)準(zhǔn)差為0.6℃，年平均氣溫最高為13.6℃（1989年和2007年），最低為11℃（1956年和1969年）.塘沽站點(diǎn)的平均氣溫為12.6℃，標(biāo)準(zhǔn)差為0.7℃，年平均氣溫最高為13.9℃（2007年），最低為10.6℃（1969年）.渤海A平臺(tái)和寶坻站點(diǎn)的平均年氣溫分別為12.9℃和11.6℃.

基于長(zhǎng)時(shí)間序列的天津和塘沽站點(diǎn)的年平均氣溫?cái)?shù)據(jù)，通過(guò)線性傾向估計(jì)可知，天津地區(qū)年平均氣溫呈升高趨勢(shì)，趨勢(shì)性達(dá)到極顯著水平.天津站點(diǎn)年平均氣溫的線性傾向率為0.181℃/（10 a）（n= 60，r=0.53）；塘沽站點(diǎn)的線性傾向率為0.309℃/（10 a）（n=60，r=0.72），均達(dá)到極顯著水平.這與天津?yàn)I海新區(qū)每10 a平均氣溫上升0.29℃的已有研究成果［22］基本一致.

天津地區(qū)年平均氣溫M-K檢驗(yàn)的統(tǒng)計(jì)曲線如圖5所示.天津站點(diǎn)1954年—1958年的UF＜0，表明氣溫下降，氣候變冷；1958年—1968年的UF＞0，表明氣溫升高，氣候增暖；1968年—1982年，氣溫下降；1982年后，氣溫升高；1991年后氣溫升高的趨勢(shì)達(dá)到0.05的顯著水平，表明氣候顯著增暖.在塘沽站點(diǎn)，氣溫的變化趨勢(shì)與天津站點(diǎn)基本相同，年平均氣溫也在1991年后顯著增加.

通過(guò)突變分析可知，天津地區(qū)年平均氣溫的變化趨勢(shì)與降水量變化趨勢(shì)不同，年平均氣溫的UF和UB統(tǒng)計(jì)曲線交點(diǎn)唯一，表明僅有1個(gè)突變時(shí)間.天津站點(diǎn)和塘沽站點(diǎn)年平均氣溫在1954年—1982年波動(dòng)變化，而在1982年后，呈現(xiàn)氣溫升高趨勢(shì)，并在1991年后顯著增加.天津站點(diǎn)氣溫突變的發(fā)生時(shí)間早于塘沽站點(diǎn)，年平均氣溫突變時(shí)間為1986年，而塘沽站點(diǎn)的突變時(shí)間是1988年，表明天津地區(qū)的年平均氣溫在20世紀(jì)80年代末發(fā)生突變，氣候變?yōu)樵雠厔?shì).

圖4　天津地區(qū)氣象站點(diǎn)年平均氣溫線性趨勢(shì)Fig.4　Linear trend of annual temperature at the meteorological observation sites in Tianjin area

圖5　天津地區(qū)年平均氣溫Mann-Kendall檢驗(yàn)Fig.5　Mann-Kendall test of annual temperature in Tianjin area

3.3降水和氣溫的小波分析

采用Morlet小波函數(shù)分別對(duì)天津地區(qū)的年降水量和平均氣溫進(jìn)行小波變換，得到1954年—2013年天津站點(diǎn)和塘沽站點(diǎn)的年降水量和年均氣溫的連續(xù)小波譜，結(jié)果如圖6所示.圖6中紅色和藍(lán)色分別表示能量密度的峰值和谷值，反映出主導(dǎo)波動(dòng)組分時(shí)頻變換的局部性和動(dòng)態(tài)性特征，顏色深淺表示能量密度的相對(duì)變化.黑色粗實(shí)線為95%置信區(qū)間邊界，通過(guò)了紅噪聲檢驗(yàn)；黑色細(xì)實(shí)線為小波影響錐邊界，為受連續(xù)小波變換數(shù)據(jù)邊緣效應(yīng)影響較大的區(qū)域［14，20］.

運(yùn)用連續(xù)小波分析天津地區(qū)年降水量和年均氣溫的周期性.由圖6（a）可知，天津站點(diǎn)的年降水量存在3個(gè)周期，1965年—1970年存在1～2 a的周期，1970年—1978年的周期為4 a左右，1980年—1990年則存在10 a左右的周期.由圖6（b）可知，塘沽站點(diǎn)的年降水量有2個(gè)周期，1960年—1970年存在1～2 a的周期，1980年代末有3 a左右的周期.天津站點(diǎn)的年均氣溫也有3個(gè)周期，1975年—1978年的周期為1～2 a左右，1980年—1990年存在6～8 a左右的周期，1995年—2000年的周期為4 a左右（圖6（c））.塘沽站點(diǎn)的年均氣溫存在2個(gè)周期，其中1975年—1980年存在1～2 a的周期，而1985年左右的周期特征不明顯（圖6（d））.在同一站點(diǎn)，降水量與氣溫的周期性具有相似性，年降水量的周期比年均氣溫的周期在時(shí)間上提前約10 a左右.

圖6　天津地區(qū)降水和氣溫的連續(xù)小波分析Fig.6　Continuous wavelet transform of annual precipitation and annual temperature

通過(guò)對(duì)連續(xù)小波變換后的系數(shù)進(jìn)行交叉小波變換和小波相關(guān)變換，對(duì)天津站點(diǎn)和塘沽站點(diǎn)降水量和氣溫相互間的小波凝聚譜進(jìn)行分析，并運(yùn)用標(biāo)準(zhǔn)譜進(jìn)行顯著性檢驗(yàn)，從多時(shí)間尺度的角度探討兩者在時(shí)頻域中的相關(guān)性［15］，結(jié)果如圖7所示.圖7中箭頭方向反映降水量和氣溫的相位關(guān)系，其中由左向右的箭頭（→）表示兩者變化同相位，呈正相關(guān)關(guān)系；由右指向左的箭頭（←）表示反相位，呈負(fù)相關(guān)關(guān)系；而垂直向下（↓）和垂直向上（↑）分別表示降水量的小波變換提前和落后氣溫1/4個(gè)周期［20］，呈非線性相關(guān)［23］.

圖7　天津地區(qū)降水和氣溫的交叉小波分析Fig.7　Cross wavelet transform of annual precipitation and annual temperature

由圖7可知，天津地區(qū)降水量和氣溫的相關(guān)關(guān)系存在周期性波動(dòng).1960年—1970年，天津站點(diǎn)和塘沽站點(diǎn)的降水量和氣溫均存在負(fù)相關(guān)關(guān)系，降水量和氣溫的共振周期高能量區(qū)主要分布在1965年前后2 a. 1990年—2000年，天津站點(diǎn)降水量和氣溫存在非線性關(guān)系，降水量和氣溫的共振周期高能區(qū)主要分布在1995年前后的5 a和9 a.1980年—1990年，塘沽站點(diǎn)的降水量和氣溫存在非線性關(guān)系，其共振周期高能量區(qū)主要在1987年前后2 a.綜上所述，天津地區(qū)的降水量和氣溫在多時(shí)間尺度上主要呈非線性關(guān)系.

4　結(jié)論

（1）以天津地區(qū)1954年—2013年的降水量和氣溫為研究對(duì)象，線性傾向估計(jì)表明天津地區(qū)降水量呈減少趨勢(shì)，氣溫呈升高趨勢(shì).天津站點(diǎn)的降水傾向率為-16.92 m/（10a），年平均氣溫的線性傾向率為0.181℃/ （10 a）；塘沽站點(diǎn)的降水傾向率為-14.11 mm/（10 a），年平均氣溫為0.309℃/（10 a）.

（2）Mann-Kendall檢驗(yàn)表明天津地區(qū)降水量波動(dòng)變化，氣溫在20世紀(jì)80年代末發(fā)生突變，氣候呈變暖趨勢(shì).天津站點(diǎn)降水量在1990年后呈減少趨勢(shì)；塘沽站點(diǎn)在1960年發(fā)生突變，降水量出現(xiàn)先減少后增加的變化趨勢(shì).天津地區(qū)氣溫在1954年—1982年波動(dòng)變化，在1982年后呈增加趨勢(shì)，且在1991年后顯著增加.

（3）根據(jù)連續(xù)小波和交叉小波分析，天津地區(qū)降水量和氣溫的周期性具有相似性，年降水量的周期比年均氣溫的周期在時(shí)間上提前約10 a.降水量和氣溫在1960年—1970年呈負(fù)相關(guān)關(guān)系，其共振周期高能量區(qū)主要分布在1965年前后的2 a，但天津地區(qū)的降水量和氣溫在多時(shí)間尺度上主要為非線性關(guān)系.

（4）天津地區(qū)降水量減少而氣溫升高，總體呈現(xiàn)干旱化趨勢(shì).降水量和氣溫是主要的氣象因素，但也受到其他氣象因素和環(huán)境因素的影響，降水量和氣溫變化與各種氣象因子間的綜合關(guān)系尚待進(jìn)一步研究.

參考文獻(xiàn)：

［1］田萌，李萬(wàn)慶.天津地區(qū)水資源問(wèn)題戰(zhàn)略研究［J］.環(huán)境保護(hù)，2008 （6）：67-69. TIAN M，LI W Q.The research on water resources strategy in Tianjin［J］. Environmental Protection，2008（6）：67-69（in Chinese）.

［2］中華人民共和國(guó)國(guó)家統(tǒng)計(jì)局.中國(guó)統(tǒng)計(jì)年鑒—2014［M］：中國(guó)統(tǒng)計(jì)出版社，2014. National Bureau of Statistics of the People′s Republic of China.China Statistical Yearbook-2014［M］.Beijing：China Statistics Press，2014（in Chinese）.

［3］于贏東，楊志勇，劉永攀，等.變化環(huán)境下海河流域降水演變研究綜述［J］.水文，2010，30（4）：32-35. YU Y D，YANG Z Y，LIU Y P，et al.Review of study on precipitation of Haihe river basin under changing environment［J］.Journal of China Hydrology，2010，30（4）：32-35（in Chinese）.

［4］張兵，韓靜艷，王中良，等.海河流域極端降水事件時(shí)空變化特征分析［J］.水電能源科學(xué)，2014，32（2）：15-18，34. ZHANG B，HAN J Y，WANG Z L，et al.Temporal and spatial variation characteristics of extreme precipitation events in Haihe basin［J］. Water Resources and Power，2014，32（2）：15-18，34（in Chinese）.

［5］劉小莽，鄭紅星，劉昌明，等.海河流域潛在蒸散發(fā)的氣候敏感性分析［J］.資源科學(xué)，2009，31（9）：1470-1476. LIU X M，ZHENG H X，LIU C M，et al.Sensitivity of the potential evapotranspiration to key climatic variables in the Haihe River basin［J］. Resources Science，2009，31（9）：1470-1476（in Chinese）.

［6］夏軍.華北地區(qū)水循環(huán)與水資源安全：?jiǎn)栴}與挑戰(zhàn)［J］.地理科學(xué)進(jìn)展，2002，21（6）：517-526. XIA J.A perspective on hydrological base of water security problem and its application study in north China［J］.Progress in Geography，2002，21（6）：517-526（in Chinese）.

［7］袁再健，沈彥俊，褚英敏，等.海河流域近40年來(lái)降水和氣溫變化趨勢(shì)及其空間分布特征［J］.水土保持研究，2009，16（3）：24-26. YUAN Z J，SHEN Y J，CHU Y M，et al.Variations and distribution of temperature and precipitation of Haihe river basin in recent 40 years［J］. Research of Soil and Water Conservation，2009，16（3）：24-26（in Chinese）.

［8］夏范燕，吳鞏勝，李麗，等.近50年內(nèi)滇西北極端氣候變化［J］.云南師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2014，34（3）：68-73. XIA F Y，WU G S，LI L，et al.The extreme climate change of northwest Yunnan in recent 50 year［J］.Journal of Yunnan Normal University：Natural Sciences Edition，2014，34（3）：68-73.

［9］覃順萍，吳鞏勝，李麗，等.1961—2010年云南省極端降水時(shí)空變化特征［J］.云南師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2015，35（1）：71-76. XIA F Y，WU G S，LI L，et al.Temporal-spatial variations of extreme precipitation events from 1961 to 2010 in Ynnan province［J］.Journal of Yunnan Normal University：Natural Sciences Edition，2015，35（1）：71-76.

［10］牛存穩(wěn)，張利平，夏軍.華北地區(qū)降水量的小波分析［J］.干旱區(qū)地理，2004，27（1）：66-70. NIU C W，ZHANG L P，XIA J.Wavelet analysis on the precipitation in northChina［J］.AridLandGeography，2004，27（1）：66-70（inChinese）.

［11］王永財(cái)，孫艷玲，張靜，等.近5 1年海河流域氣候變化特征分析［J］.天津師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2014，34（4）：58-63. WANG Y C，SUN Y L，ZHANG J，et al.Climate change characteristics of Haihe River Basin in recent 51 years［J］.Journal of Tianjin Normal University：Natural Science Edition，2014，34（4）：58-63（in Chinese）.

［12］張健，章新平，王曉云，等.近47年來(lái)京津冀地區(qū)降水的變化［J］.干旱區(qū)資源與環(huán)境，2010，24（2）：74-80. ZHANG J，ZHANG X P，WANG X Y，et al.Precipitation change in Beijing-Tianjin-Hebei province［J］.Journal of Arid Land Resources and Environment，2010，24（2）：74-80（in Chinese）.

［13］孫衛(wèi)國(guó)，程炳巖.交叉小波變換在區(qū)域氣候分析中的應(yīng)用［J］.應(yīng)用氣象學(xué)報(bào)，2008，19（4）：479-487. SUN W G，CHENG B Y.Application of cross wavelet transformation to analysis on regional climate variations［J］.Journal of Applied Meteorological Science，2008，19（4）：479-487（in Chinese）.

［14］GRINSTED A，MOORE J C，JEVREJEVA S.Application of the cross wavelet transform and wavelet coherence to geophysical time series［J］. Nonlinear Processes in Geophysics，2004，11：561-566.

［15］劉志方，劉友存，郝永紅，等.黑河出山徑流過(guò)程與氣象要素多尺度交叉小波分析［J］.干旱區(qū)地理，2014，37（6）：1137-1146. LIU Z F，LIU Y C，HAO Y H，et al.Multi-time scale cross-wavelet transformation between runoff and climate factors in the upstream of Heihe River［J］.Arid Land Geography，2014，37（6）：1137-1146（in Chinse）.

［16］劉東云，王琦，杜林芳，等.1999—2007年天津濕地景觀格局及水鳥(niǎo)生境變化［J］.濕地科學(xué)，2012，10（3）：350-358. LIU D Y，WANG Q，DU L F，et al.Changes of wetland landscape pattern and water birds’habitate in Tianjin during 1999-2007［J］.Wetland Science，2012，10（3）：350-358（in Chinese）.

［17］高軍，王中良.1976—2009年間天津濕地面積變化特征及驅(qū)動(dòng)因素分析［J］.天津師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，33（4）：32-38. GAO J，WANG Z L.Analysis on variation characteristics of Tianjin wetland area and its main controlling facts during 1976-2009［J］.Journal of Tianjin Normal University：Natural Science Edition，2013，33（4）：32-38（in Chinese）.

［18］周園園，師長(zhǎng)興，范小黎，等.國(guó)內(nèi)水文序列變異點(diǎn)分析方法及在各流域應(yīng)用研究進(jìn)展［J］.地理科學(xué)進(jìn)展，2011，30（11）：1361-1369. ZHOU Y Y，SHI C X，F(xiàn)AN X L，et al.Advances in the research methods of abrupt changes of hydrologic sequences and their applications in drainage basins in China［J］.Progress in Geography，2011，30（11）：1361-1369（in Chinese）.

［19］魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)［M］.2版.北京：氣象出版社，2007. WEI F Y.Climatic Statistical Diagnosis and Forecasting Methods［M］. SecondEdition.Beijing：ChinaMeteorologicalPress，1999（inChinse）.

［20］祁曉凡，楊麗芝，韓曄，等.濟(jì)南泉域地下水位動(dòng)態(tài)及其對(duì)降水響應(yīng)的交叉小波分析［J］.地球科學(xué)進(jìn)展，2012，27（9）：969-978. QI X F，YANG L Z，HAN Y，et al.Cross wavelet analysis of groundwater level regimes and precipitation groundwater level regime in Ji’nan spring region［J］.Advances in Earth Science，2012，27（9）：969-978 （in Chinese）.

［21］LABAT D.Cross wavelet analyses of annual continental freshwater discharge and selected climate indices［J］.Journal of Hydrology，2010，385（1/2/3/4）：269-278.

［22］趙玉潔，宋國(guó)輝，徐明娥，等.天津?yàn)I海區(qū)50年局地氣候變化特征［J］.氣象科技，2004，32（2）：86-89，96. ZHAO Y J，SONG G H，XU M E，et al.Analysis of local climate changes in Tianjin coastal area in last 50 years［J］.Meteorological Science and Technology，2004，32（2）：86-89，96（in Chinese）.

［23］LIU J，F(xiàn)U G，SONG X，et al.Stable isotopic compositions in Australian precipitation［J］.Journal of Geophysical Research，2010，115：23.

（責(zé)任編?？涸颍?/p>

第一作者：張兵（1983—），男，助理研究員，主要從事水循環(huán)與水環(huán)境方面的研究.

文章編號(hào)：1671-1114（2016）01-0032-08

中圖分類號(hào)：P458

文獻(xiàn)標(biāo)志碼：A

收稿日期：2015-10-09

基金項(xiàng)目：天津市科學(xué)基金聯(lián)合資助項(xiàng)目（15JCQNJC44200）；天津市高等學(xué)?！皠?chuàng)新團(tuán)隊(duì)培養(yǎng)計(jì)劃”資助項(xiàng)目（TD12-5037）；天津師范大學(xué)博士基金資助項(xiàng)目（52XB1401）.

通信作者：王中良（1970—），男，研究員，主要從事環(huán)境地理與環(huán)境地球化學(xué)方面的研究.

Trend and multi-time scale cross-wavelet analysis on precipitation and temperature in Tianjin area

ZHANG Bing，WANG Zhongliang
（Tianjin Key Laboratory of Water Resources and Environment，Tianjin Normal University，Tianjin 300387，China）

Abstract：Because of the analysis on the change of precipitation and temperature is important to the water cycle research and water resources planning and utilization，the annual precipitation and annual average temperature at the main weather stations in Tianjin area from 1954 to 2013 were analyzed using linear tendency estimation，Mann-Kendall test and cross wavelet transform methods.The results show that the tendency rate of precipitation is-16.92 mm/（10 a）to-14.11mm/（10 a），and the linear trend rate of annual average temperature is 0.181℃/（10 a）to 0.309℃/（10 a）.The temperature of Tianjin area changes from 1954 to 1982，the mutation occurs in 1982 showing a trend of increase，and there is a significantly rise after 1991.The periodic variations of precipitation and air temperature are similar，the precipitation and temperature show a negative correlation relationship from 1960 to 1970，and the high-energy area of resonant period is mainly distributed around 1965 as 2 a period.The relationship between precipitation and temperature is mainly nonlinear relationship in Tianjin area.Tianjin area is facing a tendency of drying because of the reducing precipitation and increasing temperature.

Keywords：annual precipitation；annual average temperature；mutation test；cross wavelet transform；Tianjin area