利用3维变分方法同化多普勒天气雷达资料的试验研究

以美国新近研发的天气研究预报模式(WRF)配置的3维变分同化系统WRF 3D-Var为平台,针对2003年7月4-5日的一次暴雨过程进行了一组同化多普勒天气雷达径向风和回波强度资料的试验研究.其中同化径向风资料的方法是先用两步变分法由多普勒资料反演出水平风,再同化反演风场;同化雷达回波强度资料的方法是通过雷达回波强度与雨水混合比的经验关系以暖雨过程为约束的同化方法.试验结果表明同化雷达资料后对降水预报有显著的正作用.其中单独同化雷达回波强度资料的正作用在前3～6 h效果非常显著;单独同化反演的水平风场对降水预报的正作用能持续12 h;二者同时同化对提前9 h的降水预报改进更加明显.     

同化多普勒雷达风资料的两种方法比较

以美国新近研发的天气研究预报模式(WRF)配置的三维变分同化系统WRF 3D-Var为平台,比较了两种不同的同化多普勒雷达径向风资料的方法。一种是WRF 3D-Var系统现有的径向风资料直接同化方法;另一种是首先用两步变分法由多普勒资料反演出水平风,再同化反演风场。针对2003年7月4~5日的一次淮河暴雨过程进行的同化试验结果表明,同化了雷达风资料后得到的水平风场包含了更多的中尺度特征;从降水预报评分和预报的雷达回波来看,两种方法都能够明显改进降水预报,这种正作用能维持6 h左右;相对而言,同化反演的水平风场的效果略优于直接同化雷达径向风的效果。     

不同雷达观测资料同化对一次罕见飑线天气模拟的影响

选用WRF(weather research and forecasting)模式及其3D-Var(three-dimensional variation)同化系统,针对2018年3月4日发生在江西的一次罕见强飑线天气,探讨同化多普勒雷达不同观测资料对极端雷暴大风天气模拟预报的影响。结果表明:仅同化由雷达反射率反演的雨水、雪和霰粒子以及由其估算的水汽不能稳定改善模式对飑线雷达反射率的预报效果,尤其对地面大风和降水的预报起反效果;当联合同化雷达反射率与雷达径向风资料后,显著改进了模式对飑线发展演变过程中雷达反射率、地面降水和地面大风的预报效果,雷达反射率的同化呈现显著正效果。原因是仅同化雷达反射率对初始水成物及热力场影响较大,而对动力场调整微弱,随着积分时间增加,热力场对动力场的反馈作用不真实,高层出现虚假辐散风场,飑线前侧模拟出虚假层状云区,且未能改进飑线系统低层垂直风切变、冷池以及对流层中下层后侧入流的模拟,模拟的飑线移动和演变过程与实况有很大差距;当联合同化雷达反射率与雷达径向风资料后明显调整了初始动力、热力和水成物场,物理配置更符合实际,形成更有利于强飑线发生的垂直风切变和风场结构,产生与实况接近的强冷池,模拟结果与实况的吻合度明显高于仅同化雷达径向风资料的试验。     

雷达资料的3DVAR和EnSRF直接同化方法对一次暴雨预报的影响分析

利用合肥站雷达资料和NCEP再分析资料,采用WRF模式、WRF-3DVAR和WRF-En SRF同化系统设计了直接同化雷达径向风和反射率的试验,对2003年7月4—5日发生在江淮流域的一次梅雨锋暴雨过程进行数值模拟分析。结果表明:(1)同化雷达径向风资料可以增加暴雨区低层的扰动,增强对流运动,同时可以减少模拟的过量降水;(2)3DVAR同化方法对水平风场调整明显,同化反射率改进了模式中的微物理量和动力场,同时同化径向风和反射率可以更好地调整风场,改善模拟雷达回波带的强度和位置;(3)En SRF同化方法可以将更多小尺度的风场信息同化进入模式,单独同化反射率可以较准确预报出降水中心,同时同化径向风和反射率可以改进模拟雷达回波带的位置;(4)使用3DVAR和En SRF同化雷达径向风和反射率可以提高降水预报,改善雨带的范围,对降水中心位置和强度的预报也有改进。     

多普勒天气雷达资料对中尺度模式短时预报的影响

利用中尺度模式ARPS(The Advanced Regional Prediction System)及其资料分析系统ADAS(ARPS Data Analysis System),将国内新一代多普勒雷达(CINRAD)反射率及径向风资料直接用于中尺度数值模拟,通过一次华北地区暴雨过程的模拟对比试验,分析了雷达资料对初始场的改进效果及其对模拟结果的影响,结果表明:(1)利用雷达径向风资料对初始风场进行调整后,自近地面到对流层顶的u,v,w都发生了变化,调整后的初始风场在对流层中层变化最大.(2)利用雷达反射率进行微物理调整和云分析能调整初始场中的云水信息,使得雷达回波附近3 km以下的水汽混合比(qv)增加,4 km以下的雨水混合比(qr)增加,对流层(约10 km以下)的云水混合比(qc)增加,4～9 km的对流层上部云冰混合比(qi)和雪混合比(qs)增加.ADAS通过非绝热初始化调整温度场,从而得到了一个动力和热力上平衡的初始场.(3)模拟的1 h雨量与实况的对比表明,同时利用雷达反射率和径向风改进过的初始场能明显增强3 h内的降水强度和落区预报,改善中尺度数值模式短时定量降水预报.模拟的1 h流场对比分析表明,经雷达径向风调整后,能够在初始场中增加气旋性涡旋等中小尺度风信息,明显减少模式的spin-up时间.(4)通过对雷达径向风和反射率对模式初始场和模拟结果影响的对比分析发现,雷达径向风主要是改进初始风场,而雷达反射率主要是改进初始场中的湿度参数,增加初始场中云水等的含量,调整温度场.通过模拟的6 h降水对比发现,利用雷达径向风调整初始场后,对降水模拟有一定的改进,但效果不甚明显,而雷达反射率资料对定量降水预报改进效果明显,同时使用雷达径向风和反射率资料改进初始场后对降水的模拟效果最明显.     

基于GSI的华南地区对流尺度快速循环同化预报试验

针对对流尺度快速循环同化系统多次循环同化带来的预报效果改进和资料应用问题,利用GSI同化技术和WRFARW区域模式,设计了华南地区对流尺度快速循环同化方案,对2016年4月17一18日华南地区的飑线天气强降水过程进行模拟试验,分析不同循环同化方案和雷达径向风资料同化对雷达回波、相对湿度、降水量级等的预报效果,以期提高华南地区飑线强降水过程预报技巧。检验结果表明:尽管只同化常规资料对预报效果的改进有局限性,但是多次循环同化对于模式预报的降水有一定改善作用;同时同化雷达径向风资料与常规资料对湿度和降水等模拟技巧均有所提高,大雨以上量级的ETS评分改进尤为明显;尽管模式模拟降水峰值小于真实观测值,但同化雷达径向风资料有效改善了飑线最强时段内的垂直上升速度,使得强降水发生时间和强度更接近真实观测。     

雷达热动力反演资料的高分辨率同化模拟试验

目前,北京地区的天气预报系统对局地对流性定量降水预报能力较弱,远不能满足人们生产、生活和防灾、减灾工作的需要。针对北京地区对提高0-12 h短时临近天气,尤其是夏季局地对流性降水预报能力的需求,基于中国气象局北京城市气象研究所变分多普勒雷达分析系统（VDRAS）的雷达热动力反演资料,建立了WRF模式初始化模块,采用四维资料同化（FDDA）方法,将VDRAS系统高时空分辨率三维热动力结构分析场资料同化到WRF模式中,实现了北京地区VDRAS分析场资料在WRF中尺度模式系统中的应用。通过降水个例的高分辨率同化模拟试验分析了雷达热动力反演资料同化对模式预报结果的影响。研究结果表明:雷达热动力反演资料的同化能够提高模式系统对近地面温、湿、风大气要素和降水过程的模拟能力,改善2 m比湿、降水落区、降水量级、降水时间的预报效果,减少降水漏报的现象。温度和比湿的同化比风的同化对模拟降水结果的改善更重要。虽然研究表明雷达热动力反演资料在WRF模式中的同化能够明显改善模式对选取降水个例的模拟效果,但其对模式尤其是数值业务模式系统预报效果的影响需要进一步更全面、更系统的检验,为业务化应用奠定更坚实的基础。      

基于雷达资料同化的飑线过程数值模拟试验研究

利用中尺度WRF模式及其3DVAR同化系统对2014年3月30日发生在我国华南地区的一次飑线过程展开多普勒天气雷达资料的同化效果试验研究。首先对雷达资料进行去地物杂波、退速度模糊等预处理,后设计了基于不同雷达观测量的同化试验及同化频次的敏感性试验。结果表明:直接循环同化雷达径向风资料和雷达反射率因子能够增加数值模式中的中小尺度信息,提供可靠的水汽分布;不同的同化频次对同化结果影响显著,每12 min同化间隔的结果略优于30 min、60 min同化间隔;同化雷达反射率因子和径向风资料分别对模式的总水场和风场有显著调整,联合同化雷达反射率因子和径向风资料比单独同化反射率因子或径向风更能改善飑线垂直结构配置,促使地面冷池和雷暴高压配合,进一步改善模式对大雨和暴雨量级降水预报效果。     

用多普勒雷达反射率调整模式大气的云微物理变量

一种简单云分析方案, 用于由多普勒天气雷达反射率反演中尺度大气模式初值分析中的云微物理变量（云水混合比和雨水混合比）和空气湿度变量（比湿），使模式积分初始场反映出观测空间的云微物理特征以及哪些空间位置上的大气处于饱和状态。应用于2002年6月梅雨期安徽省马鞍山市一次降水过程的临近数值预报试验结果表明，模式预报的大气综合反射率与雷达观测的回波图像相近，由云微物理变量变化表示的模式云系演变与雷达观测的回波图像一致, 伴随模拟的中小尺度云系, 模式大气能很快调整出合理的中小尺度流场辐散、辐合结构；它们明显好于模式初始场不引入雷达反射率时的结果，即这种方法对改进临近数值天气预报准确率是有效的。     

SVD-En3DVar方法同化多普勒雷达速度观测资料Ⅰ.模拟资料试验

利用WRF模式及模式模拟的资料,开展了利用SVD-En3DVar(基于集合和SVD技术的三维变分同化方法)方法同化雷达径向速度资料的试验.由于雷达观测经常出现大面积空缺,同化时引入了一种局地化方法避免远距离虚假相关的影响.试验着重研究了不同的初始扰动样本产生方法以及不同的样本积分时间对同化结果的影响.提出了一种为预报集合提供初始扰动场的新方法,这一方法将温度和比湿的伪随机扰动场当作观测增量,通过3DVar (three-dimensional variational technique)系统生成所有变量的初始扰动场.试验表明,用这种方法给出的初始扰动样本各个变量间有较好的协调性,积分后扰动不会快速衰减,可以减少模式调整的时间,达到缩短同化循环时间窗的目的.同化雷达径向风资料后对12小时内的温度,湿度和水平风的预报都有所改进,对降水的预报也有一定改进.     

The WRF 3DVar System Combined with Physical Initialization for Assimilation of Doppler Radar Data

The three-dimensional variational data assimilation (3DVar) system of the Weather Research and Forecasting (WRF) model (WRF-Var) is further developed with a physical initialization (PI) procedure to assimilate Doppler radar radial velocity and reflectivity observations. In this updated 3DVar system, specific humidity, cloud water content, and vertical velocity are first derived from reflectivity with PI, then the model fields of specific humidity and cloud water content are replaced with the modified ones, and finally, the estimated vertical velocity is added to the cost-function of the existing WRF-Var (version 2.0) as a new observation type, and radial velocity observations are assimilated directly by the method afforded by WRF-Var. The new assimilation scheme is tested with a heavy convective precipitation event in the middle reaches of Yangtze River on 19 June 2002 and a Meiyu front torrential rain event in the Huaihe River Basin on 5 July 2003. Assimilation results show that the increments of analyzed variables correspond well with the horizontal distribution of the observed reflectivity. There are positive increments of cloud water content, specific humidity, and vertical velocity in echo region and negative increments of vertical velocity in echo-free region where the increments of horizontal winds present a clockwise transition. Results of forecast experiments show that the effects of adjusting cloud water content or vertical velocity directly with PI on forecast are not obvious. Adjusting specific humidity shows better performance in forecasting the precipitation than directly adjusting cloud water content or vertical velocity. Significant improvement in predicting precipitation as well as in reducing the model's spin-up time are achieved when radial velocity and reflectivity observations are assimilated with the new scheme.     

雷达反射率资料的三维变分同化研究

应用天气研究和预报模式(WRF)三维变分系统中一种新的雷达反射率资料间接同化方法来进行反射率资料的三维变分同化研究,评估雷达反射率资料对夏季短时定量降水预报的作用.该方法不直接同化雷达反射率资料,而是同化由反射率资料反演出的雨水和估计的水汽.以2009年夏季北京地区发生的4次强降水过程为例,考察了北京市气象局业务运行的快速更新循环同化预报系统对京津冀地区雷达网的雷达反射率资料的同化性能以及雷达反射率资料和径向风资料同时同化的效果.数值试验结果表明:(1)同化反演雨水或水汽都能改善降水预报,但同化反演水汽对降水预报效果的改善起了更重要的作用;(2)同化反射率资料能极大地提高短时降水预报的效果,其稳定的正面效果可以延伸到6h的预报时效,而同化径向风资料不能得到稳定的正效果;(3)同化雷达资料时,应用快速更新循环同化预报系统是提高短时定量降水预报的一个有效途径.     

The impact of radar data assimilation on atmosphere state analysis in the Moscow region

The results of numerical experiments with the data assimilation system including the WRF-ARW mesoscale atmospheric model and WRFDA analysis package in the 3D-Var mode are considered. The focus is on the impact of Doppler weather radar data on the quality of short-range weather forecasting. The maps of weather events and cloud top constructed from the modeling by GIMET-2010 software package are analyzed using the web-GIS METEORAD. The experiments in the Moscow region demonstrated that if only radial wind is assimilated, the intensity of simulated cloud and precipitation formation is the closest to the observed one; however, the spatial mismatches of simulated and observed cloud systems are possible. If Doppler radar data on radial wind and reflectivity are assimilated, the general pattern and location of clouds and precipitation are simulated more accurately, but there is a certain overestimation of convection intensity which leads to the overestimation of the number of thunderstorms and rainfall rate.     

The implementation of reverse Kessler warm rain scheme for radar reflectivity assimilation using a nudging approach in New Zealand

Reverse Kessler warm rain processes were implemented within the Weather Research and Forecasting Model (WRF) and coupled with a Newtonian relaxation, or nudging technique designed to improve quantitative precipitation forecasting (QPF) in New Zealand by making use of observed radar reflectivity and modest computing facilities. One of the reasons for developing such a scheme, rather than using 4D-Var for example, is that radar VAR scheme in general, and 4D-Var in particular, requires computational resources beyond the capability of most university groups and indeed some national forecasting centres of small countries like New Zealand. The new scheme adjusts the model water vapor mixing ratio profiles based on observed reflectivity at each time step within an assimilation time window. The whole scheme can be divided into following steps: (i) The radar reflectivity is firstly converted to rain water, and (ii) then the rain water is used to derive cloud water content according to the reverse Kessler scheme; (iii) The cloud water content associated water vapor mixing ratio is then calculated based on the saturation adjustment processes; (iv) Finally the adjusted water vapor is nudged into the model and the model background is updated. 13 rainfall cases which occurred in the summer of 2011/2012 in New Zealand were used to evaluate the new scheme, different forecast scores were calculated and showed that the new scheme was able to improve precipitation forecasts on average up to around 7 hours ahead depending on different verification thresholds.     

Application of a Three-dimensional Variational Method for Radar Reflectivity Data Correction in a Mudslide-inducing Rainstorm Simulation

Various types of radars with different horizontal and vertical detection ranges are deployed in China, particularly over complex terrain where radar blind zones are common. In this study, a new variational method is developed to correct threedimensional radar reflectivity data based on hourly ground precipitation observations. The aim of this method is to improve the quality of observations of various types of radar and effectively assimilate operational Doppler radar observations. A mudslide-inducing local rainstorm is simulated by the WRF model with assimilation of radar reflectivity and radial velocity data using LAPS(Local Analysis and Prediction System). Experiments with different radar data assimilated by LAPS are performed. It is found that when radar reflectivity data are corrected using this variational method and assimilated by LAPS,the atmospheric conditions and cloud physics processes are reasonably described. The temporal evolution of radar reflectivity corrected by the variational method corresponds well to observed rainfall. It can better describe the cloud water distribution over the rainfall area and improve the cloud water analysis results over the central rainfall region. The LAPS cloud analysis system can update cloud microphysical variables and represent the hydrometeors associated with strong convective activities over the rainfall area well. Model performance is improved and the simulation of the dynamical processes and moisture transport is more consistent with observation.     

雷达资料同化对2016年6月23日阜宁龙卷模拟的改进

基于GSI（Community Gridpoint Statistical Interpolation）同化系统和WRF（Weather Research and Forecasting）模式,探讨了多部多普勒雷达的反射率因子和径向速度资料同化对2016年6月23日江苏阜宁龙卷模拟的改进效果和影响过程。结果表明:（1）仅同化雷达反射率因子和仅同化径向速度均能在一定程度上改进模式对阜宁龙卷及其环境场的模拟,且雷达径向速度同化的改进作用更大;同时同化两种资料改进效果最佳。（2）雷达反射率因子同化是利用复杂云分析技术,直接修正了水凝物含量,增加了潜热释放,对低层大气热力场进行了正温度扰动调整,从而主要改进了初始场的水汽条件和热力条件;而雷达径向速度同化通过三维变分技术直接修正了风场,进而引起水汽输送变化影响水凝物的调整和大气热力场的变化,对初始场的动力条件和热力条件修正较大;同时同化两种资料修正了初始场的动力和热力结构,保证了两者物理上的协调,综合了两者的改进作用,从而取得最佳模拟效果。（3）同时同化雷达反射率因子和径向速度后,模式在阜宁附近模拟出了明显的涡旋结构,尽管涡旋强度和龙卷结构与实况仍有一定差距,但涡旋发生发展过程、路径、地面小时极大风和降水等模拟与实况吻合度均明显高于对照试验。      

Preliminary Study on Direct Assimilation of Cloud-affected Satellite Microwave Brightness Temperatures

Direct assimilation of cloud-affected microwave brightness temperatures from AMSU-A into the GSI three-dimensional variational(3D-Var) assimilation system is preliminarily studied in this paper. A combination of cloud microphysics parameters retrieved by the 1D-Var algorithm(including vertical profiles of cloud liquid water content, ice water content, and rain water content) and atmospheric state parameters from objective analysis fields of an NWP model are used as background fields. Three cloud microphysics parameters(cloud liquid water content, ice water content, and rain water content) are applied to the control variable. Typhoon Halong(2014) is selected as an example. The results show that direct assimilation of cloud-affected AMSU-A observations can effectively adjust the structure of large-scale temperature, humidity and wind analysis fields due to the assimilation of more AMSU-A observations in typhoon cloudy areas, especially typhoon spiral cloud belts. These adjustments, with temperatures increasing and humidities decreasing in the movement direction of the typhoon,bring the forecasted typhoon moving direction closer to its real path. The assimilation of cloud-affected satellite microwave brightness temperatures can provide better analysis fields that are more similar to the actual situation. Furthermore, typhoon prediction accuracy is improved using these assimilation analysis fields as the initial forecast fields in NWP models.     

Evaluation of All-Sky Assimilation of FY-3C/MWHS-2 on Mei-yuPrecipitation Forecasts over the Yangtze-Huaihe River Basin

正1School of Atmospheric Sciences, Chengdu University of Information Technology, Chengdu 610025, China2International Center for Climate and Environment Sciences, Institute of Atmospheric Physics,Chinese Academy of Sciences, Beijing 100029, China3University of Chinese Academy of Sciences, Beijing 100049, China     

集合卡尔曼滤波同化多普勒雷达资料的数值试验

利用集合卡尔曼滤波(EnKF)在云数值模式中同化模拟多普勒雷达资料,并考察了不同条件下EnKF同化方法的性能.结果显示,经过几个同化周期后,EnKF分析结果非常接近真值.单多普勒雷达资料EnKF同化对雷达位置不太敏感,双雷达资料同化结果在同化的初期阶段比单雷达资料同化结果准确.同化由反射率导出的雨水比直接同化反射率资料更有效,联合同化径向速度和雨水有利于提高同化分析效果.协方差对EnKF同化效果起着非常重要的作用,考虑模式全部预报变量与径向速度协方差的同化效果比仅考虑速度场与径向速度协方差的同化效果好.雷达资料缺值降低了同化效果,此时增加地面常规观测资料的同化可以明显提高同化分析效果.EnKF同化技术对雷达观测资料误差不太敏感.初始集合对同化分析有较大影响.EnKF同化受集合大小和观测资料影响半径.同化对模式误差较敏感.利用EnKF同化双多普勒雷达资料,分析了一次梅雨锋暴雨过程的中尺度结构.结果表明,EnKF同化技术能够从双多普勒雷达资料反演暴雨中尺度系统的动力场、热力场和微物理场,反演的风场是较准确的,反演的热力场和微物理场分布也是基本合理的.中低层切变线是此次暴雨的主要动力特征,对流云表现为低层辐合、高层辐散并有垂直上升运动伴随,其热力特征表现为低层是低压区,高层为高压区,中部为暖区而上、下部为冷区,水汽、云水和雨水分别集中在对流云体内、上升气流区和强回波区.