雷达资料同化对一次飑线过程的模拟影响

雷达资料同化能够改善强对流天气的预报,但是不同的模式方案配置会得到不同的结果。本文针对中国南部2018年3月4日一次飑线过程,以全球预报模式GFS分析场和预报场为背景场,采用中尺度区域气象预报模式ARPS 3DVAR系统同化多普勒雷达径向速度,用云分析处理反射率数据,考虑同化间隔、频次、云分析中不同参数调整,采用1 h同化窗口,设计不同同化方案,最后用WRF模式进行预报,研究雷达资料同化对飑线系统触发及发展机制的影响。结果表明,同化间隔过短时,由于模式热动力变量没有平衡产生虚假回波,同化间隔过长时,系统触发和发展的特征普遍偏弱;采用12 min间隔同化得到了最好的初始场,并且同化频次越高得到的降水预报结果越好。此外,ARPS云分析能大大改善初始场,减少模式自调整时间,其中湿度调整、温度调整、雨水调整及水汽调整对系统动力过程和水凝物初始场分布都有较大的影响,而垂直速度相关参数调整影响较小。     

EnSRF雷达资料同化在一次飑线过程中的应用研究

本文利用包含复杂冰相微物理过程的WRF（Weather Research and Forecasting）模式,针对2007年4月23日发生在我国华南地区的一次典型飑线天气过程,分别进行了确定性预报和集合预报试验,发现确定性预报能大致捕捉到飑线系统的发生发展过程,但对飑线后部的层云区模拟效果较差。集合预报能够有效地减少模式的不确定性,大部分集合成员对飑线的模拟效果优于确定性预报。进一步将集合预报得到的40个成员作为背景场,采用EnSRF（Ensemble Square Root Filter）同化多普勒天气雷达资料,并将分析得到的集合作为初始场进行集合预报,通过与未同化雷达资料的集合对比,考察了EnSRF同化多部雷达资料对飑线系统的影响。结果表明:EnSRF雷达资料同化增加了模式初始场的中小尺度信息,大部分集合成员的分析场能够较准确地再现飑线的热力场、动力场和微物理场的细致特征,并且模拟出飑线后部的层云结构。通过对EnSRF分析的集合进行模拟发现,大部分集合成员较未同化雷达资料时模拟效果有明显改善。同化后的集合预报ETS（Equitable Threat Score）评分最高,其次是未同化的集合预报,确定性预报的最低。     

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

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

Improving the Analyses and Forecasts of a Tropical Squall Line Using Upper Tropospheric Infrared Satellite Observations

The advent of modern geostationary satellite infrared radiance observations has noticeably improved numerical weather forecasts and analyses.However,compared to midlatitude weather systems and tropical cyclones,research into using infrared radiance observations for numerically predicting and analyzing tropical mesoscale convective systems remain mostly fallow.Since tropical mesoscale convective systems play a crucial role in regional and global weather,this deficit should be addressed.This study is the first of its kind to examine the potential impacts of assimilating all-sky upper tropospheric infrared radiance observations on the prediction of a tropical squall line.Even though these all-sky infrared radiance observations are not directly affected by lower-tropospheric winds,the high-frequency assimilation of these all-sky infrared radiance observations improved the analyses of the tropical squall line’s outflow position.Aside from that,the assimilation of all-sky infrared radiance observations improved the analyses and prediction of the squall line’s cloud field.Finally,reducing the frequency of assimilating these all-sky infrared radiance observations weakened these improvements to the analyzed outflow position,as well as the analyses and predictions of cloud fields.     

循环同化雷达资料对一次飑线系统临近预报的改进作用

利用WRF（Weather Research Forecast）模式及其3D-Var（Three-Dimensional Variational assimilation）变分系统,针对2017年7月7日一次飑线进行了雷达资料的循环同化敏感性试验。结果表明：以循环同化雷达资料至飑线成熟期时刻的试验预报效果最好,主要原因在于预报的低层西北冷空气较强,从而导致较强的低层切变,再配合强的热力不稳定条件从而使飑线的发展得以维持。通过七组试验对比表明,对于单次雷达资料,同化的时机更为重要。同化飑线成熟阶段的雷达反射率,对低层热力层结有改善作用,为飑线发展提供了不稳定能量;对于循环同化,通过观测的影响和模式自身的热动力调整,对流场也有较好的修正作用,为对流系统的持续发展提供了充分的动力条件。     

Analyses and forecasts of a tornadic supercell outbreak using a 3DVAR system ensemble

As part of NOAA’s "Warn-On-Forecast" initiative, a convective-scale data assimilation and prediction system was developed using the WRF-ARW model and ARPS 3DVAR data assimilation technique. The system was then evaluated using retrospective short-range ensemble analyses and probabilistic forecasts of the tornadic supercell outbreak event that occurred on 24 May 2011 in Oklahoma, USA. A 36-member multi-physics ensemble system provided the initial and boundary conditions for a 3-km convective-scale ensemble system. Radial velocity and reflectivity observations from four WSR-88 Ds were assimilated into the ensemble using the ARPS 3DVAR technique. Five data assimilation and forecast experiments were conducted to evaluate the sensitivity of the system to data assimilation frequencies, in-cloud temperature adjustment schemes, and fixed- and mixed-microphysics ensembles. The results indicated that the experiment with 5-min assimilation frequency quickly built up the storm and produced a more accurate analysis compared with the 10-min assimilation frequency experiment. The predicted vertical vorticity from the moist-adiabatic in-cloud temperature adjustment scheme was larger in magnitude than that from the latent heat scheme. Cycled data assimilation yielded good forecasts, where the ensemble probability of high vertical vorticity matched reasonably well with the observed tornado damage path. Overall, the results of the study suggest that the 3DVAR analysis and forecast system can provide reasonable forecasts of tornadic supercell storms.     

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

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

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.     

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

利用中尺度模式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 (Weather Research and Forecasting) 模式及其Hybrid ETKF-3DVAR同化系统,考察不同生成方案的样本对同化地面观测的影响。集合样本创建方式包括3类:扰动初始背景场的方案 (RCV)、使用不同的物理参数化方案 (PPMP) 以及前两者集成方案 (BLE)。基于增量场分析,同化地面观测主要调整850 hPa以下水平风和水汽混合比的空间结构,其中RCV方案侧重于改变水平风的空间分布,PPMP方案侧重于改变水汽混合比的空间结构,BLE方案兼具二者特征。同化地面观测可以间接改善6 h降水预报,其中PPMP试验的降水预报最好,尤其是对降水位置和强度的预报。对比雷达回波观测,RCV试验和BLE试验对弓状回波模拟得较好,BLE试验的模拟较多体现RCV特征。PPMP试验和RCV试验还可改变冷池的位置和强度,同时影响飑线出现和消亡时间,相对而言,PPMP试验影响更大。     

Assimilation of Doppler Radar Observations with an Ensemble Square Root Filter：A Squall Line Case Study

The effectiveness of using an Ensemble Square Root Filter（EnSRF） to assimilate real Doppler radar observations on convective scale is investigated by applying the technique to a case of squall line on 12July 2005 in midwest Shandong Province using the Weather Research and Forecasting（WRF） model.The experimental results show that：（1） The EnSRF system has the potential to initiate a squall line accurately by assimilation of real Doppler radar data.The convective-scale information has been added into the WRF model through radar data assimilation and thus the analyzed fields are improved noticeably.The model spin-up time has been shortened,and the precipitation forecast is improved accordingly.（2） Compared with the control run,the deterministic forecast initiated with the ensemble mean analysis of EnSRF produces more accurate prediction of microphysical fields.The predicted wind and thermal fields are reasonable and in accordance with the characteristics of convective storms.（3） The propagation direction of the squall line from the ensemble mean analysis is consistent with that of the observation,but the propagation speed is larger than the observed.The effective forecast period for this squall line is about 5-6 h,probably because of the nonlinear development of the convective storm.     

Assimilating Best Track Minimum Sea Level Pressure Data Together with Doppler Radar Data Using an Ensemble Kalman Filter for Hurricane Ike (2008) at a Cloud-Resolving Resolution

Extending an earlier study, the best track minimum sea level pressure (MSLP) data are assimilated for landfalling Hurricane Ike (2008) using an ensemble Kalman filter (EnKF), in addition to data from two coastal ground-based Doppler radars, at a 4-km grid spacing. Treated as a sea level pressure observation, the MSLP assimilation by the EnKF enhances the hurricane warm core structure and results in a stronger and deeper analyzed vortex than that in the GFS (Global Forecast System) analysis; it also improves the subsequent 18-h hurricane intensity and track forecasts. With a 2-h total assimilation window length, the assimilation of MSLP data interpolated to 10-min intervals results in more balanced analyses with smaller subsequent forecast error growth and better intensity and track forecasts than when the data are assimilated every 60 minutes. Radar data are always assimilated at 10-min intervals. For both intensity and track forecasts, assimilating MSLP only outperforms assimilating radar reflectivity (Z) only. For intensity forecast, assimilating MSLP at 10-min intervals outperforms radar radial wind (Vr) data (assimilated at 10-min intervals), but assimilating MSLP at 60-min intervals fails to beat Vr data. For track forecast, MSLP assimilation has a slightly (noticeably) larger positive impact than Vr(Z) data. When Vr or Z is combined with MSLP, both intensity and track forecasts are improved more than the assimilation of individual observation type. When the total assimilation window length is reduced to 1h or less, the assimilation of MSLP alone even at 10-min intervals produces poorer 18-h intensity forecasts than assimilating Vr only, indicating that many assimilation cycles are needed to establish balanced analyses when MSLP data alone are assimilated; this is due to the very limited pieces of information that MSLP data provide.     

PRELIMINARY STUDY ON THE INFLUENCE OF FY-4 LIGHTNING DATA ASSIMILATION ON PRECIPITATION PREDICTIONS

Data from the lightning mapping imager on board the Fengyun-4 meteorological satellite (FY-4) were used to study the assimilation of lightning data and its influence on precipitation predictions. We first conducted a quality control check on the events observed by the first Fengyun-4 satellite (FY-4A) lightning mapping imager, after which the noise points were removed from the lightning distribution. The subsequent distribution was more consistent with the spatial distribution and range of ground-based observations and precipitation. We selected the radar reflectivity, which was closely related to the lightning frequency, as the parameter to assimilate the lightning data and utilized a large sample of lightning frequency and radar reflectivity data from the eastern United States provided by Vaisala. Based on statistical analysis, we found the empirical relationship between the lightning frequency and radar reflectivity and established a look-up table between them. We converted the lightning event data into radar reflectivity data and found that the converted reflectivity and composite reflectivity of ground-based radar observations showed high consistency. We further assimilated the lightning data into the model, adjusted the model cloud analysis process and adjusted the model hydrometeor field by using the lightning data. A rainstorm weather process that occurred on August 8, 2017, in south China was used for the numerical forecast experiment, and three experiments were designed for comparison and analysis: a control experiment, an experiment without the assimilation of FY-4 lightning data (NoLig), and an experiment with the assimilation of FY-4 lightning data (Lig). The results show that after assimilating the FY-4A lightning data, the accuracies of the intensity, central location and range of the precipitation predicted by the Lig experiment were obviously superior to those predicted by the control and NoLig experiments, and the effect was especially obvious in the short-term (1–2 hour) forecast. The studies in this paper highlight the application value and potential of FY-4 lightning data in precipitation predictions.     

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

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

雷达反射率资料的三维变分同化研究

应用天气研究和预报模式(WRF)三维变分系统中一种新的雷达反射率资料间接同化方法来进行反射率资料的三维变分同化研究,评估雷达反射率资料对夏季短时定量降水预报的作用.该方法不直接同化雷达反射率资料,而是同化由反射率资料反演出的雨水和估计的水汽.以2009年夏季北京地区发生的4次强降水过程为例,考察了北京市气象局业务运行的快速更新循环同化预报系统对京津冀地区雷达网的雷达反射率资料的同化性能以及雷达反射率资料和径向风资料同时同化的效果.数值试验结果表明:(1)同化反演雨水或水汽都能改善降水预报,但同化反演水汽对降水预报效果的改善起了更重要的作用;(2)同化反射率资料能极大地提高短时降水预报的效果,其稳定的正面效果可以延伸到6h的预报时效,而同化径向风资料不能得到稳定的正效果;(3)同化雷达资料时,应用快速更新循环同化预报系统是提高短时定量降水预报的一个有效途径.     

逐时云迹风资料同化对暴雨预报的模拟试验

文中基于不同的云迹风同化方案,用GRAPES模式对2005年7月11-12日长江中下游一次暴雨强降水过程进行了云迹风资料同化试验及数值模拟,通过对比分析不同方案所得的分析场及预报场的差异,研究逐时云迹风资料三维变分同化对分析场及暴雨预报的影响.首先,根据连续性原理及双通道各层次云迹风资料的误差分析,分3个步骤对7月11日00:00-12:00 UTC共12个时次的双通道云迹风资料进行了初步的质量控制;然后,把经过质量控制的云迹风资料放入基于GRAPES 3D-VAR三维变分同化方案开发的逐小时循环同化系统中进行同化,将得出的分析场与单一时次未经质量控制的云迹风资料同化得出的分析场进行对比,探讨了逐时云迹风资料同化对数值预报分析场的影响;最后,把同化后的分析场作为初始场,用 GRAPES模式对 2005年7月11-12日长江流域暴雨过程做24 h降水预报试验,分析两个同化方案所模拟得到的预报场的差异.结果表明,经资料的筛选、同经纬度单点通道的选择及资料的稀疏化3个步骤控制后,各层次云迹风资料的误差有明显减小;加入经质量控制的逐时云迹风资料其三维变分同化可以提高分析场中风压场及水汽场的质量;而且在暴雨预报试验中可以相对更准确地预报暴雨落区及雨强.     

多普勒雷达资料同化在“7.21”北京特大暴雨个例中的应用

基于WRF(Weather Research Forecast)模式和GSI(Gridpoint Statistical Interpolation)同化系统,研究了同化4部多普勒雷达探测资料对"7.21"北京特大暴雨过程中降水预报的改善作用。GSI系统直接同化径向风,而采用云分析的方式间接同化反射率。2012年7月20日21时—21日00时(世界时)雷达探测资料同化试验采用30 min循环同化径向风和反射率资料。结果表明,循环同化雷达探测资料改善了短时(0—6 h)和短期(0—24 h)降水预报,ETS评分提高了约0.2。同化反射率资料增加了初始场的水凝物,改善了温度场分布,直接影响了降水的形成,同时还使650—250 hPa位势高度的均方根误差平均降低了8 gpm。直接同化径向风资料对中尺度风场产生了一定影响。ETS评分结果表明:同化反射率资料的效果要优于同化径向风。     

雷达资料同化对2016年6月23日阜宁龙卷模拟的改进

基于GSI（Community Gridpoint Statistical Interpolation）同化系统和WRF（Weather Research and Forecasting）模式,探讨了多部多普勒雷达的反射率因子和径向速度资料同化对2016年6月23日江苏阜宁龙卷模拟的改进效果和影响过程。结果表明:（1）仅同化雷达反射率因子和仅同化径向速度均能在一定程度上改进模式对阜宁龙卷及其环境场的模拟,且雷达径向速度同化的改进作用更大;同时同化两种资料改进效果最佳。（2）雷达反射率因子同化是利用复杂云分析技术,直接修正了水凝物含量,增加了潜热释放,对低层大气热力场进行了正温度扰动调整,从而主要改进了初始场的水汽条件和热力条件;而雷达径向速度同化通过三维变分技术直接修正了风场,进而引起水汽输送变化影响水凝物的调整和大气热力场的变化,对初始场的动力条件和热力条件修正较大;同时同化两种资料修正了初始场的动力和热力结构,保证了两者物理上的协调,综合了两者的改进作用,从而取得最佳模拟效果。（3）同时同化雷达反射率因子和径向速度后,模式在阜宁附近模拟出了明显的涡旋结构,尽管涡旋强度和龙卷结构与实况仍有一定差距,但涡旋发生发展过程、路径、地面小时极大风和降水等模拟与实况吻合度均明显高于对照试验。      

GNSS反演资料在GRAPES_Meso三维变分中的应用

为了进一步提高GRAPES_Meso的分析和预报效果,该文在GRAPES_Meso三维变分同化系统中建立了同化GNSS/RO反演的大气资料的观测算子,实现了对GNSS/RO反演的大气资料的同化应用,并通过2013年7月1个月的同化和预报试验分析了GNSS/RO反演大气资料对GRAPES_Meso模式系统分析和预报的影响。结果表明:增加了GNSS/RO反演大气资料的同化后,GRAPES_Meso位势高度场的分析误差明显减小,平均分析误差减小约8%,预报误差略有减小,平均预报误差减小约1%;湿度场的分析误差和预报误差变化不明显,常规观测资料稀少的青藏高原地区的降水预报技巧有所提高,小雨到大雨的ETS (equitable threat score) 评分提高约0.01,对全国及其他分区的降水预报技巧总体上有正效果。     

Impacts of Multigrid NLS-4DVar-based Doppler Radar Observation Assimilation on Numerical Simulations of Landfalling Typhoon Haikui (2012)

We applied the multigrid nonlinear least-squares four-dimensional variational assimilation(MG-NLS4DVar) method in data assimilation and prediction experiments for Typhoon Haikui(2012) using the Weather Research and Forecasting(WRF) model. Observation data included radial velocity(V_r) and reflectivity(Z) data from a single Doppler radar, quality controlled prior to assimilation. Typhoon prediction results were evaluated and compared between the NLS-4DVar and MG-NLS4DVar methods. Compared with a forecast that began with NCEP analysis data, our radar data assimilation results were clearly improved in terms of structure, intensity, track, and precipitation prediction for Typhoon Haikui(2012). The results showed that the assimilation accuracy of the NLS-4DVar method was similar to that of the MG-NLS4DVar method,but that the latter was more efficient. The assimilation of V_r alone and Z alone each improved predictions of typhoon intensity, track, and precipitation; however, the impacts of V_r data were significantly greater that those of Z data.Assimilation window-length sensitivity experiments showed that a 6-h assimilation window with 30-min assimilation intervals produced slightly better results than either a 3-h assimilation window with 15-min assimilation intervals or a 1-h assimilation window with 6-min assimilation intervals.