陳祥義，肖文發(fā)，黃志霖，曾立雄

中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所//國家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室，北京 100091

1951─2012年三峽庫區(qū)降水時(shí)空變化研究

陳祥義，肖文發(fā)，黃志霖*，曾立雄

中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護(hù)研究所//國家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室，北京 100091

對(duì)區(qū)域降水時(shí)空分布變化趨勢(shì)的分析可以為區(qū)域水資源的合理管理與利用提供參考，利用三峽庫區(qū)及周邊共27個(gè)氣象站點(diǎn)的逐月降水?dāng)?shù)據(jù)，通過線性傾向統(tǒng)計(jì)、滑動(dòng)平均兩種方法對(duì)三峽庫區(qū)1951─2012年年均、雨季、旱季降水與降水日數(shù)的變化趨勢(shì)進(jìn)行了研究，并利用多元回歸分析對(duì)三峽庫區(qū)內(nèi)年均降水量和降水日數(shù)的空間分布情況進(jìn)行了分析。結(jié)果表明，三峽庫區(qū)近62 a來的年均、年均雨季、年均旱季降水量和年均降水日數(shù)都有減少趨勢(shì)，年均降水日數(shù)的減少趨勢(shì)比降水量的減少趨勢(shì)更加顯著，并且降水日數(shù)明顯減少的時(shí)間要早于降水量，這表明平均每個(gè)降水日內(nèi)的降水量呈增加的趨勢(shì)，也即短歷時(shí)強(qiáng)降雨事件呈增加趨勢(shì)。多年平均降水量與經(jīng)度（P=0.081）、緯度（P=0.367）的相關(guān)性均不顯著，而與海拔高度（P=1.90E-4）達(dá)到極顯著相關(guān)水平。多年平均降水日數(shù)與經(jīng)度（P=0.539）相關(guān)系不顯著，而與緯度（P=8.77E-4）和海拔高度（P=1.82E-12）均達(dá)到極顯著相關(guān)。降水量和降水日數(shù)均與海拔高度達(dá)到極顯著正相關(guān)，可以利用海拔高度相對(duì)準(zhǔn)確預(yù)測(cè)研究區(qū)內(nèi)年降水和降水日數(shù)的空間分布特征，海拔每升高100 m，對(duì)應(yīng)降水量與降水日數(shù)分別約增加30 mm、4.5 d。

三峽庫區(qū)；降水量； 降水日數(shù)；海拔；線性趨勢(shì)；多元回歸分析

全球氣候變化受到越來越多的關(guān)注，降水變化和氣溫變化是氣候變化的兩種主要表現(xiàn)形式。水資源作為人類社會(huì)生存發(fā)展的主要限制，其管理與利用成為人類可持續(xù)發(fā)展成敗的關(guān)鍵。水資源的合理利用與配置必須建立在對(duì)區(qū)域水資源的時(shí)空分布和變化趨勢(shì)的準(zhǔn)確把握基礎(chǔ)上。降水作為水資源的直接補(bǔ)給形式以及整個(gè)流域水循環(huán)的輸入來源，準(zhǔn)確預(yù)測(cè)其時(shí)空分布和變化趨勢(shì)可為水資源的合理利用配置提供參考。

三峽庫區(qū)作為三峽水庫的淹沒區(qū)和主要集水區(qū)，其區(qū)域內(nèi)的降水變化將直接影響到三峽水庫及下游的防汛抗旱調(diào)度。針對(duì)三峽庫區(qū)不同降水特征前人已進(jìn)行了一些研究，發(fā)現(xiàn)三峽庫區(qū)最大連續(xù)5 d降水量線性變化趨勢(shì)不顯著，但是在未來氣候變化背景下呈增加趨勢(shì)（張?zhí)煊畹龋?010），夏季（6─8月）降水也呈增加趨勢(shì)（廖要明等，2008；鄭剛等，2009），但三峽大壩建成前后的年降水和四季降水沒有明顯變化（張樹奎等，2013）。已有研究主要集中在降水量的變化，而對(duì)表征降水年內(nèi)分布均勻程度的降水日數(shù)變化還研究較少。對(duì)年內(nèi)降水變化的分析主要以春、夏、秋、冬4季為時(shí)間尺度，而對(duì)作為流域水文調(diào)控主要依據(jù)的雨季、旱季的降水變化研究較少。

我國氣象監(jiān)測(cè)規(guī)定日降水量（或夾有霧、露、霜量）大于或等于0.1 mm為一個(gè)降水日（顧鈞禧等，1994）。降水日數(shù)可以在一定程度上表征降水在年內(nèi)分布均勻情況，如年降水量不變的情況下降水日數(shù)越多說明降水在年內(nèi)分布越均勻。本文在對(duì)三峽庫區(qū)多年平均降水、降水日數(shù)分析的基礎(chǔ)上，分析雨季（4─9月）、旱季（10月─次年3月）的降水以及降水日數(shù)的變化，同時(shí)分析了庫區(qū)多年平均降水量、降水日數(shù)與站點(diǎn)經(jīng)度、緯度和海拔高度之間的相關(guān)性，探究庫區(qū)降水時(shí)空分布規(guī)律，可為三峽水庫年際或中長期水庫庫容調(diào)度和水資源的合理利用提供參考，保證長江中下游居民的生活生產(chǎn)用水安全。

1 研究區(qū)概況

三峽庫區(qū)是指三峽大壩建成后蓄水淹沒所涉及的地區(qū)，西起重慶江津、東至湖北宜昌市，包括重慶市所轄22個(gè)縣（區(qū)）和湖北省所轄4個(gè)縣（圖1）。研究區(qū)地處亞熱帶的北緣，屬亞熱帶大陸性季風(fēng)氣候，溫暖濕潤，四季分明，光照水熱充足。年均溫18.0 ℃，年均空氣相對(duì)濕度77%，無霜期307 d，全年平均日照時(shí)數(shù)1200～1600 h，多年平均日照率25%，年降水量950～1590 mm，且75%左右集中在4─9月，5─9月多有暴雨出現(xiàn)。受地形的影響，局地氣候明顯，年均溫和年降雨量高于同緯度的中國東部地區(qū)（潘磊等，2012）。

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

2.1 數(shù)據(jù)來源

本文所用數(shù)據(jù)源于國家科技基礎(chǔ)條件平臺(tái)中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)（http://cdc.cma.gov.cn/）中三峽庫區(qū)內(nèi)及周邊共27個(gè)氣象站點(diǎn)（圖1）1951─2012年的月降水?dāng)?shù)據(jù)，氣象站點(diǎn)基本信息見表1。在數(shù)據(jù)統(tǒng)計(jì)過程中，對(duì)個(gè)別站點(diǎn)某年之中月降水?dāng)?shù)據(jù)統(tǒng)計(jì)不全的年份進(jìn)行了剔除。

圖1 研究區(qū)氣象站點(diǎn)分布Fig. 1 The distribution of weather stations in the study area

表1 站點(diǎn)基本信息統(tǒng)計(jì)表Table 1 Basic information and statistic characteristics of the weather stations

2.2 研究方法

常用的氣象趨勢(shì)統(tǒng)計(jì)方法主要有線性傾向估計(jì)、滑動(dòng)平均、累計(jì)距平、3次樣條函數(shù)等方法（魏鳳英，2007）。線性傾向估計(jì)、滑動(dòng)平均方法因簡單易懂、直觀明了，在當(dāng)今氣象趨勢(shì)統(tǒng)計(jì)中應(yīng)用較多（馮亞文等，2013；李珊珊等，2012）。所以本文主要應(yīng)用這兩種方法對(duì)長江流域三峽庫區(qū)的多年平均、雨季以及旱季的平均降水量和降水日數(shù)的變化趨勢(shì)進(jìn)行分析，同時(shí)利用線性回歸方法對(duì)不同站點(diǎn)的降水量、降水日數(shù)與經(jīng)度、緯度和海拔之間的關(guān)系進(jìn)行了分析。

3 結(jié)果分析

3.1 降水、降水日數(shù)時(shí)間變化趨勢(shì)分析

3.1.1 降水、降水日數(shù)年際變化趨勢(shì)分析

三峽庫區(qū)近 62年來的平均年降水量為 1167 mm，最高降水量為1954年的1538 mm，最低降水量為2001年的853 mm，最高降水量約為最低降水量的1.8倍，這說明三峽庫區(qū)年際降水差異較大，存在明顯的豐水年、枯水年。從1950年到2000年平均年降水量在均線附近波動(dòng)，2000年以后降水量呈減少趨勢(shì)，總的線性傾向率為-11.8 mm/10a（圖2a）。自由度為60所對(duì)應(yīng)的相關(guān)系數(shù)顯著性水平分別為r0.05=0.2500，r0.01=0.3248，而r=0.1382

三峽庫區(qū)近 62年來的平均年降水日數(shù)為 140 d，最高降水日數(shù)為1954年的175 d，最低降水日數(shù)為2011年的117 d，最高降水日數(shù)為最低降水日數(shù)的1.50倍。從20世紀(jì)50年代到90年代初期年降水日數(shù)在均線附近波動(dòng)，且大部分在均線之上，而從20世紀(jì)90年代初期至今年降水日數(shù)出現(xiàn)了明顯的下降趨勢(shì)，總的線性傾向率為-3.04d/10a（圖2b）。r=0.4715>r0.01=0.3248，說明在0.01顯著性水平下降水日數(shù)的減少趨勢(shì)達(dá)到極顯著水平。

圖2 三峽庫區(qū)年均降水量Fig. 2 Change trend of average annual precipitation and average annual rainy days in TGRA

綜上分析得知，降水量減少趨勢(shì)弱于降水日數(shù)減少趨勢(shì)，這表明每個(gè)降水日內(nèi)的降水量有增加趨勢(shì)，這會(huì)導(dǎo)致短歷時(shí)強(qiáng)降水的出現(xiàn)，容易引發(fā)山洪、泥石流等地質(zhì)災(zāi)害。

3.1.2 雨季、旱季降水、降水日數(shù)變化趨勢(shì)分析

三峽庫區(qū)近 60年平均雨季降水量約為 877 mm，平均旱季降水量為292 mm，雨季降水量約為旱季降水量的3倍，雨季降水量約占到全年平均降水量的75%。從圖3a和4a可以看出，雨季、旱季降水量的線性傾向率分別為-8.92、-2.0 mm/10a，線性相關(guān)系數(shù)r分別為0.1058、0.0632，均小于r0.05的值，所以雨季、旱季降水量的減少趨勢(shì)與年均降水量的減少趨勢(shì)一樣也是不顯著的。

圖3 三峽庫區(qū)雨季降水量Fig. 3 Change trend of precipitation and rainy days of rainy season in TGRA

三峽庫區(qū)多年平均雨季降水日數(shù)約為78 d，平均旱季降水日數(shù)為62 d，平均雨季降水日數(shù)約為平均旱季降水日數(shù)的1.26倍，明顯小于雨季、旱季降水量之比的3倍，這表明雨季、旱季降水日數(shù)相差不大，而雨季每個(gè)降水日內(nèi)的降水量顯著高于旱季每個(gè)降水日的降水量，從而導(dǎo)致了降水量在年內(nèi)的分布不均。雨季、旱季降水日數(shù)的線性傾向率分別為-1.648、-1.329 d/10a（圖3b、圖4b），線性相關(guān)系數(shù)r分別為0.3521、0.3286，均大于r0.05，小于r0.01，所以雨季、旱季降水日數(shù)的減少趨勢(shì)在 0.05顯著性水平上是顯著的。

3.2 降水、降水日數(shù)空間變化分析

3.2.1 降水空間變化分析

對(duì)三峽庫區(qū) 27個(gè)氣象站點(diǎn)的多年平均降水量與經(jīng)度、緯度和海拔高度進(jìn)行回歸分析，得到回歸方程（1）：

式中：Q——多年平均降水量（mm）；x——經(jīng)度（°）；y——緯度（°）；z——站點(diǎn)海拔高度（m）。

圖4 三峽庫區(qū)旱季降水量Fig. 4 Change trend of precipitation and rainy days of dry season in TGRA

表2 多年平均降水量與經(jīng)度、緯度、海拔回歸分析統(tǒng)計(jì)表Table 2 Statistical table of regression analysis between average annual rainfall and longitude, latitude, altitude

分析得知，復(fù)相關(guān)系數(shù)r=0.721，P=6.338E-4< 0.001，這表明回歸方程的線性關(guān)系是極顯著的。多年平均降水量與經(jīng)度（P=0.081）、緯度（P=0.367）之間的相關(guān)性是不顯著的，多年平均降水量與海拔高度（P=1.899E-4）呈極顯著正相關(guān)（表2），這表明多年平均降水的空間差異主要受海拔高度的影響。對(duì)多年平均降水量與海拔進(jìn)行回歸分析，得到回歸方程（2），結(jié)果表明海拔每升高100 m，對(duì)應(yīng)降水量約增加30 mm。利用Arcgis 9.2進(jìn)行柵格計(jì)算得到三峽庫區(qū)多年平均降水量空間分布如圖 5所示。

圖5 三峽庫區(qū)多年平均降水量空間分布Fig. 5 Spatial distribution of average annual precipitation in TGRA

3.2.2 降水日數(shù)空間變化分析

對(duì)三峽庫區(qū) 27個(gè)氣象站點(diǎn)的多年平均降水日數(shù)與經(jīng)度、緯度和海拔高度進(jìn)行回歸分析，得到回歸方程（3）：

式中：D——多年平均降水日數(shù)（d）；x——經(jīng)度（°）；y——緯度（°）；z——站點(diǎn)海拔高度（m）。

分析得知，復(fù)相關(guān)系數(shù)r=0.958，P=1.386E-12，這表明回歸方程的線性關(guān)系是極顯著的。多年平均降水日數(shù)與經(jīng)度之間的相關(guān)性是不顯著的（P=0.539），多年平均降水日數(shù)與緯度（P=8.77E-4）和海拔高度（P=1.82E-12）呈極顯著正相關(guān)（表3），多年平均降水日數(shù)的空間差異主要受海拔高度的影響。僅對(duì)多年平均降水日數(shù)與海拔進(jìn)行回歸分析，得到回歸方程（4），結(jié)果表明僅用海拔對(duì)降水日數(shù)的回歸模擬可以達(dá)到精度要求，海拔每升高100 m，對(duì)應(yīng)降水日數(shù)約增加4.5 d。利用Arcgis 9.2進(jìn)行柵格計(jì)算得到三峽庫區(qū)多年平均降水日數(shù)空間分布如圖6所示。

表3 多年平均降水日數(shù)與經(jīng)度、緯度、海拔回歸分析統(tǒng)計(jì)表Table 3 Statistical table of regression analysis between average annual rainy days and longitude, latitude, altitude

圖6 三峽庫區(qū)多年平均降水日數(shù)空間分布Fig. 6 Spatial distribution of average annual rainy days in TGRA

4 結(jié)論與討論

4.1 結(jié)論

（1）三峽庫區(qū)年降水量在 2000年之前在均線附近波動(dòng)，之后呈減少趨勢(shì)，但減少趨勢(shì)不是很明顯。年降水日數(shù)在 1995年之前在均線附近波動(dòng)，之后呈現(xiàn)顯著的減少趨勢(shì)。降水日數(shù)的減少趨勢(shì)大于降水量的減少趨勢(shì)，這表明降水日內(nèi)平均降水量有增加的趨勢(shì)，降水日內(nèi)容易出現(xiàn)短時(shí)強(qiáng)降雨。

（2）雨季、旱季降水量在多年均線附近波動(dòng)，減少趨勢(shì)均不顯著。而雨季、旱季降水日數(shù)在1995年之前在多年均線附近波動(dòng)，之后呈顯著減少趨勢(shì)，降水日數(shù)在研究時(shí)段內(nèi)整體呈現(xiàn)顯著的減少趨勢(shì)。雨季降水量的減少可以緩解流域內(nèi)防洪調(diào)控壓力，但同時(shí)也要注意短歷時(shí)極端降水事件可能造成的洪峰。旱季降水量的減少則會(huì)對(duì)下游旱季工農(nóng)業(yè)用水以及航運(yùn)等造成很大的壓力。

（3）多年平均降水量與經(jīng)度（P=0.081）、緯度（P=0.367）的相關(guān)性均不顯著，而與海拔高度（P=1.90E-4）達(dá)到極顯著相關(guān)水平。多年平均降水日數(shù)與經(jīng)度（P=0.539）相關(guān)系不顯著，而與緯度（P=8.77E-4）和海拔高度（P=1.82E-12）均達(dá)到極顯著相關(guān)。降水量和降水日數(shù)均與海拔高度極顯著正相關(guān)，可以利用海拔高度預(yù)測(cè)研究區(qū)內(nèi)年降水和降水日數(shù)的空間分布特征，海拔每升高100 m，對(duì)應(yīng)降水量與降水日數(shù)分別約增加 30 mm、4.5 d。

4.2 討論

（1）在本研究空間尺度范圍內(nèi)，對(duì)降水特性空間分布特征影響的顯著性排序?yàn)楹０?緯度>經(jīng)度，也即局地范圍內(nèi)的降水空間分布主要受地形的影響，但在區(qū)域或全球尺度范圍內(nèi)，降水的空間分布主要受到緯度或經(jīng)度的影響。有研究發(fā)現(xiàn)在全國尺度上夏季降水量與海拔高度之間是負(fù)相關(guān)的，并且降水量在海拔高度上的分布還存在分段現(xiàn)象（盧愛剛等，2008）。

（2）地形是影響降水的重要因子。山地對(duì)水汽具有阻擋抬升作用，容易形成降雨。隨著海拔的升高，降水量與降水日數(shù)均有增加趨勢(shì)，由此構(gòu)成了山地的水塔功能（魯春霞等，2007）。在水資源嚴(yán)重缺乏的今天，山地的水塔功能受到越來越多的關(guān)注（Immerzeel et al.，2010；Messerli et al.，2004），對(duì)降水量與海拔之間關(guān)系的研究可以幫助更加全面的掌握區(qū)域降水的空間分布，為區(qū)域水資源綜合利用以及低地河流防洪調(diào)控提供參考。三峽庫區(qū)監(jiān)測(cè)站點(diǎn)降水量、降水日數(shù)與站點(diǎn)海拔高度之間關(guān)系非常密切，在分析水資源空間分布時(shí)要充分考慮地形的影響。

（3）對(duì)降水空間分布特性進(jìn)行分析時(shí)，由于能夠獲取的高海拔氣象站點(diǎn)較少，所以研究結(jié)論可能會(huì)存在一定的不確定性。

IMMERZEEL W W, VAN BEEK L P, BIERKENS M F. 2010. Climate change will affect the Asian water towers. Science [J]. 328(5984): 1382-1385.

MESSERLI B, VIVIROLI D, WEINGARTNER R. 2004. Mountains of the world: vulnerable water towers for the 21st century. Ambio [J]. 13(7): 29-34.

馮亞文, 任國玉, 劉志雨, 等. 2013. 長江上游降水變化及其對(duì)徑流的影響[J]. 資源科學(xué), 35(6): 1268-1276.

顧鈞禧, 章基嘉, 巢紀(jì)平. 1994. 大氣科學(xué)詞典[M]. 北京: 氣象出版社.

李珊珊, 張明軍, 汪寶龍, 等. 2012. 近 51 年來三江源區(qū)降水變化的空間差異[J]. 生態(tài)學(xué)雜志, 31(10): 2635-2643.

廖要明, 張強(qiáng), 陳德亮. 2008. 1951—2006 年三峽庫區(qū)夏季氣候特征[J].氣候變化研究進(jìn)展, 3(6): 368-372.

盧愛剛, 康世昌, 龐德謙, 等. 2008. 中國降水量高度效應(yīng)及全球升溫對(duì)它的影響[J]. 生態(tài)環(huán)境, 17(5): 1875-1878.

魯春霞, 王菱, 謝高地, 等. 2007. 青藏高原降水的梯度效應(yīng)及其空間分布模擬[J]. 山地學(xué)報(bào), 25(6): 655-663.

潘磊, 唐萬鵬, 肖文發(fā), 等. 2012. 三峽庫區(qū)不同退耕還林模式林地水文效應(yīng)[J]. 水土保持通報(bào), 32(5): 103-106.

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

張樹奎, 魯子愛, 張楠. 2013. 三峽水庫蓄水對(duì)庫區(qū)降水量的影響分析[J]. 水電能源科學(xué), 31(5): 21-23.

張?zhí)煊? 范莉, 程炳巖, 等. 2010. 1961─2007年三峽庫區(qū)極端降水指數(shù)R5d 的變化規(guī)律及其未來情景預(yù)估[J]. 中國農(nóng)業(yè)氣象, 31(2): 176-182.

鄭剛, 劉祥梅, 郭志華, 等. 2009. 基于GIS的三峽庫區(qū)生態(tài)環(huán)境綜合評(píng)價(jià)——Ⅳ. 降水量變化 (1951-2004) [J]. 中國農(nóng)業(yè)氣象, 30(4): 486-491.

Spatiotemporal Distribution of Precipitation in Three Gorges Reservoir Area (TGRA) during 1951─2012

CHEN Xiangyi, XIAO Wenfa, HUANG Zhilin, ZENG Lixiong
Key Laboratory of Forest Ecology and Environment of State Forestry Administration//Research Institute of Forest Ecology, Environment and Protection, Chinese Academy of Forestry, Beijing 100091, China

Analysis on the spatial and temporal distribution of precipitation could provide useful reference for the management and utilization of regional water resources. By using monthly precipitation data of 27 weather stations within and around Three Gorges Reservoir Area (TGRA), the changing trend of precipitation in TGRA from 1951 to 2012 (62 years) were analyzed through linear tendency estimation and moving average methods. The spatial distribution of precipitation and total number of precipitation days were also analyzed through multivariate regression method. The results showed that the average annual precipitation, average annual precipitation in the rainy season, average annual precipitation of the dry season and average annual rainy days all presented a decreasing tendency. The average annual rainy days decreased more significantly than the average annual precipitation, and the start time of average annual rainy days was earlier than that of average annual precipitation, which indicated that the precipitation depth in a rainy day would increase, which also means that there will be more intensified rainstorms. The average annual precipitation was found has no significant correlations with longitude (P=0.081) and latitude (P=0.367), while has a strong correlation with the altitude (P=1.90E-4). Furthermore, the average annual rainy days were also found to have a strong positive correlation with both the latitude (P=8.77E-4) and altitude (P=1.82E-12), while has no significant correlation with longitude (P=0.539). Owing to the strong positive correlations between altitude, precipitation and the rainy days, the altitude could be possibly used for predicting the spatial distribution of precipitation and rainy days. According to the forecast model, precipitation and rainy days would increase about 30mm and 4.5 days respectively, with each increase of 100m in the altitude

Three Gorges Reservoir Area; precipitation; rainy days; altitude; linear trend; multivariate regression analysis

10.16258/j.cnki.1674-5906.2015.08.008

X16；P426

A

1674-5906（2015）08-1310-06

陳祥義，肖文發(fā)，黃志霖，曾立雄. 1951─2012年三峽庫區(qū)降水時(shí)空變化研究[J]. 生態(tài)環(huán)境學(xué)報(bào), 2015, 24(8): 1310-1315.

CHEN Xiangyi, XIAO Wenfa, HUANG Zhilin, ZENG Lixiong. Spatiotemporal Distribution of Precipitation in Three Gorges Reservoir Area (TGRA) during 1951─2012 [J]. Ecology and Environmental Sciences, 2015, 24(8): 1310-1315.

“十二五”科技支撐計(jì)劃課題（2015BAD07B04）

陳祥義（1986年生），男，博士研究生，主要從事流域水文方面研究。E-mail: chenxiangyichen@163.com *通信作者：黃志霖，男，副研究員，主要從事景觀生態(tài)、水土保持研究。Email: hzlin66@163.com

2015-03-23