基于VDRAS的快速更新雷达四维变分分析系统

基于雷达资料快速更新四维变分同化 (RR4DVar) 技术和三维数值云模式,初步研发了一个针对对流尺度数值模拟的快速更新雷达四维变分分析系统。系统通过对京津冀6部多普勒天气雷达资料进行RR4DVar同化,并融合5 min自动气象站观测和中尺度数值模式结果,可快速分析得到12~18 min更新的低层大气三维动力、热力场的对流尺度结构特征。针对2009年7月22日发生在京津冀的一次强风暴个例,通过一系列敏感性试验,并利用局地加密资料进行检验对比,表明有效的雷达资料RR4DVar同化及自动气象站和中尺度模式资料融合方案、恰当的中尺度背景场设置和动力约束方法是获得合理结果的关键。研究也表明:恰当的系统配置能够模拟出与对流生消发展密切相关的近风暴环境特征,包括低层入流、垂直风切变、低层辐合上升和暖舌等,以及风暴自身形成的冷池、出流等与风暴演变密切相关的对流尺度结构。     

Dynamical and Microphysical Retrieval from Simulated Doppler Radar Observations Using the 4DVAR Assimilation Technique

Based on a cloud model and the four-dimensional variational (4DVAR) data assimilation method developed by Sun and Crook (1997), simulated experiments of dynamical and microphysical retrieval from Doppler radar data were performed. The 4DVAR data assimilation technique was applied to a cloud scale model with a warm rain parameterization scheme. The 3D wind, thermodynamical, and microphysical fields were determined by minimizing a cost function, defined by the difference between both radar observed radial velocities and reflectivities and their model predictions. The adjoint of the numerical model was used to provide the gradient of the cost function with respect to the control variables. Experiments have demonstrated that the 4DVAR assimilation method is able to retrieve the detailed structure of wind, thermodynamics, and microphysics by using either dual-Doppler or single-Doppler information. The quality of retrieval depends strongly on the magnitude of constraint with respect to the variables. Retrieving the temperature field, cloud water and water vapor is more difficult than the recovery of the wind field and rainwater. Accurate thermodynamic retrieval requires a longer assimilation period. The inclusion of a background term, even mean fields from a single sounding, helped reduce the retrieval errors. Less accurate velocity fields were obtained when single-Doppler data were used. It was found that the retrieved velocity is sensitive to the location of the retrieval domain relative to the radars while the other fields have very little changes. Two radar volumetric scans are generally adequate for providing the evolution, although the use of additional volumes improves the retrieval. As the amount of the observations decreases, the performance of the retrieval is degraded. However, the missing observations can be compensated by adding a background term to the cost function. The technique is robust to random errors in radial velocity and calibration errors in reflectivity. The boundary conditions from the dual-Doppler synthesized winds are sufficient for the retrieval. When the retrieval is mainly controlled by the observations in the regions away from the boundaries, the simple boundary conditions from velocity azimuth display (VAD) analysis are also available. The microphysical retrieval is sensitive to model errors.     

多普勒雷达四维变分分析系统概述

《热带气象学报》是中国气象局广州热带海洋气象研究所主办的学术期刊,由气象出版社出版;其前身是《热带气象》,创刊于1984年,已经连续出版20年了。《热带气象学报》主要刊登热带大气科学的基础理论研究及应用研究方面的论文、热带地区天气预报的先进方法和成功的经验、推动学科发展的综合评述和学术动态。当前一段时间的刊登重点是:海气相互作用,中低纬大气环流相互作用,低频振荡及遥相关,低纬大气环流异常及其机制,热带大气环流异常对全球天气、气候的影响以及季风动力学、热带气旋动力学与运动学等的学术成果。由于国家重视气象科学研究,…     

多普勒雷达资料4DVAR同化反演的模拟研究

利用Sun等建立的同化模式和四维变分同化方法对多普勒雷达资料反演大气风场、热力场和微物理场进行了模拟试验研究.反演的基本思路是将4DVAR同化方法应用到三维云模式,定义价值函数表征雷达资料与模式预报结果之间的差别,通过极小化价值函数得到反演场,价值函数相对模式控制变量的梯度由伴随模式求取.试验结果表明,4DVAR同化技术能够从单(双)多普勒雷达资料反演大气三维风场、热力场和微物理场.各个变量反演精度高低与同化过程中变量受约束的大小程度呈正相关.速度场和雨水场反演精度较高,温度场、云水和水汽的反演精度次之,温度场的准确反演需要较长的同化时间.价值函数中加入背景场,哪怕是单点探空给出的平均场信息也有利于提高反演精度.在采用单部多普勒雷达资料进行反演时,速度场的反演误差较大.反演区相对雷达站的位置变化对速度场反演结果有一定的影响,而对其他变量的反演影响很小.两个时次的雷达观测资料基本足够提供反演所需的时间演变信息,同化更多时次的雷达资料,反演效果改进很小.雷达观测资料的缺值会显著降低同化效果,甚至可能导致同化失败,引入背景场可以改善这一状况.4DVAR同化技术对于雷达观测资料误差不太敏感.利用双多普勒雷达合成风场提供水平风场边界条件是比较准确可靠的.在反演主体离边界较远时,VAD风场也基本可用作水平风场边界条件.微物理场的反演对模式中的微物理参数化方案较敏感.     

基于云依赖背景场误差协方差的雷达资料同化及对降雨预报影响研究

传统变分同化方法中使用各向同性和均质的背景场误差协方差,忽略了背景场误差协方差的天气系统依赖性,而在变分框架下引入集合流依赖的背景场误差协方差还需要额外的集合预报.为在变分同化中引入更合理的背景场误差协方差,通过引入云指数构建"云依赖"背景场误差协方差,提出了一种云依赖背景场误差协方差的同化方案,并应用于雷达等多源观测...     

Microphysical Retrieval from Doppler Radar Reflectivity Using Variational Data Assimilation

One of the microphysical variables, the rainwater mixing ratio qr, is retrieved from the observed reflectivity of Doppler radar by a 3D variational data assimilation system. The qr as an analysis variable is obtained by minimizing a cost function defined as the difference between observed radar reflectivity and its retrieval from qr, plus the difference between qr and its background field from a mesoscale model's prediction. Covariance matrix of the background field's error is determined by the so-called NMC method. A method called the second-order auto-regression (SOAR) is used to calculate the coefficients of regressive filtering to fit in with small spatial scale such as cumulus in the process of spatial transformation. An ideal experiment demonstrates the correctness of this system and a sensitivity experiment proves that the random error of observed reflectivity has effect on the analyzed results. At last an experiment with observed data from the Doppler radar at Ma'anshan City in Anhui Province on 19 June 2002 was performed. The retrieved analysis variable qr in this test shows structures in detail, which coincide with the distribution of the echo picture observed by the radar.     

A two-step variational method for three-dimensional wind retrieval from single Doppler radar

Summary A two-step variational method is developed to retrieve the three-dimension wind field from single-Doppler radar observation for operational applications. The first step produces a proxy background wind field, and the second step retrieves the detailed wind structures in data dense areas on a high-resolution grid. The background field is called proxy because it is obtained in the first step by a spectral analysis of the same radial-wind observations that will be used in the second step and thus the background is not independent of the observations. In the second step, the wind retrieval is constrained by the proxy background and the advection equation of liquid water content in the cost function. Tested with a dataset collected during Meiyu front heavy rainfall on 22 July 2002 in China, the two-step approach is shown to yield the best performance in comparison with a single-step approach that uses either only the proxy background field or the retrieval without the proxy background.     

Three-Dimensional Wind Field Retrieved from Dual-Doppler Radar Based on a Variational Method: Refinement of Vertical Velocity Estimates

In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.     

风暴的多普勒雷达自动识别

3种基于雷达的风暴自动识别方法:(1)美国WSR-88D Build 7.0风暴算法,它利用多个预设阈值来检验回波的强度和连续性,以构造具有三维连续结构的风暴,该方法在风暴合并、分裂以及多个单体相距较近时误差较大。(2)为美国WSR-88D Biuld 9.0风暴算法(B9SI),它用7个反射率因子识别阈值替代此前唯一的一个反射率因子阈值,增加了特征核抽取和相近单体处理技术,并保留远距离上的强的2D分量。该方法在面对成串或成簇多单体时,能够识别出多个单体核,并准确定位。B9SI没有考虑反射率因子纹理结构和空间梯度的变化,也没有利用径向速度资料,因此无法描述风暴对流的发展状况。(3)CSI方法,它在降低B9SI反射率因子识别阈值的基础上,利用模糊逻辑技术对B9SI输出结果和雷达基资料做进一步的处理,以计算描述风暴对流发展强弱的对流指数。CSI首先提取一组最能描述风暴对流性特征的物理量,包括反射率因子纹理结构、反射率因子空间变化率、垂直积分含水量和径向速度标准方差,并分配权重;其次,利用每一个物理量的统计结果,结合其物理意义,设计出相应的隶属函数,以计算风暴与该物理量描述的对流性特征相匹配的概率;最后对多个概率值进行加权平均即得对流指数。此外,计算了2004年8月11日发生在广州的超级单体演变过程中的对流指数,分析表明:对流指数两次加大对应了超级单体的合并增长和辐合增长过程;风暴最强盛时对流指数为0.744;随后对流指数减小,雷达观测到的最大反射率因子对应高度明显降低,地面上开始出现大范围的强降水。     

An improved method for doppler wind and thermodynamic retrievals

A variational method is developed to retrieve winds in the first step and then thermodynamic fields in the second step from Doppler radar observations. In the first step, wind fields are retrieved at two time levels: the beginning and ending times of the data assimilation period, simultaneously from two successive volume scans by using the weak form constraints provided by the mass continuity and vorticity equations.As the retrieved wind fields are expressed by Legendre polynomial expansions at the beginning and ending times, the time tendency term in the vorticity equation can be conveniently formulated, and the retrieved winds can be compared with the radar observed radial winds in the cost function at the precise time and position of each radar beam. In the second step, the perturbation pressure and temperature fields at the middle time are then derived from the retrieved wind fields and the velocity time tendency by using the weak form constraints provided by the three momentum equations. The merits of the new method are demonstrated by numerical experiments with simulated radar observations and compared with the traditional least squares methods which consider neither the precise observation times and positions nor the velocity time tendency. The new method is also applied to real radar data for a heavy rainfall event during the 2001 Meiyu season in China.     

Assimilating surface observations in a four-dimensional variational Doppler radar data assimilation system to improve the analysis and forecast of a squall line case

This paper examines how assimilating surface observations can improve the analysis and forecast ability of a fourdimensional Variational Doppler Radar Analysis System(VDRAS).Observed surface temperature and winds are assimilated together with radar radial velocity and reflectivity into a convection-permitting model using the VDRAS four-dimensional variational(4DVAR) data assimilation system.A squall-line case observed during a field campaign is selected to investigate the performance of the technique.A single observation experiment shows that assimilating surface observations can influence the analyzed fields in both the horizontal and vertical directions.The surface-based cold pool,divergence and gust front of the squall line are all strengthened through the assimilation of the single surface observation.Three experiments—assimilating radar data only,assimilating radar data with surface data blended in a mesoscale background,and assimilating both radar and surface observations with a 4DVAR cost function—are conducted to examine the impact of the surface data assimilation.Independent surface and wind profiler observations are used for verification.The result shows that the analysis and forecast are improved when surface observations are assimilated in addition to radar observations.It is also shown that the additional surface data can help improve the analysis and forecast at low levels.Surface and low-level features of the squall line—including the surface warm inflow,cold pool,gust front,and low-level wind—are much closer to the observations after assimilating the surface data in VDRAS.     

Synchronization of Radar Observations with Multi-Scale Storm Tracking

The 3-D radar reflectivity data has become increasingly important for use in data assimilation towards convective scale numerical weather prediction as well as next generation precipitation estimation. Typically, reflectivity data from multiple radars are objectively analyzed and mosaiced onto a regional 3-D Cartesian grid prior to being assimilated into the models. One of the scientific issues associated with the mosaic of multi-radar observations is the synchronization of all the observations. Since radar data is usually rapidly updated (～every 5--10 min), it is common in current multi-radar mosaic techniques to combine multiple radar' observations within a time window by assuming that the storms are steady within the window. The assumption holds well for slow evolving precipitation systems, but for fast evolving convective storms, this assumption may be violated and the mosaic of radar observations at different times may result in inaccurate storm structure depictions. This study investigates the impact of synchronization on storm structures in multiple radar data analyses using a multi-scale storm tracking algorithm.     

应用雷达回波强度资料反演大气云微物理量

文中利用三维变分原理 ,由多普勒天气雷达回波强度反演模式大气的雨水混合比 qr,分析变量 qr 的确定取决于对目标函数的极小化计算 ,目标函数由分析变量 qr 与 qr 的背景场之差 ,以及由分析变量演算的大气回波强度与雷达观测的回波强度之差构成。采用所谓的NMC方法求取模式变量的背景场误差协方差矩阵 ,采用二阶自回归型相关函数确定空间水平变换中的递归滤波系数以适用于反映积云结构的小尺度空间。通过理想数值试验 ,验证了变分同化系统和结果的合理性 ,由敏感性试验说明了雷达观测的回波强度随机误差对于反演的qr 值有一定的影响。对 2 0 0 2年 6月 19日安徽省马鞍山市一次多普勒雷达实际观测数据的个例试验表明 ,经过变分同化反演后的积云尺度模式大气的 qr 场与观测的雷达回波图像的细节相吻合。     

Doppler Radar Data Assimilation with a Local SVD-En3DVar Method

An observation localization scheme is introduced into an ensemble-based three-dimensional variational (3DVar) assimilation method based on the singular value decomposition technique (SVD-En3DVar) to improve assimilation skill. A point-by-point analysis technique is adopted in which the weight of each observation decreases with increasing distance between the analysis point and the observation point. A set of numerical experiments, in which simulated Doppler radar data are assimilated into the Weather Research and Forecasting (WRF) model, is designed to test the scheme. The results are compared with those obtained using the original global and local patch schemes in SVD-En3DVar, neither of which includes this type of observation localization. The observation localization scheme not only eliminates spurious analysis increments in areas of missing data, but also avoids the discontinuous analysis fields that arise from the local patch scheme. The new scheme provides better analysis fields and a more reasonable short-range rainfall forecast than the original schemes. Additional forecast experiments that assimilate real data from 10 radars indicate that the short-term precipitation forecast skill can be improved by assimilating radar data and the observation localization scheme provides a better forecast than the other two schemes.     

2004年两次对流风暴的多普勒雷达资料对比分析研究

误码用常规天气资料、卫星云图、多普勒雷达资料对2004年6月20日、2004年6月24日在河北境内发生的两次风暴天气进行了对比分析,分析表明：（1）速度场的涡旋分析与M产品的结合,可以从雷达资料中获得局地强风和冰雹预报的参考信息;（2）云图中圆形的云团强弱和雷达垂直积分液态含水量的大小,能够反映风暴的强弱。     

Establishment of a multiple-Doppler radar wind retrieval system and its application to the analysis of mesoscale convective systems

A multiple-Doppler radar wind retrieval system is established by using a three dimensional variation method. The system consists of two parts. One is to interpolate reflectivities and radial velocities of nine Doppler radars of the Korea Meteorological Administration (KMA) into one analysis domain. The other is to retrieve three dimensional winds by minimizing a cost function that includes the following costs and constraints: the observed radial velocity cost, background wind cost, continuity constraint and smoothness constraint. In order to verify the performance of the system, retrieved winds are compared with observed winds obtained from five wind profilers of KMA. The performance of the system depends on the relative position to the baselines between Doppler radars. However, the performance of the system is enhanced when the number of overlaps among the radial velocities increases. The system is applied to the analysis of the evolution of a mesoscale convective system (MCS) on the Changma front on 1 July 2005. The analysis result shows that a new convective cell is developed by the convergence of the low troposphere winds at the organizing stage. The analysis of the vertical vorticity reveals that, among the two vorticity generation terms to be calculated utilizing the retrieved winds, tilting or twisting source dominates the divergence source in most convective regions. The strong downdrafts associated with the storm are produced on the downdraft branch of a meridional direction secondary circulation across the Changma front.     

RETRIEVAL OF HORIZONTAL WIND PERTURBATION FIELDS FROM SIMULATED SINGLE DOPPLER RADAR OBSERVATIONS*

The application of the single Doppler radar dataset analysis is usually confined to the assumption that the actual wind is linearly distributed or uniform locally.Following some dynamic features of convective weather,a conceptual model of moderate complexity is constructed,wherewith a horizontal wind perturbation field is retrieved directly from the single Doppler radar measurements.The numerical experiments are based on a 3-D cloud model-generated convective cell,whose radial velocity component is taken as the radar observations that are put into the closed equations based on the conceptual model to retrieve the horizontal wind perturbation field.After the initial field is properly treated,the retrieval equation is solved in terms of the 2-D FFT technique and the sensitivity to noise is examined.Finally,contrast analysis is done of the retrieved and the cloud model output wind fields,indicating the usefulness of the approach proposed in this paper.     

RETRIEVAL OF HORIZONTAL WIND PERTURBATION FIELDS FROM SIMULATED SINGLE DOPPLER RADAR OBSERVATIONS

The application of the single Doppler radar dataset analysis is usually confined to the assumption that the actualwind is linearly distributed or uniform locally.Following some dynamic features of convective weather,a conceptualmodel of moderate complexity is constructed,wherewith a horizontal wind perturbation field is retrieved directly fromthe single Doppler radar measurements.The numerical experiments are based on a 3-D cloud model-generatedconvective cell,whose radial velocity component is taken as the radar observations that are put into the closed equationsbased on the conceptual model to retrieve the horizontal wind perturbation field.After the initial field is properlytreated,the retrieval equation is solved in terms of the 2-D FFT technique and the sensitivity to noise is examined.Finally,contrast analysis is done of the retrieved and the cloud model output wind fields,indicating the usefulness of theapproach proposed in this paper.     

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

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

新疆地区一次对流性降水的三维中尺度风场研究

利用2004年外场试验获得的乌鲁木齐和五家渠C波段双多普勒雷达资料, 分析了双多普勒雷达风场反演方法和资料的可靠性, 研究了2004年8月8日发生在乌鲁木齐和五家渠的一次强对流性降水的回波和风场中尺度结构及演变过程。结果表明:这两部雷达观测的回波强度相关很好, 雷达基线上的径向速度基本一致, 资料可靠, 适合进行双多普勒雷达观测; Cressman插值的影响半径的变化对风场的中尺度结构基本没有影响, 径向速度误差引起的风场反演误差与该点所处的位置有关, 1 m/s径向速度误差也不会改变风场的中尺度结构。该过程为对流单体发展为对流带状回波的过程, 在对流单体的左侧生成新的对流单体, 逐步发展为长度约90 km范围的带状对流系统, 该系统恰与较强的东北风和较弱的西风形成的辐合相对应, 上升气流与强对流回波相对应, 不同对流单体有各自相独立的风场结构。用双多普勒雷达观测得到对流系统的内部风场有利于了解对流系统的内部动力过程, 从而探讨降水的形成和演变机理。