笛卡儿坐标系的双多普勒天气雷达三维风场反演技术

文中研究了笛卡尔坐标系下双多普勒天气雷达三维风场反演技术 ,提出了包括雷达原始资料的预处理 ,空间插值以及可靠性检验的新方法 ,提高了反演结果的可信度和精确度。使用模拟的双多普勒雷达体扫资料进行了反演试验 ,结果表明 :本文的方法能够比较真实地反映风场的三维结构 ,可以用于真实风场的反演。     

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.     

双多普勒雷达反演强风暴三维风场的数值试验

两部雷达测得的径向速度再加上质量连续方程可以求解大气风场，利用这一原理，以1996年6月29日发生在北京东北部京冀交界处的一次强风暴过程为例，用模式反演出风暴体的三维风场结构，误差分析表明和实况基本一致。利用模式输出的三维风场，先模拟两部多普勒雷达扫描的数据，在此基础上，进行三维风场反演，结果表明，反演得到的流场和模拟的流场总体趋势完全一致，风暴的中心位置吻合较好，上升气流和下沉气流的分布也很接近，水平速度分量反演的误差很小，基本上可以反映三维风场的真实情况。垂直速度的平均离差在各个高度上都较小，反演结果较好。     

ERROR COMPARISON OF WIND FIELD RETRIEVAL FROM SINGLE AND DUAL-DOPPLER RADAR OBSERVATIONS

The error distributions of the wind fields retrieved from single and dual-Doppler radar observations are given inthis paper.The results indicate that the error of dual-Doppler retrieval depends on the position in the scan region of thedual-Doppler radar.The error of single-Doppler retrieval by using velocity azimuth processing(VAP)technique de-pends on the angle between the directions of wind and the radar beam.Generally,the winds retrieved from single Doppl-er radar are close to those retrieved from dual-Doppler radar.But,the error distribution of the single-Doppler retrievalis different from the dual-Doppler retrieval.We simulate the retrievals of single Doppler observation by the use of theoutput wind data from a 3-D numerical model of severe convection.The comparison of the simulated single-anddual-Doppler retrievals shows that the VAP may be a suitable technique for the operational analysis of mesoscale windfields.It can also be used as a supplement to wind field retrieval in the field experiment.     

双多普勒雷达风场反演误差和资料的质量控制

提出了用双多普勒雷达观测资料进行回波强度、径向速度和方位定位的质量控制方法，以及利用概率分布法订正配对的双多普勒雷达回波强度的方法。以2001年973“中国暴雨”外场试验期间获取的双多普勒雷达观测资料为例，比较了合肥、马鞍山和宜昌、荆州的两对双多普勒雷达同步观测的回波强度、径向速度和方位的改变对观测数据对比的影响；分析了这两对双多普勒雷达径向速度测量误差引起的风场反演的误差。结果表明:两对双多普勒雷达观测的回波强度和径向速度的空间位置和变化趋势比较一致，合肥与马鞍山雷达的回波强度有一定差异，径向速度也有1～2 m/s的差异。风场反演的误差与风场的方向、大小、空间位置等有关。在两个径向速度夹角在40°～140°范围内，共面上的风场的反演误差在1～2倍的雷达探测径向速度的误差范围内。     

一次强烈雹暴的三维结构和形成机制的单、双多普勒雷达分析

利用单多普勒天气雷达的反射率因子、径向速度数据和双多普勒雷达反演的三维风场,分析了一次强烈雹暴的产生、发展和维持机制.主要结果为:该雹暴是产生于中等偏上垂直风切变和较大的对流有效位能条件下的右移风暴,它在旧无序多单体风暴的右侧产生.旧的多单体风暴和雹暴初始阶段主要雷达回波区域负的径向速度(向着雷达的速度)占主导地位,随...     

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

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

Radar data assimilation of the GRAPES model and experimental results in a typhoon case

Constructing β-mesoscale weather systems in initial fields remains a challenging problem in a mesoscale numerical weather prediction (NWP) model. Without vertical velocity matching the β-mesoscale weather system, convection activities would be suppressed by downdraft and cooling caused by precipitating hydrometeors. In this study, a method, basing on the three-dimensional variational (3DVAR) assimilation technique, was developed to obtain reasonable structures of β-mesoscale weather systems by assimilating radar data in a next-generation NWP system named GRAPES (the Global and Regional Assimilation and Prediction System) of China. Single-point testing indicated that assimilating radial wind significantly improved the horizontal wind but had little effect on the vertical velocity, while assimilating the retrieved vertical velocity (taking Richardson's equation as the observational operator) can greatly improve the vertical motion. Experiments on a typhoon show that assimilation of the radial wind data can greatly improve the prediction of the typhoon track, and can ameliorate precipitation to some extent. Assimilating the retrieved vertical velocity and rainwater mixing ratio, and adjusting water vapor and cloud water mixing ratio in the initial fields simultaneously, can significantly improve the tropical cyclone rainfall forecast but has little effect on typhoon path. Joint assimilating these three kinds of radar data gets the best results. Taking into account the scale of different weather systems and representation of observational data, data quality control, error setting of background field and observation data are still requiring further in-depth study.     

多普勒雷达数据三维格点化的方法研究

天气雷达数据三维格点化是天气雷达组网拼图的基础,本文研究了几种把球坐标系下的分辨率不均匀的雷达反射率资料（dBZ）插值到统一的分辨率均匀的笛卡尔坐标系下的方法,并采用模拟的天气雷达回波数据对这几种方法进行了比较分析,主要对空间反射率强度值的连续性及是否保留了雷达资料中原始的反射率结构特征进行对比分析。      

ANALYSIS OF THE EFFECT OF 3DVAR AND ENSRF DIRECT ASSIMILATION OF RADAR DATA ON THE FORECAST OF A HEAVY RAINFALL EVENT

The present study designs experiments on the direct assimilation of radial velocity and reflectivity data collected by an S-band Doppler weather radar (CINRAD WSR-98D) at the Hefei Station and the reanalysis data produced by the United States National Centers for Environmental Prediction using the Weather Research and Forecasting (WRF) model, the WRF model with a three-dimensional variational (3DVAR) data assimilation system and the WRF model with an ensemble square root filter (EnSRF) data assimilation system. In addition, the present study analyzes a Meiyu front heavy rainfall process that occurred in the Yangtze -Huaihe River Basin from July 4 to July 5, 2003, through numerical simulation. The results show the following. (1) The assimilation of the radar radial velocity data can increase the perturbations in the low-altitude atmosphere over the heavy rainfall region, enhance the convective activities and reduce excessive simulated precipitation. (2) The 3DVAR assimilation method significantly adjusts the horizontal wind field. The assimilation of the reflectivity data improves the microphysical quantities and dynamic fields in the model. In addition, the assimilation of the radial velocity and reflectivity data can better adjust the wind fields and improve the intensity and location of the simulated radar echo bands. (3) The EnSRF assimilation method can assimilate more small-scale wind field information into the model. The assimilation of the reflectivity data alone can relatively accurately forecast the rainfall centers. In addition, the assimilation of the radial velocity and reflectivity data can improve the location of the simulated radar echo bands. (4) The use of the 3DVAR and EnSRF assimilation methods to assimilate the radar radial velocity and reflectivity data can improve the forecast of precipitation, rain-band areal coverage and the center location and intensity of precipitation.     

用EVAP方法反演冷流暴雪的风场及检验

EVAP(Extended VAP)雷达反演方法是在忽略空气垂直运动速度而考虑云雨粒子下落速度影响的条件下,导出从雷达多仰角多普勒速度资料产生不同高度上的水平径向速度计算公式,然后利用VAP方法来反演出等高面上的二维水平风场.利用烟台多普勒天气雷达资料,采用EVAP方法反演了2005年12月山东半岛两个风场背景条件和暴雪落区均有明显差异的冷流暴雪个例水平风场.将雷达实测的径向速度、反射率因子插值到各等高面上,并与反演的同等高面风场叠加显示,以此分析反演的风场是否与原始雷达数据相符.此外,还与地面自动站风场、高空风场等实测资料进行对比,综合检验反演结果的合理性和可靠性,并进一步分析了冷流暴雪的中尺度风场结构.结果表明,反演风场与实测资料相符;EVAP方法对强风和弱风的冷流暴雪个例都可反演出不同高度的中尺度风场结构,反演的西南风与东北风、西北风与东北风之间的切变线和强回波带相对应,揭示了暴雪产生的动力机制,说明EVAP方法对不同风场条件的暴雪中尺度系统都有一定的反演能力.高时空分辨率的EVAP反演风场弥补了天气图上仅能分析出山东半岛西北风的缺憾,且反演过程中一个体扫的计算时间不足1 min,因此EVAP反演的中尺度风场信息适合在冷流降雪的短时临近预报业务中应用.     

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.     

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.     

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.     

On the application of MUSCAT to a ground-baseddual-Doppler radar system

Summary The multiple-Doppler synthesis and continuity adjustment technique (MUSCAT), developed for at least a pair of Doppler observations, is based on both a simultaneous resolution of the equations relating the wind components to the measured radial velocities and of the continuity equation, and the use of second-order derivatives in the cost function. It has been developed to solve the geometrical limitations that result from a dual-Doppler radar configuration, and has been successfully applied to airborne radar equipped for a fore-aft scanning. This paper intends to show the performance of this method in the case of a differing geometry from a ground-based dual-Doppler radar system, using simulated Doppler data deduced from an analytical wind field. It is shown that the constraint on the second-order derivatives allows to also obtain reliable wind vectors outside the Doppler lobes, i.e., in regions along the radar baseline where ill-conditioned analysis occurs. Received February 14, 2000/Revised April 4, 2001     

强热带风暴碧利斯 (0604) 引发的特大暴雨中尺度结构多普勒雷达资料分析

强热带风暴 “碧利斯” (Bilis) 于2006年7月14日12:50在福建省霞浦县北壁镇再次登陆, 与南海季风相互作用, 在福建省引发特大暴雨。作者利用双多普勒雷达三维风场反演技术对厦门和龙岩新一代多普勒雷达时间同步探测资料进行了风场反演, 综合利用雷达回波强度资料, 对造成长泰、 漳州特大暴雨的中尺度对流系统的三维结构及其演变特征进行了详细分析。结果表明: (1) 此次特大暴雨主要是由中低层西南[CD*2]东北走向的β中尺度辐合线引发的, 辐合线对于水汽输送以及暴雨的形成、 触发、 维持具有重要作用, 辐合线在暴雨的整个生命史过程中经历了由弱变强、 由强变弱的演变过程, 变化过程与降水的强弱演变过程基本同步。 (2) 由于丰富的水汽供应和中低层辐合线的动力抬升作用, 西南[CD*2]东北走向的β中尺度回波带的西南不断有新的γ中尺度对流单体生成, 在沿着辐合线向东北移动过程中进一步发展、 合并形成β中尺度对流线, 造成持续的强降水。最后, 还给出了此次特大暴雨的三维云系结构模型。     

一次强雹云结构的双多普勒雷达观测分析

利用石家庄SA多普勒天气雷达和饶阳X波段双偏振雷达探测资料，结合常规综合观测资料，对2018年6月13日下午发生于太行山东麓的雹云的天气背景、降雹特征、雷达回波演变特征及回波三维立体结构进行了综合分析，重点利用双多普勒雷达径向速度资料反演出格点的三维风速（流场），并结合回波特征分析了雹云云体结构。结果显示，高空强劲的偏北风急流促使涡后横槽转竖，槽后冷空气沿偏北气流南下，形成上干冷、下暖湿的不稳定层结，在低空低涡附近及地面辐合线上发展造成这次风雹天气；双多普勒雷达反演风场表明，雹云的中层强回波中心呈明显的“S”形水平流场和悬挂回波配置特征，构成了具有成雹的“0线”结构；不仅佐证了雹云降雹“0线”结构的存在，而且证明其呈现形式具有多样性。      

A Model Study of Three-Dimensional Wind Field Analysis from Dual-Doppler Radar Data

ＡＭｏｄｅｌＳｔｕｄｙｏｆＴｈｒｅｅ－ＤｉｍｅｎｓｉｏｎａｌＷｉｎｄＦｉｅｌｄＡｎａｌｙｓｉｓｆｒｏｍＤｕａｌ－ＤｏｐｐｌｅｒＲａｄａｒＤａｔａＫｏｎｇＦａｎｙｏｕ（孔凡铀）；ＭａｏＪｉｅｔａｉ（毛节泰）（ＤｅｐａｒｔｍｅｎｔｏｆＧｅｏｐｈｙｓｉｃｓ...     

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.     

利用单多普勒雷达资料做冷流暴雪的中尺度分析

基于单多普勒天气雷达资料采用EVAP方法反演风场, 并结合径向速度、 反射率因子、 自动气象站和探空风场等观测资料, 对2005年12月6～7日山东半岛一次冷流暴雪过程的中尺度特征进行了深入分析, 结果表明： (1)雷达回波呈狭窄带状, 移动缓慢近乎停滞, 径向速度上存在风向切变线, 烟台和威海的暴雪不同步是冷流暴雪的典型特征; (2)首次通过雷达反演证实了逆风区实际就是风场切变在径向速度图上的反映, 垂直各层水平风场存在中尺度切变线, 且与强回波带相对应, 切变线的位置决定暴雪的落区; (3)通过雷达反演风场和风廓线共同揭示出强降雪产生时对流层中层有西北风、 西南风和东北风三股气流, 明显的西南气流位于850～700 hPa, 表明冷流降雪过程并非传统认为的仅有西北冷平流, 而是不同气流辐合的结果; (4)对流层中层的西南暖平流为云的播种和反馈机制提供了有利的天气背景条件, 使得冷流降雪增强, 这在常规观测资料中无法看到。