摘要
利用多普勒雷达资料反演的风场及WRF中尺度模式对2006年登陆我国的台风过程“碧利斯”的结构进行了分析研究,并利用反演风场进行同化,试图改进对台风的预报。
1)基于两步变分法,利用两部多普勒雷达资料反演的风场对0604号台风过程结构进行了分析研究。结果表明,该方法可以较好的反演出台风的水平、垂直风场特征,水平风场呈现有不对称性且有明显的偏心结构,流场随高度呈现出漏斗形特征,垂直风速与水平流场对应较好,台风中心有下沉气流,外围有上升气流。通过对台风物理量的分析发现,此次台风过程存在以最大风速半径随高度向外倾斜的主环流圈和低层为向中心流入,高层向外流出的次环流圈,并且利用WRF模式模拟结果检验了这次台风过程反演结果的合理性。
2)利用反演的水平风场对常用的对称风场涡旋模型进行了验证,发现在最大风速圈内取Rankine模型,最大风速圈外取Chen 3模型对此次台风过程拟合较好。
3)对比分析了单、双雷达资料对于此次台风过程的反演结果,将经验对称涡旋模型作为弱约束引入两步变分方法,提高了两步变分法对此次台风过程的反演能力,改进了单雷达资料反演的效果。
4)利用中尺度模式WRF及其3D-Var分别直接同化径向风、间接同化径向风(同化反演风)资料,以期能够改进预报。结果显示:同化反演风对风场的调整较大,台风呈现非对称性,有明显偏心结构。对台风路径的预报表明,在预报的前9个小时,直接同化和间接同化径向风资料都能提高对台风路径的预报能力,而所有试验预报的台风中心气压都偏强,间接同化反演风的同化预报效果好于直接同化预报效果。
The two-step variational method was used to analyse the 3D wind structure of typhoon Bilis by using two Doppler radars. WRF model and WRF 3D-Var were employed to estimate the effect of radar data assimilation on the forecasts of this typhoon process.
Firstly, The 3D wind fields of this typhoon process was analyzed by using the two-step variational wind retrieval method. The results show that the method has ability to retrieve the wind characteristics of typhoon both in the horizontal and vertical sections. The asymmetrical structure can be obviously seen from the horizontal wind fields. The wind field is presented the funnel-shape characteristic with the increase of height. The downdraft occurred in the inner core and the updraft in the external part of the typhoon. By analyzing on the variables, it is found that there are a principal circulation in which the radius of the maximum wind speed dipped outward with the height and a secondary circulation with the inflow at the lower layers and outflow at the higher layers, then we utilize the simulation result by wrf model to test the retrieval result of this typhoon process.
Subsequently, the retrieval wind fields were used to verify the common vortex models. The results show that the Rankine model fits well within the radius of the maximum tangential velocity, and the Chen 3 model fits well outside the radius of the maximum tangential velocity.
Compared the retrieval wind from single radar with dual-doppler radar data, a constraint on this typhoon pattern was introduced into the cost function in the two-step variational method. It improves the retrieval ability of two-step variational method from single-Doppler radar data. The retrieved winds were consistent with that from two Doppler radar data.
Finally, the radial velocities were assimilated into WRF 3D-Var system and the forecasts were performed by WRF model to test the effect of radial velocity assimilaton. Three experiments was designed: 1) No data was assimilated; 2) radial velocity was preprocessed to the size (about) model horizontal resolution and assimilated directly into 3D-Var system; 3) The retrieval winds by two-step variational method was assimilated into 3D-Var system. The results show that the wind is greatly adjusted and the horizontal wind is asymmetrical and eccentrical while assimilating retrieval wind. In the first nine-hour forecasts, radar data assimilation can improve the skill of the typhoon tracks forecasts. However, for all the experiments the forecasted central sea level pressures of this typhoon deepen than the observations. The forecasts of assimilating retrieval wind are better than assimilating radial velocity directly.
引文
[1]Lhermitte RM,Atlas D,Precipitation motion by Pulse Doppler Radar,Preprints 19~(th) Weather Radar Conf[J].American Meteor Society,1961,218-223.
[2]Browning K A,Wexler R A.The determination of kinematic properties of a wind field using Doppler radar[J].Journal of applied Meteorology,1968,7:105-113.
[3]Persson POG,Andersson T.A real-time system for automatic single-Doppler wind field analysis[M].Proc.Syrup.Mesoscale Analysis & Forecasting,Vancouver,ESA Publication SP-282,1987,61-66.
[4]陶祖钰.单Doppler速度场反演风矢量的VAP方法[J].气象学报,1992,50(I):81-90.
[5]刘晓阳,刘桂馥.方差反演多普勒速度水平风场[J].南京气象学院学报,2000,23(4):549-554.
[6]郎需兴,魏鸣,葛文忠等.一种新的单多普勒雷达风场反演方法[J].气象科学,2001,21(4):417-424.
[7]Zhao K,Liu GQ,Ge WZ.Retrieval of single-Doppler radar wind field by nonlinear approximation[J].Advances in Atmospheric Sciences,2003,20(2):195-204.
[8]Waldteufel P,Corbin H.On the analysis of single-Doppler radar Data.Joumal of Applied Meteorology[J],1979,18:532-542.
[9]Tulle JD.Foote GB.Determination of the boundary layer airflow from a Single Doppler Radar[J].Journal of Atmospheric and Oceanic Technology,1990,7(2):218-232.
[10]Sun J,Ficker DW,Lilly DK.Recovery of Three-Dimensional Wind and Temperature Fields from Simulated Single-Doppler Radar Data[J].Journal of the Atmospheric Sciences,1991,48:876-890.
[11]Sun J,Crook A.Wind and thermodynamic retrievals from single-Doppler measurements of a gust front observed during Phoenix Ⅱ[J].Monthly Weather Review,1994,122:1075-1091.
[12]Qiu CJ,Xu Q.A Simple Adjoint Method of Wind Analysis for Single-Doppler Data[J].Journal of Atmospheric and Oceanic Technology,1992,9:588-598.
[13]Xu Q,Gu H,Yang S.Simple adjoint method for three-dimensional wind retrievals from single-Doppler radar[J].Quarterly Journal of the Royal Meteorological Society,2001,127:1053-1067.
[14]Qiu CJ,Shao AM,Liu S,Xu Q.A Two-step Variational Method for Three-Dimensional Wind Retrieval from Single Doppler Radar[J],Meteorology and Atmospheric Physics.2006,91:1-8.
[15]姜海燕,葛润生.一种新的单多普勒雷达反演技术[J].应用气象学报,1997,8(2):219-223.
[16】 吴绍荣,徐宝祥,周秀骥.同化方法反演单多普勒雷达水平风场[M].大气遥感论文集,1997,22-27
[17]Nissen R,Hudak D,LAROCHE S,Elia R,et al.3D Wind Field Retrieval Applied to Snow Events Using Doppler Radar.Joumal of Atmospheric and Oceanic Technology,2001,18:348-362.
[18]Xiao QN,Kuo YH,Sun JZ,Lee WC.Assimilation of Doppler Radar Observation with a regional 3DVAR system:Impact of Doppler velocities on forecasts of a heavy rainfall case [J].Journal of Applied Meteorology,2005,44:768-788.
[19]Gao J,Xue M,Shapiro A,Droegemeier KK.A variational method for the analysis of three-dimensional wind fields from two Doppler radars[J].Monthly Weather Review,1999,127:2128-2142.
[20]Tong ML.Xue M.Ensemble Kalman Filter Assimilation of Doppler Radar Data with a Compressible Nomhydrostatic ModeI:OSS Experiments[J].Monthly Weather Review,2005,133:1789-1805.
[21】 郭霞,党人庆,葛文忠.雷达资料在江淮暴雨数值模拟中的应用[J].热带气象学报,1999,15(4):356-362.
[22]马清云,李泽椿等.单多普勒天气雷达风场反演及其在数值预报中的应用试验[J].应用气象学报,2001,12(4):488-493.
[23]徐枝芳,徐玉貌,葛文忠.雷达和卫星资料在中尺度模式中的初步应用[J].气象科学,2002,22(2):167-174.
[24】 朱玉祥,多普勒雷达资料在MMS模式中的4维变分同化试验[D],南京气象学院硕士论文,2004.
[25】李柏,多普勒天气雷达资料分析及同化在暴雨中尺度天气系统数值模拟中的应用研究[D].中国气象科学院博士论文,2005.
[26]李萍阳,雷达反演资料在中尺度模式中的应用[A】.奥运气象保障技术研究2004年度论 文汇编气象出版社,2004.
[27]许小永,郑国光,刘黎平.多普勒雷达资料4DVAR同化反演的模拟研究[J].气象学报,2004,62(4):410-422.
[28]万齐林,薛纪善,陈子通等.雷达TREC风的三维变分同化应用与试验[J】.热带气象学报,2005,21(5):449-457.
[29]盛春岩,薛德强,雷霆等.雷达资料同化与提高模式水平分辨率对短时预报影响的数值对比试验[J].气象学报,2006,64(3):293-308.
[30]顾建峰.多普勒雷达资料三维变分直接同化方法研究[D].中国气象科学研究所博士论文,2006.
[31]杨毅.Doppler雷达资料同化技术研究[D】.兰州大学博士论文,2007:48-49.
[32]闵锦忠,彭霞云,赖安伟等.反演同化和直接同化多普勒雷达径向风的对比试验[J].南京气象学院学报,2007,30(6):745-754.
[33】 吴陈锋,吴智辉.“艾利”热带风暴多普勒雷达定位方法的探讨[J].广西气象,2005,26:125-127.
[34]张培昌,杜秉玉,戴铁丕.雷达气象学[M].北京,气象出版社,2001,511pp.
[35]许映龙,矫梅燕,毕宝贵,张涛.近海台风雷达定位方法应用研究[J].大气科学,2006,30(6):1119-1126.
[36】 陈绍伟.使用多普勒速度进行台风拟合定位的试验[J].气象研究与应用,2007,8:95-97.
[37]Shea D J,Gray W M.The Hurricane's Inner core Region.I.Symmetric and Asymmetric structure[J].Journal of the Atmospheric Sciences,1973,30:1544-1564.
[38]Liu Y,Zhang DE,Yau MK.A Multiscale numerical study of Hurricane Andrew(1992).Part Ⅱ:kinematics and Inner-core structures[J].Monthly Weather Review,1999,127:2597-2616.
[39]Marks FD Jr,Houze RA Jr,Gamache JE Dual-Aircraft Investigation of the Inner core of Hurricane Norbert.Part Ⅰ:kinematic structure[J].Journal of the Atmospheric Sciences,1992,49(11):919-942.
[40]Marks FD,Houze RA.Inner core structure of Hurricane Alicia from Airborne Doppler Radar observations[J].Journal of the Atmospheric Sciences,1986,44(9):1296-1317.
[41]Roux F,Viltard N.Structure and evolution of Hurricane Claudette on 7 September 1991 from Airborne Doppler Radar observations.part 1:kinematics[J].Monthly Weather Review,1995,123:2611-2639.
[42]Dodge P,Burpee RW,Marks FD.The Kinematic Structure of a Hurricane with Sea Level Pressure Less Than 900 mb[J].Monthly Weather Review,1999,127:987-1004.
[43]Corbosiero K L,J Molinari,A R Aiyyer et al.The structure and Evolution of Hurricane Elenal(1985).Part Ⅱ:convective Asymmetries and Evidence for vortex Rossby waves[J].Monthly Weather Review,2006,134:3073-3091.
[44]Lee W C,Jou BJD,Chang PL,et al.Tropical cyclone kinematic structure retrieved from single Doppler radar observations.Part Ⅲ:Evolution and structures of Typhoon Alex(1987)[J].Monthly Weather Review,2000,128:3982-4001.
[45]Knaff JA,Seseske SA,Demaria JL,et al.On the influences of vertical wind shear on symmetric tropical cyclone structure derived from AMSU[J].Monthly Weather Review,2004,132:2503-2510.
[46]Chen Y,Yau MK.Asymmetric structures in a simulated landfalling Hurricane[J].Journal of the Atmospheric Sciences,2003,60:2294-2312.
[47]Black PG,Anthes RA.On the asymmetric structure of the tropical cyclone outflow layer[J].Journal of the Atmospheric Sciences,1971,28:1348-1366.
[48]魏应植,汤达章,许键民等.多普勒雷达探测”艾利”台风风场不对称结构[J].应用气象学报.2007,18(3):285-293.
[49】 彭霞云,多普勒雷达径向风反演同化和直接同化试验研究[D].南京信息工程大学硕士学位论文,2006,33-53.
[50]闽锦忠,彭霞云,赖安伟等.反演同化和直接同化多普勒雷达径向风的对比试验[J].南京气象学院学报.2007,30(6):475-454.
[51]傅洁,多普勒雷达径向速度三维变分同化方法的试验研究[D].南京信息工程大学硕士学位论文,2006,39-45.
[52]Xiao Q,Xiaoyan Zhang1,Christopher Davis et al.Hurricane Initialization Using Airborne Doppler Radar Data for WRF(ERAD).
[53]Xiao Q,Lim E,Won DJ,et al.Doppler Radar Data Assimilation in KMA's Operational Forecasting[J].Bulletin of the American Meteorological Society,2008:39-43.
[54]Xiao Q,Kuo YH,Sun J,et al.Assimilation of Doppler radar observations and its impacts on forecasting of the landfalling typhoon Rusa(2002).Proceedings of the Third European Conference on Radar in Meteorology and Hydrology(ERAD),2004:178-182.
[55]减重清.气象研究所研究报告[R],1977,28(2):83-95.
[56]高桥浩一郎.气象杂志[J].1938,17(11):417-421.
[57]Fujita T.Pressure distribution within a typhoon[J].Geophysical Magazine.1952,23:437-451.
[58]Mayers VA.Characteristics of United States hurricanes pertinent to level design for Lake Okeechobee Florida[R].U.S.Weather Bureau Hydro meteorological Report,1954,32.
[59]Jelesnianski CR A numerical calculation of storm tides induced by a tropical storm impinging on a continental shelf[J].Monthly Weather Review,1965,93(6):343-358
[60]Jelesnianski CE Numerical computation of storm surges without bottom strees[J].Monthly Weather Review,1966,94(6):374-379.
[61]陈孔沫.一种计算台风风场的方法[J].热带海洋,1994,13(2):41-48.
[62]陈联寿.热带气旋运动的研究和业务预报的现状和发展[A】.台风会议文集.北京气象出版社,1985:6-30.
[63]朱首贤,沙文钰,丁平兴等.近岸台风非对称风场模型[J].华东师范大学学报,2002,(3):66-71.
[64]杨支中,沙文钰,朱首贤等.一种新型的非对称台风海面气压场和风场模型[J].海洋通报,2005,24(1):62-68.
[65]黄小刚,费建芳,张根生等.一种台风海面非对称风场的构造方法[J】.热带气象学报,2004,20(2):129-136.
[66]杨毅,邱崇践.利用多普勒雷达资料分析一次强降水的中尺度流场[J].高原气象,2006,25(51:925-931.
[67]Zawadzki I,Ostiguy L,Laprise JPR.Retrieval of the microphysical properties in a CASP storm by integration of a numerical kinemetic model[J].Atmosphere Ocean,1993,31:201-233.
[68]蒋小平,刘春霞,费志宾等.南海夏季风对强热带风暴Bilis(0604)引发暴雨的影响[J].热带气象学报,2008,24(4):379-384.
[69】 许美玲,段旭,张秀年.登陆孟湾风暴结构演变特征数值模拟研究[J].高原气象,2006,25(6):1139-1146
[70]徐国强,梁旭东,余晖等.不同云降水方案对一次登陆台风的降水模拟[J].高原气象,2007,26(5):891-900
[71]康志明,陈涛,钱传海等.0604号强热带风暴“碧利斯”特大暴雨的诊断研究[J].高原气象,2008,27(3):596-607.
[72]张东,蔡安安,林良勋.强热带风暴“碧利斯”致洪暴雨的特征及成因[J].广东气象,2007,29(1):22-24
[73]郭达烽,郑婧,许爱华.台风“碧利斯”的结构与江西暴雨诊断分析[J].气象与减灾研究,2006,29(3):52-59.
[74]张海霞,崔晓鹏,康凤琴等.邯郸地区一次登陆台风大暴雨过程观测分析[J].高原气象,2007,26(5):980-991
[75]郭荣芬,鲁亚斌,李燕等.“伊布都”台风影响云南的暴雨过程分析[J].高原气象,200524(5):784-791
[76]毛连海,尹洁,金米娜.“碧利斯”的特点及动力成因分析[J].气象与减灾研究.2006,29(3):47-51.
[77]支树林.台风“碧利斯”影响庐山降水的探测资料分析[J].气象与减灾研究,2006,29(3):60-64.
[78]孙瑞,龙春华,郑芬.强热带风暴“碧利斯”的多普勒雷达回波分析[J].广西气象,2006,27(3):36-37.
[79]Barnes SL.A technique for maximizing details in numerical weather map analysis.Journal of Applied Meteorology,1964,3:396-409.
[80]陈孔沫.新的台风风场计算方法[J].海洋预报,1992,9(3):60-65.
[81]Kessler E.On the distribution and continuity of water substance in atmospheric circulation [J].American Meteor Society,1969,32:84pp.
[82]Lin YL,Farley RD,Orville HD.Bulk parameterization of the snow field in a cloud model[J].Journal of Applied Meteorology,1983,22:1065-1092.
[83]Skamarock WC,Klemp JB,Dudhia J,Gill DO,Barker DM,Wang W,Powers JG.A description of the Advanced Research WRF Version 2.NCAR Tech.Note,ncar/TN-468+STR,2005.
[84]Lorenc AC.Analysis methods for numerical weather prediction[J].Quarterly Journal of the Royal Meteorological Society,1986,112:1177-1194.
[85]Courtier P,Thepaut JN,Hollingsworth A.A strategy for operational implementation of 4D-Var,using an incremental approach[J].Quarterly Journal of the Royal Meteorological Society, 1994,120: 1367-1387.
[86]Parrish DF, Derber JC.The National Meteorological Center's Spectral Statistical Interpolation analysis system [J].Monthly Weather Review, 1992: 120: 1747-1763.
[87]Sun J, Crook NA.Dynamical and microphysical retrieval from Doppler Radar observations using a cloud model and its adjoint.Part 1: Model development and simulated data experiments [J].Journal of the Atmospheric Sciences, 1997, 54: 1642-1661.
[88]Sun J, Crook NA.Dynamical and microphysical retrieval from Doppler Radar observations using a cloud model and its adjoint.Part Ⅱ: Retrieval experiments of an observed Florida convective storm [J].Journal of the Atmospheric Sciences, 1998, 55: 835-852.