scholarly journals Comparison of Single-Parameter and Multiparameter Ensembles for Assimilation of Radar Observations Using the Ensemble Kalman Filter

2012 ◽  
Vol 140 (2) ◽  
pp. 562-586 ◽  
Author(s):  
Nusrat Yussouf ◽  
David J. Stensrud

Observational studies indicate that the densities and intercept parameters of hydrometeor distributions can vary widely among storms and even within a single storm. Therefore, assuming a fixed set of microphysical parameters within a given microphysics scheme can lead to significant errors in the forecasts of storms. To explore the impact of variations in microphysical parameters, Observing System Simulation Experiments are conducted based on both perfect- and imperfect-model assumptions. Two sets of ensembles are designed using either fixed or variable parameters within the same single-moment microphysics scheme. The synthetic radar observations of a splitting supercell thunderstorm are assimilated into the ensembles over a 30-min period using an ensemble Kalman filter data assimilation technique followed by 1-h ensemble forecasts. Results indicate that in the presence of a model error, a multiparameter ensemble with a combination of different hydrometeor density and intercept parameters leads to improved analyses and forecasts and better captures the truth within the forecast envelope compared to single-parameter ensemble experiments with a single, constant, inaccurate hydrometeor intercept and density parameters. This conclusion holds when examining the general storm structure, the intensity of midlevel rotation, surface cold pool strength, and the extreme values of the model fields that are most helpful in determining and identifying potential hazards. Under a perfect-model assumption, the single- and multiparameter ensembles perform similarly as model error does not play a role in these experiments. This study highlights the potential for using a variety of realistic microphysical parameters across the ensemble members in improving the analyses and very short-range forecasts of severe weather events.

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.


2009 ◽  
Vol 137 (7) ◽  
pp. 2126-2143 ◽  
Author(s):  
P. L. Houtekamer ◽  
Herschel L. Mitchell ◽  
Xingxiu Deng

Since 12 January 2005, an ensemble Kalman filter (EnKF) has been used operationally at the Meteorological Service of Canada to provide the initial conditions for the medium-range forecasts of the ensemble prediction system. One issue in EnKF development is how to best account for model error. It is shown that in a perfect-model environment, without any model error or model error simulation, the EnKF spread remains representative of the ensemble mean error with respect to a truth integration. Consequently, the EnKF can be used to quantify the impact of the various error sources in a data-assimilation cycle on the quality of the ensemble mean. Using real rather than simulated observations, but still not simulating model error in any manner, the rms ensemble spread is found to be too small by approximately a factor of 2. It is then attempted to account for model error by using various combinations of the following four different approaches: (i) additive isotropic model error perturbations; (ii) different versions of the model for different ensemble members; (iii) stochastic perturbations to physical tendencies; and (iv) stochastic kinetic energy backscatter. The addition of isotropic model error perturbations is found to have the biggest impact. The identification of model error sources could lead to a more realistic, likely anisotropic, parameterization. Using different versions of the model has a small but clearly positive impact and consequently both (i) and (ii) are used in the operational EnKF. The use of approaches (iii) and (iv) did not lead to further improvements.


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.


2012 ◽  
Vol 140 (3) ◽  
pp. 841-859 ◽  
Author(s):  
Yonghui Weng ◽  
Fuqing Zhang

Abstract Through a Weather Research and Forecasting model (WRF)-based ensemble Kalman filter (EnKF) data assimilation system, the impact of assimilating airborne radar observations for the convection-permitting analysis and prediction of Hurricane Katrina (2005) is examined in this study. A forecast initialized from EnKF analyses of airborne radar observations had substantially smaller hurricane track forecast errors than NOAA’s operational forecasts and a control forecast initialized from NCEP analysis data for lead times up to 120 h. Verifications against independent in situ and remotely sensed observations show that EnKF analyses successfully depict the inner-core structure of the hurricane vortex in terms of both dynamic (wind) and thermodynamic (temperature and moisture) fields. In addition to the improved analyses and deterministic forecast, an ensemble of forecasts initiated from the EnKF analyses also provided forecast uncertainty estimates for the hurricane track and intensity. Also documented here are the details of a series of data thinning and quality control procedures that were developed to generate superobservations from large volumes of airborne radial velocity measurements. These procedures have since been implemented operationally on the NOAA hurricane reconnaissance aircraft that allows for more efficient real-time transmission of airborne radar observations to the ground.


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.


2015 ◽  
Vol 143 (8) ◽  
pp. 3067-3086 ◽  
Author(s):  
Ryan A. Sobash ◽  
Louis J. Wicker

Abstract Storm-scale ensemble Kalman filter (EnKF) studies routinely use methods to accelerate the spinup of convective structures when assimilating convective-scale radar observations. This typically involves adding coherent perturbations into analyses at regular intervals in regions where radar observations indicate convection is ongoing. Significant uncertainty remains as to the most effective use of these perturbations, including appropriate perturbation magnitudes, spatial scales, fields, and smoothing kernels, as well as flexible strategies that can be applied across a spectrum of convective events with negligible a priori tuning. Here, several idealized experiments were performed to elucidate the impact and sensitivity of adding coherent perturbations into storm-scale analyses of convection. Through the use of toy experiments, it is demonstrated that various factors exhibit substantial influence on the postsmoothed perturbation magnitudes, making tuning challenging. Several OSSEs were performed to document the impact of these perturbations on the analyses, particularly thermodynamic analyses within convection. The repeated addition of coherent perturbations produced temperature and moisture biases that are most pronounced in analyses of the surface cold pool and aloft near the tropopause, and eventually lead to biases in the dynamic fields. In an attempt to reduce these biases and make the noise procedure more adaptive, reflectivity innovations were used to restrict the addition of noise to areas where these innovations are large. This produced analyses with reduced thermodynamic biases and RMSE values comparable to the best-performing experiment where the noise magnitudes were manually adjusted. The impact of these findings on previous and future convective-scale EnKF analyses and forecasts are discussed.


2012 ◽  
Vol 69 (11) ◽  
pp. 3147-3171 ◽  
Author(s):  
Humberto C. Godinez ◽  
Jon M. Reisner ◽  
Alexandre O. Fierro ◽  
Stephen R. Guimond ◽  
Jim Kao

Abstract In this work the authors determine key model parameters for rapidly intensifying Hurricane Guillermo (1997) using the ensemble Kalman filter (EnKF). The approach is to utilize the EnKF as a tool only to estimate the parameter values of the model for a particular dataset. The assimilation is performed using dual-Doppler radar observations obtained during the period of rapid intensification of Hurricane Guillermo. A unique aspect of Guillermo was that during the period of radar observations strong convective bursts, attributable to wind shear, formed primarily within the eastern semicircle of the eyewall. To reproduce this observed structure within a hurricane model, background wind shear of some magnitude must be specified and turbulence and surface parameters appropriately specified so that the impact of the shear on the simulated hurricane vortex can be realized. To identify the complex nonlinear interactions induced by changes in these parameters, an ensemble of model simulations have been conducted in which individual members were formulated by sampling the parameters within a certain range via a Latin hypercube approach. The ensemble and the data, derived latent heat and horizontal winds from the dual-Doppler radar observations, are utilized in the EnKF to obtain varying estimates of the model parameters. The parameters are estimated at each time instance, and a final parameter value is obtained by computing the average over time. Individual simulations were conducted using the estimates, with the simulation using latent heat parameter estimates producing the lowest overall model forecast error.


2009 ◽  
Vol 6 (4) ◽  
pp. 8279-8309 ◽  
Author(s):  
W. Ju ◽  
S. Wang ◽  
G. Yu ◽  
Y. Zhou ◽  
H. Wang

Abstract. Soil and atmospheric water deficits have significant influences on CO2 and energy exchanges between the atmosphere and terrestrial ecosystems. Model parameterization significantly affects the ability of a model to simulate carbon, water, and energy fluxes. In this study, an ensemble Kalman filter (EnKF) and observations of gross primary productivity (GPP) and latent heat (LE) fluxes were used to optimize model parameters significantly affecting the calculation of these fluxes for a subtropical coniferous plantation in southeastern China. The optimized parameters include the maximum carboxylation rate (Vcmax), the Ball-Berry coefficient (m) and the coefficient determining the sensitivity of stomatal conductance to atmospheric water vapor deficit D0). Optimized Vcmax and m showed larger seasonal and interannual variations than D0. Seasonal variations of Vcmax and m are more pronounced than the interannual variations. Vcmax and m are associated with soil water content (SWC). During dry periods, SWC at the 20 cm depth can explain 61% and 64% of variations of Vcmax and m, respectively. EnKF parameter optimization improves the simulations of GPP, LE and sensible heat (SH), mainly during dry periods. After parameter optimization using EnKF, the variations of GPP, LE and SH explained by the model increased by 1% to 4% at half-hourly steps and by 3% to 5% at daily time steps. Efforts are needed to develop algorithms that can properly describe the variations of these parameters under different environmental conditions.


2008 ◽  
Vol 8 (11) ◽  
pp. 2975-2983 ◽  
Author(s):  
C. Lin ◽  
Z. Wang ◽  
J. Zhu

Abstract. An Ensemble Kalman Filter (EnKF) data assimilation system was developed for a regional dust transport model. This paper applied the EnKF method to investigate modeling of severe dust storm episodes occurring in March 2002 over China based on surface observations of dust concentrations to explore the impact of the EnKF data assimilation systems on forecast improvement. A series of sensitivity experiments using our system demonstrates the ability of the advanced EnKF assimilation method using surface observed PM10 in North China to correct initial conditions, which leads to improved forecasts of dust storms. However, large errors in the forecast may arise from model errors (uncertainties in meteorological fields, dust emissions, dry deposition velocity, etc.). This result illustrates that the EnKF requires identification and correction model errors during the assimilation procedure in order to significantly improve forecasts. Results also show that the EnKF should use a large inflation parameter to obtain better model performance and forecast potential. Furthermore, the ensemble perturbations generated at the initial time should include enough ensemble spreads to represent the background error after several assimilation cycles.


2020 ◽  
Vol 148 (6) ◽  
pp. 2365-2389
Author(s):  
Jonathan Labriola ◽  
Nathan Snook ◽  
Youngsun Jung ◽  
Ming Xue

Abstract Ensemble Kalman filter (EnKF) analyses of the storms associated with the 8 May 2017 Colorado severe hail event using either the Milbrandt and Yau (MY) or the NSSL double-moment bulk microphysics scheme in the forecast model are evaluated. With each scheme, two experiments are conducted in which the reflectivity (Z) observations update in addition to dynamic and thermodynamic variables: 1) only the hydrometeor mixing ratios or 2) all microphysical variables. With fewer microphysical variables directly constrained by the Z observations, only updating hydrometeor mixing ratios causes the forecast error covariance structure to become unreliable, and results in larger errors in the analysis. Experiments that update all microphysical variables produce analyses with the lowest Z root-mean-square innovations; however, comparing the estimated hail size against hydrometeor classification algorithm output suggests that further constraint from observations is needed to more accurately estimate surface hail size. Ensemble correlation analyses are performed to determine the impact of hail growth assumptions in the MY and NSSL schemes on the forecast error covariance between microphysical and thermodynamic variables. In the MY scheme, Z is negatively correlated with updraft intensity because the strong updrafts produce abundant small hail aloft. The NSSL scheme predicts the growth of large hail aloft; consequently, Z is positively correlated with storm updraft intensity and hail state variables. Hail production processes are also shown to alter the background error covariance for liquid and frozen hydrometeor species. Results in this study suggest that EnKF analyses are sensitive to the choice of MP scheme (e.g., the treatment of hail growth processes).


Sign in / Sign up

Export Citation Format

Share Document