Assimilating polarimetric radar data with an ensemble Kalman filter: OSSEs with a tornadic supercell storm simulated with a two‐moment microphysics scheme

2020 ◽  
Vol 146 (729) ◽  
pp. 1880-1900
Author(s):  
Kefeng Zhu ◽  
Ming Xue ◽  
Kun Ouyang ◽  
Youngsun Jung
2014 ◽  
Vol 142 (1) ◽  
pp. 141-162 ◽  
Author(s):  
Bryan J. Putnam ◽  
Ming Xue ◽  
Youngsun Jung ◽  
Nathan Snook ◽  
Guifu Zhang

Abstract Doppler radar data are assimilated with an ensemble Kalman Filter (EnKF) in combination with a double-moment (DM) microphysics scheme in order to improve the analysis and forecast of microphysical states and precipitation structures within a mesoscale convective system (MCS) that passed over western Oklahoma on 8–9 May 2007. Reflectivity and radial velocity data from five operational Weather Surveillance Radar-1988 Doppler (WSR-88D) S-band radars as well as four experimental Collaborative and Adaptive Sensing of the Atmosphere (CASA) X-band radars are assimilated over a 1-h period using either single-moment (SM) or DM microphysics schemes within the forecast ensemble. Three-hour deterministic forecasts are initialized from the final ensemble mean analyses using a SM or DM scheme, respectively. Polarimetric radar variables are simulated from the analyses and compared with polarimetric WSR-88D observations for verification. EnKF assimilation of radar data using a multimoment microphysics scheme for an MCS case has not previously been documented in the literature. The use of DM microphysics during data assimilation improves simulated polarimetric variables through differentiation of particle size distributions (PSDs) within the stratiform and convective regions. The DM forecast initiated from the DM analysis shows significant qualitative improvement over the assimilation and forecast using SM microphysics in terms of the location and structure of the MCS precipitation. Quantitative precipitation forecasting skills are also improved in the DM forecast. Better handling of the PSDs by the DM scheme is believed to be responsible for the improved prediction of the surface cold pool, a stronger leading convective line, and improved areal extent of stratiform precipitation.


2012 ◽  
Vol 140 (5) ◽  
pp. 1457-1475 ◽  
Author(s):  
Youngsun Jung ◽  
Ming Xue ◽  
Mingjing Tong

Abstract The performance of ensemble Kalman filter (EnKF) analysis is investigated for the tornadic supercell on 29–30 May 2004 in Oklahoma using a dual-moment (DM) bulk microphysics scheme in the Advanced Regional Prediction System (ARPS) model. The comparison of results using single-moment (SM) and DM microphysics schemes evaluates the benefits of using one over the other during storm analysis. Observations from a single operational Weather Surveillance Radar-1988 Doppler (WSR-88D) are assimilated and a polarimetric WSR-88D in Norman, Oklahoma (KOUN), is used to assess the quality of the analysis. Analyzed reflectivity and radial velocity in the SM and DM experiments compare favorably with independent radar observations in general. However, simulated polarimetric signatures obtained from analyses using a DM scheme agree significantly better with polarimetric signatures observed by KOUN in terms of the general structure, location, and intensity of the signatures than those generated from analyses using an SM scheme. These results demonstrate for the first time for a real supercell storm that EnKF data assimilation using a numerical model with an adequate microphysics scheme (i.e., a scheme that predicts at least two moments of the hydrometeor size distributions) is capable of producing polarimetric radar signatures similar to those seen in observations without directly assimilating polarimetric data. In such cases, the polarimetric data also serve as completely independent observations for the verification purposes.


2019 ◽  
Vol 147 (7) ◽  
pp. 2511-2533 ◽  
Author(s):  
Bryan Putnam ◽  
Ming Xue ◽  
Youngsun Jung ◽  
Nathan Snook ◽  
Guifu Zhang

Abstract Real polarimetric radar observations are directly assimilated for the first time using the ensemble Kalman filter (EnKF) for a supercell case from 20 May 2013 in Oklahoma. A double-moment microphysics scheme and advanced polarimetric radar observation operators are used together to estimate the model states. Lookup tables for the observation operators are developed based on T-matrix scattering amplitudes for all hydrometeor categories, which improve upon previous curved-fitted approximations of T-matrix scattering amplitudes or the Rayleigh approximation. Two experiments are conducted: one assimilates reflectivity (Z) and radial velocity (Vr) (EXPZ), and one assimilates in addition differential reflectivity (ZDR) below the observed melting level at ~2-km height (EXPZZDR). In the EnKF analyses, EXPZZDR exhibits a ZDR arc that better matches observations than EXPZ. EXPZZDR also has higher ZDR above 2 km, consistent with the observed ZDR column. Additionally, EXPZZDR has an improved estimate of the model microphysical states. Specifically, the rain mean mass diameter (Dnr) in EXPZZDR is higher in the ZDR arc region and the total rain number concentration (Ntr) is lower downshear in the forward flank than EXPZ when compared to values retrieved from the polarimetric observations. Finally, a negative gradient of hail mean mass diameter (Dnh) is found in the right-forward flank of the EXPZZDR analysis, which supports previous findings indicating that size sorting of hail, as opposed to rain, has a more significant impact on low-level polarimetric signatures. This paper represents a proof-of-concept study demonstrating the value of assimilating polarimetric radar data in improving the analysis of features and states related to microphysics in supercell storms.


2013 ◽  
Vol 2013 ◽  
pp. 1-13 ◽  
Author(s):  
Daniel T. Dawson II ◽  
Louis J. Wicker ◽  
Edward R. Mansell ◽  
Youngsun Jung ◽  
Ming Xue

The impact of increasing the number of predicted moments in a multimoment bulk microphysics scheme is investigated using ensemble Kalman filter analyses and forecasts of the May 8, 2003 Oklahoma City tornadic supercell storm and the analyses are validated using dual-polarization radar observations. The triple-moment version of the microphysics scheme exhibits the best performance, relative to the single- and double-moment versions, in reproducing the low-ZDRhail core and high-ZDRarc, as well as an improved probabilistic track forecast of the mesocyclone. A comparison of the impact of the improved microphysical scheme on probabilistic forecasts of the mesocyclone track with the observed tornado track is also discussed.


2011 ◽  
Vol 139 (11) ◽  
pp. 3446-3468 ◽  
Author(s):  
Nathan Snook ◽  
Ming Xue ◽  
Youngsun Jung

Abstract One of the goals of the National Science Foundation Engineering Research Center (ERC) for Collaborative Adaptive Sensing of the Atmosphere (CASA) is to improve storm-scale numerical weather prediction (NWP) by collecting data with a dense X-band radar network that provides high-resolution low-level coverage, and by assimilating such data into NWP models. During the first spring storm season after the deployment of four radars in the CASA Integrated Project-1 (IP-1) network in southwest Oklahoma, a tornadic mesoscale convective system (MCS) was captured by CASA and surrounding Weather Surveillance Radars-1988 Doppler (WSR-88Ds) on 8–9 May 2007. The MCS moved across northwest Texas and western and central Oklahoma; two tornadoes rated as category 1 on the enhanced Fujita scale (EF-1) and one tornado of EF-0 intensity were reported during the event, just to the north of the IP-1 network. This was the first tornadic convective system observed by CASA. To quantify the impacts of CASA radar data in storm-scale NWP, a set of data assimilation experiments were performed using the Advanced Regional Prediction System (ARPS) ensemble Kalman filter (EnKF) system configured with full model physics and high-resolution terrain. Data from four CASA IP-1 radars and five WSR-88Ds were assimilated in some of the experiments. The ensemble contained 40 members, and radar data were assimilated every 5 min for 1 h. While the assimilation of WSR-88D data alone was able to produce a reasonably accurate analysis of the convective system, assimilating CASA data in addition to WSR-88D data is found to improve the representation of storm-scale circulations, particularly in the lowest few kilometers of the atmosphere, as evidenced by analyses of gust front position and comparison of simulated Vr with observations. Assimilating CASA data decreased RMS innovation of the resulting ensemble mean analyses of Z, particularly in early assimilation cycles, suggesting that the addition of CASA data allowed the EnKF system to more quickly achieve a good result. Use of multiple microphysics schemes in the forecast ensemble was found to alleviate underdispersion by increasing the ensemble spread. This work is the first assimilating real CASA data into an NWP model using EnKF.


2008 ◽  
Vol 136 (6) ◽  
pp. 2246-2260 ◽  
Author(s):  
Youngsun Jung ◽  
Ming Xue ◽  
Guifu Zhang ◽  
Jerry M. Straka

Abstract A data assimilation system based on the ensemble square-root Kalman filter (EnSRF) is extended to include the additional capability of assimilating polarimetric radar variables. It is used to assess the impact of assimilating additional polarimetric observations on convective storm analysis in the Observing System Simulation Experiment (OSSE) framework. The polarimetric variables considered include differential reflectivity ZDR, reflectivity difference Zdp, and specific differential phase KDP. To simulate the observational data more realistically, a new error model is introduced for characterizing the errors of the nonpolarimetric and polarimetric radar variables. The error model includes both correlated and uncorrelated error components for reflectivities at horizontal and vertical polarizations (ZH and ZV, respectively). It is shown that the storm analysis is improved when polarimetric variables are assimilated in addition to ZH or in addition to both ZH and radial velocity Vr. Positive impact is largest when ZDR, Zdp, and KDP are assimilated all together. Improvement is generally larger in vertical velocity, water vapor, and rainwater mixing ratios. The rainwater field benefits the most while the impacts on horizontal wind components and snow mixing ratio are smaller. Improvement is found at all model levels even though the polarimetric data, after the application of thresholds, are mostly limited to the lower levels. Among ZDR, Zdp, and KDP, ZDR is found to produce the largest positive impact on the analysis. It is suggested that ZDR provides more independent information than the other variables. The impact of polarimetric data is also expected to be larger when they are used to retrieve drop size distribution parameters. The polarimetric radar data thresholding prior to assimilation is found to be necessary to minimize the impact of noise. This study is believed to be the first to directly assimilate (simulated) polarimetric data into a numerical model.


Sign in / Sign up

Export Citation Format

Share Document