Study on the Variational Assimilation Technique for the Retrieval of Wind Fields from Doppler Radar Data

This paper introduces a variational assimilation technique for the retrieval of wind fields from Doppler radar data. The assimilated information included both the radial velocity (RV) and the movement of radar echo. In this assimilation technique, the key is transforming the movement of radar echo to a new radar measuring variable- "apparent velocity" (AV). Thus, the information of wind is added, and the indeterminacy of recovering two-dimensional wind only by AV was overcome effectively by combining RV with AV. By means of CMA GRAPES-3Dvar and CINRAD data, some experiments were performed. The results show that the method of retrieval of wind fields is useful in obtaining the construction of the weather system.     

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.     

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

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

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

利用中尺度WRF模式及其3DVAR 同化系统对2014年3 月30 日发生在我国华南地区的一次飑线过程展开多普勒天气雷达资料的同化效果试验研究.首先对雷达资料进行去地物杂波、退速度模糊等预处理,后设计了基于不同雷达观测量的同化试验及同化频次的敏感性试验.结果表明:直接循环同化雷达径向风资料和雷达反射率因子能够增加数值模...     

雷达资料同化对一次浙江初春罕见飑线过程数值模拟的影响分析

针对2019年6月5日重庆地区发生的一次飑线天气过程,利用中尺度WRF (Weather Research and Forecasting)模式和ARPS (The Advanced Regional Prediction System)的3DVAR三维变分同化系统及其ADAS (ARPS Data Assimilation System)云分析系统,探究了不同模式水平分辨率下雷达资料同化对该飑线系统的模拟改进效果。结果表明：(1) 未同化雷达资料时,模式水平分辨率从900 m提高到300 m,模式模拟结果无明显改进。(2) 不同模式水平分辨率下,同化试验的模拟效果相比同化前,对于雷达回波的形状、强度和落区都有一定改善。(3) 使用雷达资料同化的情况下,同时提高模式水平分辨率,能进一步优化调整模式的动力、热力及水汽条件,使得对本次飑线系统的发生发展和组织结构特征的模拟与实况更接近。

     

基于物理约束扰动的EnSRF雷达资料同化

在完善WRF集合均方根滤波同化系统的双向嵌套预报方案、雷达资料质量控制方案和多微物理过程积分方案的基础上,引入了一种扰动势函数、流函数的尺度较小、物理意义更明显的集合初始扰动新方案.通过台风“桑美”登陆过程中的多普勒雷达资料直接同化试验,对比新扰动方案与目前应用较多的背景误差协方差随机扰动和全场随机扰动方案,对集合滤波同化效果的影响.试验结果表明,3种初始扰动方案都有效改善了台风内部结构,但新方案由于尺度更小、结构更为合理,其同化效果也更为显著,提高了对台风路径、中心最低气压、最大风速及降水的预报质量.     

Evaluation of two momentum control variable schemes and their impact on the variational assimilation of radarwind data: Case study of a squall line

Different choices of control variables in variational assimilation can bring about different influences on the analyzed atmospheric state. Based on the WRF model’s three-dimensional variational assimilation system, this study compares the behavior of two momentum control variable options—streamfunction velocity potential (ψ–χ) and horizontal wind components (U–V)—in radar wind data assimilation for a squall line case that occurred in Jiangsu Province on 24 August 2014. The wind increment from the single observation test shows that the ψ–χ control variable scheme produces negative increments in the neighborhood around the observation point because streamfunction and velocity potential preserve integrals of velocity. On the contrary, the U–V control variable scheme objectively reflects the information of the observation itself. Furthermore, radial velocity data from 17 Doppler radars in eastern China are assimilated. As compared to the impact of conventional observation, the assimilation of radar radial velocity based on the U–V control variable scheme significantly improves the mesoscale dynamic field in the initial condition. The enhanced low-level jet stream, water vapor convergence and low-level wind shear result in better squall line forecasting. However, the ψ–χ control variable scheme generates a discontinuous wind field and unrealistic convergence/divergence in the analyzed field, which lead to a degraded precipitation forecast.     

集合卡尔曼滤波同化多普勒雷达资料的观测系统模拟试验

本文将集合卡尔曼滤波同化技术应用到对流尺度系统中,实施了基于WRF模式的同化单部多普勒雷达径向风和反射率因子的观测系统模拟试验,验证了其在对流尺度中应用的可行性和有效性,并对同化系统的特性进行了探讨。试验表明:WRF-EnKF雷达资料同化系统能较准确分析模式风暴的流场、热力场、微物理量场的细致特征;几乎所有变量的预报和分析误差经过同化循环后都能显著下降,同化分析基本上能使预报场在各层上都有所改进,对预报场误差较大层次的更正更为显著;约8个同化循环后,EnKF能在雷达反射率、径向风观测与背景场间建立较可靠的相关关系,使模式各变量场能被准确分析更新,背景场误差协方差在水平方向和垂直方向都有着复杂的结构,是高度非均匀、各项异性和流依赖的;集合平均分析场做的确定性预报在短时间内能较好保持真值场风暴的细节结构,但预报误差增长较快。     

Estimating the Soil Moisture Profile by Assimilating Near-Surface Observations with the Ensemble Kalman Filter （EnKF）

The paper investigates the ability to retrieve the true soil moisture profile by assimilating near-surface soil moisture into a soil moisture model with an ensemble Kalman filter （EnKF） assimilation scheme, including the effect of ensemble size, update interval and nonlinearities in the profile retrieval, the required time for full retrieval of the soil moisture profiles, and the possible influence of the depth of the soil moisture observation. These questions are addressed by a desktop study using synthetic data. The ＂true＂ soil moisture profiles are generated from the soil moisture model under the boundary condition of 0.5 cm d^-1 evaporation. To test the assimilation schemes, the model is initialized with a poor initial guess of the soil moisture profile, and different ensemble sizes are tested showing that an ensemble of 40 members is enough to represent the covariance of the model forecasts. Also compared are the results with those from the direct insertion assimilation scheme, showing that the EnKF is superior to the direct insertion assimilation scheme, for hourly observations, with retrieval of the soil moisture profile being achieved in 16 h as compared to 12 days or more. For daily observations, the true soil moisture profile is achieved in about 15 days with the EnKF, but it is impossible to approximate the true moisture within 18 days by using direct insertion. It is also found that observation depth does not have a significant effect on profile retrieval time for the EnKF. The nonlinearities have some negative influence on the optimal estimates of soil moisture profile but not very seriously.     

多普勒雷达资料同化对江苏一次飑线过程的数值模拟

应用新一代中尺度预报模式WRF模式及其3DVar同化系统, 针对江苏地区2009年6月14日飑线过程进行了多普勒雷达资料的同化试验研究, 在对雷达资料进行严格质量控制的基础上, 设计一系列尺度化因子优化调整及同化频率的敏感性试验。试验结果表明:同化后初始场得到不同程度改善, 适当的尺度化因子设定, 能够有效改进对模式初始场中700 hPa风场和850 hPa温度场以及组合反射率因子等要素的分析, 进而改善短时降水预报和风暴的垂直结构配置;并且同化频率越高, 对初始场的组合反射率因子分布与观测更为接近, 短时降水预报越准确。     

Lightning and Doppler radar observations of a squall line system

A typical squall line with damaging wind and hailstones occurred on 28 April 2006 in Shandong Province, middle eastern China, and caused great economic loss. The characteristics of cloud-to-ground lightning (CG) in the squall line were studied in detail by combining the data from the ground-based CG location network, two Doppler radars and the Lightning Imaging Sensor on the TRMM satellite. Results show that positive CG flashes accounted for 54.7% of the total CG flashes. During the initial developing stage, the CG flash rate was lower than 0.5fl min^{− 1} and most of the CG flashes were positive. It increased significantly, up to 4.5fl min^{− 1}, along with the rapid development of the squall line, and the percentage of positive CG was more than 75% during this period. The CG flash rate began to decrease but the percentage of negative CG flash increased gradually and exceeded that of positive CG during the mature and dissipating stages. Positive CG flashes tended to occur on the right flank and negative ones on the left flank. Strong wind at the surface occurred in or near the regions with dense positive CG flashes. Almost all positive CG flashes occurred near the strong radar echo regions, in the front parts of the squall line. However, the negative CG flashes almost exclusively occurred in the regions with weak and uniform radar echoes. The total flash rate in the storm was very high, up to 136fl min^{− 1}, and its ratio of intracloud flashes (IC) to CG flashes was 35:1. Dense positive CG flashes corresponded to updraft regions, they did not occur in the core of the updraft, but just behind and close to the main updraft instead. The rear inflow jet, between 3 and 6 km, played an important role in the formation of the bow echo and very strong wind at surface. The CG distribution features in the squall line were obviously different from that of an ordinary MCS. The charge structure could be roughly described as an inverted charge structure.