双线偏振雷达资料在数值模式中的应用: 模拟器的构建

基于Rayleigh-Gans散射原理构造一个S波段双线偏振雷达模拟器。模拟器考虑了云冰、雪、雨和雹4种水凝物,以水凝物的混合比和数浓度以及水凝物粒子的轴比、相对介电常数、下落倾角为输入量,计算得到水平/垂直偏振反射率因子、差分反射率、比差分相位等偏振量。通过二维理想飑线模拟的试验结果表明,模拟器合理地再现了二维理想飑线系统成熟期的主要偏振特征:雹的反射率较高,差分反射率较低(又称Zdr洞);对流云区的Zdr柱;层云区的反射率和比差分相位的0℃层亮带特征;雨滴反射率与差分反射率因子的“雨线”统计特征。利用该模拟器建立了模式变量和偏振雷达观测的联系,有助于未来将偏振雷达观测资料应用于模式预报效果评估及对流尺度资料同化等方面的研究和应用。      

Parameterized Forward Operators for Simulation and Assimilation of Polarimetric Radar Data with Numerical Weather Predictions

Many weather radar networks in the world have now provided polarimetric radar data(PRD)that have the potential to improve our understanding of cloud and precipitation microphysics,and numerical weather prediction(NWP).To realize this potential,an accurate and efficient set of polarimetric observation operators are needed to simulate and assimilate the PRD with an NWP model for an accurate analysis of the model state variables.For this purpose,a set of parameterized observation operators are developed to simulate and assimilate polarimetric radar data from NWP model-predicted hydrometeor mixing ratios and number concentrations of rain,snow,hail,and graupel.The polarimetric radar variables are calculated based on the T-matrix calculation of wave scattering and integrations of the scattering weighted by the particle size distribution.The calculated polarimetric variables are then fitted to simple functions of water content and volumeweighted mean diameter of the hydrometeor particle size distribution.The parameterized PRD operators are applied to an ideal case and a real case predicted by the Weather Research and Forecasting(WRF)model to have simulated PRD,which are compared with existing operators and real observations to show their validity and applicability.The new PRD operators use less than one percent of the computing time of the old operators to complete the same simulations,making it efficient in PRD simulation and assimilation usage.     

WRF模式云参数化方案对一次深对流系统模拟的验证和改进

使用中尺度数值模式WRF中的双参数云微物理方案WDM6针对2008年台风“凤凰”登陆过程中造成的强降水进行数值模拟,通过卫星模拟器利用MTSAT-1R和TRMM卫星观测的红外云顶黑体亮温TBB、PR雷达反射率资料使用统计方法验证模拟结果。通过修改云水向雨水自动转化过程、冰晶核化过程、雪和霰的下落末速度、雪和霰的截距进行敏感性试验,减小模拟结果和卫星观测结果的差异。研究结果表明:WDM6方案模拟的台风“凤凰”登陆后的降水,强对流云系及对流柱状雷达回波基本符合实况,但模拟结果局部偏强。WDM6方案模拟产生了较多的浅对流云,低估了对流云系的出现频率。不同云类型模拟的雷达回波均偏强,对流云系雷达回波垂直分布接近观测。敏感性试验结果说明修改WDM6方案中云水向雨水自动转化率有效地改善了模拟效果。同时发现云滴初始数浓度影响云水向雨水自动转化率并最终影响云系结构和雷达反射率的模拟结果,过高的云滴初始数浓度会使模拟结果变差。     

三维变分和物理初始化方法相结合同化多普勒雷达资料的试验研究

以美国新近研发的天气研究预报模式(WRF)配置的三维变分(3D-Var)同化系统WRF 3D-Var为平台,结合物理初始化方法(Physical Initialization,简称PI)来同化多普勒雷达径向风和回波强度观测资料.其基本做法是首先用物理初始化方法由雷达回波资料估计出比湿、云水混合比和垂直速度,然后用估计的比湿和云水混合比对模式的相应变量进行调整,最后再将估计出的垂直速度作为一种新的观测类型添加到现有的WRF 3D-Var目标函数中,同时以WRF 3D-Var提供的方法直接同化径向风.针对2002年6月19日的一次强对流性降水过程和2003年7月5日的一次梅雨锋暴雨过程进行了一组同化多普勒雷达径向风和回波资料的试验研究.同化结果表明:分析变量的增量场和观测的雷达回波有很好的对应关系.在雷达回波区,有正的比湿增量、云水含量增量和垂直速度增量,并且水平风增量在此辐合;在没有雷达回波的地方有负的垂直速度增量.预报结果表明,调整云水含量对降水预报改善不明显,调整比湿对降水预报改进明显,直接用物理初始化估计出的垂直速度替代模式的初始垂直速度,对降水预报改进不明显,但以新的方案同化雷达资料能有效地缩短模式的起转时间(spin-up time),明显改进短时降水预报.     

A STUDY ON VARIABLE QUANTITATIVE PRECIPITATION ESTIMATION USING DOPPLER RADAR DATA

With the pros and cons of the traditional optimization and probability pairing methods thoroughly considered, an improved optimal pairing window probability technique is developed using a dynamic relationship between the base reflectivity Z observed by radar and real time precipitation I by rain gauge. Then, the Doppler radar observations of base reflectivity for typhoons Haitang and Matsa in Wenzhou are employed to establish various Z-I relationships, which are subsequently used to estimate hourly precipitation of the two typhoons. Such estimations are calibrated by variational techniques. The results show that there exist significant differences in the Z-I relationships for the typhoons, leading to different typhoon precipitation efficiencies. The typhoon precipitation estimated by applying radar base reflectivity is capable of exhibiting clearly the spiral rain belts and mesoscale cells, and well matches the observed rainfall. Error statistical analyses indicate that the estimated typhoon precipitation is better with variational calibration than the one without. The variational calibration technique is able to maintain the characteristics of the distribution of radar-estimated typhoon precipitation, and to significantly reduce the error of the estimated precipitation in comparison with the observed rainfall.     

QVP方法在双偏振雷达冬季降水观测中的应用

双偏振雷达是研究降水微物理过程的重要探测设备,为研究我国东部沿海地区冬季降水的微物理特性,选取2019—2020年不同天气背景下（包含相态转换）影响上海的冬季降水过程,基于上海南汇WSR-88D双偏振雷达资料,采用准垂直廓线（QVP）方法反演3次降水过程的反射率因子Z_H、差分反射率Z_DR和相关系数ρ_hv的高度-时间廓线。基于QVP,结合探空、再分析资料、地面自动气象站和雨滴谱数据,分析过程时段内融化层高度变化、云中粒子微物理特征,同时借助云雷达观测比对QVP方法的实际效果。结果表明:QVP方法反演的廓线可以体现贝吉龙过程的发生以及融化层高度的变化,并能够有效区分过程时段内的凇附和聚合过程。与此同时,对于非连续性或非均匀性质的降水系统,有针对性地选择方位利用QVP方法进行处理可准确获取重点关注区域降水云系中的微物理过程变化。综上所述,QVP方法反演的高度-时间廓线能够用于降水相态的快速诊断及降水粒子从空中到地面演变的微物理过程分析,为冬季降水相态短时临近预报提供依据。     

面向资料同化的S波段双偏振雷达质量控制

双偏振雷达是强对流天气分析和云微物理研究的重要探测设备,将多普勒天气雷达升级为双偏振雷达是我国未来几年强对流天气监测发展建设的重点计划,而双偏振雷达资料质量控制是其中的关键技术问题。针对国内首批业务布建的广东省S波段双偏振雷达网研发建立了一套完整的质量控制技术方案,除了重点解决非气象回波外,还考虑了非标准波束遮挡和径向高频脉动问题。在飑线、暴雨和台风3类华南季风区强天气个例的应用表明,基于模糊逻辑的水凝物分类、偏振量（零滞后互相关系数、信噪比和比差分相位）阈值检查和杂波剔除能有效剔除非气象回波,抑制异常传播波束导致的虚假回波;线性内插较好地弥补了非标准波束遮挡带来的观测缝隙;中值滤波和滑动平均既过滤了偏振观测量在雷达径向的高频脉动,又保留了主要的偏振雷达观测特征。质量控制后的气象回波约占有效观测（反射率因子大于-30 dBZ）的40%,其偏振量取值分别为反射率因子大于5 dBZ、零滞后互相关系数大于0.8和差分反射率为-0.2~4 dB。     

登陆台风“罗莎”中云物理特征的数值模拟研究

选取2007年16号超强台风“罗莎”为个例,利用耦合了双参数混合相云降水物理方案的 GRAPES中尺度模式对其进行模拟,并结合实况雨量、雷达及卫星资料分析了本次台风暴雨的结构及云微物理特征。结果表明,低空急流和不稳定的温度层结为本次降水过程提供了有利的天气条件。强降水中心主要出现在台风云系的北部和东部,在模拟时段内6小时累积降水最大值超过100 mm。积分固态水含量与液态水相当,说明冰相粒子在形成降水的过程中起重要作用。从CloudSat卫星探测结果来看,强对流发生在台风眼周围的云墙中。根据模拟结果分析,台风东北部的强降水属于对流云降水,最强上升气流出现在0 ℃等温线之下的暖区;西南部的降水强度比较均匀,强上升气流出现在高层冷云中,有利于冰粒子的形成和生长。      

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

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

多普勒雷达资料动态定量估测台风小时降水量的研究

在充分考虑传统的最优化方法和概率配对法优缺点的基础上,使用一种改进的最佳窗概率配对法计算Z-I关系中的系数A和b,得到了雷达测得的基本反射率因子Z和雨量计实时测到的小时降水量I的动态关系.利用温州多普勒雷达体扫资料和浙江省自动雨量站资料,使用该方法对"海棠"(Haitang)和"麦莎"(Matsa)两个台风分别进行了动态计算,得到了不同系数的Z-I关系,进而对两个台风的小时降水量进行了定量估测.使用变分技术对估测的小时降水量进行了校准.结果表明,不同台风Z-I关系的系数差别较大,因而造成台风小时降水量的很大不同.使用雷达基本反射率来估测台风小时降水量,能够清楚表现出台风的螺旋雨带和其中的中小尺度雨团,估测的台风小时降水量与实况基本接近.经过变分校准的估测降水量可以较好地表现出台风雨带与地面中尺度流场动力结构的对应关系.误差统计分析表明,变分校准后的估测台风小时降水量要明显好于变分校准前的估测台风小时降水量.变分校准法既保留了雷达估测台风小时降水量的分布特征,又使估测的台风小时降水雨量与实况的误差明显减小.     

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.     

X波段双线性偏振雷达电磁波衰减对冰雹识别的影响

利用X波段双线性偏振雷达2009年4-10月降雨的ZH-ZDR统计分布,得到基于ZH-ZDR分布的冰雹识别参量(HDR)表达式。雷达发射的电磁波经过雨区导致的衰减会使HDR识别冰雹不准确,需要进行衰减订正。文中利用差分相移(DP)订正了HDR表达式中经过雨区衰减的反射率(ZH)和差分反射率(ZDR)。个例分析表明,经雨区衰减订正后X波段双线性偏振雷达回波与北京南郊观象台CIN-RAD/SA多普勒天气雷达回波具有较高的一致性,在零度层之下,HDR>0的区域能有效地反映出地面降雹情况;在零度层之上,HDR>0的区域可能是冰雹或湿霰,订正后HDR判别冰雹的虚警率有可能提高。     

CINRAD/SAD双偏振雷达非降水回波识别技术

双偏振雷达观测特征参量（如相关系数、差分反射率等）能有效抑制地物、超折射、电磁干扰、海浪和晴空等非降水回波。在上海南汇WSR-88D双偏振雷达非降水回波识别算法基础上,对我国升级布网且纳入业务运行的CINRAD/SAD双偏振雷达数据进行算法测试、算法模块适应性改进,利用江苏、广东的双偏振雷达观测冰雹、融化层、台风降水及各种杂波个例进行算法检验评估,并在组网拼图中展示质量控制效果。结果表明:总体上算法对非降水回波的识别准确率达到95.2%,降水回波的误判率为2.6%。但对夏秋季节夜晚的大面积晴空回波算法识别准确率低于90%,有待尝试利用深度学习方法改进。     

多普勒雷达资料同化对北京“7.21”大暴雨过程模拟的影响

利用新一代中尺度预报模式WRFV3.6及其三维变分同化系统(WRF-3DVAR),对2012年7月21日北京地区的一次暴雨过程进行多普勒天气雷达径向风和反射率的同化试验研究,检验和探讨高时空分辨率多普勒天气雷达资料在改进模式初始场及提高对暴雨过程预报的准确率等方面的应用效果及意义。结果发现雷达资料同化能在初始场中加入反映产生降水的低层风场辐合的动力和锋前暖区充足的水汽条件的物理信息,可以在模式积分开始后改善初始场中水汽和风的分布,较快地模拟出局地对流系统的发生、发展,改善由于中尺度观测资料不足而造成的模式初始场里中尺度信息缺乏的问题。径向速度的同化增加了中尺度信息,对初始流场的调整较为显著,侧重于改进风场。而雷达反射率资料的同化对初始温、湿度场和强回波位置的调整更明显,侧重于改进湿度场。累计降水的预报结果显示,同化径向风资料对雨带的位置、范围有较好的改进,同化雷达反射率资料对暴雨强度的预报有明显的改善。通过降水ETS评分发现,同化常规观测试验相对于控制实验,对于5、15 mm和25 mm降水评分能增加0.1左右,径向风同化试验能增加0.2左右,反射率同化试验能增加0.3左右,而径向风加反射率试验增加的评分介于0.2～0.3。雷达资料对于提高定量降水预报的精确度有着重要作用。     

Current Status and Future Challenges of Weather Radar Polarimetry: Bridging the Gap between Radar Meteorology/Hydrology/Engineering and Numerical Weather Prediction

After decades of research and development, the WSR-88 D(NEXRAD) network in the United States was upgraded with dual-polarization capability, providing polarimetric radar data(PRD) that have the potential to improve weather observations,quantification, forecasting, and warnings. The weather radar networks in China and other countries are also being upgraded with dual-polarization capability. Now, with radar polarimetry technology having matured, and PRD available both nationally and globally, it is important to understand the current status and future challenges and opportunities. The potential impact of PRD has been limited by their oftentimes subjective and empirical use. More importantly, the community has not begun to regularly derive from PRD the state parameters, such as water mixing ratios and number concentrations, used in numerical weather prediction(NWP) models.In this review, we summarize the current status of weather radar polarimetry, discuss the issues and limitations of PRD usage, and explore potential approaches to more efficiently use PRD for quantitative precipitation estimation and forecasting based on statistical retrieval with physical constraints where prior information is used and observation error is included. This approach aligns the observation-based retrievals favored by the radar meteorology community with the model-based analysis of the NWP community. We also examine the challenges and opportunities of polarimetric phased array radar research and development for future weather observation.     

Raindrop Size Distribution Parameters Retrieved from Guangzhou S-band Polarimetric Radar Observations

According to the statistical shape–slope (μ–Λ) relationship observed for the first time by several 2D-Video-Distrometers (2DVD) in southern China, a constrained gamma (C-G) model was proposed for the retrieval of rain drop size distributions (DSDs) from Guangzhou S-band polarimetric radar observations. Two typical precipitation processes were selected to verify the accuracy of the retrieval scheme. The μ–Λ relationship: Λ = 0.0241μ² + 0.867μ + 2.453 was obtained based on the 2DVD observation results from at Huizhou Longmen station, which is a very representative location in the area. Relying on the Guangzhou polarimetric radar measurements of radar reflectivity (Z_HH) and differential reflectivity (Z_DR), the gamma (Γ) size distribution parameters (N₀, μ, and Λ) can be retrieved by the C-G model retrieval scheme. The results show that the Guangzhou polarimetric radar retrievals of DSDs were close to the 2DVD observations at Guangzhou Maofengshan station. The rain rate, mass mean diameter, and normalized intercept parameter of radar retrievals were in good agreement with the 2DVD observations, and the relative errors were less than 10%. The overall accuracy of the retrieval scheme was high. The retrieval scheme has established the relationship between the polarimetric radar measurements and gamma size distribution parameters. It will be helpful to in-depth research and application of the dual-polarization radar data in microphysical precipitation processes analysis, as well as convection-resolved numerical model data assimilation and prediction effect evaluation.     

Raindrop Size Distribution Parameters Retrieved from Xinfeng C-Band Polarimetric Radar Observations

The objective of this research was to acquire a raindrop size distribution (DSDs) retrieved from C-band polarimetric radar observations scheme for the first time in south China.An observation period of the precipitation process was selected,and the shape-slope (μ-Λ) relationship of this region was statistically analyzed using the raindrop sample observations from the two-dimensional video disdrometer (2DVD) at Xinfeng Station,Guangdong Province.Simulated data of the C-band polarimetric radar reflectivity Z_(HH)and differential reflectivity Z_(DR)were obtained through scattering simulation.The simulation data were combined with DSD fitting to determine the Z_(DR)-Λand log_(10)(Z_(HH)/N_0)-Λrelationships.Using Xinfeng C-band polarimetric radar observations Z_(DR)and Z_(HH),the raindrop Gamma size distribution parametersμ,Λ,and N_0were retrieved.A scheme for using C-band polarimetric radar to retrieve the DSDs was developed.This research revealed that during precipitation process,the DSDs obtained using the C-band polarimetric radar retrieval scheme are similar to the 2DVD observations,the precipitation characteristics of rainfall intensity (R),mass-weighted mean diameter (D_m) and intercept parameter (N_w) with time obtained by radar retrieval are basically consistent with the observational results of the 2DVD.This scheme establishes the relationship between the observations of the C-band polarimetric radar and the physical quantities of the numerical model.This method not only can test the prediction of the model data assimilation system on the convective scale and determine error sources,but also can improve the microphysical precipitation processes analysis and radar quantitative precipitation estimation.The present research will facilitate radar data assimilation in the future.     

A Two-Moment Bulk Microphysics Coupled with a Mesoscale Model WRF: Model Description and First Results

The Chinese Academy of Meteorological Sciences (CAMS) two-moment bulk microphysics scheme was adopted in this study to investigate the representation of cloud and precipitation processes under different environmental conditions.The scheme predicts the mixing ratio of water vapor as well as the mixing ratios and number concentrations of cloud droplets,rain,ice,snow,and graupel.A new parameterization approach to simulate heterogeneous droplet activation was developed in this scheme.Furthermore,the improved CAMS scheme was coupled with the Weather Research and Forecasting model (WRF v3.1),which made it possible to simulate the microphysics of clouds and precipitation as well as the cloud-aerosol interactions in selected atmospheric condition.The rain event occurring on 27-28 December 2008 in eastern China was simulated using the CAMS scheme and three sophisticated microphysics schemes in the WRF model.Results showed that the simulated 36-h accumulated precipitations were generally agreed with observation data,and the CAMS scheme performed well in the southern area of the nested domain.The radar reflectivity,the averaged precipitation intensity,and the hydrometeor mixing ratios simulated by the CAMS scheme were generally consistent with those from other microphysics schemes.The hydrometeor number concentrations simulated by the CAMS scheme were also close to the experiential values in stratus clouds.The model results suggest that the CAMS scheme performs reasonably well in describing the microphysics of clouds and precipitation in the mesoscale WRF model.     

石家庄地区反射率因子垂直廓线特征分析

利用自动雨量计数据整理成的10 min一次的雨量资料和s波段多普勒天气雷达体积扫描强度数据,对石家庄地区2004～2007年4次天气过程的实时雷达反射率因子垂直廓线的特征进行了分析.结果表明:层状云和混合性降水反射率因子垂直廓线有明显的零度层亮带;短时强降水过程的反射率因子垂直廓线不存在零度层亮带.冰雹过程中反射率凼子垂直廓线变化较大,降雹前反射率因子的极大值在中上层,降雹发生时反射率因子的极大值高度下降,降雹后反射率因子的极大值减弱.降雪过程的反射率因子垂直廓线零度层亮带不明显.在石家庄西部山区,由于零度层亮带的影响.对层状云和混合性降水回波强度和降水量估计偏高.对短时强降水过程的地面降水估计用反射率因子垂直廓线的方法比最低仰角法更加准确,在均匀性降水中可较好地改善地面雨量估算结果,有利于在山区和无雨量计的地区判断强对流天气的发生、发展和估算降水量的大小.     

X 波段双极化雷达对云中水凝物粒子的相态识别

人工影响天气研究需对云中降水粒子的相态和分布结构进行准确识别,以便提高人工影响天气作业效率.中国科学院大气物理研究所的车载X波段双极化雷达可提供与云中降水粒子大小、形状、相态等特征密切相关的4个极化参数:反射率因子、差分反射率、差分相移率、水平和垂直极化相关系数.利用这4个极化参数加上环境温度作为5个输入参量,建立了降水粒子相态模糊逻辑识别算法,识别的降水粒子有10种:毛毛雨、雨、湿霰、干霰、小雹、大雹、雨加雹、湿雪、干雪、冰晶.利用此雷达的实际观测资料,并与地面和飞机空中实测资料对照,对我国南、北方地区观测的降水天气过程进行分析,结果表明:建立的模糊逻辑算法对云内水凝物粒子的相态识别分类合理.