楊 杰 王啟光 支 蓉 封國林2)?

1)(蘭州大學(xué)大氣科學(xué)學(xué)院，蘭州 730000)

2)(中國科學(xué)院大氣物理研究所東亞區(qū)域氣候-環(huán)境重點(diǎn)實(shí)驗(yàn)室，北京100029)

3)(國家氣候中心，中國氣象局氣候研究開放實(shí)驗(yàn)室，北京100081)

(2010年6月22日收到;2010年7月15日收到修改稿)

動(dòng)態(tài)最優(yōu)多因子組合的華北汛期降水模式誤差估計(jì)及預(yù)報(bào)*

楊 杰1)2)王啟光1)支 蓉3)封國林2)3)?

1)(蘭州大學(xué)大氣科學(xué)學(xué)院，蘭州 730000)

2)(中國科學(xué)院大氣物理研究所東亞區(qū)域氣候-環(huán)境重點(diǎn)實(shí)驗(yàn)室，北京100029)

3)(國家氣候中心，中國氣象局氣候研究開放實(shí)驗(yàn)室，北京100081)

(2010年6月22日收到;2010年7月15日收到修改稿)

利用國家氣候中心季節(jié)預(yù)報(bào)1983—2009年27年模式預(yù)報(bào)結(jié)果，結(jié)合74項(xiàng)環(huán)流指數(shù)及美國國家海洋局和大氣管理局提供的40個(gè)氣候指數(shù)和美國氣候預(yù)報(bào)中心實(shí)際降水分析資料，采用資料診斷分析和數(shù)值模擬實(shí)驗(yàn)相結(jié)合的方法，通過多因子的歷史相似信息提取預(yù)報(bào)相似年，獲得預(yù)報(bào)場(chǎng)的誤差訂正項(xiàng).在這一訂正思路的基礎(chǔ)上，考慮前期關(guān)鍵影響因子的選取、多因子組合的優(yōu)化配置，構(gòu)建適用于不同預(yù)報(bào)年的區(qū)域動(dòng)力-統(tǒng)計(jì)模式預(yù)報(bào)誤差訂正方案.以華北為例，探索多因子最優(yōu)組合的多元客觀相似判據(jù)，發(fā)展基于多因子動(dòng)力-統(tǒng)計(jì)模式預(yù)報(bào)誤差的動(dòng)態(tài)訂正新技術(shù)，改善華北夏季降水預(yù)報(bào)效果，提高預(yù)報(bào)技巧.通過2005—2009年獨(dú)立樣本回報(bào)結(jié)果表明，動(dòng)態(tài)最優(yōu)多因子組合相似訂正方法距平相關(guān)系數(shù)評(píng)分相對(duì)于系統(tǒng)訂正方法有著顯著的提高，該訂正方案對(duì)華北地區(qū)的夏季降水預(yù)測(cè)有著很好的業(yè)務(wù)前景，具有重要的應(yīng)用價(jià)值，即將投入業(yè)務(wù)運(yùn)行.

關(guān)鍵因子集，模式誤差估計(jì)，汛期降水，預(yù)測(cè)

PACS:92.60.Wc

1.引 言

洪澇和干旱是影響我國最大的兩種主要?dú)夂驗(yàn)?zāi)害，50% 以上的自然災(zāi)害都是由旱、澇引起［1—3］.隨著氣候系統(tǒng)的不確定性演變，包括異常干旱、破紀(jì)錄洪災(zāi)、超低溫雨雪冰凍等極端氣候事件發(fā)生的頻率和強(qiáng)度不斷加劇，對(duì)國民經(jīng)濟(jì)影響巨大［4—7］.如2006，2007年夏季我國重慶地區(qū)分別遭受大旱和大澇，2009年東北夏季低溫，2009—2010年西南五省夏冬春三季持續(xù)干旱，造成了前所未有的損失.因此，提高短期氣候預(yù)測(cè)、預(yù)警以及應(yīng)急決策能力是保障國家安全、發(fā)展經(jīng)濟(jì)和防災(zāi)減災(zāi)的現(xiàn)實(shí)需要［8—10］.近代以來華北地區(qū)水資源貧乏問題變得越來越嚴(yán)重，如何有效利用和規(guī)劃該地區(qū)的有限的水資源是目前十分迫切的一個(gè)科學(xué)問題［11—13］.華北地區(qū)旱澇時(shí)空分布的不均勻性帶來的水資源利用的不確定性嚴(yán)重制約了華北地區(qū)社會(huì)和經(jīng)濟(jì)發(fā)展，自然生態(tài)系統(tǒng)和環(huán)境也在不斷惡化.解決這些問題的前提仍然是如何準(zhǔn)確預(yù)測(cè)夏季旱澇的分布趨勢(shì)，只有掌握了空中水資源活動(dòng)規(guī)律，才能從根本上緩解華北地區(qū)的水資源匱乏，確保該地區(qū)的可持續(xù)發(fā)展，這也是氣象科技工作者歷來關(guān)注的科學(xué)問題［14，15］.

主導(dǎo)華北地區(qū)夏季降水的氣候背景十分復(fù)雜，一方面由于華北處于季風(fēng)活動(dòng)的北緣，汛期降水主要受夏季風(fēng)活動(dòng)等影響;另一方面，該地區(qū)還受中高緯西風(fēng)環(huán)流系統(tǒng)影響，這些復(fù)雜因素導(dǎo)致預(yù)測(cè)每年華北夏季的旱澇趨勢(shì)難度很大［16，17］.因此，單純的統(tǒng)計(jì)學(xué)方法和動(dòng)力學(xué)方法都不能有效開展該地區(qū)的短期氣候預(yù)測(cè).丑紀(jì)范［18，19］認(rèn)為，統(tǒng)計(jì)學(xué)方法雖善于利用大量歷史資料，但其自身局限性使它并不能成為主要的預(yù)測(cè)工具;而動(dòng)力學(xué)方法雖然善于運(yùn)用物理規(guī)律，但同樣面臨著模式物理過程不夠準(zhǔn)確和初值不精確的問題.解決的關(guān)鍵是將兩者有機(jī)結(jié)合，取長(zhǎng)補(bǔ)短，通過充分利用歷史統(tǒng)計(jì)信息達(dá)到修正動(dòng)力學(xué)方法積分過程中的系統(tǒng)性誤差.如此一來，可以從很大程度上避免直接進(jìn)行數(shù)值模式動(dòng)力框架的大幅修改，減少對(duì)初值同化方案和物理參數(shù)的復(fù)雜處理，以相對(duì)較小的代價(jià)達(dá)到相同的模式預(yù)測(cè)效果.這種思想的核心是在動(dòng)力模式預(yù)報(bào)基礎(chǔ)上結(jié)合物理統(tǒng)計(jì)方法充分利用歷史資料對(duì)模式預(yù)報(bào)結(jié)果加以修正，目前已成為備受關(guān)注的研究方向.在國內(nèi)外該研究已取得了一定的進(jìn)展［20—23］，包括相似動(dòng)力方法［24—27］、預(yù)報(bào)誤差訂正以及模式后處理［28］、基于大氣自記憶原理的方法［29—31］等，使得短期氣候預(yù)測(cè)發(fā)展很快并已初步建立了適合中國條件的短期氣候預(yù)測(cè)業(yè)務(wù)動(dòng)力模式系統(tǒng)［32，33］，但業(yè)務(wù)應(yīng)用水平仍然不高，仍需進(jìn)一步提升［34—36］.

本文從動(dòng)力和統(tǒng)計(jì)相結(jié)合角度出發(fā)，利用1983—2009年美國氣候預(yù)報(bào)中心組合降雨分析資料(Climale Prediction Center Mergod Analysis of Precipitation，CMAP)和國家氣候中心(National Climate Center，NCC)季節(jié)/年際預(yù)測(cè)業(yè)務(wù)系統(tǒng)模式回報(bào)的資料挖掘歷史相似信息，以單因子交叉檢驗(yàn)的距平相關(guān)系數(shù)(anomaly correlation coefficient，ACC)確定關(guān)鍵影響因子集，建立適用于華北夏季降水的動(dòng)態(tài)相似因子庫，通過數(shù)值模擬與ACC檢驗(yàn)獲取相似因子配置方案.由于影響因子具有周期性動(dòng)態(tài)變化，不同因子影響季節(jié)不同，為了兼顧不同預(yù)報(bào)時(shí)段動(dòng)態(tài)因子的最優(yōu)組合，以4個(gè)歷史最相似年的信息為主訂正預(yù)報(bào)結(jié)果，形成最終的動(dòng)力-相似預(yù)報(bào)修正方案，達(dá)到改善華北地區(qū)夏季降水預(yù)報(bào)效果的目的.

2.資料和方法

2.1.資料

本文的研究資料來自三個(gè)方面:

1)美國氣候預(yù)報(bào)中心組合降雨分析資料的1983—2009年夏季總降水量，水平分辨率為2.5°× 2.5°;

2)國家氣候中心提供的季節(jié)/年際預(yù)測(cè)業(yè)務(wù)系統(tǒng)模式預(yù)報(bào)的歷史回報(bào)產(chǎn)品中的逐月降水場(chǎng)資料，經(jīng)緯度格點(diǎn)2.5°×2.5°，共144×73個(gè)格點(diǎn);

3)國家氣候中心氣候系統(tǒng)診斷預(yù)測(cè)室提供的74項(xiàng)環(huán)流特征量和美國國家海洋局和大氣管理局(National Oceanic and Atmospheric Administration，NOAA)提供的40項(xiàng)氣候指數(shù).

2.2. 模式預(yù)報(bào)誤差估計(jì)原理［37，38］

實(shí)際大氣的準(zhǔn)確數(shù)值模式表示為

其中ψ(x，t)為模式預(yù)報(bào)變量，x，t分別表示空間坐標(biāo)向量和時(shí)間，L是ψ的微分算子.t0為初始時(shí)刻，ψ0為初值.t>t0時(shí)刻的值可由初值進(jìn)行數(shù)值積分得到ψ或者其泛函M(ψ).E為模式的誤差算子，反映模式中未知的總誤差項(xiàng)，即模式預(yù)報(bào)誤差.將預(yù)報(bào)場(chǎng)ψ看作歷史相似場(chǎng)與疊加在相似場(chǎng)上的小擾動(dòng)ψ′之和，即ψ=+ψ′，并且將上述方程進(jìn)行積分推導(dǎo)可得到一個(gè)動(dòng)力相似預(yù)報(bào)方程

該方程引入歷史相似對(duì)應(yīng)的預(yù)報(bào)誤差信息估計(jì)當(dāng)前的模式預(yù)報(bào)誤差，即在模式預(yù)報(bào)結(jié)果的基礎(chǔ)上添加了一個(gè)歷史相似誤差訂正項(xiàng)

該式將動(dòng)力模式預(yù)報(bào)改進(jìn)問題轉(zhuǎn)化為模式誤差的估計(jì)問題，模式的預(yù)報(bào)即包含誤差預(yù)報(bào).

3.模式誤差訂正方法

3.1.華北地區(qū)氣候背景及系統(tǒng)誤差訂正方法

華北地區(qū)屬于半濕潤(rùn)、半干旱大陸性季風(fēng)氣候區(qū)，受東亞季風(fēng)的影響，降水量年內(nèi)分配很不均勻，多年平均的逐月降水量呈單峰型分布，1—7月份呈逐步遞增趨勢(shì)，8月以后逐漸減少，80%左右集中在6—9月，且往往集中在幾次強(qiáng)降雨過程，降雨相對(duì)集中使流域洪水具有峰值高、峰型陡、一次性洪量大的特點(diǎn)［39］.華北雨季主要集中在7月和8月，所以夏季降水的多寡程度將直接影響著華北地區(qū)全年的旱澇分布.降水量年際變化很大，豐水年可達(dá)1000 mm以上，枯水年僅200 mm左右.1951年以來，1979年以前為多雨時(shí)段，1979年以后處于少雨的時(shí)段.總體來說，目前仍處于少雨的氣候背景.圖1給出了1983—2009年華北地區(qū)夏季降水量的年際變化，對(duì)這27年的實(shí)際降水量進(jìn)行聚類分析，結(jié)果顯示，相對(duì)27年夏季降水情況而言，降水偏少年為 1983，1985，2002 年;降水偏多年為 1988，1995，1996，1998，2003，2004 年;降水異常偏少年為 1997，1999年;降水異常偏多年為1990年.對(duì)比1983—2009年華北地區(qū)的夏季降水量變化趨勢(shì)與模式預(yù)報(bào)結(jié)果，可以發(fā)現(xiàn)在整個(gè)華北地區(qū)處于少雨的氣候背景下，該地區(qū)每年的實(shí)際降水量都要比模式原始的預(yù)報(bào)結(jié)果小很多，且實(shí)際降水量的年際變化很明顯，最大的年降水量為370 mm，最小的年降水量為192 mm，二者之間相差達(dá)到約180 mm，而模式預(yù)報(bào)結(jié)果年際變化很小，基本上結(jié)果都是在470 mm上下波動(dòng).可見模式的原始預(yù)報(bào)結(jié)果與實(shí)況之間的誤差較大，若無適當(dāng)訂正，將失去應(yīng)有的參考價(jià)值.

模式原始結(jié)果訂正最直接的是預(yù)報(bào)誤差訂正，如能夠通過統(tǒng)計(jì)方法較為準(zhǔn)確地得到模式在預(yù)報(bào)年的預(yù)報(bào)誤差場(chǎng)，將預(yù)報(bào)誤差場(chǎng)與模式原始結(jié)果進(jìn)行疊加即可對(duì)預(yù)報(bào)結(jié)果有很好的改進(jìn).如前所述，降水預(yù)報(bào)問題就轉(zhuǎn)換成模式預(yù)報(bào)誤差的估計(jì)問題.最簡(jiǎn)單的預(yù)報(bào)誤差訂正方法是對(duì)統(tǒng)計(jì)歷史資料中已知的歷年模式預(yù)報(bào)誤差取平均，以多年預(yù)報(bào)誤差的平均值作為預(yù)報(bào)年的模式原始預(yù)測(cè)結(jié)果與實(shí)際降水的誤差來對(duì)預(yù)測(cè)結(jié)果進(jìn)行訂正，即系統(tǒng)誤差訂正方法.圖2(a)，(b)分別給出了1983—2009年交叉檢驗(yàn)的系統(tǒng)誤差訂正后的夏季降水的預(yù)報(bào)結(jié)果以及相應(yīng)的ACC評(píng)分情況，可以看出經(jīng)過系統(tǒng)訂正后的預(yù)報(bào)結(jié)果較模式原始結(jié)果有了較大的提高，預(yù)報(bào)誤差也有很大的改進(jìn)，交叉檢驗(yàn)系統(tǒng)訂正后預(yù)報(bào)結(jié)果的所有年份的距平相關(guān)系數(shù)平均值為0.15.結(jié)合圖1可以看出，系統(tǒng)誤差訂正結(jié)果與實(shí)際觀測(cè)誤差異常大的年份為:1985，1992，1997，1999，2002年，而其中 1992，1997，1999，2002 年都對(duì)應(yīng)了華北地區(qū)的降水偏少或異常偏少.我們可以看出系統(tǒng)訂正后的預(yù)報(bào)效果仍然不是十分理想，特別是對(duì)于一些夏季降水關(guān)鍵的異常年份，例如1997，1999這兩個(gè)降水異常偏少年，其預(yù)報(bào)結(jié)果的 ACC都只有－0.7左右，幾乎是給出了完全相反的預(yù)報(bào)結(jié)果.因此，模式的系統(tǒng)訂正結(jié)果對(duì)于降水偏少年的預(yù)報(bào)效果較差，對(duì)降水偏多的預(yù)報(bào)效果則相對(duì)較好.綜上所述，系統(tǒng)訂正模式預(yù)報(bào)結(jié)果的方法是比較簡(jiǎn)單、粗糙的，對(duì)每個(gè)預(yù)報(bào)年份缺乏一定的針對(duì)性，使得預(yù)報(bào)準(zhǔn)確率仍處于較低的水平;提高預(yù)報(bào)準(zhǔn)確率，需要用具有一定技巧性的統(tǒng)計(jì)方法結(jié)合模式，來對(duì)預(yù)報(bào)年的預(yù)報(bào)誤差進(jìn)行有針對(duì)性的、較準(zhǔn)確的訂正.

3.2.動(dòng)力-相似訂正

短期氣候預(yù)測(cè)中應(yīng)用相似的概念由來已久，這是由于繼天氣過程的相似發(fā)展之后，未來一段時(shí)間內(nèi)也往往繼續(xù)出現(xiàn)相似的天氣過程.早在20世紀(jì)40年代，前蘇聯(lián)就在長(zhǎng)期天氣預(yù)報(bào)方面提出過兩個(gè)自然天氣季節(jié)的相似條件，這些相似條件在理論上是嚴(yán)謹(jǐn)?shù)?，但在?shí)際工作中很難完全滿足［40］.相似預(yù)報(bào)的成效在很大程度上取決于相似條件選取得是否客觀合理，因此相似預(yù)報(bào)的關(guān)鍵問題在于歷史相似的選取.由于數(shù)值模式的初值條件的自由度非常大，要在歷史資料中找到與其合適的相似，在現(xiàn)有資料條件下幾乎是不可能的，但對(duì)某一特定區(qū)域而言，并不是所有前期環(huán)境變量都有同等重要的作用，都存在著前期關(guān)鍵變量，但是要想通過統(tǒng)計(jì)方法從初始時(shí)刻場(chǎng)中確定關(guān)鍵變量也是不太現(xiàn)實(shí)的.由于大氣是不斷地向外強(qiáng)迫適應(yīng)的，所以前期外強(qiáng)迫的異常信號(hào)會(huì)導(dǎo)致數(shù)值模式初邊值的異常.統(tǒng)計(jì)預(yù)報(bào)方法正是通過歷史資料建立起前期外強(qiáng)迫異常與未來異常的聯(lián)系，那么就可用統(tǒng)計(jì)預(yù)報(bào)中前期的環(huán)境變量選取歷史相似，在縮小了自由度的同時(shí)又抓住了關(guān)鍵的歷史相似信息.

3.3.區(qū)域關(guān)鍵因子的選取

汛期降水是由多種影響因子綜合作用的結(jié)果，對(duì)中國各個(gè)區(qū)域都存在各自占主導(dǎo)作用并直接影響區(qū)域夏季降水多寡的前期關(guān)鍵因子.例如，副高是直接影響中國夏季降水的天氣系統(tǒng).研究表明，在副高、海溫、積雪、阻高、季風(fēng)等影響長(zhǎng)江流域夏季降水的因子中，副高與降雨的對(duì)應(yīng)關(guān)系最好，作用最為顯著.姚愚和嚴(yán)華生［41］研究了副高的5種指數(shù)對(duì)中國降雨的影響，發(fā)現(xiàn)副高北界和脊線位置主要影響汛期雨帶南北位置的變化，副高加強(qiáng)西伸時(shí)則有利于中國大部分地區(qū)的降水，強(qiáng)度指數(shù)和面積指數(shù)與局部地區(qū)降水的多寡有關(guān).其他研究也表明，副高的形狀、強(qiáng)弱和位置變動(dòng)與我國大部地區(qū)主汛期旱澇趨勢(shì)有密切聯(lián)系.西太平洋副熱帶高壓和阻塞高壓形勢(shì)、西風(fēng)急流、大陸副熱帶高壓等因子對(duì)于華北夏季降水量的多少具有很好的指示意義，也為選擇關(guān)鍵的前期相似因子提供了可靠的依據(jù).

大量研究表明，前冬海溫對(duì)夏季副高的形狀、位置和強(qiáng)弱等特征有顯著影響，是具有前兆意義的預(yù)測(cè)信號(hào)，可以為相似預(yù)報(bào)方法提供參考，檢驗(yàn)在影響區(qū)域降水的海溫關(guān)鍵區(qū)域是否真的存在相似之處.我們給出了華北夏季降水與前冬平均海溫的相關(guān)，如圖3所示.從圖3中可以看出，華北的夏季降水與前冬印度洋、赤道中東太平洋海溫以及加勒比海域的海溫呈現(xiàn)出顯著的正相關(guān)，即前期赤道印度洋及阿拉伯海域和孟加拉灣的海溫如果處于異常偏冷狀態(tài)，而印度洋南部是偏暖分布時(shí)，則華北地區(qū)夏季降水具有整體上偏少的傾向，反之亦然.冬季印度洋海溫異常對(duì)預(yù)測(cè)未來華北地區(qū)降水具有較明確的指示意義，所以可以用海溫場(chǎng)來檢驗(yàn)利用關(guān)鍵因子選擇的歷史相似年與預(yù)報(bào)年之間在特定的關(guān)鍵區(qū)域是否真的相似.

利用相似預(yù)報(bào)方法對(duì)預(yù)報(bào)年選取相似誤差場(chǎng)時(shí)，從國家氣候中心氣候系統(tǒng)診斷預(yù)測(cè)室提供的74項(xiàng)環(huán)流特征量以及NOAA提供的40項(xiàng)氣候指數(shù)(由于環(huán)流特征量是逐月數(shù)據(jù)，對(duì)于同一個(gè)指數(shù)12個(gè)月的數(shù)值看作12個(gè)獨(dú)立的前期影響因子，因此我們一共可以得到1368個(gè)前期因子)中篩選出一部分對(duì)華北地區(qū)有物理意義并且可以對(duì)模式的預(yù)報(bào)誤差做很好訂正的因子作為關(guān)鍵因子.首先，對(duì)前期的預(yù)報(bào)因子與華北汛期降水量進(jìn)行相關(guān)性檢驗(yàn)，初步提取出對(duì)華北夏季降水影響較大的因子，在相關(guān)系數(shù)都通過95%可信度檢驗(yàn)的環(huán)流特征量中尋找前期關(guān)鍵因子，對(duì)相關(guān)性高的因子還需進(jìn)行進(jìn)一步的篩選，因?yàn)橛械囊蜃与m然與華北地區(qū)的夏季降水有很好的相關(guān)，但模式結(jié)果對(duì)其沒有很好的響應(yīng)，從而不能對(duì)預(yù)報(bào)誤差進(jìn)行估計(jì).因此，還要通過預(yù)報(bào)誤差訂正回報(bào)結(jié)果的評(píng)分，提取出評(píng)分較高的因子，將這些因子作為模式預(yù)報(bào)的前期關(guān)鍵因子，在進(jìn)行相關(guān)性檢驗(yàn)之后，對(duì)通過相關(guān)性檢驗(yàn)的因子進(jìn)行單因子相似誤差訂正交叉檢驗(yàn)回報(bào)評(píng)分，即以每個(gè)因子作為選相似的依據(jù)，將歐氏距離作為選取相似年的條件，計(jì)算預(yù)報(bào)年中因子強(qiáng)度與之前每年強(qiáng)度的歐氏距離選取相似，距離越短則越相似，為了保證相似質(zhì)量，選取最接近的4個(gè)相似年進(jìn)行誤差線性估計(jì)，以此確定預(yù)報(bào)年的相似誤差，通過該方法可篩選出能夠提高模式預(yù)報(bào)結(jié)果準(zhǔn)確率的關(guān)鍵因子.

選擇合適的因子進(jìn)行相似誤差的訂正對(duì)于提高預(yù)報(bào)水平的效果是比較明顯的.我們這里以交叉檢驗(yàn)平均ACC最高的東太平洋濤動(dòng)指數(shù)為例，給出了預(yù)報(bào)年前期8月東太平洋濤動(dòng)指數(shù)進(jìn)行華北汛期降水相似誤差訂正的27年ACC及均方根誤差，如圖4所示.由圖4可得，總體而言單因子(8月東太平洋濤動(dòng)指數(shù))相似訂正的方法較系統(tǒng)訂正方法有一定的改進(jìn)，提高較為明顯的年份是1992，2001—2005年，其中2002年為降水偏少年，2003，2004年為降水偏多年，其他為降水正常年，這幾個(gè)年份中的8月東太平洋濤動(dòng)指數(shù)都是處于偏小的狀態(tài)，通過對(duì)ACC有一定提高的年份以及圖4(c)中8月東太平洋濤動(dòng)指數(shù)的年際變化進(jìn)行分析研究，可以發(fā)現(xiàn)當(dāng)8月東太平洋濤動(dòng)指數(shù)偏小時(shí)，利用該因子相似訂正的方法效果較好，對(duì)于8月東太平洋濤動(dòng)指數(shù)強(qiáng)度較大時(shí)訂正效果則并不明顯，反而會(huì)起到相反的作用.可以看到，由于近幾年該因子變小的趨勢(shì)，所以近幾年的相似訂正效果相對(duì)系統(tǒng)訂正基本都有一定的提高，對(duì)于降水異常年的預(yù)報(bào)結(jié)果仍然不理想，對(duì)關(guān)鍵的降水年份例如1997和1999年兩個(gè)異常年，訂正結(jié)果相對(duì)系統(tǒng)訂正提高非常有限，沒有給出很好的預(yù)報(bào)結(jié)果.此外，由于其他環(huán)境要素的影響以及因子對(duì)于區(qū)域降水影響的年際變化，單個(gè)因子不能保持每年都和區(qū)域的降水有著很好的對(duì)應(yīng)關(guān)系，不能保證預(yù)報(bào)結(jié)果都要優(yōu)于系統(tǒng)訂正的結(jié)果，且從圖4(b)均方根誤差的改進(jìn)也是非常有限的，在部分年份的預(yù)報(bào)結(jié)果還要低于系統(tǒng)訂正方法的結(jié)果.綜上所述，單因子相似誤差訂正也存在著一定的局限性:由于汛期降水是多種因素綜合作用的結(jié)果，只用單因子選擇相似就忽略了其他因子對(duì)于華北地區(qū)的影響，單個(gè)因子不可能完全反映出預(yù)報(bào)年前期環(huán)境變量場(chǎng)的整體信息，因此利用單個(gè)影響因子的前期信號(hào)作為相似年選取依據(jù)有時(shí)候會(huì)比較片面，在因子對(duì)華北地區(qū)降水影響較強(qiáng)的年份選取的相似年與預(yù)報(bào)年的相似度較好，但在該因子作用較弱時(shí)，選取的相似年可能會(huì)出現(xiàn)較大的偏差.

圖4 系統(tǒng)訂正與單因子(8月東太平洋濤動(dòng))相似訂正的近27年變化趨勢(shì) (a)ACC評(píng)分，(b)均方根誤差，(c)8月東太平洋濤動(dòng)指數(shù)

由此可見，單因子選取相似歷史的方法對(duì)預(yù)報(bào)效果的訂正和提高仍是非常有限的，且因子對(duì)于華北地區(qū)的降水影響可能存在著年代際的變化，使得單因子相似訂正不能夠保證對(duì)每年的預(yù)報(bào)效果都有提高，只對(duì)該因子信號(hào)影響較大的年份有著較好的訂正，對(duì)于因子影響較弱的年份則效果并不理想.因此，在利用影響因子的前期信號(hào)選取歷史相似年時(shí)，必須同時(shí)考慮多個(gè)影響因子的作用，將這些因子進(jìn)行組合以后得到的前期綜合作用信號(hào)作為選取相似年的最佳判據(jù).

3.4.最優(yōu)多因子組合相似判據(jù)的配置

對(duì)于多因子組合相似判據(jù)的配置，由于前期環(huán)境場(chǎng)變量的自由度很大，導(dǎo)致計(jì)算量非常大，另外，多因子進(jìn)行組合配置時(shí)的噪聲問題以及因子間可能存在的自相關(guān)會(huì)使得相似訂正誤差出現(xiàn)非線性增長(zhǎng).研究表明，通過正交經(jīng)驗(yàn)函數(shù)(empirical orthogonal function，EOF)分解方法能有效壓縮自由度同時(shí)消除部分噪聲影響，從而降低前期環(huán)境場(chǎng)的自由度.因此，對(duì)多因子組合構(gòu)成的變量場(chǎng) Xm×n進(jìn)行EOF分解，即

其中，Vm×n空間函數(shù)矩陣，Tm×n為時(shí)間系數(shù)矩陣.

通過計(jì)算得到每個(gè)特征向量λk的方差貢獻(xiàn)，即

由于前幾個(gè)主分量的方差貢獻(xiàn)已經(jīng)包含了原變量場(chǎng)的大部分信息，為了進(jìn)一步壓縮相似判據(jù)的自由度，我們?nèi)±鄯e方差貢獻(xiàn)達(dá)到80%，即

此時(shí)，將得到的主分量時(shí)間系數(shù)作為相似年的選取依據(jù)，進(jìn)行交叉檢驗(yàn)以及多因子獨(dú)立樣本檢驗(yàn)的雙重檢驗(yàn)，以此得到最優(yōu)的多因子組合配置.

多因子組合中元素個(gè)數(shù)的變化以及不同因子間的搭配對(duì)預(yù)報(bào)效果也會(huì)有很大的影響.為得到最佳的組合個(gè)數(shù)以及最佳的因子搭配，首先研究組合中因子個(gè)數(shù)的變化對(duì)于ACC的影響，如圖5所示.圖 5 給出了 1983—2006，1983—2007，1983—2008，1983—2009年4個(gè)時(shí)段內(nèi)交叉檢驗(yàn)的平均ACC隨影響因子個(gè)數(shù)的變化趨勢(shì).圖中1983—2006指的是對(duì)1983—2006年時(shí)段進(jìn)行多因子組合交叉檢驗(yàn)，1983—2007以此類推.這里分不同時(shí)段進(jìn)行多因子組合交叉檢驗(yàn)是因?yàn)榭紤]到影響因子的實(shí)測(cè)數(shù)據(jù)是不斷更新的，并且部分影響因子對(duì)于華北地區(qū)汛期降水的影響作用可能存在著年代際的變化，會(huì)導(dǎo)致對(duì)不同時(shí)段進(jìn)行交叉檢驗(yàn)篩選出最佳的因子組合搭配有所差別，因此，有必要考慮多因子組合的配置關(guān)系隨不同時(shí)段的變化，從而可對(duì)不同的預(yù)報(bào)時(shí)期因子影響作用的變化進(jìn)行有針對(duì)性的動(dòng)態(tài)最優(yōu)多因子組合配置.

從圖5中可以看出，4個(gè)時(shí)段的多因子配置交叉檢驗(yàn)ACC評(píng)分隨影響因子個(gè)數(shù)的變化規(guī)律是非常一致的:當(dāng)組合影響因子個(gè)數(shù)從1個(gè)增加至10個(gè)的過程中，ACC的提高非常明顯，當(dāng)增加到10個(gè)以上時(shí)，ACC增幅逐漸減小并趨于穩(wěn)定，當(dāng)因子個(gè)數(shù)進(jìn)一步增大至60個(gè)以上時(shí)，平均ACC都出現(xiàn)了不同程度的減小，但減幅較小.同時(shí)我們也可看到，不同時(shí)段內(nèi)交叉檢驗(yàn)得到最終平均ACC穩(wěn)定值之間差別不大，基本在0.5—0.6之間，說明在現(xiàn)有的資料條件下，通過相似誤差訂正后的華北地區(qū)模式預(yù)報(bào)結(jié)果可以達(dá)到的平均 ACC上限為0.6左右.因此，采用多因子組合配置是提高汛期降水預(yù)測(cè)水平的一個(gè)有效途徑.通過上述變化我們也可得到，因子組合個(gè)數(shù)在增加至15個(gè)以上時(shí)對(duì)于ACC的提高已經(jīng)很小，此時(shí)ACC已經(jīng)基本達(dá)到了上限值，同時(shí)考慮到計(jì)算量的問題，因此在最優(yōu)因子組合篩選的時(shí)候因子個(gè)數(shù)在15個(gè)左右是最優(yōu)化的.表1分別給出了 1983—2004，1983—2005，1983—2006，1983—2007，1983—2008年5個(gè)時(shí)段配置出的 15個(gè)因子的動(dòng)態(tài)最優(yōu)組合.考慮到業(yè)務(wù)應(yīng)用以及影響因子更新等問題，配置的最優(yōu)組合中1月份因子視為預(yù)報(bào)年當(dāng)年前期因子，而2—12月份的因子則采用預(yù)報(bào)年前年因子.

表1 不同時(shí)段配置出的動(dòng)態(tài)最優(yōu)多因子組合

從表1中可看出，在最優(yōu)因子組合配置的因子中，表示極渦和副高的因子占了絕大部分，這是因?yàn)楦备吲c極渦是影響我國天氣氣候變化的兩個(gè)主要的大氣環(huán)流實(shí)體，并且二者對(duì)氣候的影響并不是孤立的，而是相互作用、相互配合.極渦主要體現(xiàn)的是中高緯環(huán)流特征，而副高很大程度上反映了中低緯大氣環(huán)流的變化.東亞大氣環(huán)流的年際和年代際變化對(duì)華北氣候異常具有重要的影響.氣候系統(tǒng)中大氣環(huán)流因子對(duì)大氣中水分的輸送起到了非常重要的作用，大氣中的水汽含量不僅與大氣環(huán)流關(guān)系密切，而且作為能量和水循環(huán)過程的主要因子對(duì)區(qū)域水分平衡有很大的影響［42］.從前面的相關(guān)計(jì)算可知，副高脊線、極渦指數(shù)與華北汛期降水之間以正相關(guān)為主，而副高面積、強(qiáng)度與華北降水基本呈負(fù)相關(guān).當(dāng)亞洲和歐洲區(qū)極渦異常南擴(kuò)，北非、大西洋、北美副高顯著收縮減弱，西太平洋和南海副高明顯北抬時(shí)，華北降水易偏多.除此之外，可以看到在動(dòng)態(tài)最優(yōu)多因子組合中出現(xiàn)較為頻繁、起到主導(dǎo)作用的影響因子還有太平洋暖池、加勒比海的海表溫度(sea surface temperature，SST)等.太平洋暖池是全球空氣對(duì)流最強(qiáng)烈的地區(qū)，且活動(dòng)持久，是氣候異常的源地之一，太平洋暖池的變化制約著亞洲、太平洋區(qū)域，甚至全球氣候變化和某些重大自然災(zāi)害的形成與變化.而對(duì)于加勒比海 SST，在前面圖3我們已經(jīng)得到加勒比海SST與華北地區(qū)的汛期降水具有非常好的正相關(guān)，該海域的海溫對(duì)于華北地區(qū)的汛期降水具有很好的指示意義.

圖6給出了1983—2009年27年的模式誤差組合a與組合e的EOF分解第一主分量時(shí)間系數(shù)與之間的相關(guān)關(guān)系，對(duì)比圖6(a)和(b)我們可以發(fā)現(xiàn)，組合a的EOF第一主分量的時(shí)間系數(shù)在前面大部分年份里與模式誤差保持著較好的正相關(guān)關(guān)系，而從2004年往后則表現(xiàn)為負(fù)相關(guān)關(guān)系，而組合e EOF第一主分量時(shí)間系數(shù)在前半段時(shí)期內(nèi)與模式誤差的相關(guān)關(guān)系并不非常明顯，但在后期特別是近10年里與模式誤差之間表現(xiàn)出了較強(qiáng)的正相關(guān)關(guān)系，正是由于兩個(gè)組合對(duì)于模式誤差的相關(guān)關(guān)系的年代際變化使得兩個(gè)多因子組合對(duì)于不同時(shí)段的相似預(yù)報(bào)能力也不一致，因此，隨著時(shí)段的推移以及觀測(cè)資料的不斷更新，最優(yōu)多因子組合配置也是動(dòng)態(tài)的、不斷更新的，應(yīng)針對(duì)各時(shí)段內(nèi)具體情況配置最優(yōu)多因子組合.

綜合表1中的多因子組合配置以及圖6中多因子組合EOF主分量時(shí)間系數(shù)與模式誤差的相關(guān)關(guān)系可知，影響因子對(duì)于華北地區(qū)汛期降水的影響作用存在著年際或年代際的變化，導(dǎo)致對(duì)不同時(shí)段內(nèi)進(jìn)行交叉檢驗(yàn)配置出的多因子組合有所差別，多因子組合是動(dòng)態(tài)變化的.由于我們的最終目的是利用多因子組合相似判據(jù)對(duì)未來的汛期降水情況進(jìn)行較準(zhǔn)確的訂正預(yù)報(bào)，有必要研究多因子組合的配置關(guān)系隨時(shí)間的變化對(duì)比其實(shí)際預(yù)報(bào)效果的變化，從而找到適合近階段內(nèi)汛期降水的最優(yōu)多因子組合配置.因此，為了對(duì)比檢驗(yàn)多因子組合的真正預(yù)報(bào)能力，對(duì)不同時(shí)段配置出的最優(yōu)多因子組合進(jìn)行2005—2009年5年的獨(dú)立樣本回報(bào)檢驗(yàn)，如圖7所示.圖7(a)，(b)分別給出了5個(gè)不同時(shí)段交叉檢驗(yàn)配置的最優(yōu)多因子組合的5年回報(bào)的ACC和均方根誤差結(jié)果.圖中組合 a至組合 e分別代表了利用1983—2004年，1983—2005年，1983—2006年，1983—2007年，1983—2008年5個(gè)時(shí)段配置的最優(yōu)多因子組合.

由于獨(dú)立樣本回報(bào)不能用到預(yù)報(bào)年當(dāng)年的降水?dāng)?shù)據(jù)，因此1983—2004時(shí)段配置的最優(yōu)多因子組合，即組合a可以對(duì)2005年以后的年份進(jìn)行獨(dú)立樣本回報(bào)，而1983—2005時(shí)段資料配置的最優(yōu)組合則只能對(duì)2006年以后的年份進(jìn)行回報(bào)，以此類推.從圖7中可以看出，利用1983—2004年資料配置出的組合a可以對(duì)2005年汛期降水進(jìn)行很好的預(yù)測(cè)，相較系統(tǒng)訂正提高非常明顯，而對(duì)于2006的汛期降水組合a的相似訂正預(yù)報(bào)結(jié)果則不太理想，訂正結(jié)果要低于系統(tǒng)訂正，但1983—2005時(shí)段配置的組合b的相似訂正結(jié)果則對(duì)2006年的汛期降水則顯示出良好的預(yù)報(bào)能力，預(yù)報(bào)效果要明顯高于組合a與系統(tǒng)訂正，類似地，對(duì)于2007年汛期獨(dú)立樣本回報(bào)，組合c的相似訂正要好于組合a與組合b.因此，通過對(duì)比各個(gè)組合對(duì)每個(gè)預(yù)報(bào)年的獨(dú)立樣本回報(bào)結(jié)果分析可以看出，交叉檢驗(yàn)時(shí)段越長(zhǎng)配置出來的最優(yōu)多因子組合對(duì)于該預(yù)報(bào)年汛期的預(yù)報(bào)效果相對(duì)較好，一方面因?yàn)闀r(shí)段越長(zhǎng)對(duì)應(yīng)的相似年的樣本量也就越多，另一方面，汛期降水有2—4年周期變化，時(shí)段越長(zhǎng)配置的因子所包含的預(yù)報(bào)年的歷史關(guān)鍵相似信息也越多，配置出來的最優(yōu)多因子組合也就越有針對(duì)性，能夠?qū)υ擃A(yù)報(bào)年的模式誤差進(jìn)行較好的訂正.由于樣本量較少使得以上結(jié)論仍然存在一定的不確定性，但從近5年回報(bào)結(jié)果可以看出，該方法是現(xiàn)階段樣本量較少的情況下提高汛期預(yù)測(cè)能力較為有效的方法.因此，對(duì)不同預(yù)報(bào)年采用不同的最優(yōu)多因子組合即動(dòng)態(tài)最優(yōu)多因子組合可以有效地提高預(yù)報(bào)結(jié)果.

綜上所述，不同組合在不同時(shí)間段內(nèi)顯示出的預(yù)報(bào)能力也是不一樣的，并且一般來說，交叉檢驗(yàn)所選取的時(shí)段越長(zhǎng)其相應(yīng)配置出的多因子組合對(duì)預(yù)測(cè)年的效果越好.因此，隨著資料的不斷更新，在不同年份進(jìn)行實(shí)際的汛期降水預(yù)報(bào)時(shí)，由于可使用的樣本量長(zhǎng)度不同以及前期因子對(duì)降水影響的年際或年代際變化，配置出最優(yōu)多因子組合也應(yīng)當(dāng)是動(dòng)態(tài)變化的.

3.5.預(yù)報(bào)個(gè)例分析

為了較直觀地了解最優(yōu)因子組合的相似選取方法對(duì)于夏季降水預(yù)報(bào)的改進(jìn)能力，選擇了2005—2009年夏季降水作為獨(dú)立樣本進(jìn)行研究，限于篇幅，以2007年為典型預(yù)報(bào)個(gè)例，給出了實(shí)況與預(yù)報(bào)圖，其他年份只給出了相應(yīng)的評(píng)分結(jié)果.之所以選擇2007年，一方面為了采用盡可能多的資料選取較為準(zhǔn)確的歷史相似，另一方面，華北的北部地區(qū)在2007年出現(xiàn)較強(qiáng)的負(fù)異常，同時(shí)在華北南部靠近淮河流域的地區(qū)則是呈現(xiàn)很強(qiáng)的降水正異常，以此可以檢驗(yàn)該動(dòng)力-相似方法對(duì)降水異常的區(qū)域是否有異常信號(hào)的預(yù)報(bào)能力，能否針對(duì)其給出降水異常的預(yù)測(cè).

圖8(a)—(c)分別給出了2007年汛期降水的實(shí)況、系統(tǒng)訂正和相似訂正的預(yù)報(bào)結(jié)果.總體來說，2007年汛期西太平洋副熱帶高壓持續(xù)偏強(qiáng)，亞洲地區(qū)西風(fēng)帶以緯向環(huán)流為主，北半球極渦偏西，7，8月份極渦龜縮于極地，導(dǎo)致入侵我國的冷空氣勢(shì)力偏弱，主要雨帶在黃河和長(zhǎng)江之間，而其中淮河流域汛期出現(xiàn)了僅次于1954年的大洪水.從2007年華北汛期降水實(shí)況圖可以看出，華北2007年降水量分布呈現(xiàn)較明顯的緯向分布，北部地區(qū)整體呈現(xiàn)降水偏少的態(tài)勢(shì)，其中北京、河北中部以及山西北部有較強(qiáng)的降水負(fù)異常，西太平洋副高持續(xù)偏強(qiáng)以及北方冷空氣勢(shì)力偏弱是造成華北北部地區(qū)汛期降水量偏少的重要原因之一;而華北南部地區(qū)則呈現(xiàn)不同程度的降水偏多，其中山東南部地區(qū)降水偏多達(dá)到5成左右.從圖8(b)中可以看到系統(tǒng)誤差訂正給出的結(jié)果是整個(gè)華北都是降水偏多，對(duì)于南部地區(qū)降水預(yù)報(bào)較為準(zhǔn)確，而對(duì)整個(gè)華北北部地區(qū)的預(yù)報(bào)效果并不理想，沒有預(yù)報(bào)出北部地區(qū)的降水偏少，系統(tǒng)誤差訂正ACC為0.54.圖8(c)的相似誤差訂正結(jié)果則相對(duì)系統(tǒng)訂正結(jié)果有了一定的提高，雖然相似誤差訂正也沒有很好地給出華北北部的降水偏少，但對(duì)華北南部地區(qū)的降水偏多則給出了較準(zhǔn)確的預(yù)報(bào)，并且給出了非常強(qiáng)的降水正異常預(yù)報(bào)，強(qiáng)度與實(shí)況大致相當(dāng)，范圍略小于實(shí)況，相似誤差訂正ACC為0.69，較系統(tǒng)誤差訂正 ACC提高了0.15，相似訂正方法具有很好的可行性.由此也可看出，多因子優(yōu)化組合的相似選取方法對(duì)前期的異常信號(hào)會(huì)有很好的響應(yīng)并且能夠?qū)⒃摦惓Ｐ盘?hào)反映到異常降水區(qū)域的預(yù)報(bào)中.

表2給出了各個(gè)預(yù)報(bào)年進(jìn)行相似訂正時(shí)所選取的相似年以及預(yù)報(bào)年的實(shí)況與相似年實(shí)況之間的距平相關(guān)系數(shù)，通過計(jì)算距平相關(guān)系數(shù)可以判斷預(yù)報(bào)年與相似年之間的相似度的大小，進(jìn)行判斷選取的相似年的準(zhǔn)確性.結(jié)果顯示，各預(yù)報(bào)年的實(shí)際降水與所選取相似年降水之間具有較好的相似性，只有2006和2009年選取的4個(gè)相似年中有一個(gè)相似年選的不夠準(zhǔn)確，差異較大，距平相關(guān)系數(shù)達(dá)到了負(fù)值，對(duì)最終的預(yù)報(bào)結(jié)果造成了一定的影響，使得訂正效果有一定下降，但由于預(yù)報(bào)年選取的其他相似年都較為準(zhǔn)確，對(duì)最終的結(jié)果影響并不大，預(yù)報(bào)年的實(shí)況與選取的相似年實(shí)況相似度較高，最終的相似訂正結(jié)果也較為理想，預(yù)報(bào)結(jié)果與實(shí)況比較接近，基本預(yù)報(bào)出了預(yù)報(bào)年的降水主要分布形式.

圖8 2007年獨(dú)立樣本回報(bào) 色標(biāo)為距平百分率P (a)實(shí)況，(b)系統(tǒng)訂正，(c)相似訂正

表2 獨(dú)立樣本回報(bào)2005—2009年選取的相似年及預(yù)報(bào)年與相似年的距平相關(guān)系數(shù)

為了直觀地了解動(dòng)態(tài)最優(yōu)多因子組合相似誤差訂正對(duì)于模式誤差的訂正效果，圖9(a)，(b)分別給出了獨(dú)立樣本回報(bào)的2005—2009年汛期降水系統(tǒng)誤差訂正和相似誤差訂正的距平相關(guān)系數(shù)以及均方根誤差對(duì)比.從圖中可以看出，相似誤差訂正的結(jié)果都要優(yōu)于系統(tǒng)誤差訂正結(jié)果，將系統(tǒng)訂正較差的年份的預(yù)報(bào)ACC評(píng)分提高到較高的水平，近5年獨(dú)立樣本的系統(tǒng)訂正平均ACC為0.40，而相似訂正的平均ACC約為0.61，相比提高了0.21，并且始終保持較高的準(zhǔn)確率.對(duì)于圖9(b)中的均方根誤差而言，相似訂正效果較系統(tǒng)訂正也有一定的提高，但并不十分明顯，對(duì)于系統(tǒng)訂正的均方根誤差較大的2005和2007年都有明顯的減小，顯示出明顯的優(yōu)越性.因此，相似訂正結(jié)果能夠穩(wěn)定地提高模式預(yù)報(bào)的準(zhǔn)確率，由此可見，基于動(dòng)態(tài)最優(yōu)多因子集組合的相似誤差訂正方法能有效提高模式的預(yù)報(bào)技巧.

4.動(dòng)態(tài)最優(yōu)多因子組合相似預(yù)報(bào)方案

通過前面的研究可知，多因子組合相似預(yù)報(bào)方法可以有效地改善華北地區(qū)汛期預(yù)報(bào)效果，提高預(yù)報(bào)水平，具有很好的應(yīng)用前景.為了與業(yè)務(wù)應(yīng)用相結(jié)合，給出了最終的動(dòng)態(tài)最優(yōu)多因子組合的相似預(yù)報(bào)方案，如圖10所示.從圖10中可以得到該預(yù)報(bào)方案的框架流程:首先，利用汛期降水 CMAP資料和國家氣候中心季節(jié)預(yù)報(bào)模式對(duì)每年汛期降水預(yù)報(bào)結(jié)果得到汛期降水模式預(yù)報(bào)誤差集，將114項(xiàng)指數(shù)作為影響因子，對(duì)所有影響因子進(jìn)行單因子交叉檢驗(yàn)預(yù)報(bào)實(shí)驗(yàn)，通過交叉檢驗(yàn)的ACC排序篩選出對(duì)華北地區(qū)模式誤差訂正較好的前期關(guān)鍵因子，建立關(guān)鍵因子集;針對(duì)關(guān)鍵因子集中的各項(xiàng)因子，對(duì)其進(jìn)行組合配置實(shí)驗(yàn)，通過交叉檢驗(yàn)ACC敏感性實(shí)驗(yàn)確定最優(yōu)的多因子組合成員個(gè)數(shù)，通過得到區(qū)域預(yù)報(bào)年前期最優(yōu)多因子組合，并針對(duì)影響因子對(duì)區(qū)域汛期降水的年際或年代際變化導(dǎo)致不同時(shí)段內(nèi)最優(yōu)多因子配置的變化，結(jié)合歷史近期最優(yōu)多因子組合得到適合不同預(yù)報(bào)時(shí)段的動(dòng)態(tài)最優(yōu)多因子配置，最后進(jìn)行模式誤差預(yù)報(bào)并對(duì)模式預(yù)報(bào)結(jié)果進(jìn)行訂正并給出訂正后模式預(yù)報(bào)結(jié)果.

5.結(jié) 論

基于動(dòng)力-相似的基本原理，利用已有的物理統(tǒng)計(jì)氣候?qū)W方面的研究成果，通過對(duì)前期因子進(jìn)行單因子交叉檢驗(yàn)篩選，建立適用于華北夏季降水的區(qū)域特點(diǎn)的普適性固定關(guān)鍵相似因子集，對(duì)關(guān)鍵因子集中的影響因子進(jìn)行單因子交叉檢驗(yàn)回報(bào)實(shí)驗(yàn)，給出單因子交叉檢驗(yàn)ACC排序;針對(duì)單因子交叉檢驗(yàn)ACC排序，篩選出能反映模式預(yù)報(bào)誤差分布特征的關(guān)鍵預(yù)報(bào)因子;對(duì)關(guān)鍵預(yù)報(bào)因子進(jìn)行配置試驗(yàn)，通過獨(dú)立樣本檢驗(yàn)ACC得到區(qū)域最優(yōu)多因子配置，考慮到影響因子對(duì)區(qū)域汛期降水的年際或年代際變化導(dǎo)致在不同的研究時(shí)段內(nèi)的最優(yōu)多因子配置的動(dòng)態(tài)變化，結(jié)合歷史近期最優(yōu)多因子組合得到適合不同預(yù)報(bào)時(shí)段的動(dòng)態(tài)最優(yōu)多因子配置，結(jié)合歷史近期最優(yōu)因子配置得到預(yù)報(bào)時(shí)段內(nèi)穩(wěn)定的最優(yōu)關(guān)鍵因子組合;針對(duì)預(yù)報(bào)年，通過最優(yōu)多因子組合分別選取最優(yōu)的4個(gè)相似年估計(jì)當(dāng)前的預(yù)報(bào)誤差，利用計(jì)算出的預(yù)報(bào)誤差訂正模式的預(yù)報(bào)結(jié)果，改善華北地區(qū)夏季降水預(yù)報(bào)效果.經(jīng)獨(dú)立樣本回報(bào)檢驗(yàn)顯示，動(dòng)態(tài)最優(yōu)關(guān)鍵因子組合的相似誤差訂正方法可以顯著提高華北地區(qū)汛期降水的預(yù)報(bào)技巧，對(duì)華北地區(qū)夏季降水具有很好的預(yù)報(bào)可行性，并最終形成了動(dòng)態(tài)最優(yōu)多因子組合的相似預(yù)報(bào)方案.

［1］Feng G L，Gong Z Q，Zhi R，Zhang D Q 2008Chin.Phys.B 17 2745

［2］Feng G L，Dong W J，Li J P 2004Chin.Phys.13 1582

［3］Zhao S R，Song Z S 1999Plateau Meteorol.18 535(in Chinese)［趙聲蓉、宋正山1999高原氣象 18 535］

［4］Hou W，Yang P，F(xiàn)eng G L 2008Acta Phys.Sin.57 3932(in Chinese)［侯 威、楊 萍、封國林 2008物理學(xué)報(bào) 57 3932］

［5］Barnett T P 1990Nature343 743

［6］Feng G L，Dong W J，Jia X J，Cao H X 2002Acta Phys.Sin.51 1181(in Chinese)［封國林、董文杰、賈曉靜、曹鴻興 2002物理學(xué)報(bào) 51 1181］

［7］Hou W，Yang P，Zheng Z H，Gong Z Q 2009Acta Phys.Sin.58 2106(in Chinese)［侯 威、楊 萍、鄭志海、龔志強(qiáng) 2009物理學(xué)報(bào) 58 2106］

［8］Gong Z Q，F(xiàn)eng G L，Wan S Q，Li J P 2006Acta Phys.Sin.55 477(in Chinese)［龔志強(qiáng)、封國林、萬仕全、李建平 2006物理學(xué)報(bào) 55 477］

［9］Zhang D Q，Zhang L，Yang J，F(xiàn)eng G L 2010Acta Phys.Sin.59 655(in Chinese)［章大全、張 璐、楊 杰、封國林 2010物理學(xué)報(bào) 59 655］

［10］Yang J，Hou W，F(xiàn)eng G L 2010Acta Phys.Sin.59 664(in Chinese)［楊 杰、侯 威、封國林 2010物理學(xué)報(bào)59 664］

［11］Zhi R，Lian Y，F(xiàn)eng G L 2007Acta Phys.Sin.56 1837(in Chinese)［支 蓉、廉 毅、封國林2007物理學(xué)報(bào) 56 1837］

［12］Xie K，Ren X J 2008Sci.Meteorol.Sin.28 508(in Chinese)［謝 坤、任雪娟2008氣象科學(xué) 28 508］

［13］Zhang L J，Xia J，Hu Z F 2008Adv.Clim.Chang.Res.4 140(in Chinese)［張利平、夏 軍、胡志芳 2008氣候變化研究 4 140］

［14］Lu R Y 2005Chin.Sci.Bull.50 1131(in Chinese)［陸日宇2005科學(xué)通報(bào) 50 1131］

［15］Gong Z Q，F(xiàn)eng G L 2008Acta Phys.Sin.57 3920(in Chinese)［龔志強(qiáng)、封國林2008物理學(xué)報(bào) 57 3920］

［16］Fan K，Lin M J，Gao Y Z 2008Sci.Chin.38 1452(in Chinese)［范 可、林美靜、高煜中 2008中國科學(xué)38 1452］

［17］Feng G L，Dai X G，Wang A H，Chou J F 2001Acta Phys.Sin.50 606(in Chinese)［封國林、戴新剛、王愛慧、丑紀(jì)范2001物理學(xué)報(bào)50 606］

［18］Chou J F 2003Bimonth.Xinjiang Meteorol.26 1(in Chinese)［丑紀(jì)范2003新疆氣象 26 1］

［19］Chou J F 1986Plateau Meteorol.5 367(in Chinese)［丑紀(jì)范1986高原氣象5 367］

［20］Thomas A G 1970J.Appl.Meteorol.8 333

［21］Michael K T，Lisa G，Anthony G B 2005Notes Corresp.18 1831

［22］Gu Z C 1958Acta Meteorol.Sin.29 176(in Chinese)［顧震潮1958氣象學(xué)報(bào)29 176］

［23］Chou J F 1974Sci.Chin.Ser.A 6 635(in Chinese)［丑紀(jì)范1974中國科學(xué)A輯 6 635］

［24］Qiu C J，Chou J F 1989Chin.J.Atmos.Sci.13 22(in Chinese)［邱崇踐、丑紀(jì)范1989大氣科學(xué)13 22］

［25］Huang J P，Wang S W 1991Sci.Chin.Ser.B 2 216(in Chinese)［黃建平、王紹武 1991中國科學(xué)B輯2 216］

［26］Huang J P，Yi Y H，Wang S W，Chou J F 1993Quart J.Roy.Meteor.Soc.119 547

［27］Ren H L，Zhang P Q，Li W J，Chou J F 2006Acta Phys.Sin.55 4388(in Chinese)［任宏利、張培群、李維京、丑紀(jì)范 2006物理學(xué)報(bào)55 4388］

［28］Chen H，Lin ZH 2006Adv.Atmos.Sci.23 425

［29］Gu X Q 1998Chin.Sci.Bull.43 1(in Chinese)［谷湘潛1998科學(xué)通報(bào)43 1］

［30］Cao H X 1993Sci.Chin.Ser.B 23 104(in Chinese)［曹鴻興1993中國科學(xué)B輯23 104］

［31］Feng G L，Cao H X，Wei F Y，Chou J F 2001Acta Meteorol.Sin.59 206(in Chinese)［封國林、曹鴻興、魏鳳英、丑紀(jì)范2001氣象學(xué)報(bào) 59 206］

［32］Ding Y H，Liu Y M，Song Y J，Li Q Q 2002Clim.Envir.7 236(in Chinese)［丁一匯、劉一鳴、宋永加、李清泉 2002氣候與環(huán)境研究7 236］

［33］Li W J，Zhang P Q，Li Q Q，Wang L N，Liu Y M，Shi X L，Zhang Z Q，Liu Y M，Hu G Q，Dang H Y，Zhang F，Chen L J，Sun C R，Zhao Q G，Dong M 2005Quart.J.Appl.Meteorol.16 1(in Chinese)［李維京、張培群、李清泉、王蘭寧、劉益民、史學(xué)麗、張祖強(qiáng)、劉一鳴、胡國權(quán)、黨鴻雁、張 芳、陳麗娟、孫除榮、趙其慶、董 敏 2005應(yīng)用氣象學(xué)報(bào)16 1］

［34］Ren H L，Chou J F 2007Sci.Chin.Ser.D 37 988(in Chinese)［任宏利、丑紀(jì)范2007中國科學(xué)D輯37 988］

［35］Ren H L，Chou J F 2007Adv.Earth Sci.22 376(in Chinese)［任宏利、丑紀(jì)范2007地球科學(xué)進(jìn)展22 376］

［36］Feng G L，Dong W J 2003Acta Phys.Sin.52 2347(in Chinese)［封國林、董文杰 2003物理學(xué)報(bào) 52 2347］

［37］Zheng Z H，Ren H L，Huang J P 2009Acta Phys.Sin.58 7359(in Chinese)［鄭志海、任宏利、黃建平 2009物理學(xué)報(bào) 58 7359］

［38］Bao M，Ni Y Q，Chou J F 2004Chin.Sci.Bull.49 1112(in Chinese)［鮑 名、倪允琪、丑紀(jì)范 2004科學(xué)通報(bào)49 1112］

［39］Shi Y S，Yao X X，Yang X L，Li Z T 2008Sci.Meteorol.Sin.28 377(in Chinese)［史印山、姚學(xué)祥、楊曉亮、李宗濤 2008氣象科學(xué)28 377］

［40］Mao Y J 2005Bull.Sci.Tech.5 533(in Chinese)［毛燕軍2005科技通報(bào) 5 533］

［41］Yao Y，Yan H S 2008J.Trop.Meteorol.24 483(in Chinese)［姚 愚、嚴(yán)華生2008熱帶氣象學(xué)報(bào)24 483］

［42］Zhang H D，Jin R H，Zhang Y S 2008J.Trop.Meteorol.24 417(in Chinese)［張恒德、金榮花、張友姝 2008熱帶氣象學(xué)報(bào)24 417］

PACS:92.60.Wc

Dynamic optimal multi-indexes configuration for estimating the prediction errors of dynamical climate model in North China*

Yang Jie1)2)Wang Qi-Guang1)Zhi Rong3)Feng Guo-Lin2)3)?
1)(College of Atmospheric Scicences，Lanzhou University，Lanzhou 730000，China)
2)(Key Laboratory of Regional Climate-Environment Research for Temperate East Asia of the Institute of Atmospheric Physics，
Chinese Academy of Sciences，Beijing 100029，China)
3)(Laboratory for Climate Studies of China Meteorological Administration，National Climate Center，Beijing 100081，China)
(Received 22 June 2010;revised manuscript received 15 July 2010)

Based on the idea of using the historical-analogue information to revise the prediction errors of National Climate Centre numerical business model，for North China，based on analysis data of the CMAP from 1983 to 2009，40 pieces of climate indices from NOAA，27 years of the season prediction model results from 1983 to 2009 and 74 pieces of circulation characteristics materials provided by Weather Diagnostic Forecasting Room of National Climate Center，using the method of combining data analysis and numerical simulation of diagnostic tests，taking the advantage of the prediction error of the key information of similar model from the historical data，by identifying key factors，optimizing allocation of the different factors of different forecasting years，we established specific multi-factor dynamic optimal portfolios to revise prediction errors in different periods of the power-statistical model in North China，and constructed early environmental factors similar to field multiple objective criteria，to develop new technology of revising prediction errors from the power-statistical model based on dynamic optimal combination of multi-factor，and improved the prediction effect in the summer precipitation in North China and the forecasting skills.Results of independent sample return of 2005—2009 shows that，the score of similarity revised method has improved significantly compared with the score of systematic revised method.The method has a good prospect for summer precipitation forecast in North China，and is going to be put into operation.

set of key factors，estimation of model errors，precipitation，prediction

*公益性行業(yè)(氣象)科研專項(xiàng)基金(批準(zhǔn)號(hào):GYHY200806005)、國家自然科學(xué)基金(批準(zhǔn)號(hào):40875040，40930952)和國家科技支撐計(jì)劃(批準(zhǔn)號(hào):2007BAC29B01，2009BAC51B04)資助的課題.

?通訊聯(lián)系人.E-mail:fenggl@cma.gov.cn

*Project supported by the Special Scientific Research Fund of Meteorological Public Welfare Profession of China(Grant No.GYHY200806005)，the National Natural Science Foundation of China(Grant Nos.40875040，40930952)and the National Key Technology Research and Development Program of the Ministry of Science and Technology of China(Grant Nos.2007BAC29B01，2009BAC51B04).

?Corresponding author.E-mail:ferggl@cma.gov.cn