Ensemble Probabilistic Forecasts of a Tornadic Mesoscale Convective System from Ensemble Kalman Filter Analyses Using WSR-88D and CASA Radar Data

2012 ◽  
Vol 140 (7) ◽  
pp. 2126-2146 ◽  
Author(s):  
Nathan Snook ◽  
Ming Xue ◽  
Youngsun Jung

Abstract This study examines the ability of a storm-scale numerical weather prediction (NWP) model to predict precipitation and mesovortices within a tornadic mesoscale convective system that occurred over Oklahoma on 8–9 May 2007, when the model is initialized from ensemble Kalman filter (EnKF) analyses including data from four Engineering Research Center for Collaborative Adaptive Sensing of the Atmosphere (CASA) X-band and five Weather Surveillance Radar-1988 Doppler (WSR-88D) S-band radars. Ensemble forecasts are performed and probabilistic forecast products generated, focusing on prediction of radar reflectivity (a proxy of quantitative precipitation) and mesovortices (an indication of tornado potential). Assimilating data from both the CASA and WSR-88D radars for the ensemble and using a mixed-microphysics ensemble during data assimilation produces the best probabilistic mesovortex forecast. The use of multiple microphysics schemes within the ensemble aims to address at least partially the model physics uncertainty and effectively plays a role of flow-dependent inflation (in precipitation regions) during EnKF data assimilation. The ensemble predicts with high probability (approximately 0.65) the near-surface mesovortex associated with the first of three reported tornadoes. Though a bias toward stronger precipitation is noted in the ensemble forecasts, all experiments produce skillful probabilistic forecasts of radar reflectivity on a 0–3-h time scale as evaluated by multiple probabilistic verification metrics. These results suggest that both the inclusion of CASA radar data and use of a mixed-microphysics ensemble during EnKF data assimilation positively impact the skill of 2–3-h ensemble forecasts of mesovortices, despite having little impact on the quality of precipitation forecasts (analyzed in terms of predicted radar reflectivity), and are important steps toward an operational EnKF-based ensemble analysis and probabilistic forecast system to support convective-scale warn-on-forecast operations.

2017 ◽  
Vol 145 (6) ◽  
pp. 2257-2279 ◽  
Author(s):  
Bryan J. Putnam ◽  
Ming Xue ◽  
Youngsun Jung ◽  
Nathan A. Snook ◽  
Guifu Zhang

Abstract Ensemble-based probabilistic forecasts are performed for a mesoscale convective system (MCS) that occurred over Oklahoma on 8–9 May 2007, initialized from ensemble Kalman filter analyses using multinetwork radar data and different microphysics schemes. Two experiments are conducted, using either a single-moment or double-moment microphysics scheme during the 1-h-long assimilation period and in subsequent 3-h ensemble forecasts. Qualitative and quantitative verifications are performed on the ensemble forecasts, including probabilistic skill scores. The predicted dual-polarization (dual-pol) radar variables and their probabilistic forecasts are also evaluated against available dual-pol radar observations, and discussed in relation to predicted microphysical states and structures. Evaluation of predicted reflectivity (Z) fields shows that the double-moment ensemble predicts the precipitation coverage of the leading convective line and stratiform precipitation regions of the MCS with higher probabilities throughout the forecast period compared to the single-moment ensemble. In terms of the simulated differential reflectivity (ZDR) and specific differential phase (KDP) fields, the double-moment ensemble compares more realistically to the observations and better distinguishes the stratiform and convective precipitation regions. The ZDR from individual ensemble members indicates better raindrop size sorting along the leading convective line in the double-moment ensemble. Various commonly used ensemble forecast verification methods are examined for the prediction of dual-pol variables. The results demonstrate the challenges associated with verifying predicted dual-pol fields that can vary significantly in value over small distances. Several microphysics biases are noted with the help of simulated dual-pol variables, such as substantial overprediction of KDP values in the single-moment ensemble.


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.


2014 ◽  
Vol 142 (9) ◽  
pp. 3243-3263 ◽  
Author(s):  
Dustan M. Wheatley ◽  
Nusrat Yussouf ◽  
David J. Stensrud

A Weather Research and Forecasting Model (WRF)-based ensemble data assimilation system is used to produce storm-scale analyses and forecasts of the 4–5 July 2003 severe mesoscale convective system (MCS) over Indiana and Ohio, which produced numerous high wind reports across the two states. Single-Doppler observations are assimilated into a 36-member, storm-scale ensemble during the developing stage of the MCS with the ensemble Kalman filter (EnKF) approach encoded in the Data Assimilation Research Testbed (DART). The storm-scale ensemble is constructed from mesoscale EnKF analyses produced from the assimilation of routinely available observations from land and marine stations, rawinsondes, and aircraft, in an attempt to better represent the complex mesoscale environment for this event. Three EnKF simulations were performed using the National Severe Storms Laboratory (NSSL) one- and two-moment and Thompson microphysical schemes. All three experiments produce a linear convective segment at the final analysis time, similar to the observed system at 2300 UTC 4 July 2003. The higher-order schemes—in particular, the Thompson scheme—are better able to produce short-range forecasts of both the convective and stratiform components of the observed bowing MCS, and produce the smallest temperature errors when comparing surface observations and dropsonde data to corresponding model data. Only the higher-order microphysical schemes produce any appreciable rear-to-front flow in the stratiform precipitation region that trailed the simulated systems. Forecast performance by the three microphysics schemes is discussed in context of differences in microphysical composition produced in the stratiform precipitation regions of the rearward expanding MCSs.


2018 ◽  
Vol 2018 ◽  
pp. 1-18 ◽  
Author(s):  
Shibo Gao ◽  
Jinzhong Min

Using radar observations, the performances of the ensemble square root filter (EnSRF) and an indirect three-dimensional variational (3DVar) data assimilation method were compared for a mesoscale convective system (MCS) that occurred in the Front Range of the Rocky Mountains, Colorado (USA). The results showed that the root mean square innovations (RMSIs) of EnSRF were lower than 3DVar for radar reflectivity and radial velocity and that the spread of EnSRF was generally consistent with its RMSIs. EnSRF substantially improved the analysis of the MCS compared with an experiment without radar data assimilation, and it produced a slight but noticeable improvement over 3DVar in terms of both coverage and intensity. Forecast results initiated from the final analysis revealed that EnSRF generally produced the best prediction of the MCS, with improved quantitative reflectivity and precipitation forecast skills. EnSRF also demonstrated better performance than 3DVar in the prediction of neighborhood probability for reflectivity at thresholds of 20 and 35 dBZ, which better matched the observed radar reflectivity in terms of both shape and extension. Additionally, the humidity, temperature, and wind fields were also improved by EnSRF; the largest error reduction was found in the water vapor field near the surface and at upper levels.


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.


2010 ◽  
Vol 138 (4) ◽  
pp. 1273-1292 ◽  
Author(s):  
Altuğ Aksoy ◽  
David C. Dowell ◽  
Chris Snyder

Abstract The quality of convective-scale ensemble forecasts, initialized from analysis ensembles obtained through the assimilation of radar observations using an ensemble Kalman filter (EnKF), is investigated for cases whose behaviors span supercellular, linear, and multicellular organization. This work is the companion to Part I, which focused on the quality of analyses during the 60-min analysis period. Here, the focus is on 30-min ensemble forecasts initialized at the end of that period. As in Part I, the Weather Research and Forecasting (WRF) model is employed as a simplified cloud model at 2-km horizontal grid spacing. Various observation-space and state-space verification metrics, computed both for ensemble means and individual ensemble members, are employed to assess the quality of ensemble forecasts comparatively across cases. While the cases exhibit noticeable differences in predictability, the forecast skill in each case, as measured by various metrics, decays on a time scale of tens of minutes. The ensemble spread also increases rapidly but significant outlier members or clustering among members are not encountered. Forecast quality is seen to be influenced to varying degrees by the respective initial soundings. While radar data assimilation is able to partially mitigate some of the negative effects in some situations, the supercell case, in particular, remains difficult to predict even after 60 min of data assimilation.


2014 ◽  
Vol 142 (1) ◽  
pp. 183-202 ◽  
Author(s):  
Alexandre O. Fierro ◽  
Jidong Gao ◽  
Conrad L. Ziegler ◽  
Edward R. Mansell ◽  
Donald R. MacGorman ◽  
...  

Abstract This work evaluates the short-term forecast (≤6 h) of the 29–30 June 2012 derecho event from the Advanced Research core of the Weather Research and Forecasting Model (WRF-ARW) when using two distinct data assimilation techniques at cloud-resolving scales (3-km horizontal grid). The first technique assimilates total lightning data using a smooth nudging function. The second method is a three-dimensional variational technique (3DVAR) that assimilates radar reflectivity and radial velocity data. A suite of sensitivity experiments revealed that the lightning assimilation was better able to capture the placement and intensity of the derecho up to 6 h of the forecast. All the simulations employing 3DVAR, however, best represented the storm’s radar reflectivity structure at the analysis time. Detailed analysis revealed that a small feature in the velocity field from one of the six selected radars in the original 3DVAR experiment led to the development of spurious convection ahead of the parent mesoscale convective system, which significantly degraded the forecast. Thus, the relatively simple nudging scheme using lightning data complements the more complex variational technique. The much lower computational cost of the lightning scheme may permit its use alongside variational techniques in improving severe weather forecasts on days favorable for the development of outflow-dominated mesoscale convective systems.


Atmosphere ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 95
Author(s):  
Xinyao Qian ◽  
Haoliang Wang

Lightning simulation is important for a variety of applications, including lightning forecast, atmospheric chemical simulation, and lightning data assimilation. In this study, the potential of five storm parameters (graupel volume, precipitation ice mass, radar echo volume, maximum updraft, and updraft volume) to be used as the proxy for the diagnosis of gridded total lightning flash rates has been investigated in a convection-allowing model. A mesoscale convective system occurred in the Guangdong province of China was selected as the test case. Radar data assimilation was used to improve the simulation accuracy of the convective clouds, hence providing strong instantaneous correlations between observed and simulated storm signatures. The areal coverage and magnitude of the simulated lightning flash rates were evaluated by comparing to those of the total lightning observations. Subjective and the Fractions Skill Score (FSS) evaluations suggest that all the five proxies tested in this study are useful to indicate general tendencies for the occurrence, region, and time of lightning at convection-allowing scale (FSS statistics for the threshold of 1 flash per 9 km2 per hour were around 0.7 for each scheme). The FSS values were decreasing as the lightning flash rate thresholds used for FSS computation increased for all the lightning diagnostic schemes with different proxies. For thresholds from 1 to 3 and 16 to 20 flashes per 9 km2 per hour, the graupel contents related schemes achieved higher FSS values compared to the other three schemes. For thresholds from 5 to 15 flashes per 9 km2 per hour, the updraft volume related scheme yielded the largest FSS. When the thresholds of lightning flash rates were greater than 13 flashes per 9 km2 per hour, the FSS values were below 0.5 for all the lightning diagnostic schemes with different proxies.


Sign in / Sign up

Export Citation Format

Share Document