條件非線性最優(yōu)擾動在可預(yù)報(bào)性問題研究中的應(yīng)用

穆穆 段晚鎖

1中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室，北京100029

2中國科學(xué)院海洋研究所，青島266071

1 引言

可預(yù)報(bào)性是數(shù)值天氣預(yù)報(bào)和氣候預(yù)測研究的基本問題之一。在可預(yù)報(bào)性問題研究中，關(guān)于預(yù)報(bào)結(jié)果不確定性的估計(jì)是可預(yù)報(bào)性研究的重要問題之一。Tennekes（1991）指出，對于一個預(yù)報(bào)產(chǎn)品，如果沒有估計(jì)其可能的預(yù)報(bào)誤差，那么這個預(yù)報(bào)產(chǎn)品是不能夠被使用的。事實(shí)上，早在 1957年，Thompson（1957）已經(jīng)意識到了估計(jì)預(yù)報(bào)結(jié)果不確定性的重要性，業(yè)務(wù)天氣預(yù)報(bào)也從此開始嘗試在每一次預(yù)報(bào)中定量估計(jì)初始場不確定性的發(fā)展。這種觀點(diǎn)后來也被海氣耦合系統(tǒng)的預(yù)報(bào)所采納（Moore and Kleeman, 1996; Samelson and Tziperman, 2002;Mu et al., 2007a; Duan et al., 2009; Yu et al., 2012a,2012b）。通過對預(yù)報(bào)誤差演變規(guī)律的研究以及對預(yù)報(bào)誤差的估計(jì)，人們對熱帶海氣耦合系統(tǒng)的可預(yù)報(bào)性，獲得了更為清晰的認(rèn)識和理解。

可預(yù)報(bào)性研究的另一個關(guān)鍵問題是探尋所關(guān)注的天氣和氣候事件發(fā)生的前期征兆，即何種初始異常會發(fā)展成所關(guān)注的天氣或氣候事件（Moore and Kleeman, 1996; Duan et al., 2004; Duan et al., 2012）。這也是動力氣象學(xué)與氣候動力學(xué)中的中心問題之一。在一定條件下，最容易發(fā)展成某一天氣或氣候事件的初始異常，我們稱其為該事件的最優(yōu)前期征兆（Duan et al., 2004）。如何通過觀測（或同化）捕捉到初始場（分析場）中的這種前期征兆信號是天氣和氣候預(yù)報(bào)中的關(guān)鍵問題。國際上有研究使用線性奇異向量（LSV）方法研究天氣和氣候事件的最優(yōu)前兆信號。如，Moore and Kleeman（1996）將LSV初始距平描述的海溫狀態(tài)作為El Ni?o事件發(fā)生的前期征兆。

LSV是線性模式中增長最快的一類初始擾動，該方法最早由 Lorenz（1965）將其引入氣象學(xué)研究。它不僅被用于可預(yù)報(bào)性最優(yōu)前期征兆問題的研究，同時也被廣泛地用于研究預(yù)報(bào)誤差增長的可預(yù)報(bào)性問題。LSV基于線性理論，其使用的前提條件是初始擾動充分小，其發(fā)展能夠近似地由切線性模式（TLM）來描述。然而，大氣和海洋的運(yùn)動是非線性的，天氣和氣候的可預(yù)報(bào)性不可避免地會受到非線性物理過程的影響。所以，基于線性模式的LSV不能充分反映非線性物理過程的影響，不能用于描述非線性模式中增長最快的有限振幅的初始擾動（Mu et al., 2003）。在可預(yù)報(bào)性研究中，LSV不能有效代表天氣和氣候事件的最優(yōu)前期征兆和增長最快初始誤差（Duan et al., 2004; Mu et al.,2007a, 2007b）。

為了考察非線性的影響, Mu et al.(2003) 提出了條件非線性最優(yōu)擾動（CNOP）方法。CNOP代表了具有一定物理約束條件且在預(yù)報(bào)時刻具有最大非線性發(fā)展的一類初始擾動；LSV代表了切線性模式中具有最大增長率的一類初始擾動，CNOP是LSV在非線性系統(tǒng)中的自然推廣。CNOP在預(yù)報(bào)時刻的增長率可能沒有LSV大，但其在預(yù)報(bào)時刻的非線性發(fā)展是最大的，從而代表了在預(yù)報(bào)時刻對預(yù)報(bào)結(jié)果有最大影響的一類初始誤差，也代表了最容易發(fā)展成一類天氣或氣候事件的前期征兆。因此，在可預(yù)報(bào)性研究中，CNOP比LSV具有更重要的作用。尤其在非線性系統(tǒng)中，CNOP可以刻畫有限振幅擾動大小的發(fā)展，LSV只能描述充分小初始擾動的演變，而這正是非線性影響的結(jié)果。所以，通過比較CNOP和LSV，可以揭示非線性物理過程對可預(yù)報(bào)性的影響。

CNOP方法首先被應(yīng)用于ENSO可預(yù)報(bào)性的研究，用來尋找ENSO事件的最優(yōu)前期征兆和探討預(yù)報(bào)結(jié)果不確定性問題（尤其是春季可預(yù)報(bào)性障礙(SPB) 問題）（Duan et al., 2004; Mu et al., 2007a,2007b; Duan et al., 2009; Duan et al., 2012; Yu et al.,2009; Yu et al., 2012a, 2012b）。繼而又將CNOP方法應(yīng)用于熱鹽環(huán)流（THC）穩(wěn)定性和敏感性問題的研究中，探討非線性過程對 THC穩(wěn)定性和敏感性的影響（Mu et al., 2004; Wu and Mu, 2009）。所有這些工作表明，在考慮有限振幅初始擾動或/和較長預(yù)報(bào)時間的情況下，CNOP和LSV存在明顯區(qū)別。對于非線性大氣和海洋系統(tǒng)，我們相信CNOP的結(jié)果更能刻畫大氣和海洋的非線性運(yùn)動特征，因此在估計(jì)天氣和氣候的預(yù)報(bào)結(jié)果不確定性的研究中，我們應(yīng)該使用CNOP方法。Riviere et al.（2008）也將CNOP方法用于一個二層準(zhǔn)地轉(zhuǎn)模式研究非線性物理過程對斜壓不穩(wěn)定流的影響。Terwisscha van Scheltinga （2007）通過一種隱式4D-Var方法計(jì)算了描述Double Gyre海洋環(huán)流模式的CNOP，研究了有限振幅double gyre環(huán)流的穩(wěn)定性，揭示了非線性物理過程的影響。CNOP方法還被應(yīng)用于集合預(yù)報(bào)和目標(biāo)觀測敏感區(qū)的研究，取得了滿意的結(jié)果。這些研究拓展了CNOP方法的應(yīng)用，揭示了CNOP更豐富的物理意義，為改進(jìn)天氣預(yù)報(bào)和氣候預(yù)測的預(yù)報(bào)技巧提供了新思路（Duan and Mu, 2009）。

近年來，CNOP方法得到了進(jìn)一步發(fā)展，而且在對我國天氣和氣候異常有重要影響的 ENSO事件、黑潮路徑變異、臺風(fēng)和阻塞事件的可預(yù)報(bào)性研究中的應(yīng)用取得了重要進(jìn)展。Mu et al.（2010）考慮到實(shí)際天氣預(yù)報(bào)和氣候預(yù)測初始誤差和模式誤差同時存在的情形，將CNOP方法拓展到既考慮初始擾動又考慮數(shù)值模式參數(shù)擾動。拓展的CNOP方法要求求解初始擾動和模式參數(shù)擾動的組合模態(tài)的最優(yōu)模態(tài)，尋找導(dǎo)致最大預(yù)報(bào)誤差的初始誤差和模式參數(shù)誤差的組合模態(tài)。應(yīng)用拓展的CNOP方法，Yu et al.（2009）和 Yu et al.（2012a）考察了影響ENSO事件SPB的誤差因子，強(qiáng)調(diào)了初始誤差是導(dǎo)致SPB的主要誤差來源，并且Yu et al.（2012b）和Duan and Wei（2012）先后在中等復(fù)雜程度模式和完全耦合的GCM模式的回報(bào)試驗(yàn)中發(fā)現(xiàn)了類似的初始誤差模態(tài)，從而闡明在實(shí)際預(yù)報(bào)中，如果通過資料同化或目標(biāo)觀測方法濾掉CNOP類型的初始誤差，ENSO預(yù)報(bào)技巧可能會被大大提高；拓展的CNOP方法也被Wang et al.（2012）用于黑潮路徑變異的可預(yù)報(bào)性研究，揭示了初始誤差是導(dǎo)致黑潮路徑變異預(yù)報(bào)不確定性的主要誤差因子，其伴隨的局地性特征為黑潮路徑變異預(yù)報(bào)的目標(biāo)觀測敏感區(qū)提供了有用的信息。另外，CNOP方法也被應(yīng)用于阻塞的可預(yù)報(bào)性，探討了阻塞事件最優(yōu)前期征兆（OPR）和最快增長初始誤差（OGR）的相似性及其伴隨的局地性特征，從而為確定阻塞預(yù)報(bào)目標(biāo)觀測敏感區(qū)提供了新方法，為提高數(shù)值天氣預(yù)報(bào)和氣候預(yù)測的預(yù)報(bào)技巧開拓了新思路（Jiang and Wang,2010; Mu and Jiang, 2011）。用CNOP確定目標(biāo)觀測敏感區(qū)的思路被較成功地應(yīng)用于臺風(fēng)的目標(biāo)觀測研究，為確定臺風(fēng)預(yù)報(bào)的初值敏感區(qū)、進(jìn)而提高臺風(fēng)預(yù)報(bào)技巧提供了重要指導(dǎo)意義（Wang and Tan,2009; Qin and Mu, 2011a, 2011b; Zhou and Mu, 2011;Chen and Mu, 2012; Zhou and Mu, 2012a, 2012b）?？梢?，CNOP方法在原有應(yīng)用研究的基礎(chǔ)上，學(xué)者們近年來進(jìn)一步將其用于實(shí)際天氣預(yù)報(bào)或氣候預(yù)測直接面對的一些困難問題（如目標(biāo)觀測敏感區(qū)）的研究，獲得了有特色的成果，為提高數(shù)值天氣預(yù)報(bào)和氣候預(yù)測的預(yù)報(bào)技巧提供了新思路。

該文將簡要回顧上述CNOP方法的拓展及其應(yīng)用的研究進(jìn)展。第二節(jié)將介紹拓展的CNOP方法及其物理意義；第三節(jié)敘述CNOP方法在ENSO可預(yù)報(bào)性問題研究中的應(yīng)用；第四節(jié)回顧C(jī)NOP在黑潮路徑變異可預(yù)報(bào)性研究中的應(yīng)用；第五節(jié)闡述阻塞前期征兆和最快增長初始誤差的相似性及其對目標(biāo)觀測的指示意義；第六節(jié)將說明CNOP在臺風(fēng)目標(biāo)觀測中的應(yīng)用。最后，在第七節(jié)討論上述成果對于改進(jìn)天氣預(yù)報(bào)和氣候預(yù)測的預(yù)報(bào)技巧的科學(xué)指導(dǎo)意義和應(yīng)用前景。

2 條件非線性最優(yōu)擾動方法的拓展

Mu et al.(2003) 提出的CNOP是指給定物理約束條件的具有最大非線性發(fā)展的一類初始擾動，它既可以用來研究最容易導(dǎo)致某一天氣或氣候事件的初始異常，又可以研究對預(yù)報(bào)結(jié)果具有最大影響的一類初始誤差，這兩個問題是天氣和氣候可預(yù)報(bào)性研究中的兩個關(guān)鍵問題。因此，CNOP在有關(guān)初始場的可預(yù)報(bào)性問題研究中具有重要作用?？紤]到實(shí)際天氣預(yù)報(bào)和氣候預(yù)測初始誤差和模式誤差同時存在的情形，Mu et al.（2010）進(jìn)一步將上述CNOP方法拓展到初始誤差和數(shù)值模式參數(shù)誤差同時存在的情形，并提出了條件非線性最優(yōu)參數(shù)擾動（CNOP-P）的概念。拓展的CNOP（為方便敘述，以下我們也稱之為 CNOP），即是求解初始擾動和模式參數(shù)擾動的最優(yōu)組合模態(tài)。CNOP有兩種特殊情形：一種是Mu et al.（2003）提出的僅與初始擾動有關(guān)，并且在預(yù)報(bào)時刻具有最大非線性發(fā)展的CNOP初始擾動，以下將之記為 CNOP-I；另一種情形則只與數(shù)值模式參數(shù)擾動有關(guān)，即CNOP-P，CNOP-P導(dǎo)致在預(yù)報(bào)時刻，預(yù)報(bào)結(jié)果與參考態(tài)有最大差別。CNOP方法既可以用來研究由初始誤差引起的第一類可預(yù)報(bào)性問題，又可以研究由參數(shù)誤差引起的第二類可預(yù)報(bào)性問題，同時也可用于研究初始誤差和參數(shù)誤差同時存在的情形。本節(jié)將簡要回顧C(jī)NOP方法及其物理意義。

2.1 CNOP的定義

假定狀態(tài)向量w的發(fā)展由以下動力系統(tǒng)方程描述:

其中，w( x,t) = (w1( x,t), w2(x,t) , …,wl( x,t))為包含l個狀態(tài)變量（如，海表溫度、溫躍層深度、風(fēng)速等）的狀態(tài)向量；w0是初始條件；x =(x1, x2,… ,xn),x ∈ Ω（Ω是歐氏空間Rn的一個集合）；t是時間，0≤t＜∞；p =( p1, p2, … , pm)是模式參數(shù)向量；F是非線性算子。

若動力系統(tǒng)方程（1）和初始條件已知，則未來時刻狀態(tài)可以通過對方程（1）積分確定，

這里Mτ(p)代表含有參數(shù)p的動力系統(tǒng)方程（1）從0時刻到τ時刻的傳播算子。如果方程（2）對應(yīng)于初始條件U0和U0+u0的解為U(t)和U(t)+u(t)，那么其在τ時刻的解為：

其中，u0是U(t) 的初始擾動，uI(τ)則代表了初始擾動的非線性發(fā)展。

如果p'是在模式參數(shù)向量P上的參數(shù)向量擾動,方程（2）對應(yīng)于P和P+ p'的兩個解可以表示為

這里up(τ) 代表了由于參數(shù)擾動p'導(dǎo)致的預(yù)報(bào)結(jié)果和參考狀態(tài)U(τ) 的偏差。

當(dāng)同時考慮初始條件和模式參數(shù)的擾動時, 方程（2）的兩個解為

這里uI,p(τ) 代表了由初始擾動和模式參數(shù)擾動共同導(dǎo)致的預(yù)報(bào)結(jié)果和參考態(tài)U (τ)的偏差。

定義如下非線性優(yōu)化問題:

其中，J (u0, p ') = Mτ(P + p')( U0+ u0)- Mτ(P) (U0),u0和p' 分別是疊加在參考態(tài)初值U0和模式參數(shù)向量P上的擾動，u0∈ Cδ, p′ ∈Cσ是各自的約束條件。通過解優(yōu)化問題（6），即可得到最優(yōu)初始擾動和參數(shù)擾動的聯(lián)合模態(tài)(u0δ; pσ′)。在給定約束條件下，該模態(tài)導(dǎo)致預(yù)報(bào)時刻τ的預(yù)報(bào)結(jié)果和參考態(tài)Mτ(P) (U0)的偏差最大。Mu et al.（2010）將該模態(tài)定義為 CNOP。CNOP包含兩種特殊情形：一種特殊情形是 CNOP-I (即優(yōu)化問題（7）式中的u0δ,I),即Mu et al.（2003）提出的CNOP初始擾動，其在預(yù)報(bào)時刻具有最大的非線性發(fā)展。CNOP-I可通過求解下面優(yōu)化問題得到:

CNOP的另一種特殊情形是CNOP-P (即（8）式中的p′σ,p)。CNOP-P代表了導(dǎo)致預(yù)報(bào)時刻預(yù)報(bào)結(jié)果與參考態(tài)的偏差最大的一類參數(shù)擾動，可通過解下述優(yōu)化問題得到：

2.2 CNOP的物理意義

CNOP代表了滿足一定物理約束條件的初始擾動和模式參數(shù)擾動的最優(yōu)組合模態(tài)，該模態(tài)導(dǎo)致預(yù)報(bào)結(jié)果與參考態(tài)的偏差最大。在可預(yù)報(bào)性研究中，CNOP代表了導(dǎo)致最大預(yù)報(bào)誤差的一類初始誤差和參數(shù)誤差的組合模態(tài)。

CNOP-I是CNOP的一種特殊情形。作為初始擾動，CNOP-I不僅可以描述最容易發(fā)展成某一天氣或氣候事件的最優(yōu)前期征兆，而且當(dāng)CNOP-I被考慮為疊加在某一天氣或氣候事件上的初始誤差時，它可以描述在預(yù)報(bào)時刻對預(yù)報(bào)結(jié)果影響最大的一類初始誤差；CNOP-I也可用于表征非線性最不穩(wěn)定（或最敏感）的初始擾動模態(tài)，該類擾動在預(yù)報(bào)時刻具有最大的非線性發(fā)展。

CNOP-P是CNOP的另一種特殊情形。CNOP-P可以代表在預(yù)報(bào)時刻對預(yù)報(bào)結(jié)果具有最大影響的一類參數(shù)誤差；在敏感性研究中，它也可以用來描述最敏感的參數(shù)擾動，從而為參數(shù)的敏感性排序提供有用的信息；通過CNOP-P描述的敏感參數(shù)，我們也可以據(jù)此確定最需要通過觀測定量確定的參數(shù)。

如上所述，CNOP可以代表導(dǎo)致最大預(yù)報(bào)誤差的初始誤差和參數(shù)誤差的組合模態(tài)，CNOP-I誤差和CNOP-P誤差分別代表導(dǎo)致最大預(yù)報(bào)誤差的初始誤差和參數(shù)誤差。因此，通過CNOP誤差，我們不僅可以探討初始誤差和參數(shù)誤差對預(yù)報(bào)結(jié)果不確定性的影響，而且可以分辨初始誤差和參數(shù)誤差的預(yù)報(bào)結(jié)果不確定性產(chǎn)生中的相對重要性，最終確定導(dǎo)致預(yù)報(bào)結(jié)果不確定性的主要誤差因子。

2.3 CNOP的計(jì)算

已有的優(yōu)化算法主要是用于求解極小值問題，但 CNOP優(yōu)化問題（6）是一個極大值問題，所以當(dāng)計(jì)算CNOP時，需要將極大值問題轉(zhuǎn)變?yōu)闃O小值問題。即將目標(biāo)函數(shù)（7）改寫為：

其中＜·＞表示內(nèi)積，這樣極大值問題（6）就轉(zhuǎn)換為極小值問題，從而可用已有的優(yōu)化算法求解（9）的極小值，得到 CNOPs。求解 CNOPs，一個重要的步驟是計(jì)算目標(biāo)函數(shù)關(guān)于初始擾動和參數(shù)擾動的梯度。Mu et al.（2010）給出了利用伴隨模式計(jì)算目標(biāo)函數(shù)J1的梯度公式，為進(jìn)一步使用數(shù)值模式計(jì)算目標(biāo)函數(shù)梯度提供了理論指導(dǎo)。Mu et al.（2010）給出的計(jì)算梯度公式為

其中[]?·表示伴隨算子，1()tλ和2()tλ滿足，利用上述梯度信息，我們可以使用優(yōu)化算法，如SPG2 (Spectral Projected Gradient 2) (Birgin et al.,2000)，SQP (Sequential Quadratic Programming)Powell et al., 1982) 和LBFGS-B (Limited memory Broyden-Fletcher-Goldfarb-Shanno for Boundconstrained optimization) (Liu and Nocedal, 1989;Zhu et al., 1997) 等計(jì)算CNOPs。

上述已經(jīng)表明，CNOP的求解需要采用伴隨技術(shù)，然而伴隨模式的編寫非常繁雜，需要花費(fèi)大量時間，限制了該方法的推廣應(yīng)用。為了克服這一困難，Wang and Tan（2009）提出了基于集合投影技術(shù)的CNOP求解算法，避免了伴隨編寫困難。該算法使用多個相互線性獨(dú)立或者說正交的初始擾動（x′）得到相應(yīng)的預(yù)報(bào)增量（y′），通過集合方法建立x′和y′之間的統(tǒng)計(jì)相關(guān)模型，用該統(tǒng)計(jì)模型向后積分計(jì)算目標(biāo)函數(shù)的梯度，以此進(jìn)行非線性最優(yōu)化迭代，從而得到滿足初始約束條件且在預(yù)報(bào)時刻有最大非線性發(fā)展的CNOPs。Wang and Tan（2009）表明，新算法能夠很好地抓住CNOP發(fā)展的非線性特征，而且基于集合投影的新算法和基于伴隨算法產(chǎn)生的CNOPs空間結(jié)構(gòu)幾乎完全相同。也就是說，基于集合投影計(jì)算CNOP的新算法是有效可行的。

新算法的提出在一定程度上推廣了CNOP在不同模式中的應(yīng)用。例如，Wang and Tan（2009）將新算法應(yīng)用于 GRAPES模式等研究臺風(fēng)的目標(biāo)觀測問題，從而肯定了CNOP新算法的有效性。

3 CNOP在El Ni?o事件春季可預(yù)報(bào)性障礙（SPB）問題研究中的應(yīng)用

ENSO事件的“春季可預(yù)報(bào)性障礙”（SPB）問題是ENSO可預(yù)報(bào)性研究的一個重要內(nèi)容。SPB現(xiàn)象是指ENSO預(yù)報(bào)具有較大預(yù)報(bào)誤差，且預(yù)報(bào)誤差的最快發(fā)展發(fā)生在春季和夏初，即預(yù)報(bào)誤差的發(fā)展呈現(xiàn)顯著的季節(jié)依賴性特征。國際上許多研究探討了ENSO事件的SPB問題，但SPB產(chǎn)生的機(jī)制至今仍撲朔迷離，具有很大爭議。而且，在業(yè)務(wù)預(yù)報(bào)中，ENSO事件的 SPB問題仍然嚴(yán)重影響著ENSO的預(yù)報(bào)技巧，采用何種有效的預(yù)報(bào)策略減弱或消除SPB至今尚不清楚。作者們近期的工作強(qiáng)調(diào)了初始誤差的特定結(jié)構(gòu)及其局地性特征在 El Ni?o事件SPB現(xiàn)象中的重要作用，進(jìn)而提出了用目標(biāo)觀測的方法提高ENSO預(yù)報(bào)技巧的新思路。

CNOP誤差導(dǎo)致最大預(yù)報(bào)誤差，所以最有可能導(dǎo)致 El Ni?o事件的 SPB現(xiàn)象 (Mu et al., 2007a,2007b)。因此，Duan and Zhang（2010）用一個理論ENSO模式（WF96模式；Wang and Fang, 1996）研究了CNOP誤差（即初始誤差和模式參數(shù)誤差的最優(yōu)聯(lián)合模態(tài)）對SPB的影響。結(jié)果表明，CNOP誤差導(dǎo)致El Ni?o事件產(chǎn)生最大的預(yù)報(bào)誤差，尤其在春季誤差增長最快，具有明顯的季節(jié)依賴性特征。CNOP誤差導(dǎo)致了El Ni?o事件顯著的SPB現(xiàn)象。進(jìn)一步研究表明，CNOP-P誤差導(dǎo)致較小的預(yù)報(bào)誤差，預(yù)報(bào)誤差在每個季節(jié)的增長率也很小，而且沒有明顯的季節(jié)依賴性，沒有導(dǎo)致El Ni?o事件的SPB現(xiàn)象。這些結(jié)果說明，初始誤差可能是導(dǎo)致El Ni?o事件SPB現(xiàn)象的主要誤差因子。Yu et al.（2012）進(jìn)一步利用中等復(fù)雜程度的Zebiak-Cane模式，研究了初始誤差和模式參數(shù)誤差對El Ni?o事件SPB現(xiàn)象的影響，考察了它們在SPB產(chǎn)生中的相對重要性。結(jié)果表明，CNOP-P誤差導(dǎo)致的預(yù)報(bào)誤差既沒有顯著的季節(jié)依賴性發(fā)展，對預(yù)報(bào)結(jié)果不確定性也沒有大的貢獻(xiàn)，因而未能導(dǎo)致El Ni?o事件的SPB現(xiàn)象。然而對于初始誤差，CNOP-I誤差不僅導(dǎo)致顯著的預(yù)報(bào)結(jié)果不確定性，而且呈現(xiàn)明顯的季節(jié)依賴性發(fā)展，發(fā)生了顯著的 SPB現(xiàn)象。CNOP-I誤差和CNOP-P誤差的組合模態(tài)也導(dǎo)致了明顯的 SPB現(xiàn)象，而且組合模態(tài)導(dǎo)致的預(yù)報(bào)誤差與CNOP-I誤差導(dǎo)致的預(yù)報(bào)誤差相差不大。所以，CNOP-I誤差可能更容易導(dǎo)致 SPB現(xiàn)象。綜上，Duan and Zhang(2010) 和Yu et al.(2012) 一致表明，初始誤差可能比參數(shù)誤差更容易導(dǎo)致El Ni?o事件的SPB現(xiàn)象，即初始誤差在El Ni?o事件SPB現(xiàn)象的產(chǎn)生中具有更加重要的作用，可能是SPB誤差增長的主要誤差來源。

Mu et al.（2007b）研究發(fā)現(xiàn)CNOP-I誤差容易導(dǎo)致El Ni?o事件顯著的SPB現(xiàn)象，而且他們也發(fā)現(xiàn)存在一些其他類型初始誤差，該類誤差與CNOP-I誤差大小相同，只是空間結(jié)構(gòu)不同，然而它們卻沒有導(dǎo)致大的預(yù)報(bào)誤差，也沒有明顯的季節(jié)依賴性發(fā)展，最終未能導(dǎo)致SPB。該結(jié)果說明SPB的發(fā)生依賴于初始誤差的特定空間結(jié)構(gòu)，即具有特定空間結(jié)構(gòu)的初始誤差可能更容易導(dǎo)致El Ni?o事件的SPB現(xiàn)象。為進(jìn)一步研究該問題，Yu et al.（2009）針對不同強(qiáng)度的El Ni?o事件，考察了CNOP-I誤差、LSV誤差和隨機(jī)初始誤差在上述El Ni?o事件SPB產(chǎn)生中的作用（圖1和2）。結(jié)果仍然表明，CNOP-I誤差導(dǎo)致El Ni?o事件顯著的SPB現(xiàn)象，但LSV誤差導(dǎo)致的SPB現(xiàn)象較CNOP-I誤差弱，隨機(jī)初始誤差導(dǎo)致的預(yù)報(bào)誤差則在每個季節(jié)的增長率可以忽略不計(jì)，而且其最大季節(jié)增長率沒有明顯的季節(jié)依賴性，未能導(dǎo)致SPB現(xiàn)象。所以，Yu et al.（2009）進(jìn)一步表明，El Ni?o事件SPB的發(fā)生依賴于初始誤差的空間結(jié)構(gòu)，具有CNOP-I誤差空間結(jié)構(gòu)的初始誤差更容易導(dǎo)致SPB。CNOP-I誤差可以分為兩類（Yu et al., 2009）：一類誤差的 SSTA分量在赤道東太平洋呈現(xiàn)正異常，而在赤道中西太平洋呈現(xiàn)負(fù)異常，溫躍層深度則沿著赤道一致加深；另一類誤差則呈現(xiàn)與前者幾乎相反的模態(tài)。由此可見，CNOP-I誤差盡管存在不同的空間結(jié)構(gòu)，但空間模態(tài)一致地具有局地性特征，即誤差能量集中于某一特定局地區(qū)域。

圖1 Zebiak-Cane模式中8次具有代表性的El Ni?o事件（初始變暖的時間分別為1月、4月、7月和10月的相對強(qiáng)和弱的El Ni?o事件）的CNOP-I誤差和LSV誤差導(dǎo)致的SSTA（SST anomaly）預(yù)報(bào)誤差（單位：°C）的季節(jié)增長率κ.預(yù)報(bào)的起始月為 (a) July (―1), (b) October (―1), (c) January (0)和 (d) April (0) (0和―1表示El Ni?o年和前一年).LSV誤差的季節(jié)增長率幾乎在每個季節(jié)都小于CNOP-I誤差的增長率.這里的預(yù)報(bào)誤差季節(jié)增長率用一個季節(jié)內(nèi)預(yù)報(bào)誤差的增長與該季節(jié)時間跨度的比值來度量Fig.1 Ensemble mean of the seasonal growth rate κ of the SSTA prediction errors (unit: °C) caused by the CNOP-I-type errors and the LSV-type errors for eight El Ni?o events in the Zebiak-Cane model, where the seasonal growth rate of the prediction errors is measured by the ratio of the growth of the prediction errors in one season to the time length of the season.The start months of the predictions are (a) July (―1), (b) October (―1), (c) January (0), and (d) April (0).The growth rates of the LSV-type errors in almost each season are smaller than those of the CNOP-type errors

圖2 Zebiak-Cane模式中8次El Ni?o事件的（a）隨機(jī)初始誤差的季節(jié)增長率及其（b）導(dǎo)致的SSTA預(yù)報(bào)誤差（單位：°C）的集合平均，橫坐標(biāo)的數(shù)字表示8次El Nino事件。引自Yu et al.(2009)Fig.2(a) The ensemble mean of the seasonal growth rate κ of the random initial errors for eight El Ni?o events in the Zebiak-Cane model; (b) the SSTA prediction errors (unit: °C) of the eight El Ni?o events that are caused by the random initial errors (the numbers on the horizontal axis denote eight El Nino events).From Yu et al.(2009)

上述兩類 CNOP-I誤差是通過考察模式中 El Ni?o事件而得到的對預(yù)報(bào)結(jié)果具有最大影響的初始誤差。那么，這些誤差是否存在于ENSO的實(shí)際回報(bào)試驗(yàn)中？我們試想，如果在ENSO預(yù)測的初始分析場中確實(shí)存在與 CNOP-I類似的誤差（記為CNOP-I-like誤差），那么我們可以通過合適的方法濾掉CNOP-I誤差，進(jìn)而提高ENSO預(yù)報(bào)技巧。因此，首先要弄清楚的一個關(guān)鍵問題是：在 ENSO的預(yù)測試驗(yàn)中是否存在CNOP-I誤差？Yu et al.（2012b）表明，在所考察的 Zebiak-Cane模式的最新版本(LDEO5) (Chen et al., 2004) 關(guān)于1980年1月至1999年12月的回報(bào)試驗(yàn)的240個初始分析場誤差中，與CNOP-I-like誤差模態(tài)相似性較高的誤差，一般對應(yīng)于較大的預(yù)報(bào)誤差，也就是說，在LDEO5關(guān)于El Ni?o的回報(bào)試驗(yàn)中，確實(shí)存在CNOP-like-I誤差，而且該類誤差對應(yīng)了較大的預(yù)報(bào)誤差。 Duan and Wei（2012）探討了復(fù)雜GCM模式—FGOLAS-g模式關(guān)于El Ni?o回報(bào)試驗(yàn)的初始分析場誤差。結(jié)果表明，在FGOALS-g關(guān)于1982～2007年ENSO事件的回報(bào)試驗(yàn)中，導(dǎo)致El Ni?o事件顯著SPB現(xiàn)象的初始分析誤差的主模態(tài)具有CNOP-I-like誤差的空間結(jié)構(gòu)。該結(jié)果在一定程度上說明，在FGOALS-g回報(bào)試驗(yàn)中，與CNOP-I-like相似的誤差更容易導(dǎo)致El Ni?o預(yù)報(bào)結(jié)果具有較大不確定性。從Zebiak-Cane模式和FGOALS-g模式的回報(bào)試驗(yàn)看，CNOP-I-like誤差可能確實(shí)存在于ENSO預(yù)測中。因此，我們應(yīng)該研究如何根據(jù)CNOP-I誤差的特點(diǎn)想辦法濾掉ENSO預(yù)測中的CNOP-I-like誤差,進(jìn)而提高ENSO的預(yù)報(bào)技巧。

綜上，初始誤差是El Ni?o事件SPB現(xiàn)象的主要誤差來源，而具有CNOP-I誤差模態(tài)的初始誤差更容易導(dǎo)致 SPB，該類誤差不僅是所考察模式 El Ni?o事件的最優(yōu)增長初始誤差，而且存在于有關(guān)模式的實(shí)際回報(bào)試驗(yàn)中，因此，如果在預(yù)報(bào)中通過合適的方法濾掉CNOP-I-like誤差，ENSO的預(yù)報(bào)技巧可能會大大提高；另一方面，CNOP-I誤差的空間分布具有明顯的局地性特征，即誤差主要集中在某一特定區(qū)域。Yu et al.（2012b）表明，在CNOP-I誤差集中的局地性區(qū)域減小誤差，ENSO預(yù)報(bào)誤差會大大減小，CNOP-I誤差集中的區(qū)域可能代表了ENSO預(yù)測的敏感區(qū)。該結(jié)果建議我們可以通過在CNOP-I敏感區(qū)進(jìn)行目標(biāo)觀測，消除CNOP-I-like誤差，改進(jìn)初始場，進(jìn)而提高ENSO預(yù)報(bào)技巧。

4 CNOP在黑潮路徑變異可預(yù)報(bào)性研究中的應(yīng)用

圖3 兩類CNOP-I誤差。左列為SSTA分量，右列為溫躍層深度異常分量。引自Yu et al.(2009)Fig.3 The patterns for the two categories of CNOP-I-type errors.The left column is for the SSTA component, while the right column is for the thermocline depth anomaly.From Yu et al.(2009)

圖4 CNOP-I的上層厚度分量，單位：m。引自Wang et al.(2012)Fig.4 Upper-layer thickness (unit: m) component of CNOP-I.From Wang et al.(2012)

黑潮路徑變異對區(qū)域氣候產(chǎn)生重要影響（Xu et al.2010），對我國長江流域降水也具有較大影響（侍茂崇，2004）。黑潮路徑變異的可預(yù)報(bào)性研究具有重要意義。國際上關(guān)于黑潮路徑變異可預(yù)報(bào)性的研究，主要從初始誤差增長的角度進(jìn)行研究（Komori et al.2003; Ishikawa et al.2004; Fujii et al.2008）。模式誤差對黑潮路徑變異預(yù)報(bào)的影響研究則很少。近年來，作者們將CNOP方法應(yīng)用于黑潮路徑變異的可預(yù)報(bào)性研究，分辨了初始誤差和模式參數(shù)誤差對黑潮路徑變異預(yù)報(bào)影響的相對重要性，強(qiáng)調(diào)了初始誤差的作用（Wang et al., 2012）。

Wang et al.（2012）采用1.5層淺水模式計(jì)算了黑潮大彎曲路徑的 CNOP-I（圖 4）和 CNOP-P，其中CNOP-P包括參數(shù)風(fēng)應(yīng)力的最優(yōu)擾動（CNOPTAU-P）、側(cè)向摩擦系數(shù)的最優(yōu)擾動（CNOPAH-P）、界面摩擦系數(shù)的最優(yōu)擾動（CNOPRI-P）和上述三個參數(shù)擾動的最優(yōu)組合模態(tài)（CNOPAll-P）。計(jì)算結(jié)果表明，如果以預(yù)報(bào)時刻擾動的動能作為預(yù)報(bào)誤差的度量，那么CNOP-I所導(dǎo)致的預(yù)報(bào)誤差比不同參數(shù)的最優(yōu)擾動導(dǎo)致的預(yù)報(bào)誤差大（圖5），也就是說，在黑潮路徑變異的預(yù)報(bào)中，初始誤差可能比模式參數(shù)誤差對預(yù)報(bào)結(jié)果不確定性的影響更大。具體地，由CNOP-I導(dǎo)致的海洋上層厚度的均方根誤差最大（圖6），并且大約在預(yù)報(bào)時刻120天時，均方根誤差開始大于均方根變化（RMS variation），這表明在預(yù)報(bào)時刻120天時，以具有CNOP-I誤差的初始場作為模式的初始場預(yù)報(bào)黑潮路徑變異時，該時刻的預(yù)報(bào)結(jié)果不在預(yù)報(bào)技巧允許的范圍內(nèi)；然而當(dāng)模式參數(shù)誤差為CNOPAll-P時，雖然CNOPAll-P也導(dǎo)致了顯著的預(yù)報(bào)誤差，但其在預(yù)報(bào)時刻210天時（遠(yuǎn)大于120天），預(yù)報(bào)結(jié)果才開始不能被接受。也就是說，對于給定的預(yù)報(bào)時刻，初始誤差比參數(shù)誤差更容易導(dǎo)致海洋上層厚度的預(yù)報(bào)不可接受。

圖5 CNOP-I誤差和不同參數(shù)的最優(yōu)擾動導(dǎo)致的預(yù)報(bào)誤差.橫軸上的數(shù)字1～5分別代表CNOPRI-P, CNOPTAU-P, CNOPAH-P, CNOPAll-P和 CNOP-I，縱軸表示橫軸所示的最優(yōu)擾動導(dǎo)致的預(yù)報(bào)誤差（用擾動動能度量）。引自Wang et al.(2012)Fig.5 The prediction errors caused by the CNOP-I and optimal perturbations of different parameters.The numbers 1-5 in the horizontal axis denote CNOPRI-P, CNOPTAU-P, CNOPAH-P, CNOPAll-P, and CNOP-I, respectively.The vertical axis denotes the prediction errors(measured by kinetic energy of perturbations).From Wang et al.(2012)

對 1.5層淺水模式來說，初始誤差，尤其是CNOP-I誤差，比參數(shù)誤差導(dǎo)致更大的預(yù)報(bào)誤差，更易導(dǎo)致黑潮路徑變異預(yù)報(bào)結(jié)果的不確定性（圖7），因此，改進(jìn)1.5層淺水模式的初始場比改進(jìn)參數(shù)化方案更容易使得黑潮路徑變異的預(yù)報(bào)技巧顯著提高。那么如何改進(jìn)初始場？事實(shí)上，黑潮路徑變異預(yù)報(bào)的CNOP-I誤差與ENSO預(yù)報(bào)的CNOP-I誤差有類似的特征，即CNOP-I誤差的空間模態(tài)具有明顯的局地性特征，該局地性特征使我們聯(lián)想到，可以通過目標(biāo)觀測方法濾掉模式初始場中的CNOP-I誤差，進(jìn)而提高黑潮路徑變異的可預(yù)報(bào)性。

5 阻塞發(fā)生的最優(yōu)前期征兆和最快增長初始誤差

圖6 不同誤差導(dǎo)致的海洋上層厚度均方根誤差的時間序列?！癛MS variation”表示海洋上層厚度的均方根變化，作為黑潮路徑預(yù)報(bào)誤差是否可以接受的判斷標(biāo)準(zhǔn)，該標(biāo)準(zhǔn)也被Komori et al.(2003) 用于研究黑潮路徑變異的預(yù)報(bào)時效。當(dāng)均方根誤差大于“RMS variation”時，預(yù)報(bào)結(jié)果不可接受；反之，預(yù)報(bào)結(jié)果能被接受。引自Wang et al.(2012)Fig.6 Time series of root-mean-square error (RMSE) of upper-layer thickness caused by different errors.RMS variation, the root mean square variation of the upper-layer thickness, is used as a criterion determining whether the prediction results of Kuroshio path are acceptable.This criterion was also used by Komori et al.(2003) to study the time limits for the prediction of Kuroshio path variation.When RMSE is larger than RMS variation, the prediction results are unacceptable; otherwise, the predictions are acceptable.From Wang et al.(2012)

圖7 模式具有 CNOPRI-P（黃線），CNOPTAU-P（紅線），CNOPAH-P（藍(lán)線），CNOPAll-P（灰線）和CNOP-I（黑實(shí)線）誤差時，關(guān)于參考態(tài)黑潮大彎曲路徑（黑虛線）的預(yù)報(bào)結(jié)果，其中黑潮流軸由海洋上層厚度520 m等值線表征。CNOP-I誤差導(dǎo)致的黑潮路徑的變異比其他類型誤差更顯著。引自Wang et al.(2012)Fig.7 The reference Kuroshio large meander axis (black dashed line) and the predicted Kuroshio axes caused by CNOPRI-P (yellow line),CNOPTAU-P (red line), CNOPAH-P (blue line), CNOPAll-P (gray line), and CNOP-I (black solid line).The Kuroshio axis is described by the 520-m contour of the upper-layer thickness.CNOP-I results in larger prediction error than other types of errors.From Wang et al.(2012)

阻塞是一個典型的大尺度環(huán)流，其對區(qū)域天氣和氣候有著深遠(yuǎn)的影響（Rex，1950）。阻塞發(fā)生對初始場的不確定性具有極端敏感性，阻塞環(huán)流轉(zhuǎn)型預(yù)報(bào)常常因此失效。所以，阻塞預(yù)報(bào)能力的提高仍然是中長期數(shù)值天氣預(yù)報(bào)研究的核心問題之一（Tibaldi and Molteni 1990; Kimoto et al.1992;Frederiksen et al.2005）。最優(yōu)前期征兆（OPR）和增長最快初始誤差（OGR）是可預(yù)報(bào)性研究的兩個重要問題，許多研究基于線性近似方法研究阻塞的OPR 和OGR（Frederiksen, 1998；Molteni and Palmer,1993）。然而，線性近似方法不能充分刻畫非線性物理過程的影響，存在著較大的局限性。

考慮到線性方法的局限性，Jiang and Wang（2010）用CNOP方法，基于Liu（1994）提出的阻塞指數(shù)定義，用T21L3準(zhǔn)地轉(zhuǎn)模式計(jì)算了最容易激發(fā)歐洲—大西洋阻塞事件發(fā)生的 OPR（圖 8）。該OPR的空間模態(tài)呈現(xiàn)很強(qiáng)的局地性特征，擾動的主要能量位于阻塞區(qū)域的上游。該 OPR主要呈現(xiàn)為北美大陸的波列結(jié)構(gòu)，且隨高度西傾，分布于大西洋急流的北側(cè)，軸線沿著東北—西南方向。該OPR隨著時間的演變，逐漸向下游傳播并發(fā)展，最終在大西洋區(qū)域形成典型的偶極子阻塞結(jié)構(gòu)。上述OPR擾動能量的發(fā)展主要來自于基本流的水平切變，斜壓調(diào)整僅起著小部分的作用。OPR擾動引發(fā)的行星尺度渦度平流向阻塞的高壓區(qū)（低壓區(qū)）輸送反氣旋渦度（氣旋渦度）（圖 9），從而促進(jìn)了偶極子阻塞的發(fā)生。

Mu and Jiang（2011）同樣基于T21L3準(zhǔn)地轉(zhuǎn)模式，選取了歐洲—大西洋區(qū)域1985年冬至1991年冬的20個個例，進(jìn)一步探討了阻塞發(fā)生的OPR；而且他們又探討了阻塞預(yù)報(bào)的OGR。OGR分為兩類，一類為局部CNOP-I，記為type-1-OGR（圖9）；另一類為全局CNOP-I，記為type-2-OGR（圖10）。結(jié)果發(fā)現(xiàn)，OPR和OGR具有很高的相似性。以優(yōu)化時間3天為例，阻塞的OPR與type-1-OGR的相似指數(shù)達(dá)到0.95，與type-2-OGR的相似指數(shù)達(dá)到―0.87。優(yōu)化時間為5天時，OPR與type-1-OGR相似指數(shù)達(dá)到0.85，與type-2-OGR相似指數(shù)達(dá)到―0.69。而且，OPR和 type-1-OGR具有類似的演變行為，且都發(fā)展成為歐洲—大西洋區(qū)域的偶極子阻塞模態(tài)；type-2-OGR則隨時間的演變，發(fā)展成為上述偶極子阻塞的相反模態(tài)。也就是說，阻塞發(fā)生的 OPR與OGR不僅在初始結(jié)構(gòu)上存在著相似性，在發(fā)展機(jī)理上也具有很高的相似性。

綜上所述，阻塞發(fā)生的OPR與OGR的空間模態(tài)高度相似，且具有明顯的局地性特征。這種相似性及其局地性特征啟發(fā)我們：如果在 OGR局地性特征描述的誤差大值區(qū)實(shí)施目標(biāo)觀測，一方面可以消除 OGR出現(xiàn)的可能性，改進(jìn)初始場；另一方面可以改善局地觀測網(wǎng)，有助于捕捉最優(yōu)前期征兆的信號，從而提高我們對阻塞事件的預(yù)報(bào)能力。

6 CNOP在臺風(fēng)目標(biāo)觀測研究中的應(yīng)用

目標(biāo)觀測是指為了提高高影響天氣和氣候事件的預(yù)報(bào)技巧，在一些關(guān)鍵地區(qū)（敏感區(qū)）加強(qiáng)觀測（Lorenz and Emanuel, 1998）的一種觀測手段。在一系列理論和方法研究的基礎(chǔ)上，國際上針對熱帶氣旋、冬季風(fēng)暴與溫帶氣旋等，進(jìn)行了一系列外場試驗(yàn)（Snyder, 1996; Langland et al., 1999;Szunyogh et al., 2000; Wu et al., 2005; Rabier et al.,2008）。這些試驗(yàn)的結(jié)果表明，目標(biāo)觀測能夠有效提高預(yù)報(bào)技巧。目標(biāo)觀測的關(guān)鍵之一是確定敏感區(qū)。目前確定目標(biāo)觀測敏感區(qū)的主要方法有兩大類，一類是線性方法（Gelaro et al., 1998; Palmer et al., 1998; Bergot, 1999; Buizza and Montani, 1999）；另一類是集合方法（Bishop and Toth, 1999）。然而這兩類方法都是基于初始擾動發(fā)展的線性性質(zhì)這一假定，不能考慮非線性物理過程的影響。

圖8 500 hPa上優(yōu)化時間為3天的激發(fā)歐洲—大西洋阻塞發(fā)生的OPR及其演變（單位：gpm）。（a）第0天；（b）第1天；（c）第2天；（d）第3天。歐洲—大西洋阻塞的 OPR（見（a）圖所示）主要呈現(xiàn)為北美大陸的波列結(jié)構(gòu)，且隨高度西傾，分布于大西洋急流的北側(cè)，軸線沿著東北—西南方向。該 OPR隨著時間的演變，逐漸向下游傳播并發(fā)展，最終在大西洋區(qū)域形成典型的偶極子阻塞結(jié)構(gòu)（見（b），（c）和（d）圖所示）。引自Jiang and Wang (2010)Fig.8(a) The optimal precursor (OPR) triggering the North Atlantic–European blocking at 500 hPa with optimization time 3 days and its nonlinear evolution(gpm) with (b) leading time 1 day; (c) leading time 2 days; and (d) leading time 3 days.The OPR (see Fig.(a)) triggering the North Atlantic–European blocking presents wave trains over North America, westward with height.The disturbances, which mostly focus on the northward flanks of the corresponding Atlantic upper-level jet, take on northeast–southwest trend.With time development, it propagates downwards and develops into a typical dipole blocking over the North Atlantic–European area (see Figs.(b), (c), and (d)).From Jiang and Wang (2010)

為揭示非線性的影響，一些研究將CNOP方法應(yīng)用于暴雨和臺風(fēng)的目標(biāo)觀測敏感區(qū)的研究中，表明了CNOP初始擾動，即CNOP-I模態(tài)具有明顯的局地性特征，嘗試將其應(yīng)用于確定臺風(fēng)預(yù)報(bào)的敏感區(qū)，并在敏感區(qū)減小誤差，考察其對臺風(fēng)預(yù)報(bào)結(jié)果的影響。結(jié)果表明，CNOP-I敏感區(qū)較LSV敏感區(qū)更有潛力改進(jìn)預(yù)報(bào)技巧（穆穆等, 2007; 譚曉偉,2009; Mu et al., 2009; Wang and Tan, 2009; Chen,2011; Qin and Mu, 2011a, 2011b; Zhou and Mu, 2011;Chen and Mu, 2012; Zhou and Mu, 2012a, 2012b）。

圖9 500 hPa上20個個例的type-1-OGR的合成場。（a）第0天及其在（b）第1天、（c）第2天、（d）第3天的非線性演變的合成場（單位：gpm）。引自Mu and Jiang (2011)Fig.9(a) The composite field of type-1 OGR of 20 cases at 500 hPa and its nonlinear evolution with (b) lead time 1 day, (c) lead time 2 days, and (d) lead time 3 days.From Mu and Jiang (2011)

近來，作者們進(jìn)一步通過OSSE和/或OSE研究了CNOP-I確定的敏感區(qū)的有效性。Chen and Mu（2012）利用 MM5模式和 WRF模式，以臺風(fēng)Longwang（2005）和 Sinlaku（2008）為例，通過OSSE試驗(yàn)研究了不同方法確定的敏感區(qū)和不同類型初始誤差對臺風(fēng)初值敏感區(qū)的影響以及敏感區(qū)對模式的依賴性。結(jié)果表明，敏感區(qū)內(nèi)的隨機(jī)初始誤差的發(fā)展一般要大于其他隨機(jī)選擇的區(qū)域內(nèi)誤差的發(fā)展；而在所有選擇的區(qū)域中，CNOP-I所識別的敏感區(qū)內(nèi)的隨機(jī)誤差的發(fā)展是最大的；在同一個區(qū)域中，具有特定空間結(jié)構(gòu)的初始誤差的發(fā)展一般要大于隨機(jī)誤差的發(fā)展，如 CNOP-I、singular vectors（SVs）和LSV類型的初始誤差的結(jié)構(gòu)。綜合比較研究的結(jié)果，發(fā)現(xiàn)CNOP-I所識別的敏感區(qū)內(nèi)CNOP-I類型初始誤差對預(yù)報(bào)結(jié)果不確定性的影響最大。將MM5模式計(jì)算得到的臺風(fēng)預(yù)報(bào)的敏感區(qū)應(yīng)用于WRF模式，結(jié)果發(fā)現(xiàn)由MM5模式確定的敏感區(qū)在WRF模式中也能非常有效地減小預(yù)報(bào)誤差，說明臺風(fēng)預(yù)報(bào)的初值敏感區(qū)對模式（至少在MM5和WRF模式）的依賴性可能較弱。

針對上述CNOP-I確定的敏感區(qū)，Chen（2011）結(jié)合臺灣的“追風(fēng)計(jì)劃”收集的下投探空資料，通過OSE試驗(yàn)表明，在CNOP-I敏感區(qū)內(nèi)同化目標(biāo)觀測資料，各氣象要素的總擾動能量（TPE）的均方根誤差（RMSE）的減小幅度要明顯大于在隨機(jī)區(qū)域的情況。這說明在CNOP-I確定的敏感區(qū)內(nèi)進(jìn)行目標(biāo)觀測，確實(shí)能夠有效減小臺風(fēng)預(yù)報(bào)誤差。

圖10 500 hPa上20個個例的type-2-OGR的合成場。（a）第0天及其在（b）第1天、（c）第2天、（d）第3天的非線性演變的合成場（單位：gpm）。引自Mu and Jiang (2011)Fig.10(a) The composite field of type-2 OGR of 20 cases at 500 hPa and its nonlinear evolution with (b) lead time 1 day, (c) lead time 2 days, and (d) lead time 3 days.From Mu and Jiang (2011)

Qin and Mu（2011b）將其應(yīng)用于2009年西太平洋7個個例臺風(fēng)路徑預(yù)報(bào)的研究。結(jié)果表明，選擇占全場 1%的敏感區(qū)作為目標(biāo)觀測區(qū)。將同化了觀測資料后的初始場作為進(jìn)行了目標(biāo)觀測后改善的初始場，進(jìn)行每6 h一次的路徑預(yù)報(bào)，總共預(yù)報(bào)時長為72 h。試驗(yàn)結(jié)果表明：7個臺風(fēng)個例中有6個個例的路徑預(yù)報(bào)得到了明顯改善（圖 11）。在CNOP敏感區(qū)內(nèi)進(jìn)行目標(biāo)觀測，路徑預(yù)報(bào)改善程度為13%～46%，48 h路徑預(yù)報(bào)誤差平均減少89.6 km左右，驗(yàn)證區(qū)內(nèi)預(yù)報(bào)氣壓平均減小 7 hPa；在 SVs敏感區(qū)內(nèi)進(jìn)行目標(biāo)觀測，改善程度僅為14%～23%；而且上述改善不僅出現(xiàn)在優(yōu)化時間段內(nèi)（24～48 h），甚至延續(xù)到72 h，臺風(fēng)路徑72 h預(yù)報(bào)誤差平均減少25.9 km。這說明CNOP適應(yīng)性資料同化技術(shù)可以有效地改進(jìn)臺風(fēng)路徑的預(yù)報(bào)技巧。

綜上，CNOP方法比其他方法更能抓住臺風(fēng)預(yù)報(bào)初始誤差的敏感區(qū)域，OSEs (OSSEs) 試驗(yàn)驗(yàn)證了CNOP確定的敏感區(qū)在改進(jìn)臺風(fēng)預(yù)報(bào)技巧中的有效性。所以，用CNOP方法所識別的敏感區(qū)是有效的，通過在敏感區(qū)加密觀測或者進(jìn)行資料同化，可以大大提高臺風(fēng)的預(yù)報(bào)技巧。

7 討論

圖11 左圖：臺風(fēng)環(huán)高（Vamco）在CNOP（實(shí)心菱形）和SVs（空心菱形）敏感區(qū)進(jìn)行與不進(jìn)行目標(biāo)觀測預(yù)報(bào)路徑誤差（24～72 h）比較；右圖：路徑預(yù)報(bào)誤差減小程度。引自Qin and Mu (2011b)Fig.11 Scatter diagrams for cases Vamco (left).The y-axis represents the track forecast errors with dropsondes, and the x-axis represents those without dropsondes.Filled and empty diamonds denote the results of CNOPs and SVs, respectively.The colour of each diamond indicates the forecast time.Histograms on the right are relative differences corresponding to each case.From Qin and Mu (2011b)

CNOP方法已被應(yīng)用于臺風(fēng)的目標(biāo)觀測問題研究，CNOP-I在確定敏感區(qū)中發(fā)揮了重要作用。上述所考察臺風(fēng)個例的CNOP-I誤差具有明顯的局地性特征，即誤差主要分布在某一特定區(qū)域，該區(qū)域可能代表了目標(biāo)觀測敏感區(qū)，而且作者們通過OSSEs或OSEs試驗(yàn)，驗(yàn)證了這一想法，即通過在敏感區(qū)進(jìn)行資料同化，減小敏感區(qū)誤差，顯著改進(jìn)了臺風(fēng)的預(yù)報(bào)技巧。該事實(shí)告訴人們，在 CNOP-I確定的敏感區(qū)加強(qiáng)觀測，對于改進(jìn)臺風(fēng)預(yù)報(bào)技巧具有重要作用。事實(shí)上，對于本文介紹的高影響天氣和氣候事件：El Ni?o事件、黑潮路徑變異以及阻塞事件，它們的CNOP-I誤差都導(dǎo)致顯著的預(yù)報(bào)誤差，且誤差也主要分布在某一特定區(qū)域，具有明顯的局地性特征?？紤]到臺風(fēng)的目標(biāo)觀測問題，我們自然聯(lián)想到，這些事件的CNOP-I呈現(xiàn)的局地性區(qū)域是否也代表了初始誤差敏感區(qū)，在該區(qū)域加強(qiáng)觀測或進(jìn)行資料同化，上述事件的預(yù)報(bào)技巧是否也有顯著提高？等等問題都是需要進(jìn)一步解決的問題。因此，在未來的工作中，我們有必要嘗試開展關(guān)于El Ni?o、黑潮路徑變異以及阻塞事件等高影響天氣和氣候事件的目標(biāo)觀測研究。

另一方面，上述工作提到，阻塞事件的CNOP-I前期征兆和初始誤差的空間模態(tài)具有高度相似性。因此，如果阻塞事件的CNOP-I誤差局地性區(qū)域代表了阻塞預(yù)報(bào)的敏感區(qū)，那么表現(xiàn)在前期征兆中的同一局地性區(qū)域說明在 OGR敏感區(qū)增加觀測不僅可以提高預(yù)報(bào)技巧，而且有利于增強(qiáng)識別阻塞事件發(fā)生的前期信號，進(jìn)而提高阻塞事件的預(yù)報(bào)技巧。也就是說，阻塞事件前期征兆和最快增長初始誤差的相似性和局地性特征可以為目標(biāo)觀測的布控提供有用信息。既然前期征兆和初始誤差的相似性和局地性特征對于目標(biāo)觀測有重要作用，那么發(fā)生在阻塞事件中的相似性和局地性特征是否在 El Ni?o事件、黑潮路徑變異以及臺風(fēng)等高影響天氣和氣候事件中也存在？這種相似性和局地性特征是否確實(shí)能夠指導(dǎo)這些天氣和氣候事件的目標(biāo)觀測敏感區(qū)確定，提高預(yù)報(bào)技巧？這些問題也是亟待解決的重要問題。

總之，可預(yù)報(bào)性研究是一個具有挑戰(zhàn)性的課題，仍然有大量的科學(xué)問題需要解決，這不僅需要學(xué)者具備扎實(shí)的氣象學(xué)知識，而且需要具備堅(jiān)實(shí)的數(shù)理基礎(chǔ)和高性能計(jì)算機(jī)計(jì)算的能力。隨著計(jì)算機(jī)的發(fā)展以及有效的學(xué)科交叉和來自不同領(lǐng)域的學(xué)者的通力合作，大力提高天氣和氣候可預(yù)報(bào)性的目的將能很快變?yōu)楝F(xiàn)實(shí)。

（References）

Bergot T.1999.Adaptive observations during FASTEX: A systematic survey of upstream flights [J].Quart.J.Roy.Meteor.Soc., 125 (561):3271–3298.

Birgin E G, Martinez M J, Raydan M.2000.Nonmonotone spectral projected gradient methods on convex sets [J].SIAM J.Optim., 10 (4):1196–1211.

Bishop C H, Toth Z.1999.Ensemble transformation and adaptive observations [J].J.Atmos.Sci., 56: 1748–1765.

Buizza R, Montani A.1999.Targeting observations using singular vectors[J].J.Atmos.Sci., 56: 2965–2985.

Chen B Y.2011.Observation system experiments for typhoon Nida (2004)using the CNOP method and DOTSTAR data [J].Atmos.Oceanic Sci.Lett., 4 (2): 118–123.

Chen B Y, Mu M.2012.The roles of spatial locations and patterns of initial errors in the uncertainties of tropical cyclone forecasts [J].Adv.Atmos.Sci., 29 (1): 63–78.

Chen D, Cane M A, Kaplan A, et al.2004.Predictability of El Ni?o over the past 148 years [J].Nature, 428: 733–736.

Duan W S, Mu M, Wang B.2004.Conditional nonlinear optimal perturbation as the optimal precursors for ENSO events [J].J.Geophys.Res., 109 (D23): 4756–4768.

Duan W S, Mu M.2009.Conditional nonlinear optimal perturbation:Applications to stability, sensitivity, and predictability [J].Science in China (D), 52 (7): 884–906.

Duan W S, Liu X C, Zhu K Y, et al.2009.Exploring the initial errors that cause a significant “spring predictability barrier” for El Ni?o events [J].J.Geophys.Res., 114: C04022, doi:10.1029/2008JC004925.

Duan W, Zhang R.2010.Is model parameter error related to a significant spring predictability barrier for El-Ni?o event? Result from Theoretical Model [J].Adv.Atmos.Sci., 27: 1003–1013.

Duan W S, Wei C.2012.The spring predictability barrier for El Ni?o events and its possible mechanism results from a fully coupled model [J].Int.J.Climatol., doi:10.1002/joc.3513，in press.

Duan W S, Yu Y S, Xu H, et al.2012.Behaviors of nonlinearities modulating the El Ni?o events induced by optimal precursory disturbances [J].Climate Dyn., doi:10.1007/s00382–012-1557-z，in press.

Frederiksen J S.1998.Precursors to blocking anomalies: The tangent linear and inverse problems [J].J.Atmos.Sci., 55: 2419–2436.

Frederiksen J S, Collier M A, Watkins A B.2005.Dependence of ensemble prediction skill on blocking instability regimes [J].Tellus, 56A: 485–500.

Fujii Y, Tsujino H, Usui N, et al.2008.Application of singular vector analysis to the Kuroshio large meander [J].J.Geophys.Res., 113:C07026, doi:10.1029/2007JC004476.

Gelaro R, Buizza R, Palmer T N, et al.1998.Sensitivity analysis of forecast errors and the construction of optimal perturbations using singular vectors [J].J.Atmos.Sci., 55: 1012–1037.

Ishikawa Y, Awaji T, Komori N, et al.2004.Application of sensitivity analysis using an adjoint model for short-range forecasts of the Kuroshio path south of Japan [J].Journal of Oceanography, 60: 293–301.

Jiang Zhina, Wang Donghai.2010.A study on precursors to blocking anomalies in climatological flows by using conditional nonlinear optimal perturbations [J].Quart.J.Roy.Meteor.Soc., 136 (650): 1170–1180.

Kimoto M, Mukougawa H, Yoden S.1992.Medium-range forecast skill variation and blocking transition: A case study [J].Mon.Wea.Rev., 120:1616–1627.

Komori N, Awaji T, Ishikawa Y, et al.2003.Short-range forecast experiments of the Kuroshio path variabilities south of Japan using TOPEX/Poseidon altimetric data [J].J.Geophys.Res., 108 (C1): 3010,doi:10.1029/2001JC001282.

Langland R H, Toth Z, Gelaro R, et al.1999.The North Pacific Experiment(NORPEX-98).Targeted observations for improved North American weather forecasts [J].Bull.Amer.Meteor.Soc., 80: 1363–1384.

Liu D C, Nocedal J.1989.On the limited memory BFGS method for large scale optimization [J].Math.Program., 45 (3): 503–528.

Liu Q.1994.On the definition and persistence of blocking [J].Tellus, 46 (3):286–298.

Lorenz E N.1965.A study of the predictability of a 28-variable atmospheric model [J].Tellus, 17: 321–333

Lorenz E N, Emanuel K A.1998.Optimal sites for supplementary observation sites: Simulation with a small model [J].J.Atmos.Sci., 55:399–414.

Molteni F, Palmer T N.1993.Predictability and finite-time instability of the northern winter circulation [J].Quart.J.Roy.Meteor.Soc., 119: 269–198.

Moore A M, Kleeman R.1996.The dynamics of error growth and predictability in a coupled model of ENSO [J].Quart.J.Roy.Meteor.Soc., 122: 1405–1446

Mu M, Duan W S, Wang B.2003.Conditional nonlinear optimal perturbation and its applications [J].Nonlin.Processes Geophys., 10:493–501.

Mu M, Sun L, Henk D A.2004.The sensitivity and stability of the ocean’s thermocline circulation to finite amplitude freshwater perturbations [J].J.Phys.Oceanogr., 34: 2305–2315.

穆穆, 王洪利, 周菲凡.2007.條件非線性最優(yōu)擾動方法在適應(yīng)性觀測研究中的初步應(yīng)用 [J].大氣科學(xué), 31 (6): 1102–1112. Mu Mu, Wang Hongli, Zhou Feifan.2007.A preliminary application of conditional nonlinear optimal perturbation to adaptive observation [J].Chinese Journal of Atmospheric Sciences (in Chinese), 31 (6): 1102–1112.

Mu M, Duan W S, Wang B.2007a.Season-dependent dynamics of nonlinear optimal error growth and El Ni?o-Southern Oscillation predictability in a theoretical model [J].J.Geophys.Res., 112: D10113,doi:10.1029/2005JD006981.

Mu M, Xu H, Duan W S.2007b.A kind of initial errors related to “spring predictability barrier” for El Ni?o events in Zebiak-Cane model [J].Geophys.Res.Lett., 234: L03709, doi:10.1029/2006GL027412.

Mu M, Zhou F F, Wang H H.2009.A method for identifying the sensitive areas in targeted observations for tropical cyclone prediction: Conditional nonlinear optimal perturbation [J].Mon.Wea.Rev., 137: 1623–1639.

Mu M, Duan W S, Wang Q, et al.2010.An extension of conditional nonlinear optimal perturbation and its applications [J].Nonlin.Processes Geophys., 2010, 17: 211–220.

Mu M, Jiang Z N.2011.Similarities between optimal precursors that trigger the onset of blocking events and optimally growing initial errors in onset prediction [J].J.Atmos.Sci., 68: 2860–2877.

Palmer T N, Gelaro R, Barkmeijer J, et al.1998.Singular vectors, metrics,and adaptive observations [J].J.Atmos.Sci., 55: 633–653.

Powell M J D.1982.VMCWD: A FORTRAN subroutine for constrained optimization [R].DAMTP Report 1982/NA4, University of Cambridge,England.

Qin X H, Mu M.2011a.A study on the reduction of forecast error variance by three adaptive observation approaches for tropical cyclone prediction[J].Mon.Wea.Rev., 139: 2218–2232.

Qin X H, Mu M.2011b.Influence of conditional nonlinear optimal perturbations sensitivity on typhoon track forecasts [J].Quart.J.Roy.Meteor.Soc., 138: 185–197.

Rabier F, Gauthier P, Cardinali C, et al.2008.An update on THORPEX-related research in data assimilation and observing strategies [J].Nonlin.Processes Geophys., 15: 81–94.

Rex D F.1950.Blocking action in the middle troposphere and its effects upon regional climate.I: An aerological study of blocking action [J].Tellus, 2: 196–211.

Riviere O, Lapeyre G, Talagrand O.2008.Nonlinear generalization of singular vectors: Behavior in a baroclinic unstable flow [J].J.Atmos.Sci.,65: 1896–1911.

Samelson R G, Tziperman E.2001.Instability of the chaotic ENSO: The growth-phase predictability barrier [J].J.Atmos.Sci., 58: 3613–3625.

侍茂崇.2004.物理海洋學(xué) [M].濟(jì)南: 山東教育出版社.462pp.Shi Maochong.2004.Physical Oceanography (in Chinese) [M].Jinan:Shandong Education Press, 462pp.

Snyder C.1996.Summary of an informal workshop on adaptive observations and FASTEX [J].Bull.Amer.Meteor.Soc., 77: 953–961.

Szunyogh I, Toth Z, Morss R E, et al.2000.The effect of targeted dropsonde observations during the 1999 Winter Storm Reconnaissance Program [J].Mon.Wea.Rev., 128: 3520–3537.

譚曉偉.2009.CNOP新算法研究及其在目標(biāo)觀測中的應(yīng)用檢驗(yàn) [D].中國科學(xué)院大氣物理研究所博士學(xué)位論文, 99pp.Tan Xiaowei.2009.Studies on new algorithms to obtain CNOP and their applications in the targeting observations [D].Ph.D.dissertation (in Chinese), Institute of Atmospheric Physics, Chinese Academy of Sciences, 99pp.

Tennekes H.1991.Karl Popper and the accountability of numerical forecasting [C] // New Developments in Predictability.ECMWF Workshop Proceedings, ECMWF, Shinfield Park, Reading UK.

Terwisscha van Scheltinga A D.2007.Data assimilation with implicit ocean models [D].Ph.D.dissertation.Institute for Marine and Atmospheric research, Utrecht University, Utrecht, the Netherlands, 119pp.

Thompson P.1957.Uncertainty of the initial state as a factor in the predictability of large scale atmospheric flow patterns [J].Tellus, 9: 275–295

Tibaldi S, Molteni F.1990.On the operational predictability of blocking [J].Tellus, 42A: 343–365.

Wang B, Fang Z.1996.Chaotic oscillation of tropical climate: A dynamic system theory for ENSO [J].J.Atmos.Sci., 53: 2786-2802.

Wang B, Tan X.2009.A fast algorithm for solving CNOP and associated target observation tests [J].Acta Meteorologica Sinica, 23: 387–402.

Wang Q, Mu M, Dijkstra H A.2012.Application of the conditional nonlinear optimal perturbation method to the predictability study of the Kuroshio large meander [J].Adv.Atmos.Sci., 29 (1): 118–134.

Wu X G, Mu M.2009.Impact of horizontal diffusion on the nonlinear stability of thermohaline circulation in a modified box model [J].J.Phys.Oceanogr., 39: 798–805.

Wu C C, Lin P H, Aberson S, et al.2005.Dropwindsonde observations for typhoon surveillance near the Taiwan region (DOSTAR): An overview [J].Bull.Amer.Meteor.Soc., 86: 787–790.

Xu H, Tokinaga H, Xie S P.2010.Atmospheric effects of the Kuroshio large meander during 2004–05 [J].J.Climate., 23: 4704–4715.

Yu Y S, Duan M S, Xu H, et al.2009.Dynamics of nonlinear error growth and season-dependent predictability of El Ni?o events in the Zebiak-Cane model [J].Quart.J.Roy.Meteor.Soc., 135: 2146–2160.doi:10.1002/qj.526.

Yu Y S, Mu M, Duan W S.2012a.Does model parameter error cause a significant spring predictability barrier for El Ni?o events in the Zebiak-Cane model? [J].J.Climate., 25: 1263–1277.

Yu Y S, Mu M, Duan W S, et al.2012b.Contribution of the location and spatial pattern of initial error to uncertainties in El Ni?o predictions [J].J.Geophys.Res.-Ocean, 117: C06018, doi:10.1029/2011JC007758.

Zhou F F, Mu M.2011.The impact of verification area design on tropical cyclone targeted observations based on the CNOP method [J].Adv.Atmos.Sci., 28 (5): 997–1010.

Zhou F F, Mu M.2012a.The impact of horizontal resolution on the CNOP and on its identified sensitive areas for tropical cyclone predictions [J].Adv.Atmos.Sci., 29 (1): 36–46.

Zhou F F, Mu M.2012b.The time and regime dependences of sensitive areas for tropical cyclone prediction using the CNOP method [J].Adv.Atmos.Sci., doi:10.1007/s00376-012-1174-0.

Zhu C Y, Byrd R H, Lu P H, et al.1997.Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization [J].ACM T.Math.Software, 23 (4): 550–560.