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

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

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.     

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.     

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.     

Analysis and very short range forecast of cyclone “AILA” with radar data assimilation with rapid intermittent cycle using ARPS 3DVAR and cloud analysis techniques

In this study, both reflectivity and radial velocity are assimilated into the Weather Research and Forecasting (WRF) model using ARPS 3DVAR technique and cloud analysis procedure for analysis and very short range forecast of cyclone ÁILA. Doppler weather radar (DWR) data from Kolkata radar are assimilated for numerical simulation of landfalling tropical cyclone. Results show that the structure of cyclone AILA has significantly improved when radar data is assimilated. Radar reflectivity data assimilation has strong influence on hydrometeor structures of the initial vortex and precipitation pattern and relatively less influence is observed on the wind fields. Divergence/convergence conditions over cyclone inner-core area in the low-to-middle troposphere (600–900 hPa) are significantly improved when wind data are assimilated. However, less impact is observed on the moisture field. Analysed minimum sea level pressure (SLP) is improved significantly when both reflectivity and wind data assimilated simultaneously (RAD-ZVr experiment), using ARPS 3DVAR technique. In this experiment, the centre of cyclone is relocated very close to the observed position and the system maintains its intensity for longer duration. As compared to other experiments track errors are much reduced and predicted track is very much closer to the best track in RAD-ZVr experiment. Rainfall pattern and amount of rainfall are better captured in this experiment. The study also reveals that cyclone structure, intensification, direction of movement, speed and location of cyclone are significantly improved and different stages of system are best captured when both radar reflectivity and wind data are assimilated using ARPS 3DVAR technique and cloud analysis procedure. Thus optimal impact of radar data is realized in RAD-ZVr experiment. The impact of DWR data reduces after 12 h forecast and it is due to the dominance of the flow from large-scale global forecast system model. Successful coupling of data assimilation package ARPS 3DVAR with WRF model for Indian DWR data is also demonstrated.     

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

利用中尺度模式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降水对比发现,利用雷达径向风调整初始场后,对降水模拟有一定的改进,但效果不甚明显,而雷达反射率资料对定量降水预报改进效果明显,同时使用雷达径向风和反射率资料改进初始场后对降水的模拟效果最明显.     

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

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

CINRAD/CD雷达反射率因子同化对中尺度数值模式云微物理量场调整的分析

以ARPS模式及其资料三维同化系统ARPS3DVAR和复杂云分析模块为研究平台,使用贵阳多普勒雷达的体扫观测资料,分析了雷达反射率因子资料对云微物理量场的同化问题及雷达资料质量控制对同化的影响.结果表明,数值试验中两个不同的云分析方案都能同化出合理的云水场和水凝物场,同时因为水凝物场的调整,位温场和垂直速度场也得到了很好的响应;与没有使用云分析方案的试验比较,使用了云分析方案同化雷达资料的试验,在初始时刻能够调整出一个合理的云微物理量场,减少模式的热启动时间;而没有使用云分析方案的试验在2～3 h后才能调整出一个位置有较大偏差的水凝物场,雷达反射率资料的同化在强对流风暴的模拟中起到了一个非常关键的作用.     

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

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

微物理方案及其参数对EnSRF资料同化的影响研究

利用自主构建的基于风暴尺度的WRF-En SRF系统同化模拟多普勒雷达资料,讨论了微物理方案及其参数的不确定性对同化效果的影响。试验采用组合微物理方案以及扰动微物理方案中的参数的方法,结果表明,模式误差非常小甚至可以忽略时,使用单个微物理方案并扰动参数能够使真实风暴的主要特征在分析场中较未扰动参数得到更好地反映;存在模式误差时,使用单个微物理方案并扰动参数后,分析场中的各要素的分布较未扰动参数更加接近真实风暴,同化效果得到改进,且改进效果比模式误差非常小时更为明显;存在模式误差时,组合微物理方案并扰动参数后,分析场中对流云团的形态较未组合方案或未扰动参数更接近真实风暴,主要要素场的配置最能反映真实风暴的特征,同化效果最为理想。结果也表明,扰动参数时、参数扰动范围较小时,同化效果较优。     

Mesoscale Structure of Rainstorm Retrieved from Dual-Doppler Radar Observations Using the 4DVAR Assimilation Technique

The four-dimensional variational (4DVAR) data assimilation method was applied to dual-Doppler radar data about two Meiyu rainstorms observed during CHeRES (China Heavy Rain Experiment and Study). The purpose of this study is to examine the performance of the 4DVAR technique in retrieving rainstorm mesoscale structure and to reveal the feature of rainstorm mesoscale structure. Results demonstrated that the 4DVAR assimilation method was able to retrieve the detailed structure of wind, thermodynamics, and microphysics fields from dual-Doppler radar observations. The retrieved wind fields agreed with the dual-Doppler synthesized winds and were accurate. The distributions of the retrieved perturbation pressure, perturbation temperature, and microphysics fields were also reasonable through the examination of their physical consistency. Both of the two heavy rainfalls were caused by merging cloud processes. The wind shear and convergence lines at middle and lower levels were their primary dynamical characteristics. The convective system was often related to low-level convergence and upper-level divergence coupled with up-drafts. During its mature stage, the convective system was characterized by low pressure at lower level and high pressure at upper level, associated with warmer at middle level and colder at lower and upper levels than the environment. However, a region of cooling and high pressure occurred in the lower and middle levels compared to warming and low pressure in the upper level during its dissipating stage. The water vapor, cloud water, and rainwater corresponded to the convergence, the updraft and the intensive reflectivity, respectively.     

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

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

单多普勒雷达二维风场反演——Extended VPP方法

文中提出的EVPP(Extended Velocity Plan Processing)方法是根据风场连续及局部分片光滑的性质,假定在一很小的分析窗内风矢量保持线性变化关系,使得在分析窗内反演得到的风场在矢径上的投影与实际观测的多普勒雷达径向速度场之间的误差达到最小,并利用零速度线这一特点对问题的求解做出约束,使得所得解尽可能达到全局最优,即:(1)利用变分的思想求出水平风速的大小及方向;(2)利用风场局部分片光滑,对风的方向做出矫正。虽然由于风场受到线性假设这一限制,使得其使用的范围得到了局限,但是在本文中利用零速度线这一特征使得反演的结果得到了进一步的提高。文中利用模拟的理想个例和实际观测到的个例分别进行了研究分析。为了证明其理论的可行性,首先我们利用三维冰雹云模式模拟了辽宁2005年7月8日的一次冰雹云过程,考虑到个例的代表性,我们选择了辐合、辐散、旋转以及均匀风场等4种类型的流场结构进行了模拟分析,同时为了验证该种方法对于资料质量的依赖性以及噪声的敏感性,我们在理想质量中加入了随机噪声。从对模拟风场的反演,证明了该方法理论上的可行性,而且噪声对其结构的影响不大。然后针对台风、飑线、强对流等实际天气过程作了进一步分析。最后对反演结果进行了误差分析,结果表明这一方法能较好地反演出中尺度系统的水平流场结构。     

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

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

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

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

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.     

四维变分方法反演低层风场能力研究

从新一代天气雷达径向速度资料中反演出可靠的三维风场对提高新一代天气雷达的应用水平有重要的作用，将雷达直接观测的径向速度转换成台站预报员更为熟悉的风场结构，对识别中小尺度信息有很大帮助。为此该文对4DVAR同化技术在风场业务反演中应用的可能性进行了探讨，利用广州、济南新一代多普勒天气雷达观测资料，从是否加入云模式湿过程以及迭代次数、计算时间、背景场、初始场、分辨率和反演区域等方面对干模式的4DVAR系统进行了讨论，并从风场结构、均方根差别等方面对反演结果进行分析。多种试验表明，干模式的4DVAR系统与湿的云模式反演结果差异不大。模式的初始场和背景场对反演结果具有较高的敏感性，利用前一时次的反演结果作为背景场迭代15～20次的干模式结果可以很好地在业务上试运行，对台站预报员提高中小尺度天气预报的准确率有着很重要的作用。     

多普勒激光雷达风场反演方法研究

采用三维变分同化反演(3DVAR)、 四维变分同化反演(4DVAR) 对多普勒激光雷达资料反演风场的方法进行了研究, 利用车载多普勒激光雷达在2008年残奥会测试赛期间外场试验取得的数据, 反演了海面10 m高度处的风场, 并将风场反演结果与浮标资料进行了对比分析, 结果表明: 3DVAR、4DVAR风场反演方法均能实现近海面风的精细化风场反演, 并能反映出风向的变化, 反演风场与浮标数据基本一致, 在风速较大的天气情况, 3DVAR与4DVAR反演风场的一致性要好于风速较小的天气情况; 4DVAR反演方法中以浮标资料作为背景场, 使得其与浮标的符合程度要好于3DVAR方法反演风场; 反演风场的风向与浮标风向具有很好的相关关系, 反演风场的风速与浮标风速具有一定的相关关系, 反演风场的风向、风速与浮标的风向、风速之间平均均方根误差和平均绝对误差表明, 这两序列之间具有一定差别, 在风速较小的天气情况下使用时需要注意。     

华南一次强飑线过程的三维变分风场反演效果分析

利用三维变分方法对2014年3月30—31日华南一次强飑线过程进行风场反演,经与风廓线雷达探测结果、双多普勒天气雷达反演结果、原始径向速度数据等对比分析,得到如下结论：三维变分方法反演的中低层水平风场与风廓线雷达探测到的结果较为一致,且能很好地表现飑线过境时的风向切变;通过与双多普勒雷达风场反演结果对比发现,两种方法得到的风场空间分布十分相似,均能很好地表现2 km高度上系统内部强带状回波前缘的辐合线以及5 km高度上较弱的辐散;三维变分方法反演的水平风场与径向速度场有较好的一致性,2 km高度强回波带前缘阵风锋处的辐合线位置以及5 km和8 km高度上辐散区的位置均与径向速度场十分吻合;三维变分方法反演的垂直速度能较好地反映该飑线过程中气流的上升和下沉运动,平行于飑线方向的气流变化较小,而系统气流变化主要沿垂直于飑线的方向。三维变分方法反演的飑线系统的三维风场结构合理,反演结果可靠。     

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.