scholarly journals Analysis of a Tornadic Mesoscale Convective Vortex Based on Ensemble Kalman Filter Assimilation of CASA X-Band and WSR-88D Radar Data

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.

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 (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.


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.


2015 ◽  
Vol 143 (1) ◽  
pp. 165-194 ◽  
Author(s):  
Thomas A. Jones ◽  
David Stensrud ◽  
Louis Wicker ◽  
Patrick Minnis ◽  
Rabindra Palikonda

Abstract Assimilating high-resolution radar reflectivity and radial velocity into convection-permitting numerical weather prediction models has proven to be an important tool for improving forecast skill of convection. The use of satellite data for the application is much less well understood, only recently receiving significant attention. Since both radar and satellite data provide independent information, combing these two sources of data in a robust manner potentially represents the future of high-resolution data assimilation. This research combines Geostationary Operational Environmental Satellite 13 (GOES-13) cloud water path (CWP) retrievals with Weather Surveillance Radar-1988 Doppler (WSR-88D) reflectivity and radial velocity to examine the impacts of assimilating each for a severe weather event occurring in Oklahoma on 24 May 2011. Data are assimilated into a 3-km model using an ensemble adjustment Kalman filter approach with 36 members over a 2-h assimilation window between 1800 and 2000 UTC. Forecasts are then generated for 90 min at 5-min intervals starting at 1930 and 2000 UTC. Results show that both satellite and radar data are able to initiate convection, but that assimilating both spins up a storm much faster. Assimilating CWP also performs well at suppressing spurious precipitation and cloud cover in the model as well as capturing the anvil characteristics of developed storms. Radar data are most effective at resolving the 3D characteristics of the core convection. Assimilating both satellite and radar data generally resulted in the best model analysis and most skillful forecast for this event.


2014 ◽  
Vol 142 (6) ◽  
pp. 2118-2138 ◽  
Author(s):  
Weiguang Chang ◽  
Kao-Shen Chung ◽  
Luc Fillion ◽  
Seung-Jong Baek

Abstract An 80-member high-resolution ensemble Kalman filter (HREnKF) is implemented for assimilating radar observations with the Canadian Meteorological Center’s (CMC’s) Global Environmental Multiscale Limited-Area Model (GEM-LAM). This system covers the Montréal, Canada, region and assimilates radar data from the McGill Radar Observatory with 4-km data thinning. The GEM-LAM operates in fully nonhydrostatic mode with 58 hybrid vertical levels and 1-km horizontal grid spacing. As a first step toward full radar data assimilation, only radial velocities are directly assimilated in this study. The HREnKF is applied on three 2011 summer cases having different precipitation structures: squall-line structure, isolated small-scale structures, and widespread stratiform precipitation. The short-term (<2 h) accuracy of the HREnKF analyses and forecasts is examined. In HREnKF, the ensemble spread is sufficient to cover the estimated error from innovations and lead to filter convergence. It results in part from a realistic initiation of HREnKF data assimilation cycle by using a Canadian regional EnKF system (itself coupled to a global EnKF) working at meso- and synoptic scales. The filter convergence is confirmed by the HREnKF background fields gradually approaching to radar observations as the assimilation cycling proceeds. At each analysis step, it is clearly shown that unobserved variables are significantly modified through HREnKF cross correlation of errors from the ensemble. Radar reflectivity observations are used to verify the improvements in analyses and short-term forecasts achievable by assimilating only radial velocities. Further developments of the analysis system are discussed.


2011 ◽  
Vol 139 (1) ◽  
pp. 224-246 ◽  
Author(s):  
Alexander D. Schenkman ◽  
Ming Xue ◽  
Alan Shapiro ◽  
Keith Brewster ◽  
Jidong Gao

Abstract The Advanced Regional Prediction System (ARPS) model is employed to perform high-resolution numerical simulations of a mesoscale convective system and associated cyclonic line-end vortex (LEV) that spawned several tornadoes in central Oklahoma on 8–9 May 2007. The simulation uses a 1000 km × 1000 km domain with 2-km horizontal grid spacing. The ARPS three-dimensional variational data assimilation (3DVAR) is used to assimilate a variety of data types. All experiments assimilate routine surface and upper-air observations as well as wind profiler and Oklahoma Mesonet data over a 1-h assimilation window. A subset of experiments assimilates radar data. Cloud and hydrometeor fields as well as in-cloud temperature are adjusted based on radar reflectivity data through the ARPS complex cloud analysis procedure. Radar data are assimilated from the Weather Surveillance Radar-1988 Doppler (WSR-88D) network as well as from the Engineering Research Center for Collaborative and Adaptive Sensing of the Atmosphere (CASA) network of four X-band Doppler radars. Three-hour forecasts are launched at the end of the assimilation window. The structure and evolution of the forecast MCS and LEV are markedly better throughout the forecast period in experiments in which radar data are assimilated. The assimilation of CASA radar data in addition to WSR-88D data increases the structural detail of the modeled squall line and MCS at the end of the assimilation window, which appears to yield a slightly better forecast track of the LEV.


2012 ◽  
Vol 140 (12) ◽  
pp. 3992-4004 ◽  
Author(s):  
Brian C. Ancell

Abstract Mesoscale atmospheric data assimilation is becoming an integral part of numerical weather prediction. Modern computational resources now allow assimilation and subsequent forecasting experiments ranging from resolutions of tens of kilometers over regional domains to smaller grids that employ storm-scale assimilation. To assess the value of the high-resolution capabilities involved with assimilation and forecasting at different scales, analyses and forecasts must be carefully evaluated to understand 1) whether analysis benefits gained at finer scales persist into the forecast relative to downscaled runs begun from lower-resolution analyses, 2) how the positive analysis effects of bias removal evolve into the forecast, and 3) how digital filter initialization affects analyses and forecasts. This study applies a 36- and 4-km ensemble Kalman filter over 112 assimilation cycles to address these important issues, which could all be relevant to a variety of short-term, high-resolution, real-time forecasting applications. It is found that with regard to surface wind and temperature, analysis improvements gained at higher resolution persist throughout the 12-h forecast window relative to downscaled, high-resolution forecasts begun from analyses on the coarser grid. Aloft, however, no forecast improvements were found with the high-resolution analysis/forecast runs. Surface wind and temperature bias removal, while clearly improving surface analyses, degraded surface forecasts and showed little forecast influence aloft. Digital filter initialization degraded temperature analyses with or without bias removal, degraded wind analyses when bias removal was used, but improved wind analyses when bias removal was absent. No forecast improvements were found with digital filter initialization. The consequences of these results with regard to operational assimilation/forecasting systems on nested grids are discussed.


2015 ◽  
Vol 143 (2) ◽  
pp. 511-523 ◽  
Author(s):  
Sim D. Aberson ◽  
Altuğ Aksoy ◽  
Kathryn J. Sellwood ◽  
Tomislava Vukicevic ◽  
Xuejin Zhang

Abstract NOAA has been gathering high-resolution, flight-level dropwindsonde and airborne Doppler radar data in tropical cyclones for almost three decades; the U.S. Air Force routinely obtained the same type and quality of data, excepting Doppler radar, for most of that time. The data have been used for operational diagnosis and for research, and, starting in 2013, have been assimilated into operational regional tropical cyclone models. This study is an effort to quantify the impact of assimilating these data into a version of the operational Hurricane Weather Research and Forecasting model using an ensemble Kalman filter. A total of 83 cases during 2008–11 were investigated. The aircraft whose data were used in the study all provide high-density flight-level wind and thermodynamic observations as well as surface wind speed data. Forecasts initialized with these data assimilated are compared to those using the model standard initialization. Since only NOAA aircraft provide airborne Doppler radar data, these data are also tested to see their impact above the standard aircraft data. The aircraft data alone are shown to provide some statistically significant improvement to track and intensity forecasts during the critical watch and warning period before projected landfall (through 60 h), with the Doppler radar data providing some further improvement. This study shows the potential for improved forecasts with regular tropical cyclone aircraft reconnaissance and the assimilation of data obtained from them, especially airborne Doppler radar data, into the numerical guidance.


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.


Sign in / Sign up

Export Citation Format

Share Document