郭家力，郭生練，郭 靖，陳 華

（武漢大學(xué)水資源與水電工程科學(xué)國家重點實驗室，武漢 430072）

鄱陽湖流域未來降水變化預(yù)測分析

郭家力，郭生練，郭 靖，陳 華

（武漢大學(xué)水資源與水電工程科學(xué)國家重點實驗室，武漢 430072）

利用鄱陽湖流域的13個國家氣象站1961-2001年的實測降水?dāng)?shù)據(jù)和NCEP再分析數(shù)據(jù)，建立了鄱陽湖流域降水的統(tǒng)計降尺度模型；在IPCC 2000年排放情景特別報告（SRES）中的A2和B2排放情景下，應(yīng)用HadCM3的輸出數(shù)據(jù)，預(yù)測鄱陽湖流域未來3個時段（2010-2039年、2040-2069年、2070-2099年）的降水變化情況。結(jié)果表明：鄱陽湖流域大部分區(qū)域的降水量有所增加，在本世紀(jì)末最大可能增加11.15%。

氣候變化；統(tǒng)計降尺度；降水預(yù)測；鄱陽湖流域

鄱陽湖是我國最大的淡水湖泊，在我國被稱為“大陸之腎”，蘊藏著豐富的水資源、濕地資源、生物資源等自然資源。鄱陽湖還是中國最大的濕地之一，是最大的候鳥集散地、棲息地。鄱陽湖豐富的自然資源和良好的生態(tài)環(huán)境，是實現(xiàn)江西省崛起新跨越的突破點和驅(qū)動力。鄱陽湖優(yōu)越而重要的地理位置和320億m3的容積，是長江防洪體系重要的組成部分。鄱陽湖被4 300萬贛鄱人民稱為“母親湖”，其在江西省乃至整個中部地區(qū)所處的戰(zhàn)略地位可見一斑。IPCC（Intergovernmental Panel on Climate Change，政府間氣候變化專門委員會）第四次評估報告顯示，全球變暖毋庸置疑，極端天氣事件出現(xiàn)頻率增加，全球氣候系統(tǒng)正發(fā)生顯著的變化［1］。大氣環(huán)流模式（GCM）是目前研究氣候變化與氣候變異最重要也是最可行的方法。GCM預(yù)測的氣候變化情景可以通過降尺度處理得出區(qū)域的氣候變化信息［2］。本文應(yīng)用統(tǒng)計降尺度方法建立表征氣候變化信息的大尺度氣候因子與鄱陽湖流域氣象站點日降水之間的統(tǒng)計學(xué)關(guān)系，預(yù)測分析鄱陽湖流域降水特點及其變化規(guī)律，研究氣候變化影響下鄱陽湖流域的水、氣相互作用機理，以應(yīng)對氣候變化或異常形成的自然災(zāi)害，保障地區(qū)的水資源可持續(xù)利用和國民經(jīng)濟的可持續(xù)發(fā)展。

1 研究區(qū)域概況及數(shù)據(jù)資料

鄱陽湖位于長江中下游南岸的江西省北部，地理坐標(biāo)E115°49′至116°46′、N28°24′至29°46′。鄱陽湖以湖口通長江，是一個過水性吞吐型湖泊，湖口站歷史最高水位為1998年7月30日的22.59 m（吳淞基面高程），相應(yīng)湖面面積4 070 km2，總?cè)莘e為320億m3。鄱陽湖和贛江、撫河、信江、饒河、修水5大河流等組成了鄱陽湖水系，流域面積16.22萬km2，約占江西省總面積的97.2%，占長江流域面積的9%，其水系年均徑流量為1 525億m3，約占長江流域年均徑流量的16.3%。鄱陽湖水系為長江流域的多雨區(qū)，多年平均降水量約為1 620 mm；降水量年際變化大，豐水年1954年降水量達2 100 mm，枯水年1963年降水量僅1 136 mm；年內(nèi)分配不均勻，4-6月降水量占全年降水量的50%［3］。

研究需要的數(shù)據(jù)包括3類：氣象站點實測資料、大尺度氣候觀測資料及GCM模擬的當(dāng)前和未來氣候條件下的輸出數(shù)據(jù)。選取了鄱陽湖流域內(nèi)共13個國家氣象站點（見圖1），站點分布比較均勻，基本上能代表該區(qū)域的氣候特征。選取1961-2001年為研究基準(zhǔn)期，其中1961-1990年為模型率定期，1991-2001年為模型檢驗期。根據(jù)鄱陽湖流域“一湖通五水”的地理特點，將13個站點劃分到6個分區(qū)，宜春、吉安、遂川、贛州和樟樹5站位于贛江流域，南城、廣昌位于撫河流域，貴溪、玉山位于信江流域，修水位于修水流域，景德鎮(zhèn)位于饒河流域，波陽、南昌位于湖區(qū)。

在氣候變化研究中，NCEP（National Center for Environment Prediction，美國環(huán)境預(yù)報中心）再分析數(shù)據(jù)可作為診斷數(shù)據(jù)進行一系列的分析和研究，尤其在一些全球、區(qū)域氣候變化的模擬和預(yù)測中，以NCEP再分析數(shù)據(jù)作為氣候模式的初始場和驅(qū)動場數(shù)據(jù)［4］。本文采用NCEP的全球再分析日數(shù)據(jù)（該數(shù)據(jù)可以從國際互聯(lián)網(wǎng)鏈接http：／／www.esrl.noaa.gov／psd／data／gridded／reanalysis／免費下載）作為觀測的大尺度氣候數(shù)據(jù)，其空間分辨率為3.75°× 2.5°。根據(jù)已有的研究表明，大尺度氣候模式網(wǎng)格的范圍一般都比研究區(qū)域的范圍要大很多，因此選擇覆蓋鄱陽湖流域的經(jīng)緯網(wǎng)格及其周圍的網(wǎng)格共4×4個NCEP經(jīng)緯網(wǎng)格（在經(jīng)度方向上的網(wǎng)格編號是30～33，緯度方向上的網(wǎng)格編號是22～25）。站點位置和NCEP經(jīng)緯網(wǎng)格的劃分如圖1所示。徐影等的研究結(jié)果表明，英國Hadley中心的HadCM3模式對東亞地區(qū)模擬效果較好［5］。

圖1 鄱陽湖流域站點分布及NCEP網(wǎng)格Fig.1 Nationalweather stations distribution and NCEP grids for Poyang Lake basin

2 統(tǒng)計降尺度模型

統(tǒng)計降尺度方法是基于這樣一種基本假設(shè)：區(qū)域氣候是由當(dāng)?shù)氐牡乩硖匦院痛蟪叨葰庀鬆顩r所決定的。有鑒于此，區(qū)域觀測氣象變量（如降水量）和觀測的大尺度氣候因子之間的統(tǒng)計關(guān)系就可能建立。GCM模擬輸出的關(guān)于大尺度的氣候狀況再聯(lián)合這種統(tǒng)計關(guān)系，就能夠降尺度到區(qū)域?qū)哟巍；贛ATLAB環(huán)境開發(fā)的統(tǒng)計降尺度模型ASD（Automated Statistical Downscalingmodel）［6］的結(jié)構(gòu)簡圖如圖2所示。

ASD模型可以在某一事件發(fā)生的條件下（如降水）或無條件下（如氣溫）運行。因此，日降水的模擬如式（1）所示，分為2個步驟：先模擬降水發(fā)生概率，如果有降水發(fā)生再模擬降水量。

式中：Oi為日降水發(fā)生概率；Ri為日降水量；pij為預(yù)報因子；n為預(yù)報因子的數(shù)量；α，β為模型參數(shù)；ei為模型誤差，并假設(shè)其服從高斯分布。

基于回歸的降尺度方法一般采用多元線性回歸，但是當(dāng)預(yù)報因子之間的相關(guān)性很強時，會使回歸系數(shù)的最小二乘估計不穩(wěn)定。因此，ASD模型提供了2種回歸方法建立預(yù)報因子和預(yù)報量之間的統(tǒng)計關(guān)系，一般情況下采用多元線性回歸，當(dāng)預(yù)報因子之間有很強的相關(guān)性時，使用嶺回歸方法。該模型程序還具有將模擬的氣候情景結(jié)果進行統(tǒng)計分析、制圖輸出等功能。

圖2 ASD模型結(jié)構(gòu)圖Fig.2 Structural diagram of ASD model

2.1 因子選擇

選擇站點所在網(wǎng)格及其周圍的8個網(wǎng)格共9個NCEP網(wǎng)格（規(guī)定這9個網(wǎng)格為站點相關(guān)網(wǎng)格）的大氣環(huán)流因子，每個網(wǎng)格包含如海平面氣壓、地面氣溫、500 hPa和850 hPa位勢高度場、500 hPa和850 hPa濕度等26個因子（規(guī)定為備選因子），則一共有9× 26＝234個因子，由于多個網(wǎng)格和多個預(yù)報因子組成了一個多維的大尺度預(yù)報因子數(shù)據(jù)集，如果不進行降維和特征因子的選擇，必將遭遇計算災(zāi)難。在ASD模型中提供了向后逐步回歸和偏相關(guān)系數(shù)2種方法來進行因子選擇，本文采用的是后者。選擇過程主要分為3個步驟：①自動計算。分別將9個站點相關(guān)網(wǎng)格的備選因子與站點實測降雨建立統(tǒng)計關(guān)系，將各數(shù)據(jù)輸入ASD模型，進行模型自動的因子選擇。②網(wǎng)格篩選。模型會根據(jù)上一步篩選的因子，計算出站點相關(guān)網(wǎng)格內(nèi)的預(yù)報因子與預(yù)報量之間的解釋方差，該指標(biāo)表示了預(yù)報量與環(huán)流因子之間的相關(guān)性大小。根據(jù)這一指標(biāo)選擇相關(guān)性最好的3個網(wǎng)格，并選擇站點所在網(wǎng)格共4個網(wǎng)格。③因子選擇。在第②步的基礎(chǔ)上進行進一步的降維，通過比較4個網(wǎng)格內(nèi)各因子的偏相關(guān)系數(shù)，選擇偏相關(guān)系數(shù)大于0.15的因子，最終將因子總數(shù)控制在26個以內(nèi)（因為ASD模型要求輸入的最大因子數(shù)目為26）。

2.2 模型率定

模型率定需要的資料包括大尺度網(wǎng)格的NCEP再分析數(shù)據(jù)和降尺度目標(biāo)（站點或者流域）的降水量。通過建立實測降水與NCEP再分析資料之間的統(tǒng)計關(guān)系，即確定統(tǒng)計降尺度模型的參數(shù)，完成模型的率定。選擇2個評價指標(biāo)：均值相對誤差REmean和標(biāo)準(zhǔn)差相對誤差REstd來評定所建立的統(tǒng)計關(guān)系的好壞。本文的降尺度目標(biāo)為站點降水，將率定期1961-1990年的NCEP再分析數(shù)據(jù)降尺度到站點，得到各站點的降水過程。因為所建立的統(tǒng)計關(guān)系是基于站點的，所以，在這種關(guān)系建立以后，無論是以站點降水還是面降水進行結(jié)果展示，都不會影響模型的精度。為節(jié)省篇幅，將各站點的降水（包括實測和模擬）經(jīng)算術(shù)平均得到鄱陽湖流域6個分區(qū)的降水量，表1列出了鄱陽湖流域“五河”（贛江、饒河、撫河、修水和信江）和湖區(qū)實測與NCEP模擬多年平均月降水序列的均值、標(biāo)準(zhǔn)差及相對誤差。根據(jù)表1顯示，修水、饒河、信江、贛江、撫河和湖區(qū)均值相對誤差的絕對值多月平均分別為6.62%，9.19%，3.18%，4.49%，5.75%和5.62%；標(biāo)準(zhǔn)差相對誤差的絕對值多月平均分別為11.93%，14.55%，7.42%，8.01%，7.11%和8.51%，這表明由NCEP再分析資料率定的降尺度統(tǒng)計關(guān)系能較好地模擬鄱陽湖流域的降水過程。

表1 鄱陽湖流域各分區(qū)率定期實測和NCEP模擬降水對比分析Table1 The com parison of observed precipitation and estimated data w ith NCEP during calibration period for Poyang Lake basin

2.3 模型驗證

將率定好的ASD模型應(yīng)用于預(yù)報因子的觀測數(shù)據(jù)（NCEP再分析數(shù)據(jù)），根據(jù)輸入自動生成對應(yīng)的預(yù)報量序列。對比分析生成序列與實測序列的擬合程度，即使用獨立觀測資料對該統(tǒng)計降尺度模型進行驗證。同時，為對比實測降水和GCM模擬降水之間的差異，以選擇模擬未來降水的基準(zhǔn)，本節(jié)模擬了以IPCC 2000年排放情景特別報告（SRES）中的A2和B2情景下GCM輸出值作為降尺度模型輸入的降水過程。用模型模擬出各站點的降雨量之后，由算術(shù)平均法求得各分區(qū)的實測和模擬降水量。

圖3列出了檢驗期（1991-2001年）各分區(qū)的實測、NCEP模擬以及A2和B2情景下HadCM3模擬的多年平均月降水序列，從圖中可以看出，無論是NCEP再分析數(shù)據(jù)，還是A2和B2情景下HadCM3輸出數(shù)據(jù)，在枯水期（7月至次年3月）均能較好地擬合實測降水；在豐水期（4-6月），NCEP再分析數(shù)據(jù)模擬的降水同樣能較好地擬合實測降水，盡管GCM輸出數(shù)據(jù)模擬的降水較實測降水偏大，但月平均降水偏大幅度均在可接受范圍以內(nèi)，因此可以應(yīng)用HadCM3輸出數(shù)據(jù)進行降尺度分析，預(yù)測鄱陽湖流域的未來降水變化情況。

圖3 檢驗期鄱陽湖流域各分區(qū)實測數(shù)據(jù)、NCEP模擬、A2／B2情景下HadCM 3模擬降水對比Fig.3 The com parison of observed precipitations and estimated data w ith NCEP and GCM output data under A2 and B2 scenarios during validation period for the Poyang Lake basin and every section

3 未來降水預(yù)測分析

王懷清等［7］采用線性回歸方法研究了鄱陽湖五大流域1960-2006年的年降水量等要素的變化趨勢，并用Mann方法對變化趨勢做了顯著性檢驗，用Mann-Kendall方法對各時間序列進行突變檢測，認為鄱陽湖各流域的年降水量變化趨勢基本一致，呈波動上升趨勢，但并未出現(xiàn)趨勢性的突變。

為分析鄱陽湖流域未來的降水情況，以1961-2001年為基準(zhǔn)期，針對未來3個時段，即2020s（2010-2039年）、2050s（2040-2069年）、2080s（2070-2099年），應(yīng)用HadCM3輸出的未來氣候情景，輸入經(jīng)驗證過的統(tǒng)計降尺度ASD模型生成未來序列。表2列出了A2和B2情景下鄱陽湖流域及其各分區(qū)未來3個時段的多年平均降水量P（mm）及其相對于基準(zhǔn)期的變幅（%），其中基準(zhǔn)期的降水量采用GCM資料模擬，使基準(zhǔn)期和未來的降水能在相同的模型誤差下進行比較?？傮w上來說，未來3個時段中無論是在A2還是在B2情景下，距離基準(zhǔn)期的時間越長，降水的變化幅度越大。從各分區(qū)來看，贛江流域在2種情景下3個時段的降水較基準(zhǔn)期均增加；饒河流域不同于其他分區(qū)，在A2和B2情景下均遞減；修水、撫河、信江、湖區(qū)4個分區(qū)均是在A2情景下降水先減后增，在B2情景下各時段降水都增加。鄱陽湖全流域的未來降水在A2和B2情景下較基準(zhǔn)期都有所增加，且隨著時間的推移，降水的增幅也隨著變大，最大增幅可達11.15%，約為200 mm。

表2 A2和B2情景下鄱陽湖流域及各分區(qū)未來3個時段的降水變化Table2 The precipitation change under A2 and B2 scenarios during three future periods for Poyang Lake basin and every section

4 結(jié) 語

本文利用鄱陽湖流域的13個國家氣象站41年的實測降水?dāng)?shù)據(jù)和NCEP再分析數(shù)據(jù)，建立了鄱陽湖流域降水統(tǒng)計降尺度ASD模型，在此模型基礎(chǔ)之上，采用A2和B2情景下的HadCM3輸出數(shù)據(jù)，預(yù)測鄱陽湖流域的未來降水變化，結(jié)論如下：

（1）在進行因子選擇時發(fā)現(xiàn)在所有包含預(yù)報因子的9個相關(guān)網(wǎng)格中，位于站點所在網(wǎng)格的右側(cè)、下側(cè)、右下側(cè)的3個網(wǎng)格中的預(yù)報因子與站點的聯(lián)系最緊密，分析認為這可能是與鄱陽湖流域降水受西太平洋副熱帶高壓和東南季風(fēng)的影響有關(guān)。

（2）建立的統(tǒng)計降尺度模型在率定期和檢驗期均能較好地模擬鄱陽湖流域的降水變化過程，各分區(qū)未來降水有增有減，距離基準(zhǔn)期的時間越長，降水的變化幅度越大。在A2和B2排放情景下，鄱陽湖流域的未來降水較基準(zhǔn)期均有所增加，到本世紀(jì)末最大可能增加11.15%。

［1］ IPCC.Summary for Policymakers of the Synthesis Report of the IPCC Fourth Assessment Report［M］.Cambridge，UK：Cambridge University Press，2007.

［2］ 陳 華，郭 靖，郭生練，等.應(yīng)用統(tǒng)計學(xué)降尺度方法預(yù)測漢江流域降水變化［J］.人民長江，2008，39（14）：51-55.（CHEN Hua，GUO Jing，GUO Sheng-lian，et al.Application of statistical downscalingmethod in precipitation prediction for the Hanjiang River basin［J］.Yangtze River，2008，39（14）：51-55.（in Chinese））

［3］ 朱宏富，金 鋒，李榮昉.鄱陽湖調(diào)蓄功能與防災(zāi)綜合治理研究［M］.北京：氣象出版社，2003.（ZHU Hongfu，JIN Feng，LIRong-fang.Regulation and Storage Function and Comprehensive Management of Disaster Preventaion for the Poyang Lake Basin［M］.Beijing：Meterorological Press，2003.（in Chinese））

［4］ 趙天保，艾麗坤，馮錦明.NCEP再分析資料和中國站點觀測數(shù)據(jù)的分析與比較［J］.氣候與環(huán)境研究，2004，9（2）：278-294.（ZHAO Tian-bao，AI Li-kun，F(xiàn)ENG Jin-ming.An intercomparison between NCEP reanalysis and observed data over China［J］.Climatic and Environmental Research，2004，9（2）：278-294.（in Chinese））

［5］ 徐 影，丁一匯，趙宗慈.近30年人類活動對東亞地區(qū)氣候變化影響的檢測與評估［J］.應(yīng)用氣象學(xué)報，2002，13（5）：513-525.（XU Ying，DING Yi-hui，ZHAO Zong-ci.Detection and evaluation of effect of human activities on climatic change in east asia in recent 30 yesrs［J］.Quarterly Journal of Applied Meteorlolgy，2002，13（5）：513-525.（in Chinese））

［6］ HESSAMI，M，GACHONP，OUARDA T BM J，et al.Automated regression-based statistical downscaling tool［J］.Environmental Modelling and Software，2008，23：813-834.

［7］ 王懷清，趙冠男，彭 靜，等.近50年鄱陽湖五大流域降水變化特征研究［J］.長江流域資源與環(huán)境，2009，18（7）：615-619.（WANG Huai-qing，ZHAO Guan-nan，PENG Jing，et al.Precipitation characteristics over five major river systems of Poyang Lake drainage areas in recent50 years［J］.Resources and Ennironment in the Yangtze River Basin，2009，18（7）：615-619.（in Chinese） ）

（編輯：趙衛(wèi)兵）

Prediction of Precipitation Change in Poyang Lake Basin

GUO Jia-li，GUO Sheng-lian，GUO Jing，CHEN Hua
（State Key Laboratory ofWater Resources and Hydropower Engineering Science，Wuhan University，Wuhan 430072，China）

The observed precipitation series from thirteen national weather stations during 1961 to 2001 and NCEP reanalysis data were used to establish a statistical downscalingmodel in the Poyang Lake basin.The output data of HadCM3 under A2 and B2 scenarioswere inputted to automatically statistical downscalingmodel to predict the future precipitation change during 2010-2039，2040-2069，2070-2099 periods.The results show that the precipitation will increase in themost area of the Poyang Lake basin，and maximum increase rate can reach 11.15%in the end of this century.

climate change；statistical downscaling；precipitation prediction；Poyang Lake basin

P333.9

A

1001-5485（2010）08-0020-05

2009-10-05；

2010-01-27

國家自然科學(xué)基金（50809049）；教育部高等學(xué)校博士學(xué)科點專項科研基金（200804861062）

郭家力（1984-），男，湖北孝感人，碩士研究生，研究方向為流域水文模擬，（電話）027-68772066，13237190881（電子信箱）culture _me＠163.com。