青藏高原和華東地區(qū)背景誤差協(xié)方差特征的對比研究

陳耀登,趙幸,閔錦忠,范水勇,王元兵,曾臘梅

(1.南京信息工程大學(xué) 氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心,江蘇 南京 210044;

2.氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室(南京信息工程大學(xué)),江蘇 南京 210044;

3.中國氣象局 北京城市氣象研究所,北京 100089;4.湖南省氣象臺(tái),湖南 長沙 410119)

青藏高原和華東地區(qū)背景誤差協(xié)方差特征的對比研究

陳耀登1,2,趙幸1,2,閔錦忠1,2,范水勇3,王元兵1,2,曾臘梅4

(1.南京信息工程大學(xué) 氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心,江蘇 南京 210044;

2.氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室(南京信息工程大學(xué)),江蘇 南京 210044;

3.中國氣象局 北京城市氣象研究所,北京 100089;4.湖南省氣象臺(tái),湖南 長沙 410119)

摘要：背景誤差協(xié)方差特征與區(qū)域的天氣氣候特征密切相關(guān)。為了更好地理解中國華東地區(qū)和青藏高原地區(qū)的背景誤差協(xié)方差特征,利用夏季一個(gè)月的模擬結(jié)果,以最新的多元變量相關(guān)的背景誤差協(xié)方差模型為基礎(chǔ),通過提取隱含背景誤差協(xié)方差中的變量相關(guān)系數(shù)、特征值、特征向量和特征長度尺度等,對這兩個(gè)區(qū)域的背景誤差協(xié)方差特征進(jìn)行比較和分析。結(jié)果表明,相對于華東地區(qū),青藏高原地區(qū)變量之間的影響關(guān)系更顯著,背景場的誤差更大,大氣特征具有更強(qiáng)的局地性。對青藏高原地區(qū)資料同化而言,觀測資料占有更大的權(quán)重和更小的影響范圍,對青藏高原地區(qū)觀測資料提出了更高的要求。

關(guān)鍵詞：資料同化;背景誤差協(xié)方差;青藏高原;華東地區(qū)

中圖分類號：

文章編號：1674-7097(2015)05-0650-08P456.7

文獻(xiàn)標(biāo)志碼：碼：A

doi：10.13878/j.cnki.dqkxxb.20141021002

Abstract:The characteristics of background error covariances have a close relationship with the local meteorological features.In order to get a better understanding for the characteristics of background error covariances in East China and the Tibetan Plateau,the major parameters of background error covariances,such as variable correlation coefficients,eigenvalues,eigenvectors and characteristic length scales,were compared and analyzed using one month series of numerical forecast fields in summer based on the latest model of multivariate background error covariances.Results show that the relationships between different variables in the Tibetan Plateau are more notable than those in East China.In the Tibetan Plateau,the deviation of background error is more obvious and the local property of atmospheric characteristics are stronger.In terms of data assimilation in the Tibetan Plateau,the observations occupy more weight and smaller influencing area,which puts forward higher requirements on observations in the Tibetan Plateau.

收稿日期：2013-03-01；改回日期：2014-12-31

基金項(xiàng)目:國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2015CB452804);國家自然科學(xué)基金資助項(xiàng)目(41175056);上海臺(tái)風(fēng)研究基金課題(2013ST01);氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室(南京信息工程大學(xué))開放課題(KLME1407);浙江省科技計(jì)劃項(xiàng)目(2014C33056);浙江省氣象科技計(jì)劃項(xiàng)目(2013ZD01)

通信作者：許孌,博士,工程師,研究方向?yàn)榈顷懪_(tái)風(fēng)數(shù)值模擬和動(dòng)力診斷,xuluannew@126.com.

Comparative analysis of characteristics for background

error covariances in Tibetan Plateau and East China

CHEN Yao-deng1,2,ZHAO Xing1,2,MIN Jin-zhong1,2,FAN Shui-yong3,

WANG Yuan-bing1,2,ZENG La-mei4

(1.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,NUIST,Nanjing 210044,China;

2.Key Laboratory of Meteorological Disaster(NUIST),Ministry of Education,Nanjing 210044,China;

3.Institue of Urban Meteorology of CMA,Beijing 100089,China;4.Hunan Meteorological Observatory,Changsha 410119,China)

Key words:data assimilation;background error covariance;Tibetan Plateau;East China

0引言

利用資料同化的方法提取觀測資料的有效信息是改進(jìn)數(shù)值預(yù)報(bào)初始場進(jìn)而提高數(shù)值預(yù)報(bào)水平的一個(gè)重要手段(張衛(wèi)民等,2005)。目前資料同化方法有很多,主要有:變分法(三維變分和四維變分)(官元紅等,2009)、Kalman濾波(擴(kuò)展Kalman濾波和集合Kalman濾波)(閔錦忠等,2013)方法以及混合(Hybrid)同化方法等(陳耀登等,2014a;閔錦忠等,2015)。不論對于哪種資料同化方法,背景誤差協(xié)方差信息都起著關(guān)鍵作用,背景誤差協(xié)方差的好壞直接影響著同化系統(tǒng)的性能(Fisher,2003)。因此如何合理估計(jì)與構(gòu)造背景誤差協(xié)方差以及其特征的研究一直以來都是資料同化的關(guān)鍵工作和重點(diǎn)研究的問題(Berre,2000;邱曉濱,2011;王瑞春等,2012;趙延來等,2013;陳耀登等,2014b)。

在實(shí)際操作中,背景誤差協(xié)方差矩陣的計(jì)算存在兩個(gè)主要問題:首先,“真實(shí)”大氣狀態(tài)未知,在實(shí)際問題中如何準(zhǔn)確計(jì)算背景場誤差是一項(xiàng)難度很大的工作(邱崇踐,2001)。為估算出合理的背景場誤差,研究人員采用了不少方法,如更新矢量法(龔建東和趙剛,2006),NMC法(National Meteorology Cent,亦稱NCEP法)(Parrish and Derber,1992),Ensemble法(Evensen,2003)等。其次,超大規(guī)模的背景誤差協(xié)方差信息在同化系統(tǒng)中進(jìn)行直接表示和運(yùn)算都具有較大難度 景誤差協(xié)方差矩陣,研究人員提出了控制變量轉(zhuǎn)換法(Control Variable Transforms,簡寫為CVT)(Derber and Bouttier,1999)?？刂谱兞哭D(zhuǎn)換(Bannister,2008)。為構(gòu)造在同化系統(tǒng)中既可以方便操作又較為真實(shí)可靠的背通常包括物理變換、水平變換和垂直變換(張華等,2004;莊照榮等,2006)。控制變量轉(zhuǎn)換將背景誤差協(xié)方差矩陣隱含在控制變量轉(zhuǎn)換算子中,不再需要直接表示。因此控制變量轉(zhuǎn)換算子也就包含了該區(qū)域的背景誤差協(xié)方差特征。

近年來,國內(nèi)外開展了不同區(qū)域、不同天氣氣候情況下的背景誤差協(xié)方差特征的研究(Dance,2004;范水勇等,2006;曹小群等,2008;劉磊等,2009;Michel and Thomas,2010;王曼等,2011),他們的研究表明背景誤差協(xié)方差特征與區(qū)域的天氣氣候特征密切相關(guān)。青藏高原由于其復(fù)雜的地形和獨(dú)特的氣候特點(diǎn)而被稱為地球“第三極”(丁一匯和張莉,2008),對中國、亞洲甚至全球的大氣環(huán)流都有非常重要的影響(李永華等,2011)。而同處類似中緯度區(qū)域的我國華東江淮地區(qū)春夏交替時(shí)的梅雨,是我國夏季降水的重要組成部分(周曾奎,1996),梅雨期降水持續(xù)時(shí)間的長短和雨量的多寡與江淮地區(qū)的旱澇災(zāi)害和社會(huì)經(jīng)濟(jì)密切相關(guān)(丁一匯等,2007)。為對我國華東江淮地區(qū)和青藏高原地區(qū)背景誤差協(xié)方差特征有更好的理解,本文以WRFDA(Data Assimilation(DA) system for the Weather Research and Forecasting model)中多元變量相關(guān)的背景誤差協(xié)方差計(jì)算模型為研究基礎(chǔ)(Chen et al.,2013),該協(xié)方差計(jì)算模型補(bǔ)充建立了散度風(fēng)場與溫度場、表面氣壓場的相關(guān)關(guān)系,也補(bǔ)充建立了濕度場與風(fēng)場、溫度場、表面氣壓場的相關(guān)關(guān)系。研究通過提取背景誤差協(xié)方差構(gòu)造過程中,隱含背景誤差協(xié)方差信息的控制變量轉(zhuǎn)換算子的相關(guān)信息,來對這兩個(gè)區(qū)域的背景誤差協(xié)方差特征進(jìn)行比較和分析。

1背景誤差協(xié)方差矩陣的計(jì)算

如前文所述,在變分同化系統(tǒng)中B矩陣的直接描述存在著困難,目前國內(nèi)外多數(shù)資料同化系統(tǒng)采用控制變量轉(zhuǎn)換的方法來對B矩陣進(jìn)行描述?？刂谱兞哭D(zhuǎn)換算子U滿足關(guān)系式:

B=UUT。

(1)

U包括物理變換Up、垂直變換Uv和水平變換Uh。則控制變量轉(zhuǎn)換算子U和B矩陣分別表示為:

U=UpUvUh。

(2)

(3)

目前在WRFDA也采用控制變量轉(zhuǎn)換的方法來描述背景場誤差協(xié)方差信息(Baker et al.,2012):經(jīng)過控制變量轉(zhuǎn)換,背景誤差協(xié)方差矩陣的特征隱含在了物理變換的回歸系數(shù)、水平變換的特征長度尺度、垂直變換的特征值和特征向量中。

2試驗(yàn)方案介紹

研究區(qū)域分別為華東地區(qū)和青藏高原地區(qū)(圖1),兩個(gè)區(qū)域的基本設(shè)置一致:網(wǎng)格點(diǎn)為150×120,垂直層39層,頂層氣壓為50 hPa,水平分辨率12 km。預(yù)報(bào)模式使用WRF(ARW)V3.5.1版本,積云對流參數(shù)化方案為淺對流Kain-Fritcsh方法,微物理方案為WSM6方法,邊界層方案為YSU邊界層方法,Dudhia短波輻射方法和RRTM長波輻射方法。

利用NCEP(National Centers for Environmental Prediction)的FNL(Final Operational Global Analysis)資料,從2009年6月20日00時(shí)—7月20日00時(shí)(世界時(shí),下同),每天分別從00時(shí)和12時(shí)進(jìn)行冷啟動(dòng)做12和24 h預(yù)報(bào),連續(xù)進(jìn)行一個(gè)月,采用NMC方法,以模式同一時(shí)刻預(yù)報(bào)時(shí)效分別為12和24 h預(yù)報(bào)的結(jié)果之差作為預(yù)報(bào)誤差,利用WRFDA中的“gen-be”模塊計(jì)算多元變量背景誤差協(xié)方差。在WRFDA系統(tǒng)中,多元變量相關(guān)的背景誤差協(xié)方差模型(Chen et al.,2013)的選項(xiàng)為CV_OPTION等于6(簡稱CV6)。CV6的控制變量為:流函數(shù)、非平衡速度勢、非平衡地面氣壓、非平衡溫度和非平衡相對濕度。區(qū)別于傳統(tǒng)控制變量選項(xiàng)(CV_OPTION等于5),CV6在溫度場和地面氣壓場非平衡部分的計(jì)算中,除考慮二者與流函數(shù)的回歸統(tǒng)計(jì)關(guān)系,增加計(jì)算了溫度和地面氣壓場與非平衡速度勢的回歸分析項(xiàng);在相對濕度非平衡部分計(jì)算中,完整考慮了相對濕度與其他所有控制變量的相關(guān)關(guān)系,利用回歸分析的方法,建立相對濕度與其他所有控制變量的統(tǒng)計(jì)平衡關(guān)系。

圖1　研究區(qū)域設(shè)置　　a.華東地區(qū);b.青藏高原地區(qū)Fig.1　Setting of studied regions　　a.East China;b.the Tibetan Plateau

3背景誤差協(xié)方差特征分析

3.1　平衡特征

物理變換是為了消除模式變量之間的相關(guān)性而引入,物理變換通過變量間的回歸關(guān)系,將模式變量分為平衡部分和非平衡部分,平衡部分表示變量之間的相關(guān)影響部分,留下的非平衡部分則作為同化系統(tǒng)的分析變量,具體定義見(Chen et al.,2013)。圖2顯示了各控制變量的平衡部分及各變量對平衡場貢獻(xiàn)比例的垂直分布。

流函數(shù)代表了風(fēng)場中有旋無輻散部分,非平衡速度勢則代表風(fēng)場中無旋有輻散部分。從圖2a、b中可以看出,在兩個(gè)區(qū)域,流函數(shù)對速度勢的貢獻(xiàn)都很小,且在近地面和對流層中層,流函數(shù)的貢獻(xiàn)出現(xiàn)極大值,但也沒超過10%。在兩個(gè)區(qū)域相對濕度的平衡場中(圖2c、d),可以發(fā)現(xiàn)在模式的下層(13層以下)平衡場的比例較大,說明在模式中下層各個(gè)變量對濕度場有較大影響。兩個(gè)區(qū)域各個(gè)變量對相對濕度平衡部分的作用,非平衡溫度對平衡場的貢獻(xiàn)最大,接著非平衡速度勢和流函數(shù),非平衡表面氣壓影響最小。同時(shí),也能看出在高原地區(qū)各控制變量對相對濕度的影響比在華東地區(qū)略大,說明在高原地區(qū)各變量對濕度場有更大的影響。

在圖2e、f溫度的平衡場中,兩個(gè)區(qū)域非平衡速度勢的貢獻(xiàn)都比流函數(shù)的影響大,說明無旋的輻散風(fēng)場對溫度的影響比較大。溫度的平衡場中近地面層和對流層中層貢獻(xiàn)出現(xiàn)極值,說明這在近地面和對流層中部,風(fēng)場對溫度的影響比較明顯。對比圖2e和圖2f,可以發(fā)現(xiàn),在高原地區(qū)非平衡速度勢和流函數(shù)對溫度場的影響明顯比華東地區(qū)大,表明在高原地區(qū)風(fēng)場對溫度場的影響更為顯著。

圖2　華東地區(qū)(a,c,e)和青藏高原地區(qū)(b,d,f)平衡部分的比例及相關(guān)變量的貢獻(xiàn)　　a,b.溫度場;c,d.相對濕度場;e,f.平衡場Fig.2　The proportion of the balanced part and the contributions of the variables in (a,c,e)East China and (b,d,f)the Tibetan Plateau　　a,b.temperature field;c,d.relative humidity field;e,f.balanced field

總體而言,本研究中高原地區(qū)與華東地區(qū)由于在類似緯度,整體的變量平衡關(guān)系有大體類似的特征,但在高原地區(qū)變量之間的影響程度要比華東地區(qū)顯著一些。

3.2　特征值和特征向量

垂直變換通過經(jīng)驗(yàn)正交函數(shù)(EOF)的特征模分解得到特征向量和特征值,從而估計(jì)背景誤差協(xié)方差的垂直分量。圖3為在EOF空間各個(gè)模態(tài)上各控制變量的特征值,可以看出,四個(gè)控制變量隨模式層有類似的變化趨勢,在前幾個(gè)模態(tài),高原地區(qū)和華東地區(qū)的特征值均比較大,但隨模態(tài)數(shù)增加快速減小,在第20模式態(tài)后逐步趨近于0。對比兩個(gè)區(qū)域,在前幾個(gè)模態(tài),高原地區(qū)各個(gè)控制變量特征值均比華東地區(qū)的大,高原地區(qū)的流函數(shù)、非平衡溫度和相對濕度的第一模態(tài)對應(yīng)的特征值甚至是華東地區(qū)的兩倍。由于特征值主要表示誤差量級的大小,前幾個(gè)模態(tài)又代表了誤差的主要特征,所以可以得出高原地區(qū)的背景場的誤差比華東地區(qū)顯著,表明模式在高原地區(qū)的模擬效果比華東平原地區(qū)略差,這主要是由于高原地形復(fù)雜,且觀測站點(diǎn)稀少,導(dǎo)致數(shù)值模式以及資料同化等誤差相對較大。

圖3　各控制變量在垂直模態(tài)上的特征值　　a.流函數(shù);b.非平衡速度勢;c.非平衡溫度;d.非平衡相對濕度Fig.3　Eigenvalues of control variables on vertical modes　　a.stream function;b.unbalanced velocity potential;c.unbalanced temperature;d.unbalanced relative humidity

最大特征值對應(yīng)第一特征向量,代表了背景場誤差的最主要垂直結(jié)構(gòu)特征。圖4為各控制變量最大特征值對應(yīng)的第一特征向量在垂直高度上的分布。由流函數(shù)的第一特征向量,可見高原地區(qū)和華東地區(qū)在20層以下都是正的垂直誤差分量。對于非平衡速度勢的第一特征向量,在27、28層(200 hPa)附近,兩個(gè)地區(qū)均出現(xiàn)極值,由于流函數(shù)和非平衡速度勢反映風(fēng)場的運(yùn)動(dòng),因此極值的出現(xiàn),說明此處風(fēng)場的背景誤差比較大,很有可能與中緯度高空急流有關(guān),同時(shí)也表明在這個(gè)位置對高空急流的模擬能力不足,導(dǎo)致模式計(jì)算得到的背景場誤差較大。非平衡溫度的第一特征向量在兩個(gè)區(qū)域的低層出現(xiàn)極值,高原地區(qū)為極大值,華東地區(qū)為極小值。隨高度增加非平衡溫度特征向量逐步減小,在10層以上,誤差幾乎為零,說明非平衡溫度在垂直方向上傳播時(shí)衰減很快,局地性比較強(qiáng),也表明高層溫度與低層溫度的相關(guān)性很小。相對濕度的第一特征向量隨高度變化,在低層和高層,第一特征向量幾乎為0,是因?yàn)榈蛯虞^小表明在低層水汽場的垂直相關(guān)性很小,而高層小是因?yàn)樗苌?幾乎為零。由于青藏高原高海拔大地形的影響,高原地區(qū)的氣壓層分布與平原地區(qū)本身就具有較大差異,誤差特征差異進(jìn)一步的分析還需要考慮高原大地形的影響。

圖4　各控制變量的第一特征向量的垂直分布　　a.流函數(shù);b.非平衡速度勢;c.非平衡溫度;d.非平衡相對濕度Fig.4　Vertical distributions of the first eigenvectors of control variablesa.stream function;b.unbalanced velocity potential;c.unbalanced temperature;d.unbalanced relative humidity

3.3　水平特征長度尺度

特征長度尺度是水平變換遞歸濾波過程中的重要參數(shù),通過其數(shù)值大小可反映在同化過程中觀測信息的影響范圍。從圖5可以看出,兩個(gè)區(qū)域流函數(shù)和非平衡速度勢的特征長度尺度的數(shù)值大小都要遠(yuǎn)大于非平衡溫度和非平衡相對濕度,說明在同化過程中溫度和相對濕度的觀測影響較風(fēng)場觀測的影響范圍小,也說明溫度和相對濕度的水平尺度較小,局地性強(qiáng)。

圖5　控制變量在各個(gè)垂直模態(tài)上的水平特征長度尺度　　a.流函數(shù);b.非平衡速度勢;c.非平衡溫度;d.非平衡相對濕度Fig.5　Horizontal characteristic length scales of control variables on vertical modesa.stream function;b.unbalanced velocity potential;c.unbalanced temperature;d.unbalanced relative humidity

對比兩個(gè)區(qū)域的特征長度尺度也可以發(fā)現(xiàn),高原地區(qū)各個(gè)控制變量的特征長度尺度均比華東地區(qū)的小。對于流函數(shù)和非平衡速度勢,高原地區(qū)的特征長度尺度略小于華東地區(qū)。而對于非平衡溫度和非平衡相對濕度,高原地區(qū)的特征長度尺度明顯小于華東地區(qū)。表明高原地區(qū)的大氣特征相對于華東地區(qū)具有更小的水平尺度和更強(qiáng)的局地性,尤其對于本身水平尺度就較小、局地性較強(qiáng)的非平衡溫度和非平衡相對濕度,在高原地區(qū)水平尺度表現(xiàn)的更加小,局地性更強(qiáng);也表明在利用該背景場誤差協(xié)方差模型進(jìn)行同化的過程中高原的觀測資料影響范圍比華東地區(qū)影響范圍小,尤其是溫度場觀測和相對濕度場的觀測。

4結(jié)論與討論

背景誤差協(xié)方差特征與區(qū)域的天氣氣候特征密切相關(guān),對不同背景誤差協(xié)方差進(jìn)行描述是非常必要的,本文以WRFDA中最新的多元變量相關(guān)的背景誤差協(xié)方差計(jì)算模型為基礎(chǔ),對比類似緯度我國華東江淮地區(qū)和青藏高原地區(qū)的背景誤差協(xié)方差特征,得出以下結(jié)論:

1)高原地區(qū)與華東地區(qū)由于在類似緯度,整體的變量平衡關(guān)系有大體類似的特征,但由各變量之間的變量平衡關(guān)系得到的平衡部分貢獻(xiàn)在高原地區(qū)要大于華東地區(qū),表明在高原地區(qū)變量之間的相關(guān)性要比華東地區(qū)顯著一些。

2)由垂直變換中的特征值分析表明,高原地區(qū)各個(gè)控制變量前幾個(gè)模態(tài)的特征值均比華東地區(qū)的大,高原地區(qū)的流函數(shù)、非平衡溫度和相對濕度的第一模態(tài)對應(yīng)的特征值甚至是華東地區(qū)的兩倍,說明高原地區(qū)的背景場的垂直誤差比華東地區(qū)顯著。這一結(jié)果表明,在資料同化過程中,高原地區(qū)的觀測資料將比在華東地區(qū)占有更大的權(quán)重,對高原區(qū)域觀測資料提出了更高的要求。

3)對第一特征向量的研究表明,高原地區(qū)和華東地區(qū)在20層以下,流函數(shù)的特征向量均顯示為正的垂直誤差分量。對于非平衡速度勢的第一特征向量,在200 hPa附近,兩個(gè)地區(qū)均出現(xiàn)極值,說明此處風(fēng)場的背景誤差比較大,很有可能與模式對中緯度高空急流模擬能力不足有關(guān)。非平衡溫度在垂直方向上傳播時(shí)衰減很快,局地性比較強(qiáng)。

4)從特征長度尺度看,兩個(gè)區(qū)域流函數(shù)和非平衡速度勢的數(shù)值都遠(yuǎn)大于非平衡溫度和非平衡相對濕度,說明在同化過程中溫度和相對濕度的觀測影響較風(fēng)場觀測的影響范圍小,也說明在兩個(gè)區(qū)域溫度和相對濕度的水平尺度均較小,局地性較強(qiáng)。高原地區(qū)各個(gè)控制變量的特征長度尺度均比華東地區(qū)的小,表明高原地區(qū)的大氣特征相對于華東地區(qū)具有更小的水平尺度和更強(qiáng)的局地性。這與高原地區(qū)地形變化大,天氣特征具有更強(qiáng)的局地性特征是相聯(lián)系的,這一結(jié)果表明,在資料同化過程中,高原地區(qū)的觀測資料將比在華東地區(qū)有更小的影響范圍。

由于計(jì)算量和存儲(chǔ)空間的限制,本文背景誤差協(xié)方差的計(jì)算只使用了夏季一個(gè)月的模擬結(jié)果,因此得到的相關(guān)結(jié)果主要代表了該季節(jié)的背景場誤差特征,其他季節(jié)的背景誤差協(xié)方差特征分析需要利用更多或更長時(shí)間段(如3個(gè)月)的樣本資料進(jìn)行計(jì)算并分析。另外,由于青藏高原高海拔大地形的影響,高原地區(qū)的氣壓層分布與平原地區(qū)本身就具有較大差異,誤差特征差異進(jìn)一步的分析還需要考慮青藏高原大地形的影響。

致謝:美國國家大氣研究中心郭永潤博士對本文給予了幫助和指點(diǎn),謹(jǐn)致謝意!

參考文獻(xiàn)(References)：

Bannister R N.2008.A review of forecast error covariance statistics in atmospheric variational data assimilation.II:Modelling the forecast error covariance statistics[J].Quart J Roy Meteor Soc,134(637):1971-1996.

Barker D,Huang X Y,Liu Z,et al.2012.The weather research and forecasting model’s community variational/ensemble data assimilation system:WRFDA.Bull Amer Meteor Soc,93:831-843.

Berre L.2000.Estimation of synoptic and mesoscale forecast error covariances in a limited-area model[J].Mon Wea Rev,128:644-667.

曹小群,黃思訓(xùn),張衛(wèi)民,等.2008.區(qū)域三維變分同化中背景誤差協(xié)方差的模擬[J].氣象科學(xué),28(1):8-14.Cao Xiaoqun,Huang Sixun,Zhang Weimin.2008.Modeling background error covariance inregional3D-VAR[J].Scientia Meteorologica Sinica,28(1):8-14.(in Chinese).

Chen Yaodeng,Rizvi S R H,Huang Xiangyu,et al.2013.Balance characteristics of multivariate background error covariances and their impact on analyses and forecasts in tropical and Arctic regions[J].Meteor Atmos Phys,121:79-98.

陳耀登,王元兵,閔錦忠,等.2014a.一種控制變量循環(huán)擾動(dòng)和多參數(shù)化方案的混合同化方法[J].大氣科學(xué)學(xué)報(bào),37(6):732-739.Chen Yaodeng,Wang Yuanbing,Min Jinzhong,et al.2014a.A hybrid assimilation method with control variable cycle perturbation and multi-physical parameterization[J].Trans Atmos Sci,37(6):732-739.(in Chinese).

陳耀登,曾臘梅,黃向宇,等.2014b.緯度依賴背景場誤差協(xié)方差及其對臺(tái)風(fēng)同化和預(yù)報(bào)的影響研究[J].熱帶氣象學(xué)報(bào),30(4):654-662.Chen Yaodeng,Zeng Lamei,Huang Xiangyu,et al.2014b.A study of latitude-dependent background error covariance and its impacts for typhoon data assimilation and forcast[J].J Trop Meteor,30(4):654-662.(in Chinese).

Dance S L.2004.Issues in high resolution limited area data assimilation for quantitative precipitation forecasting[J].Physica D,196:1-27.

Derber J C,Bouttier F.1999.A reformulation of the background error covariance in the ECMWF global data assimilation system[J].Tellus A,51:195-221.

丁一匯,張莉.2008.青藏高原與中國其他地區(qū)氣候突變時(shí)間的比較[J].大氣科學(xué),32(4):794-805.Ding Yihui,Zhang Li.2008.Intercomparison of the time for climate abrupt change between the Tibetan Plateau and other regions in China[J].Chinese J Atmos Sci,32(4):794-805.(in Chinese).

丁一匯,柳俊杰,孫穎,等.2007.東亞梅雨系統(tǒng)的天氣—?dú)夂驅(qū)W研究[J].大氣科學(xué),31(6):1082-1101.Ding Yihui,Liu Junjie,Sun Ying,et al.2007.A study of the synoptic-climatology of the Meiyusystem in East Asia[J].Chinese J Atmos Sci,31(6):1082-1101.(in Chinese).

Evensen G.2003.The ensemble Kalmanfilter:Theoretical formulation and practical implementation[J].Ocean Dynamics,53:343-367.

范水勇,張朝林,仲躋芹.2006.MM5三維變分系統(tǒng)在北京地區(qū)冷暖季背景場誤差的對比分析[J].高原氣象,25(5):855-861.Fan Shuiyong,Zhang Chaolin,Zhong Jiqin.2006.Contrast analysis of background error of MM5 3DVARSystem in cold and warm seasons in Beijing[J].Plateau Meteorology,25(5):855-861.(in Chinese).

Fisher M.2003.Background error covariance modeling[C]//The ECMWF Seminar on recent developments in data assimilation for atmosphere and ocean.Reading:ECMWF:45-63.

龔建東,趙剛.2006.全球資料同化中誤差協(xié)方差三維結(jié)構(gòu)的準(zhǔn)確估計(jì)與應(yīng)用Ⅱ:背景誤差協(xié)方差調(diào)整與數(shù)值試驗(yàn)分析[J].氣象學(xué)報(bào),64(6):684-698.Gong Jiandong,Zhao Gang.2006.Accurate estimation and application of 3d error covariance structures in global data assimilation.partII:Background error covariance structureadjustments and numerical experiments[J].Acta Meteor Sinica,64(6):684-698.(in Chinese).

官元紅,周廣慶,陸維松.2009.一個(gè)用于氣候模式的三維變分資料同化系統(tǒng)方案設(shè)計(jì)及檢驗(yàn)[J].大氣科學(xué)學(xué)報(bào),32(3):392-404.Guan Yuanhong,Zhou Guangqing,Lu Weisong.2009.A three-dimensional variation data assimilation system for a climate model(AGCM)-basic sheme and tests[J].Trans Atmos Sci,32(3):392-404.(in Chinese).

李永華,盧楚翰,徐海明,等.2011.夏季青藏高原大氣熱源與西南地區(qū)東部旱澇的關(guān)系[J].大氣科學(xué),35(3):422-434.Li Yonghua,Lu Chuhan,Xu Haiming,et al.2011.Contemporaneous relationships between summer atmospheric heat source over the Tibetan Plateau and drought/flood in eastern Southwest China [J].Chinese J Atmos Sci,35(3):422-434.(in Chinese).

劉磊,費(fèi)建芳,程小平,等.2009.我國東部海區(qū)不同氣候背景條件下背景誤差協(xié)方差的性質(zhì)對比分析[J].海洋預(yù)報(bào),26(4):25-35.Liu Lei,Fei Jianfang,Cheng Xiaoping,et al.The comparative analysis of background error covariance in the eastern sea area under different climatic conditions[J].Marine Forecasts,26(4):25-35.(in Chinese).

Michel Y,Thomas A.2010.Inhomogeneous background error modeling and estimation over Antarctica [J].Mon Wea Rev,138:2229-2252.

閔錦忠,畢坤,陳耀登,等.2013.基于物理約束擾動(dòng)的EnSRF雷達(dá)資料同化[J].大氣科學(xué)學(xué)報(bào),36(2):129-138.Min Jinzhong,Bi Kun,Chen Yaodeng,et al.2013.Assimilation of Doppler radar data with EnSRF based onphysical constraint perturbation[J].Trans Atmos Sci,36(2):129-138.(in Chinese).

閔錦忠,劉盛玉,畢坤,等.2015.基于Hybrid EnSRF-En3DVar的雷達(dá)資料同化研究[J].大氣科學(xué)學(xué)報(bào),38(2):213-221.Min Jinzhong,LiuShengyu,Bi Kun,et al.2015.Study on the assimilation of Doppler radar data using a hybrid EnSRF-En3DVar method[J].Trans Atmos Sci,38(2):213-221.(in Chinese).

Parrish D F,Derber J C.1992.The national meteorological center’s spectral statistical interpolation analysis system[J].Mon Wea Rev,120:1747-1763.

邱崇踐.2001.變分同化中使用背景場時(shí)尺度匹配的數(shù)值研究[J].大氣科學(xué),25(1):103-110.Qiu Chongjian.2001.Maching scales of the retrieved variables with background scale in variational assimilation:Anumerical study[J].Chinese J Atmos Sci,25(1):103-110.(in Chinese).

邱曉濱.2011.流依賴的背景誤差協(xié)方差在資料同化和雷達(dá)風(fēng)場反演中的應(yīng)用[D].蘭州:蘭州大學(xué).Qiu Xiaobin.2011.Usage of flow-dependent background error covariance in data assimilation and radar wind retrieval[D].Lanzhou:Lan Zhou University.(in Chinese).

王曼,李華宏,段旭,等.2011.WRF模式三維變分中背景誤差協(xié)方差估計(jì)[J].應(yīng)用氣象學(xué)報(bào),22(4):482-492.Wang Man,Li Huahong,Duan Xu,et al.2011.Estimating background error covariance in regional 3D-VAR of WRF[J].J Appl Meteor Sci,22(4):482-492.(in Chinese).

王瑞春,龔建東,張林.2012.GRAPES變分同化系統(tǒng)中動(dòng)力平衡約束的統(tǒng)計(jì)求解[J].應(yīng)用氣象學(xué)報(bào),23(2):129-138.Wang Ruichun,Gong Jiandong,Zhang Lin.2012.Statistical estimation of dynamic balance constraints in GRAPES variational data assimilation system[J].J Appl Meteor Sci,23(2):129-138.(in Chinese).

張華,薛紀(jì)善,莊世宇.2004.GRAPES三維變分同化系統(tǒng)的理想試驗(yàn)[J].氣象學(xué)報(bào),62(1):31-41.Zhang Hua,Xue Jisan,Zhuang Shiyu,et al.2004.Experiments of GRAPES three-dimensional variational data assimilation system[J].Acta Meteor Sinica,62(1):31-41.(in Chinese).

張衛(wèi)民,朱小謙,趙軍.2005.氣象資料三維變分同化階段區(qū)域分解并行實(shí)現(xiàn)[J].計(jì)算機(jī)研究與發(fā)展,42(6):1059-1064.Zhang Weimin,Zhu Xiaoqian,Zhao Jun.2005.Implementation of phase domain decomposition parallel algorithm of three-dimensional variational data assimilation[J].Journal of Computer Research and Development,42(6):1059-1064.(in Chinese).

趙延來,黃思訓(xùn),張維峰,等.2013.三維變分同化中多變量平衡約束設(shè)計(jì)[J].大氣科學(xué)學(xué)報(bào),36(3):277-285.ZhaoYanlai,Huang Sixun,Zhang Weifeng,et al.2013.Design of multivariable balance constraints inthree-dimensional variational assimilation[J].Trans Atmos Sci,36(3):277-285.(in Chinese).

周曾奎.1996.1994年江淮地區(qū)持續(xù)高溫干旱的環(huán)流特征[J].氣象,22(7):40-42.Zhou Zengkui.1996.Thecharacteristies of the circulation causing continuous high temperature and drought in the Yangtze and Huaihe River Valley in 1994[J].Meteor Mon,22(7):40-42.(in Chinese).

莊照榮,薛紀(jì)善,朱宗申,等.2006.非線性平衡方案在三維變分同化系統(tǒng)中的應(yīng)用[J].氣象學(xué)報(bào),64(2):137-148.Zhuang Zhaorong,Xue Jishan,Zhu Zongshen,et al.2006.Application of nonlinear balance scheme in three-dimensional variational data assimilation[J].Acta Meteor Sinica,64(2):137-148.(in Chinese).

(責(zé)任編輯：張福穎)

許孌,崔曉鵬,高守亭，等.2015.臺(tái)風(fēng)“鲇魚”(1013)路徑突變過程的成因分析[J].大氣科學(xué)學(xué)報(bào),38(5):658-669.doi：10.13878/j.cnki.dqkxxb.20130301001.

Xu Luan,Cui Xiao-peng,Gao Shou-ting,et al.2015.Cause analysis of sudden track change of Typhoon Megi(1013)[J].Trans Atmos Sci,38(5):658-669.(in Chinese).