冀春曉 趙放 高守亭 劉黎平
1 浙江省氣象科學研究所,杭州 310017
2 浙江省氣象臺,杭州 310017
3 中國科學院大氣物理研究所,北京 100029
4 中國氣象科學研究院,北京 100081
登陸臺風Matsa(麥莎)中尺度擾動特征分析
冀春曉1趙放2高守亭3劉黎平4
1 浙江省氣象科學研究所,杭州 310017
2 浙江省氣象臺,杭州 310017
3 中國科學院大氣物理研究所,北京 100029
4 中國氣象科學研究院,北京 100081
地面中尺度自動站和多普勒雷達資料的分析都表明,臺風Matsa登陸后的低層螺旋云帶中活躍著中尺度氣旋性渦旋系統(tǒng)。本文使用新一代中尺度WRF模式對臺風Matsa登陸后的變化特征進行了數(shù)值模擬,使用四維變分多普勒雷達分析系統(tǒng) (4D-VDRAS)對臺風Matsa多普勒雷達徑向風進行了風場反演。在此基礎(chǔ)上對臺風Matsa登陸后中尺度擾動特性進行了初步探討;對臺風Matsa與其螺旋云帶的中尺度系統(tǒng)之間動能和渦度的相互轉(zhuǎn)換進行了診斷分析。分析表明:(1)數(shù)值模擬和雷達風場反演結(jié)果表明,登陸臺風Matsa的低層螺旋云帶中活躍著中尺度氣旋式渦旋系統(tǒng),與之相伴隨的為較強的中尺度上升區(qū),而且,中尺度垂直上升運動的強弱與雷達對流回波強度成正相關(guān),中尺度垂直上升運動越強,雷達對流回波發(fā)展越旺盛。(2)臺風Matsa與其中尺度系統(tǒng)動能轉(zhuǎn)換的診斷分析說明,低層中尺度系統(tǒng)從臺風Matsa環(huán)流中獲得動能而發(fā)展;Matsa在陸地上長久維持主要是從高層獲得動能。(3)臺風Matsa與其中尺度系統(tǒng)渦度轉(zhuǎn)換的診斷分析說明,低層中尺度系統(tǒng)向Matsa輸送正渦度主要依靠中尺度垂直運動來完成;高層正渦度的轉(zhuǎn)換通過水平輸送和垂直輸送共同來完成。所以,中尺度系統(tǒng)所產(chǎn)生的正渦度源源不斷地向Matsa輸送,使Matsa的氣旋性環(huán)流可以在陸地上長久維持。
登陸臺風 中尺度擾動 維持 動能和渦度轉(zhuǎn)換
隨著氣候的逐漸變暖,高溫熱浪、干旱、強降水等極端氣候事件發(fā)生的頻率增加、強度增大,造成洪澇、干旱等災(zāi)害更為嚴重。近幾年來,影響我國臺風的強度和頻率似乎也有增長的趨勢。例如,2007年8~10月期間,“圣帕”、“韋帕”、“羅莎”3個強臺風接連嚴重影響了浙江省,造成了大范圍的暴雨和大暴雨。登陸臺風所帶來的災(zāi)害與其強度、登陸后維持的時間等有著密切的關(guān)系。臺風強度及強度變化主要受下墊面、環(huán)境場的風垂直切變、上層槽和臺風渦旋自身結(jié)構(gòu)的影響 (Emanuel,1988;DeMaria,1996; Montgomery and Kallenbach,1997;Jones,2000)。環(huán)境因子通過與渦旋的相互作用影響渦旋內(nèi)部結(jié)構(gòu)及物理過程,進而影響臺風的強度變化。環(huán)境場、臺風及渦旋內(nèi)部多尺度之間的相互作用決定了臺風強度變化的復(fù)雜性 (Wang and Wu,2004)。Chen(1998)指出,登陸臺風在以下幾種環(huán)流條件下將維持較長時間:(1)臺風環(huán)流保持一定的水汽供給;(2)臺風環(huán)流中存在活躍的中尺度對流活動;(3)弱冷空氣侵入臺風環(huán)流引起變性;(4)登陸臺風環(huán)流移入一個高空輻散區(qū)之下。陳聯(lián)壽等 (2002)的研究表明,登陸臺風因地形或環(huán)境流場輻合作用常在其外圍形成中尺度渦旋,當這種渦旋與臺風相互作用被吸入內(nèi)區(qū)時,這時臺風渦量將增加,有利于登陸臺風在陸上的維持。李英等 (2005a)的診斷研究發(fā)現(xiàn),長久維持的熱帶氣旋在登陸后仍保持一定強度,并從外界獲得熱量和水汽補充來支持積云對流發(fā)展,而積云對流對登陸熱帶氣旋的長久維持具有正反饋作用。陳永林等 (2008)分析麥莎臺風次天氣尺度系統(tǒng)的非對稱性結(jié)構(gòu)時發(fā)現(xiàn),臺風麥莎登陸后其外圍云系中活躍著中尺度對流系統(tǒng)。段麗和陳聯(lián)壽 (2005)在對熱帶風暴“菲特”的診斷研究中發(fā)現(xiàn),當“菲特”穿越瓊州海峽時,其外圍流場和五指山山脈的輻合極有利于中尺度對流小渦系統(tǒng) (MCS)的生成,而特大暴雨就出現(xiàn)在MCS所在地區(qū)。王繼志和楊元琴(1995)研究了8807號臺風突然增強的原因后指出,臺風中中尺度系統(tǒng)與周圍環(huán)境中大尺度系統(tǒng)的相互作用對臺風的增強十分重要。而熱帶氣旋的加強則可能與其從次天氣尺度運動中獲得動能補充有關(guān)。李英等 (2004)在登陸熱帶氣旋維持的次天氣尺度環(huán)流特征中指出,臺風渦旋的強對流云團或雨帶作為次天氣尺度呈現(xiàn)出的波狀螺旋型特征,這種螺旋波結(jié)構(gòu)有利于熱帶氣旋的維持。維持和加強的登陸熱帶氣旋可以從次天氣尺度環(huán)流中獲得動能補充,所以,判斷一個熱帶氣旋登陸后能否繼續(xù)維持或再加強,除考慮大尺度環(huán)流等因素的影響外,次天氣尺度環(huán)流的作用也是一個不可忽視的因子。于玉斌等 (2008)的數(shù)值模擬試驗表明低層氣旋性渦旋并入臺風環(huán)流是“桑美”近海急劇增強的重要原因。鄧國等 (2005)對臺風數(shù)值模擬中邊界層方案進行敏感性試驗后發(fā)現(xiàn),臺風邊界層通過摩擦混合和輻射等作用與地表產(chǎn)生水汽、熱量和動量的交換,并通過湍流效應(yīng)和積云的夾卷作用將邊界層的影響擴展至整個自由大氣。勵申申等 (1992)的研究表明,登陸臺風的消亡除摩擦作用外,還可能通過支持外圍大暴雨的發(fā)展將能量傳遞出去而使本身消亡。袁金南等 (2005)發(fā)現(xiàn),邊界層摩擦對登陸熱帶氣旋強度的影響非常明顯,其中摩擦是造成登陸熱帶氣旋強度迅速減弱的一個重要因素。李英等(2005b)的研究表明,外界水汽輸送有助于熱帶氣旋雨帶中的強對流活動,使雨量加強。而對流活動過程中水汽凝結(jié)潛熱釋放是熱帶氣旋獲得能量繼續(xù)維持的重要條件。于玉斌等 (2008)在超強臺風“桑美”(2006)近海急劇增強特征及機理分析中指出,“桑美”臺風急劇增強過程中,對流層高層動能的下傳是對流層低層動能補充的重要途徑之一。這些研究均表明,登陸臺風強度的變化與中尺度系統(tǒng)活動存在著非常密切的關(guān)系。
另外,隨著科學技術(shù)的進步,新一代多普勒天氣雷達在對臺風等災(zāi)害天氣的探測中發(fā)揮了越來越重要的作用。多普勒天氣雷達不僅能夠探測臺風的螺旋云帶,而且通過風場反演技術(shù)有可能得到臺風中高時空分辨率的風場結(jié)構(gòu),這對于研究其中系統(tǒng)的動力結(jié)構(gòu)大有益處。Sun et al.(1991)、Sun and Crook(1994,1997,1998)利用四維變分方法,結(jié)合云模式及其伴隨模式,提出了多普勒雷達4DVAR同化反演方法,試驗表明該方法能夠反演風暴云團的熱、動力場和微物理場。隨著變分技術(shù)的引入,使風場反演的可靠性和精度得到了進一步提高,深化了對臺風、颮線等災(zāi)害性天氣的認識(Wu et al.,2000)。周海光和王玉彬 (2005)使用雙多普勒雷達三維風場反演技術(shù) (MUSCAT)對6月29~30日合肥和馬鞍山多普勒雷達探測到的暴雨資料進行了三維風場反演,對暴雨系統(tǒng)的中尺度三維動力結(jié)構(gòu)進行了研究,給出了此次暴雨的三維動力結(jié)構(gòu)。王俊等 (2011)雙多普勒雷達反演的三維風場分析表明,普通多單體雷暴中低層存在水平輻合風場,且雷暴的發(fā)展與水平輻合風場的強度成正比。孫建華等 (2006)的多普勒反演風場分析表明,造成滁州強降雨的系統(tǒng)中存在中尺度輻合線,且中尺度輻合線只出現(xiàn)對流層低層,輻合線隨高度向北傾斜。但利用多普勒雷達探測資料研究臺風系統(tǒng)中的中尺度系統(tǒng)三維動力結(jié)構(gòu)的國內(nèi)還不多見。
2005年9號臺風Matsa登陸后在陸地上持續(xù)維持了近75小時,沿途帶來較大的災(zāi)害,造成了巨大的損失。關(guān)于Matsa登陸后持續(xù)維持的機理冀春曉等 (2007)已經(jīng)作過分析。本文著重分析研究臺風Matsa登陸后維持期間其螺旋云帶中的中尺度擾動特性。利用多普勒雷達探測到的Matsa臺風資料進行了四維變分風場反演,探討了臺風Matsa登陸后維持期間的中尺度結(jié)構(gòu),從而進一步深化對登陸臺風中尺度三維動力結(jié)構(gòu)演變特征的認識。
本文使用美國NCAR等研制開發(fā)的新一代中尺度數(shù)值模式 WRF(V2.1),網(wǎng)格格距為15km,模式區(qū)域中心點為 (29.3°N,120.1°E),格點數(shù)為363×363,垂直方向31層σ坐標。微物理過程選取Lin等 (1983)的方案;陸面過程選取Noah(Chen and Dudhia,2001)方案,積云對流參數(shù)化方案選取Betts-Miller-Janjic(Janji ,1994,2000)方案。初始時間選取2005年8月6日00時 (協(xié)調(diào)世界時,下同),積分時間為72小時,初猜場和邊界值選用NCEP的1°×1°再分析資料,常規(guī)資料為氣象信息綜合分析處理系統(tǒng) (MICAPS)地面和高空資料,采用3DVAR技術(shù)將初始時刻的地面、高空觀測資料同化進入模式。側(cè)邊界為6小時變邊界。對臺風未經(jīng)任何人工處理。
本文使用中國氣象局氣象科學研究院災(zāi)害天氣國家重點實驗室提供的4DVAR同化反演系統(tǒng)。同化模式和變分同化方法由Sun and Crook(1997)建立,這里僅對其作一簡要介紹。同化模式為一個三維云模式,云模式建立在笛卡兒坐標系中,采用滯彈性近似,包含6個預(yù)報方程,它們是3個動量方程、熱力方程、雨水方程和總水方程。數(shù)值模式以無量綱變量形式編程,這樣可以平衡不同變量量級的差異,使得在同化過程中每個變量有相似的權(quán)重從而得到更好的收斂率。4DVAR資料同化的基本思想就是找模式變量的最優(yōu)初始場,使得模式輸出結(jié)果在一定的時間域和空間域上與相應(yīng)的觀測結(jié)果盡可能接近。為此,單部雷達價值函數(shù)J為
其中,求和針對空間區(qū)域σ、同化窗τ而言,ηv和ηz分別是徑向速度和反射率的權(quán)重系數(shù),Vobr和Zob是雷達觀測的徑向速度和反射率。Vr和Z表示模式輸出的徑向速度和反射率。Jb和Jp分別為背景場、補償項函數(shù)。具體算法見Sun and Crook(1997),在此不再贅述。
2005年第9號臺風Matsa于8月5日19:40時在浙江省玉環(huán)縣干江鎮(zhèn)沿海登陸,登陸時中心氣壓950hPa,近中心最大風速12級以上 (45m/s)。登陸后先向西北方向移動,之后Matsa一路北上,強度減弱緩慢,給所經(jīng)之地造成了重大的災(zāi)害。所以,探討臺風Matsa在登陸后強度繼續(xù)維持的中尺度系統(tǒng)特征,對業(yè)務(wù)天氣預(yù)報應(yīng)該具有十分重要的指導(dǎo)意義。
圖1 (a)臺風中心氣壓隨時間演變圖 [豎虛線為兩個階段的分界線 (下同)];(b)臺風路徑圖Fig.1 (a)Temporal evolution of the observed(■)and simulated(▲)sea level pressure(hPa)at the center of typhoon Matsa(The vertical dashed line denotes the dividing point between the two different periods of typhoon and tropical storm,the same below);(b)the track of typhoon Matsa
臺風Matsa登陸后,在陸地上持續(xù)維持了近75小時,基本上可分為兩個階段,第一階段為5日18時~6日18時,臺風中心迅速被填塞,中心氣壓從950hPa快速上升到985hPa,近中心風速從45m/s快速下降到28m/s,但強度仍保持較強;第二階段為6日18時~8日21時,這一階段里臺風強度變化很小,之后Matsa減弱為溫帶氣旋。從圖1a可以看到,模式模擬的臺風中心氣壓在開始的幾個小時里與實際觀測值相差比較大,主要原因可能是由于所采用的NCEP資料水平分辨率還不夠高,對臺風的強度達不到更精確的描述,以及模式自身的spin up等問題,導(dǎo)致模擬的臺風中心氣壓要明顯高于實際值。隨著時間的延續(xù),兩者很快變得一致,模式對Matsa強度緩慢變化過程進行了成功的模擬。從圖1b可以看到,模式對臺風Mat-sa登陸后路徑的模擬也是比較接近實況的,其中包括Matsa登陸后先向西北方向移動,而后轉(zhuǎn)向偏北方向移動。但是,隨著模擬時間的延長,模擬Matsa路徑與實況的誤差增大,這可能是由于模式模擬準確度與模擬時間成反比的緣故。
圖2 2005年8月 (a)6日15:00時和 (b)7日06:00時不同區(qū)域地面中尺度自動站2分鐘平均風向風速圖 (單位:m/s,數(shù)字為風速)▲:地面輻合中心 (下同)Fig.2 Observed 2-minute-mean surface winds(m/s)from AWS for different time and areas:(a)At 1500UTC 6Aug 2005;(b)at 0600 UTC 7Aug 2005.▲:the position of convergence center(the same below)
圖3 2005年8月不同時間和不同地點雷達組合反射率 (單位:dBZ)和地面中尺度自動站2分鐘平均風向風速圖 (單位:m/s):(a)6日15:01時溫州雷達;(b)7日06:01時南京雷達Fig.3 Composite reflectivity of Doppler radar and observed 2-minute-mean surface winds from AWS for deferent time and areas:(a)Wenzhou radar at 1501UTC 6Aug 2005;(b)Nanjing radar at 0601UTC 7Aug 2005
我國地面中尺度自動站與多普勒雷達的布網(wǎng)為我們分析中尺度對流系統(tǒng)提供了寶貴的觀測資料。圖2給出了不同時間、不同區(qū)域的地面風向風速實測分布圖。圖2a為2005年8月6日15時地面中尺度自動站2分鐘的風向風速分布圖,圖中數(shù)字表示風速值,此時 Matsa中心位置位于 (30.5°N,119.3°E)。從圖2a可以看到Matsa的氣旋式環(huán)流很明顯,強風速區(qū)位于Matsa中心的東北方。將圖2a中黑方框區(qū)域放大后 (圖3a)可以發(fā)現(xiàn),在Matsa的氣旋式環(huán)流中,即在浙江省的臺州地區(qū)存在中尺度氣旋性輻合區(qū),輻合中心 (圖3a中▲,下同)的極小值達到了-3.0×10-4s-1左右。在6日15:01時溫州多普勒雷達組合反射率圖 (圖3a)上可以看到,中尺度氣旋性輻合區(qū)域中存在著結(jié)構(gòu)密實的強對流回波,30dBZ以上的強回波呈帶狀分布,最強回波達45dBZ左右,25dBZ回波頂高達到500hPa左右 (圖5a),但輻合中心與強回波中心并不重合。2005年8月7日06時 (圖2b)臺風Matsa的中心位置已經(jīng)到達了118.3°E、31.7°N,雖然仍維持著氣旋式環(huán)流,但其強度明顯減弱。與前一時刻相同的是如圖3b(圖2b中黑方框區(qū)域)所示,在Matsa的氣旋式環(huán)流中仍然存在著中尺度氣旋性輻合區(qū),但輻合強度比6日15時要有所減弱,中心的最小值為-2.0×10-4s-1左右。相應(yīng)地在7日06:01時南京的多普勒雷達組合反射率圖(圖3b)上,存在發(fā)展旺盛的對流回波與地面中尺度氣旋性輻合區(qū)域相對應(yīng),同樣輻合中心與強回波中心不重合,這些對流回波的組合反射率最大值均達到了50dBZ左右,而25dBZ回波頂高達到350hPa左右 (圖5b),與前一時刻相比,中尺度對流回波的發(fā)展更加旺盛,可見這些中尺度對流回波都是比較深厚的對流系統(tǒng)。以上觀測事實均說明Matsa螺旋云帶中存在著比較深厚的中尺度對流系統(tǒng)。
圖4 2005年8月6日15:00時 (a、b)和7日06:00時 (c、d)950hPa(a、c)和700hPa(b、d)中尺度水平流場Fig.4 Mesoscale streamline fields at(a,c)950hPa and(b,d)700hPa at(a,b)1500UTC 6Aug 2005and(c,d)0600UTC 7Aug 2005
為了分析中尺度擾動場,濾波方案采用GrADS九點平滑算子,對積分流場進行尺度分離,得到中尺度擾動場。與觀測事實相一致,本文給出了6日15時和7日06時兩個時刻相應(yīng)區(qū)域的模擬結(jié)果。圖4為Matsa在這兩個時刻、950hPa和700hPa經(jīng)過濾波后的中尺度擾動流場。分析圖4可知,無論Matsa在第一階段還是第二階段,低層950hPa上 (圖4a、c)Matsa螺旋云帶中都活躍著中尺度氣旋式渦旋系統(tǒng)。處于第一階段時,低層950hPa上 (圖4a)Matsa螺旋云帶的中尺度氣旋式渦旋的位置與地面中尺度自動站 (圖3a)觀測到的中尺度輻合中心基本吻合;到第二階段時,Matsa螺旋云帶的中尺度氣旋式渦旋的位置(圖4b)與地面中尺度自動站 (圖3b)觀測到的中尺度輻合中心存在明細的偏離。隨著高度的增加,中尺度氣旋式渦旋逐漸轉(zhuǎn)變減弱,到700hPa高度(圖4b、d),中尺度氣旋式渦旋系統(tǒng)已經(jīng)消失并轉(zhuǎn)變?yōu)檩椛^(qū),說明中尺度氣旋式渦旋系統(tǒng)只在低層存在。
圖5給出了沿中尺度氣旋式渦旋中心 (圖4中▲)的垂直剖面圖,陰影區(qū)為對應(yīng)時刻的雷達基本反射率。圖5a為6日15時沿 (29.2°N,121.7°E)中尺度氣旋式渦旋中心的緯向垂直剖面圖,從圖5a可以看到,在中尺度氣旋式渦旋中心的上空,存在明顯的中尺度垂直上升區(qū),中尺度垂直上升速度的極大值位于750hPa附近,量值為0.6×10-1m/s。與實際雷達回波相比較可知,低層的上升運動區(qū)對應(yīng)著較強的雷達對流回波,對流回波的極大值約35dBZ,25dBZ回波頂高達到650hPa左右,但最強的雷達回波與上升運動中心不重合。在沿(29.2°N,121.7°E)中尺度氣旋式渦旋中心的經(jīng)向垂直剖面圖 (圖略)上,中尺度氣旋式渦旋中心上空同樣存在上升運動,上升速度的極大值(700hPa附近)出現(xiàn)在渦旋中心的南側(cè),與較強的雷達對流回波存在偏離。圖5b為7日06時沿位于 (33.3°N,119.7°E)中尺度氣旋式渦旋中心的經(jīng)向垂直剖面圖,同樣,中尺度氣旋式渦旋中心的上空為明顯的中尺度垂直上升區(qū),中尺度垂直上升速度的極大值中心位于700hPa附近,為4×10-1m/s。相應(yīng)的對流回波極大值約40dBZ,25dBZ回波頂高達到400hPa左右,上升區(qū)的兩側(cè)為下沉運動區(qū),同樣最強的雷達回波位于上升運動中心的南側(cè)。從(33.3°N,119.7°E)中尺度氣旋式渦旋中心的緯向垂直剖面圖 (圖略)上可知,中尺度氣旋式渦旋中心上空同樣存在上升運動,但上升速度的極大值(750hPa附近)出現(xiàn)在渦旋中心的西側(cè),同樣與較強的雷達對流回波存在偏離。以上結(jié)果說明,最強的雷達回波與最強的中尺度上升運動中心是不重合的。對比圖5a和圖5b可以發(fā)現(xiàn),7日06時中尺度擾動、對流回波極大值都比6日15時的要強很多,最強的雷達回波與中尺度上升運動中心并不完全重合,這一點在雷達風場反演結(jié)果中也得到了證實。以上分析結(jié)果說明中尺度氣旋式渦旋系統(tǒng)對于強的雷達對流回波的形成是至關(guān)重要的。而且,中尺度垂直上升運動的強弱與雷達對流回波強度成正相關(guān),中尺度垂直上升運動越強,雷達對流回波發(fā)展越旺盛。
從冀春曉等 (2007)的分析可以知道,6日15時為Matsa第一階段區(qū)域平均垂直速度出現(xiàn)極大值和低層輻合為最強的時刻,7日06時為Matsa第二階段區(qū)域平均垂直速度出現(xiàn)次極大值和低層輻合為較強的時刻,這說明中尺度對流系統(tǒng)的產(chǎn)生,不僅有利于低層輻合的加強,而且也使臺風垂直運動得到加強。以上結(jié)果說明,Matsa能夠在陸地上持續(xù)維持,與其螺旋云帶中活躍著中尺度氣旋式渦旋系統(tǒng)有非常密切的聯(lián)系,這些中尺度對流系統(tǒng)產(chǎn)生劇烈的輻合和垂直上升運動,把低層暖濕空氣向上輸送,發(fā)生凝結(jié)并釋放潛熱,對Matsa強度的維持起到了正反饋作用。
多普勒雷達資料4DVAR同化反演技術(shù)在理論上是先進的,Sun and Crook(1994)進行的一些實際運用試驗也被證明是有效的。但是到目前為止,多普勒雷達資料的4DVAR同化仍然處在試驗研究階段,特別是在國內(nèi)類似研究還剛剛起步。本文使用中國氣象局氣象科學研究院災(zāi)害天氣國家重點實驗室提供的4DVAR同化反演系統(tǒng),對臺風Matsa登陸后的對流云帶進行了局部(圖3b黑方塊區(qū)域)反演 (如圖6所示)。
從南京多普勒雷達組合反射率圖 (圖3b)中可以看到,臺風Matsa登陸后直到2005年8月7日06:01時,其外圍云帶依然很強,外圍云帶中鑲嵌著多個大于45dBZ的強回波單體。從相應(yīng)時刻仰角為1.37°的雷達徑向風場 (圖略)可以看到,在南京雷達站北側(cè)(約25km)存在一個逆風區(qū)成逆時針環(huán)流分布,說明該處存在一個中尺度氣旋式環(huán)流;在仰角為2.34°的雷達徑向風場上,該中尺度氣旋式環(huán)流的面積比低仰角的要小,且中心位置向偏北方傾斜。當仰角大于3.5°后,該逆風區(qū)變得模糊或不存在,說明這個中尺度氣旋式環(huán)流只在低層存在。圖6給出了使用2005年8月7日06:01時南京多普勒雷達探測資料反演(圖3b黑方塊區(qū)域)的不同高度的流場與基本反射率分布圖。從圖6可以看到,此時對流回波單體發(fā)展很旺盛,大于30 dBZ的強回波頂高伸展到2.5km以上。在臺風Matsa的強回波云帶中存在明顯的中尺度氣旋式渦旋,它位于3.5km以下,且中尺度氣旋式渦旋中心位置隨著高度的升高向偏北方發(fā)生傾斜,氣旋性渦旋面積隨著高度的增加而縮小。2km以下為明顯的中尺度氣旋式渦旋中心,2.5~3km之間減弱為輻合線或輻合中心,到3.5km時輻合中心基本消失。與雷達徑向風場的分析相一致。同時,從圖6可以發(fā)現(xiàn)中尺度氣旋式渦旋中心位于最強回波(≥40dBZ)的一側(cè),兩者并不完全重合,這與圖5分析的結(jié)論相一致。
圖5 2005年8月基本反射率 (彩色,單位:dBZ)、中尺度擾動垂直速度 (等值線,單位:10-1 m/s)和中尺度擾動合成流場 (矢量)的垂直剖面圖:(a)6日15:00時,矢量為 (u*,w*×10);(b)7日06:00時,矢量為 (v*,w*×10)Fig.5 The vertical cross sections of base reflectivity(shaded),mesoscale disturbance vertical velocity(isoline,units:10-1 m/s),and composite mesoscale disturbance wind fields(a)(u*,w*×10)and(b)(v*,w*×10)at(a)1500UTC 6Aug 2005and(b)0600UTC 7Aug 2005
圖6 2005年8月7日06:01時不同高度基本反射率 (彩色,單位:dBZ)和反演的水平流場分布圖:(a)1.0km;(b)1.5km;(c)2.0km;(d)2.5km;(e)3.0km;(f)3.5km。坐標數(shù)字為格點數(shù) (下同)Fig.6 Base reflectivity(shaded)and streamline fields retrieved at(a)1.0km,(b)1.5km,(c)2.0km,(d)2.5km,(e)3.0km,and(f)3.5km at 0601UTC on 7Aug 2005
上述分析表明,臺風Matsa的維持與其螺旋云帶的中尺度系統(tǒng)活動密切相關(guān)。下面從兩種尺度間動能及渦度轉(zhuǎn)換出發(fā),討論Matsa與其螺旋云帶的中尺度系統(tǒng)的相互作用。
從陳受鈞和謝安 (1981)分析可知,天氣尺度與中尺度間相互轉(zhuǎn)換項 (V*·I*)包括水平運動和垂直運動的貢獻,分別為:
其中,(V*·I*)表示中尺度系統(tǒng)與天氣尺度系統(tǒng)動能的轉(zhuǎn)換。正值表示中尺度運動向天氣尺度系統(tǒng)提供動能,負值表示中尺度系統(tǒng)從天氣尺度系統(tǒng)獲得動能。
垂直渦度方程中天氣尺度和中尺度系統(tǒng)轉(zhuǎn)換項Iζ的水平運動項和垂直運動項,分別為:
其中,Iζ>0表示兩種尺度相互作用的結(jié)果是中尺度系統(tǒng)給天氣尺度系統(tǒng)提供正渦度;Iζ<0表示天氣尺度系統(tǒng)給中尺度系統(tǒng)提供正渦度。
為了詳細討論臺風Matsa在陸地上持續(xù)維持期間與其環(huán)流中的中尺度系統(tǒng)之間的相互作用,我們選用移動坐標,取臺風為中心的 (15×15)經(jīng)緯度范圍為臺風區(qū)域,對式 (1)、(2)計算區(qū)域平均動能轉(zhuǎn)換項 [V*·I*]([]表示取區(qū)域平均),帶“*”的值由原始值減去區(qū)域平均值得到。用式(3)、(4)計算區(qū)域平均渦度轉(zhuǎn)換項 [Iζ],如圖7所示。對 [V*·I*]和 [Iζ]從地面至100hPa進行垂直積分,得到整層的 [V*·I*]和 [Iζ]值 (見圖8)。
圖7a為天氣尺度與中尺度之間區(qū)域平均動能的水平運動轉(zhuǎn)換項[V*·I*]h,圖7b為垂直運動轉(zhuǎn)換項[V*·I*]v,圖7c為水平運動和垂直運動轉(zhuǎn)換兩者之和 [V*·I*]。對比圖7a和圖7b可以看到,[V*·I*]h比 [V*·I*]v基本屬于同一個量級,最大的轉(zhuǎn)換都出現(xiàn)在高層,但兩者對動能轉(zhuǎn)換卻起著相反的作用。在第一階段 (6日00~18時),[V*·I*]h的作用在低層是Matsa從中尺度運動中獲得動能,極大值位于900hPa附近;高層是Matsa向中尺度系統(tǒng)提供動能。[V*·I*]v的作用與前者剛好相反。低層Matsa向中尺度系統(tǒng)提供動能,極大值也位于900hPa附近;高層是Matsa從中尺度系統(tǒng)中獲得動能。兩者共同的作用 [V*·I*](圖7c)表明,在800hPa以下為負值,Matsa向中尺度系統(tǒng)提供動能,極小值在900hPa附近;800 hPa以上為正值,Matsa從中尺度運動中獲得動能,極大值在650hPa附近。但 [V*·I*]負值的絕對值幾乎是其正值的兩倍,說明近地層以動能的垂直轉(zhuǎn)換為主,Matsa向中尺度系統(tǒng)提供動能。整層積分結(jié)果 (圖8a)表明,[V*·I*]均為正值,說明Matsa在這一階段從中尺度系統(tǒng)獲得了動能。所以,雖然Matsa登陸后的強度迅速減弱,但仍然可以維持熱帶風暴的強度。
到第二階段 (6日18時~8日21時),高層動能轉(zhuǎn)換明顯增加,而低層動能轉(zhuǎn)換卻有所減弱。[V*·I*]h(圖7a)的作用表現(xiàn)為兩種情況:6日18時~8日12時期間,正負分布成雙層結(jié)構(gòu),500hPa以下、250~300hPa之間為正值區(qū),300~500hPa、250hPa以上為負值區(qū),而且高空負值的絕對值要明顯大于低層的正值,說明高層動能轉(zhuǎn)換明顯要大于低層轉(zhuǎn)換。而 [V*·I*]v(圖7b)的正負結(jié)構(gòu)分布與 [V*·I*]h剛好相反,負值主要分布在600hPa以下,正值主要分布在250hPa以上,中間的垂直轉(zhuǎn)換相對較弱。8日12時~9日00時期間,[V*·I*]h與 [V*·I*]v結(jié)構(gòu)分布是一致的,即表現(xiàn)為正負分布單層結(jié)構(gòu)。[V*·I*]h在300hPa以下Matsa向中尺度系統(tǒng)提供動能,300hPa以上Matsa從中尺度系統(tǒng)獲得動能。[V*·I*]v在700hPa以下Matsa向中尺度系統(tǒng)提供動能,700hPa以上Matsa從中尺度系統(tǒng)獲得動能。兩項的綜合作用 [V*·I*](圖7c)的結(jié)果顯示,動能的
圖7 區(qū)域平均動能轉(zhuǎn)換項 (a-c,單位:10-4 W·kg-1·m-2)和區(qū)域平均渦度轉(zhuǎn)換項 (d-f,單位:10-10 s-2/m2)時間演變圖:(a)
Fig.7 Temporal evolution of(a-c)regionally averaged kinetic energy conversion(10-4W·kg-1·m-2)and(d-f)regionally averaged vorticity conversion(10-10s-2/m2):轉(zhuǎn)換隨高度的變化正負交替出現(xiàn),主要以高層的動能轉(zhuǎn)換為主。[V*·I*]低層為負,說明 Matsa失去動能,中尺度系統(tǒng)從Matsa獲得動能,這一結(jié)論說明低層Matsa環(huán)流中的中尺度系統(tǒng)得到明顯發(fā)展。從整層積分的結(jié)果 (圖8a)來看,6日18時~7日21時期間,[V*·I*]整層積分均為正值,即中尺度系統(tǒng)向Matsa提供動能。從7日21時~9日00時,[V*·I*]整層積分均為負值,即 Matsa向中尺度系統(tǒng)提供動能。以上分析結(jié)果表明,[V*·I*]低層為負,說明中尺度系統(tǒng)從臺風Matsa環(huán)流中獲得動能而發(fā)展,這應(yīng)該是Matsa環(huán)流中低層活躍著中尺度系統(tǒng)的原因之一;[V*·I*]的正值主要出現(xiàn)在高層,即高層的中尺度擾動為Matsa提供動能。所以,Matsa在陸地上長久維持主要是從高層獲得動能。
圖8 區(qū)域平均動能轉(zhuǎn)換項 (a,單位:W·kg-1·m-2·hPa-1)和區(qū)域平均渦度轉(zhuǎn)換項 (b,單位:10-6 s-2·m-2·hPa-1)從地面至200hPa垂直積分時間演變圖Fig.8 Temporal evolution of(a)regionally averaged kinetic energy conversion(W·kg-1·m-2·hPa-1)and(b)regionally averaged vorticity conversion(10-6 s-2·m-2·hPa-1)vertically integrated from surface to 100hPa
[Ihζ]是區(qū)域平均渦度水平運動轉(zhuǎn)換項(圖7d)。它的作用主要表現(xiàn)在中高層,而且在第二階段7日12時和8日00時,200hPa附近分別存在兩個極大值中心。說明中尺度系統(tǒng)向 Matsa提供正渦度。為區(qū)域平均渦度垂直運動轉(zhuǎn)換項,如圖7e所示。在低層和高層,而且低層正值要明顯大于高層正值,說明在低層中尺度系統(tǒng)主要通過垂直運動向Matsa提供正渦度。的負值區(qū)主要位于中高層,Matsa向中尺度系統(tǒng)提供正渦度。圖7f為水平運動和垂直運動轉(zhuǎn)換兩者之和 [Iζ]。[Iζ]的結(jié)果表明 (圖7f),800hPa以下和300hPa以上,[Iζ]>0,而在800~300hPa之間,幾乎 [Iζ]<0,說明Matsa所需要的正渦度主要來自于低層和高層的中尺度擾動。[Iζ]整層積分結(jié)果如圖8b所示,在第一階段,[Iζ]整層積分幾乎都為負,Matsa失去正渦度,而中尺度系統(tǒng)得到正渦度,實況表現(xiàn)為Matsa氣旋性環(huán)流強度的迅速減弱和中尺度系統(tǒng)的發(fā)展。在第二階段,除8月8日15時~9日00時為負值外,其他時間 [Iζ]整層積分結(jié)果均為正,中尺度系統(tǒng)向Matsa提供正渦度。
上面的分析結(jié)果說明,在低層,臺風Matsa通過水平運動為其螺旋云帶中的中尺度系統(tǒng)提供動能,激發(fā)中尺度系統(tǒng)的發(fā)生發(fā)展,產(chǎn)生中尺度氣旋性渦旋;而中尺度系統(tǒng)通過其垂直運動將正渦度向臺風輸送。所以,中尺度系統(tǒng)所產(chǎn)生的正渦度以及向Matsa的轉(zhuǎn)換,對Matsa在陸地的長久維持起到了非常重要的作用。
本文使用新一代中尺度WRF模式對臺風Matsa登陸后的路徑和強度變化特征進行了較成功的數(shù)值模擬,采用4DVAR風場反演技術(shù)對臺風Matsa多普勒雷達資料進行了風場反演。在此基礎(chǔ)上對臺風Matsa登陸后中尺度擾動特性進行了初步探討,對臺風Matsa與其環(huán)流中的中尺度系統(tǒng)之間相互作用進行了診斷分析。臺風Matsa登陸后強度的持續(xù)維持與其環(huán)流中活躍的中尺度擾動有著非常密切的關(guān)系。
(1)地面中尺度自動站和多普勒雷達資料的分析都表明,臺風Matsa的螺旋云帶中活躍著中尺度氣旋式渦旋系統(tǒng)。
(2)數(shù)值模擬和風場反演結(jié)果也表明,臺風Matsa的螺旋云帶中在低層活躍著中尺度氣旋式渦旋系統(tǒng),與之相伴隨的為較強的中尺度上升區(qū)域,而且,中尺度垂直上升運動的強弱與雷達對流回波強度呈正相關(guān),中尺度垂直上升運動越強,雷達對流回波發(fā)展越旺盛。但最強的雷達回波中心與中尺度上升運動中心并不完全重合。
(3)臺風Matsa與中尺度系統(tǒng)動能轉(zhuǎn)換的診斷分析說明,低層中尺度系統(tǒng)從臺風Matsa環(huán)流中獲得動能而發(fā)展,Matsa在陸地上長久維持主要是從高層獲得動能。
(4)臺風Matsa與中尺度系統(tǒng)渦度轉(zhuǎn)換的診斷分析說明,低層中尺度系統(tǒng)向Matsa輸送正渦度主要依靠中尺度垂直運動來完成,高層正渦度的轉(zhuǎn)換通過水平輸送和垂直輸送共同來完成。所以,中尺度系統(tǒng)所產(chǎn)生的正渦度源源不斷地向Matsa轉(zhuǎn)換,使Matsa的氣旋性環(huán)流可以在陸地上長久維持。
(References)
Chen F,Dudhia J.2001.Coupling an advanced land surface-h(huán)ydrology model with the Penn State-NCAR MM5modeling system.Part I:Model implementation and sensitivity[J].Mon.Wea.Rev.,129:569-585.
Chen L S.1998.Decay after landfall[C]∥WMO/TD,875:1-5.
陳聯(lián)壽,徐祥德,羅哲賢.2002.熱帶氣旋動力學引論 [M].北京:氣象出版社,317pp. Chen Lianshou,Xu Xiangde,Luo Zhexian.2002.Introduction to Tropical Cyclone Dynamics(in Chinese)[M].Beijing:China Meteorological Press,317pp.
陳受鈞,謝安.1981.次天氣尺度與天氣尺度系統(tǒng)間動能交換的診斷分析 [J].氣象學報,39(4):408-415. Chen Shoujun,Xie An.1981.The exchange of kinetic energy between synoptic and subsynoptic-scale motion[J].Acta Meteorologica Sinica(in Chinese),39(4):408-415.
陳永林,王智,楊引明.2008.0509號麥莎臺風次天氣尺度系統(tǒng)的非對稱性結(jié)構(gòu)特征分析 [J].大氣科學研究與應(yīng)用,(1):9-17.Chen Yonglin,Wang Zhi,Yang Yinming.2008.Analysis of asymmetric characteristics on the sub-synoptic scale systems of typhoon Matasa(0509)[J].Atmospheric Science Research and Application(in Chinese),(1):9-17.
DeMaria M.1996.The effect of vertical shear on tropical cyclone intensity change[J].J.Atmos.Sci.,53:2076-2087.
鄧國,周玉淑,李建通.2005.臺風數(shù)值模擬中邊界層方案的敏感性試驗I:對臺風結(jié)構(gòu)的影響 [J].大氣科學,29(3):417-428.Deng Guo,Zhou Yushu,Li Jiantong.2005.The experiments of the boundary layer schemes on simulated typhoon.Part I:The effects on the structure of typhoon [J].Chinese Journal of Atmospheric Sciences(in Chinese),29(3):417-428.
段麗,陳聯(lián)壽.2005.熱帶風暴 “菲特”(0114)特大暴雨的診斷研究 [J].大氣科學,29(3):343-353. Duan Li,Chen Lianshou.2005.Diagnostic analysis and numerical study of torrential rain associated with the tropical storm Fitow (0114)[J].Chinese Journal of Atmospheric Sciences(in Chinese),29(3):343-353.Emanuel K A.1988.The maximum intensity of hurricanes[J].J.Atmos.Sci.,45:1143-1155.
Janji Z I.1994.The step-mountain eta coordinate model:further developments of the convection,viscous sublayer,and turbulence closure schemes[J].Mon.Wea.Rev.,122:927-945.
Janji Z I.2000.Comments on“development and evaluation of a convection scheme for use in climate models”[J].J.Atmos.Sci.,57:3686.
冀春曉,陳聯(lián)壽,趙放.2007.登陸臺風Matsa維持機理的數(shù)值研究 [J].氣象學報,65(6):888-894.Ji Chunxiao,Chen Lianshou,Zhao Fang.2007.A numerical study on sustaining of landfalling typhoon Matsa [J].Acta Meteorologica Sinica (in Chinese),65(6):888-894.
Jones S C.2000.The evolution of vortices in vertical shear.PartⅢ:Baroclinic vortices[J].Quart.J.Roy.Meteor.Soc.,126:3161-3185.
勵申申,壽紹文,王信.1992.登陸臺風與其外圍暴雨的相互作用[J].氣象學報,50(1):33-41. Li Shenshen,Shou Shaowen,Wang Xin.1992.The interaction between the landed typhoon and its external heavy rain [J].Acta Meteorologica Sinica (in Chinese),50(1):33-41.
李英,陳聯(lián)壽,徐祥德.2004.登陸熱帶氣旋長久維持的次天氣尺度環(huán)流特征 [J].氣象學報,62(3):257-268. Li Ying,Chen Lianshou,Xu Xiangde.2004.The characteristics of sub-synoptic scale circulation of tropical cyclones sustaining over land[J].Acta Meteorologica Sinica(in Chinese),62(3):257-268.
李英,陳聯(lián)壽,王繼志.2005a.熱帶氣旋登陸維持和迅速消亡的診斷研究 [J].大氣科學,29(3):482-490. Li Ying,Chen Lianshou,Wang Jizhi.2005a.Diagnostic study of the sustaining and decaying of tropical cyclones after landfall[J].Chinese Journal of Atmospheric Sciences(in Chinese),29(3):482-490.
李英,陳聯(lián)壽,徐祥德.2005b.水汽輸送影響登陸熱帶氣旋維持和降水的數(shù)值試驗 [J].大氣科學,29(1):91-98. Li Ying,Chen Lianshou,Xu Xiangde.2005b.Numerical experiments of the impact of moisture transportation on sustaining of the landfalling tropical cyclone and precipitation[J].Chinese Journal of Atmospheric Sciences(in Chinese),29(1):91-98.
Lin Y L,F(xiàn)arley R D,Orville H D.1983.Bulk parameterization of the snow field in a cloud model[J].J.Climate Appl.Meteor.,22:1065-1092.
Montgomery M T,Kallenbach R J.1997.A theory for vortex Rossby-waves and its application to spiral bands and intensity changes in hurricanes[J].Quart.J.Roy.Meteor.Soc.,123:435-465.
孫建華,周海光,趙思雄.2006.2003年7月3~5日淮河流域大暴雨中尺度對流系統(tǒng)的觀測分析 [J].大氣科學,30(6):1104-1118.Sun Jianhua,Zhou Haiguang,Zhao Sixiong.2006.An observational study of mesoscale convective systems producing severe heavy rainfall in the Huaihe River basin during 3-5July 2003[J].Chinese Journal of Atmospheric Sciences(in Chinese),30(6):1104-1118.
Sun J Z,Crook A.1994.Wind and thermodynamic retrieval from single-Doppler measurements of a gust front observed during Phoenix II[J].Mon.Wea.Rev.,122:1075-1091.
Sun J Z,Crook N A.1997.Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint.Part I:Model development and simulated data experiments[J].J.Atmos.Sci.,54:1642-1661.
Sun J Z,Crook N A.1998.Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint.Part II:Retrieval experiments of an observed Florida convective storm [J].J.Atmos.Sci.,55:835-852.
Sun J Z,F(xiàn)licker D W,Lilly D K.1991.Recovery of three-dimensional wind and temperature fields from single-Doppler radar data[J].J.Atmos.Sci.,48:876-892.
王繼志,楊元琴.1995.8807號臺風突然增強與其中尺度關(guān)系的研究 [M]∥85-906-07課題組.臺風科學、業(yè)務(wù)試驗和天氣動力學理論的研究第三分冊.北京:氣象出版社,87-94. Wang Jizhi,Yang Yuanqin.1995.Numerical simulation on heavy rainfall caused by typhoon in South China Sea[M]∥Report on Typhoon Scientific Experiment and Synoptic/Dynamics Theoretical Study,Vol.3(in Chinese).Beijing:China Meteorological Press,87-94.
王俊,俞小鼎,邰慶國,等.2011.一次強烈雹暴的三維結(jié)構(gòu)和形成機制的單、雙多普勒雷達分析 [J].大氣科學,35(2):247-285.Wang Jun,Yu Xiaoding,Tai Qingguo,et al.2011.Analysis on the three-dimensional structure and formation mechanism of a severe hailstorm with single-and dual-Doppler radar data[J].Chinese Journal of Atmospheric Sciences(in Chinese),35(2):247-285.
Wang Y,Wu C C.2004.Current understanding of tropical cyclone structure and intensity changes:A review [J].Meteor.Atmos.Phys.,87:257-278.
Wu B,Verlinde J,Sun J Z.2000.Dynamical and microphysical retrievals from Doppler radar observations of a deep convective cloud[J].J.Atmos.Sci.,57(2):262-283.
于玉斌,段海霞,炎利軍,等.2008.超強臺風 “桑美”(2006)近海急劇增強過程數(shù)值模擬試驗 [J].大氣科學,32(6):1365-1378.Yu Yubin,Duan Haixia,Yan Lijun,et al.2008.Numerical simulation of rapid intensification change of super typhoon Saomai(2006)over the coastal water of China[J].Chinese Journal of Atmospheric Sciences(in Chinese),32(6):1365-1378.
袁金南,谷德軍,梁建茵.2005.地形和邊界層摩擦對登陸熱帶氣旋路徑和強度影響的研究 [J].大氣科學,29(3):429-437.Yuan Jinnan,Gu Dejun,Liang Jianyin.2005.A study of the influence of topography and boundary layer friction on landfalling tropical cyclone track and intensity[J].Chinese Journal of Atmospheric Sciences(in Chinese),29(3):429-437.
周海光,王玉彬.2005.2003年6月30日梅雨鋒大暴雨β和γ中尺度結(jié)構(gòu)的雙多普勒雷達反演 [J].氣象學報,63(3):301-312.Zhou Haiguang,Wang Yubin.2005.Structure of meso-βand-γscale on Meiyu in Huaihe River basin on 30June,2003by dual-Doppler radar[J].Acta Meteorologica Sinica(in Chinese),63(3):301-312.
Analysis of the Characteristics of Mesoscale Disturbance for Landfalling Typhoon Matsa
JI Chunxiao1,ZHAO Fang2,GAO Shouting3,and LIU Liping4
1ZhejiangInstituteofMeteorologicalSciences,Hangzhou310017
2ZhejiangProvinceMeteorologicalObservatory,Hangzhou310017
3InstituteofAtmosphericPhysics,ChineseAcademyofSciences,Beijing100029
4ChineseAcademyofMeteorologicalSciences,Beijing100081
Doppler radar reflectivity and the surface automatic weather stations(AWS)wind observations demonstrate that the mesoscale cyclonic eddies at the low level are active in the spiral-cloud bands of landfalling typhoon Matsa.Numerical simulation has been carried out to investigate the evolution characteristics of typhoon Matsa in 2005after its landfalling by using the new-generation mesoscale weather research and forecasting(WRF)model,and the four-dimensional variational Dopple radar analysis system (4D-VDRAS)is used to retrieve single-Doppler velocity.Then a preliminary study has been made on the mesoscale features of typhoon Matsa after landfalling in parallel with diagnostic analyses of the conversion of kinetic energy and vorticity between typhoon Matsa and mesoscale systems on the area mean basis.Results show that:(1)The results of numerical simulation and Doppler radial velocity retrieval indicate that the mesoscale cyclonic eddies at the low level are active in the spiral-cloud bands of landfalling typhoon Matsa,and there is severe mesoscale upward movement accompanying the mesoscale cyclonic eddies.The stronger the mesoscale upward movement is,the larger Doppler radar reflectivity is.(2)The diagnostic results of the kinetic energy conversion between Matsa and mesoscale systems demonstrate that the mesoscale cyclonic eddies at the low level obtain kinetic energy from landfalling typhoon Matsa to develop,whilst Matsa obtains kinetic energy from the mesoscale disturbances in the upper layers to be sustained after landfalling.(3)The diagnostic results of the vorticity conversion between Matsa and mesoscale systems show that Matsa obtains positive vorticity from the mesoscale cyclonic eddies at the low levels by means of the mesoscale upward movement,but in the upper layers,both the mesoscale horizontal and vertical movements play important roles.Therefore,the positive cyclonic vorticity is transported continuously to Matsa,and Matsa could be sustained for a long time after landfalling.
landfalling typhoon,mesoscale disturbance,maintenance,kinetic energy and vorticity conversion
1006-9895(2012)03-0551-13
P444
A
10.3878/j.issn.1006-9895.2011.11085
冀春曉,趙放,高守亭,等.2012.登陸臺風 Matsa(麥莎)中尺度擾動特征分析 [J].大氣科學,36(3):551-563,
10.3878/j.issn.1006-9895.2011.11085.Ji Chunxiao,Zhao Fang,Gao Shouting,et al.2012.Analysis of the characteristics of mesoscale disturbance for landfalling typhoon Matsa[J].Chinese Journal of Atmospheric Sciences(in Chinese),36(3):551-563.
2011-04-29,2011-09-16收修定稿
國家自然科學基金資助項目90815028,國家公益性行業(yè)(氣象)科研專項GYHY201106035
冀春曉,女,1963年出生,碩士,正研級高級工程師,研究方向:中尺度數(shù)值預(yù)報和雷達資料應(yīng)用。E-mail:jichunxiao@sina.com