肖鶯 陳晨 杜良敏

（武漢區(qū)域氣候中心，武漢 430074）

0 引言

長江上游五大流域分別為岷沱江流域、嘉陵江流域、烏江流域、宜賓到重慶區(qū)間流域、重慶到宜昌區(qū)間流域，其分布在四川、貴州、湖北及重慶境內(nèi)，地形、地貌十分復(fù)雜。受上游復(fù)雜地形、地貌的影響，流域強(qiáng)降水常造成長江干流水位猛漲、流量加大，在長江高水位的情況下，還可能形成致洪暴雨。此外，該流域強(qiáng)降水具有明顯的區(qū)域性特征，且降水強(qiáng)度大、過程多。國內(nèi)學(xué)者先后進(jìn)行了一系列研究。有些著重于研究流域面雨量方法，如畢寶貴等和方慈安等[1-2]；有些著重于面雨量時(shí)空分布特征，如王仁喬等研究了上游六大流域強(qiáng)降水面雨量特征，指出強(qiáng)降水面雨量出現(xiàn)頻次具有明顯的季節(jié)變化和年際變化特征；于大峰提出了長江上游流域干濕季節(jié)轉(zhuǎn)換指標(biāo)和計(jì)算方法，認(rèn)為夏季風(fēng)由東南向西北推進(jìn)緩慢，冬季風(fēng)由北向南推進(jìn)迅速；高琦等研究結(jié)果表明，不同流域強(qiáng)降水面雨量等級分布差異較大、極值分布差異明顯，重慶至萬縣、萬縣至宜昌兩個(gè)流域月極大值明顯大于其他流域[3-5]。有些對長江流域出現(xiàn)的強(qiáng)降水過程展開了分析，如施望芝等客觀分析了116次強(qiáng)降水的500 hPa環(huán)流形勢場和700 hPa天氣系統(tǒng)，并歸納出六大流域不同季節(jié)強(qiáng)降水的天氣模型及天氣系統(tǒng)；張萍萍等對30次長江上游與洞庭湖洪水遭遇過程進(jìn)行了分析，建立了5類天氣學(xué)概念模型[6-7]。這些研究都是基于氣象基本觀測事實(shí)展開“計(jì)算方法—時(shí)空分布—成因機(jī)理”的一系列研究，對不同空間、不同時(shí)間尺度降水開展分析，并建立了相關(guān)天氣模型。本文將針對長江流域水文氣象預(yù)報(bào)服務(wù)需求，著手于水庫調(diào)度所關(guān)注的強(qiáng)降水發(fā)生時(shí)間，重點(diǎn)開展長江上游五大流域之一——嘉陵江首末場強(qiáng)降水時(shí)間特征分析及預(yù)測方法研究。

1 資料與方法

采用的資料有：1961—2017年共55年嘉陵江49個(gè)氣象站點(diǎn)逐日（20時(shí)至次日20時(shí)）雨量實(shí)況資料，中國氣象局130項(xiàng)氣候系統(tǒng)監(jiān)測指數(shù)集。所選站點(diǎn)基本均勻分布（圖1）。

圖1 長江上游嘉陵江49個(gè)代表氣象站空間分布Fig. 1 Locations of 49 meteorological stations in Jialing River

面雨量的計(jì)算方法：基于逐日雨量資料，采用算術(shù)平均法，分別計(jì)算出嘉陵江的面雨量。

采用百分位法定義長江上游嘉陵江強(qiáng)降水的閾值。具體做法為：基于連續(xù)兩天平均面雨量，將面雨量按升序排列，取排序后降水序列的第95個(gè)百分位值定義為強(qiáng)降水閾值。計(jì)算得到嘉陵江強(qiáng)降水閾值為24.5 mm。當(dāng)各流域面雨量超過各自對應(yīng)的閾值時(shí)，記為一次強(qiáng)降水過程。根據(jù)強(qiáng)降水過程計(jì)算方法，統(tǒng)計(jì)歷史上強(qiáng)降水過程。選取每年的第一場強(qiáng)降水過程發(fā)生首日為首場強(qiáng)降水時(shí)間，最后一場強(qiáng)降水過程發(fā)生首日為末場強(qiáng)降水時(shí)間。

采用逐步回歸方法實(shí)現(xiàn)首末場強(qiáng)降水時(shí)間預(yù)測。具體做法：基于嘉陵江首末場強(qiáng)降水時(shí)間歷史序列和相應(yīng)時(shí)段的130項(xiàng)氣候系統(tǒng)監(jiān)測指數(shù)，利用逐步回歸方法，建立首末場強(qiáng)降水時(shí)間的回歸方程。將預(yù)報(bào)年份的氣候系統(tǒng)監(jiān)測指數(shù)代入回歸方程，得到首末場強(qiáng)降水時(shí)間的預(yù)測結(jié)論。

2 嘉陵江首末場強(qiáng)降水時(shí)間特征分析

根據(jù)定義，計(jì)算得到長江上游嘉陵江首場強(qiáng)降水時(shí)間歷史序列。從歷史平均來看（表1），嘉陵江首場強(qiáng)降水開始日發(fā)生在5月6日，末場強(qiáng)降水開始日發(fā)生在9月26日。首場強(qiáng)降水最早發(fā)生在4月2日（1973年），最晚發(fā)生在6月19日（2009年），兩者時(shí)間可相差78 d；末場強(qiáng)降水最早發(fā)生在7月1日（1997年），最晚發(fā)生在11月18日（1961年），兩者時(shí)間可相差119 d。利用功率譜方法提取首末場強(qiáng)降水時(shí)間逐年序列的顯著周期，結(jié)果顯示：嘉陵江首場強(qiáng)降水同時(shí)具有年際和年代際尺度周期，表現(xiàn)為準(zhǔn)2 a和11～18 a周期；末場強(qiáng)降水具有顯著的年代際尺度周期，分別為7～11、7～12和8～14 a周期。

表1 長江上游嘉陵江首末場強(qiáng)降水時(shí)間統(tǒng)計(jì)表Table 1 Climatic date of the first and last heavy rainfall in the Jialing River of upper Yangtze River

3 五大流域首末場強(qiáng)降水時(shí)間預(yù)測方法研究

根據(jù)三峽水庫服務(wù)需求，制作嘉陵江首場強(qiáng)降水時(shí)間預(yù)測產(chǎn)品的時(shí)間為同年的3月，制作嘉陵江末場強(qiáng)降水時(shí)間預(yù)測產(chǎn)品的時(shí)間為同年的8月。為了盡可能地囊括影響因子，采用提前一年的氣候系統(tǒng)監(jiān)測指數(shù)，即上一年9月至當(dāng)年提前1個(gè)月的指數(shù)集，進(jìn)行首末場強(qiáng)降水預(yù)測。利用1961—2005年首末場強(qiáng)降水時(shí)間序列與相應(yīng)時(shí)段的氣候系統(tǒng)監(jiān)測指數(shù)建立逐步回歸方程，2006—2017年做獨(dú)立樣本檢驗(yàn)。

圖2 給出了首末場強(qiáng)降水時(shí)間的擬合和獨(dú)立樣本檢驗(yàn)情況。從擬合情況看，首場、末場均是有42個(gè)指數(shù)進(jìn)入了逐步回歸方程。擬合結(jié)果顯示，首場和末場強(qiáng)降水時(shí)間回歸值與觀測值之間幾乎完全擬合，兩者之間的相關(guān)系數(shù)高達(dá)1.0；統(tǒng)計(jì)檢驗(yàn)表明回歸方程是顯著的，通過0.001的顯著性水平檢驗(yàn)。獨(dú)立樣本檢驗(yàn)顯示，首場強(qiáng)降水時(shí)間預(yù)測值與實(shí)況值的相關(guān)系數(shù)為0.56，通過0.05的顯著性水平檢驗(yàn)，從符號一致率看，正確率為66.7%，表明12年內(nèi)有8年的趨勢預(yù)報(bào)是正確的；末場強(qiáng)降水時(shí)間預(yù)測值與實(shí)況值的相關(guān)系數(shù)較首場有所下降，為0.4，未通過0.05顯著性水平檢驗(yàn)，從符號一致率看，正確率為66.7%，同樣的也是12年內(nèi)有8年的趨勢預(yù)報(bào)正確。

4 結(jié)論

圖2 嘉陵江首場（a）、末場（b）強(qiáng)降水時(shí)間的多元線性回歸曲線（藍(lán)線）、獨(dú)立樣本檢驗(yàn)（紅線：預(yù)測值，黑線：實(shí)況值）Fig. 2 Variations of the first (a) and last (b) heavy rainfall date (blue line) and independent sample test results(red line∶ forecasting values, black line∶ observations)

基于長江上游嘉陵江氣象站點(diǎn)雨量資料和中國氣象局130項(xiàng)氣候系統(tǒng)監(jiān)測指數(shù)集，利用功率譜和逐步回歸方法，開展了嘉陵江首末場強(qiáng)降水時(shí)間特征分析及預(yù)測方法研究。得到以下結(jié)論：

1）嘉陵江首末場強(qiáng)降水時(shí)間具有顯著的年際差異。首場強(qiáng)降水時(shí)間最早和最晚可相差兩個(gè)月以上（78 d），末場強(qiáng)降水時(shí)間最早和最晚甚至可相差3個(gè)月以上（119 d）。

2）首末場強(qiáng)降水時(shí)間周期特征顯著。嘉陵江首場強(qiáng)降水同時(shí)具有年際和年代際尺度周期，表現(xiàn)為準(zhǔn)2 a和11～18 a周期；末場強(qiáng)降水具有顯著的年代際尺度周期，分別為7～11、7～12和8～14 a周期。

3）基于氣候系統(tǒng)監(jiān)測指數(shù)建立的首末場強(qiáng)降水時(shí)間逐步回歸方程模擬效果好，復(fù)相關(guān)系數(shù)均通過了0.001顯著性水平檢驗(yàn)。2006—2017年獨(dú)立樣本檢驗(yàn)結(jié)果顯示：首場強(qiáng)降水開始時(shí)間預(yù)測值與實(shí)況值相關(guān)系數(shù)通過了0.05顯著性水平檢驗(yàn)；首場和末場預(yù)測距平符號一致率均為66.7%。

本文的預(yù)測方法研究僅初步利用中國氣象局下發(fā)的指數(shù)集建立回歸模型。從結(jié)果可以看出，本研究成果對于嘉陵江首場、末場強(qiáng)降水時(shí)間趨勢預(yù)測，有一定的預(yù)報(bào)技巧。在實(shí)際業(yè)務(wù)中，還需結(jié)合延伸期預(yù)報(bào)來一并做出最終結(jié)論。