Ensemble Transform Kalman Filter Data Research Articles

Reconstructing cardiac electrical excitations from optical mapping recordings.

The reconstruction of electrical excitation patterns through the unobserved depth of the tissue is essential to realizing the potential of computational models in cardiac medicine. We have utilized experimental optical-mapping recordings of cardiac electrical excitation on the epicardial and endocardial surfaces of a canine ventricle as observations directing a local ensemble transform Kalman filter data assimilation scheme. We demonstrate that the inclusion of explicit information about the stimulation protocol can marginally improve the confidence of the ensemble reconstruction and the reliability of the assimilation over time. Likewise, we consider the efficacy of stochastic modeling additions to the assimilation scheme in the context of experimentally derived observation sets. Approximation error is addressed at both the observation and modeling stages through the uncertainty of observations and the specification of the model used in the assimilation ensemble. We find that perturbative modifications to the observations have marginal to deleterious effects on the accuracy and robustness of the state reconstruction. Furthermore, we find that incorporating additional information from the observations into the model itself (in the case of stimulus and stochastic currents) has a marginal improvement on the reconstruction accuracy over a fully autonomous model, while complicating the model itself and thus introducing potential for new types of model errors. That the inclusion of explicit modeling information has negligible to negative effects on the reconstruction implies the need for new avenues for optimization of data assimilation schemes applied to cardiac electrical excitation.

Assimilation of radar‐derived refractivity and radar data in the context of ensemble Kalman filter: Cases study of the Southwest Monsoon Experiment

AbstractRadar‐derived refractivity provides moisture information near the surface, which plays a crucial role in convective initiation and of heavy rain. In this study, the high‐resolution Weather Research and Forecasting local ensemble transform Kalman filter data assimilation system was employed for two cases in the Southwest Monsoon Experiment to conduct two experimental sets. The first set was applied in both cases to investigate the effects of assimilating radar‐retrieved refractivity along with reflectivity and radial wind. The second set was conducted in the second case to examine the benefit of increasing the frequency of refractivity assimilation and investigate the optimal strategy to assimilate refractivity. Results of the two cases revealed that assimilating reflectivity and radial velocity modified near‐surface humidity on the basis of the flow‐dependent error correlation estimated by the ensemble, but the spatial distribution may not be fully accurate, causing underestimation of rainfall. With further refractivity assimilation, stronger convergence and more accurate low‐level moisture, temperature, and wind field corrections were obtained, leading to superior forecasts for both light and heavy rainfall during six hours. The results of the second set indicated that increasing the assimilation interval of refractivity enabled capturing the dramatic moisture variation and enhancing wind convergence, resulting in short‐term forecast improvement. The strategy that solely assimilated refractivity before the appearance of the convection system optimized the correction of environmental moisture, then accurately represented the humidity and strengthened the wind convergence when precipitation occurred. Consequently, the achieved improvement of the short‐term forecast was most noteworthy, particularly for heavy rain.

УЧЕТ СКОРРЕЛИРОВАННОСТИ ОШИБОК СПУТНИКОВЫХ ДАННЫХ НАБЛЮДЕНИЙ AMV В АНСАМБЛЕВОЙ СИСТЕМЕ УСВОЕНИЯ ДАННЫХ

The effect of considering correlated errors in AMV (Atmospheric Motion Vectors) satellite observations of wind in the local ensemble transform Kalman filter data assimilation system is studied. It is customary to use a diagonal covariance matrix of errors in observations taking part in assimilation. When assimilating satellite data with correlated errors, the data are thinned, and the values of diagonal elements of the error covariance matrix are often overestimated. This is accompanied by the loss of useful information about the correlation between the errors. The present study uses a different method: the elements of the error covariance matrix for AMV satellite observations are simulated using a second-order autoregressive function. It is shown that such approach reduces the root-mean-square error in initial data for a numerical weather prediction model, in particular on small scales, and improves the forecast quality. It is found that the application of the non-diagonal AMV observation error covariance matrix increases the accuracy of analysis and forecast fields.

Assimilating Retrieved Water Vapor and Radar Data from NCAR S-PolKa: Performance and Validation Using Real Cases

Abstract This study investigated the effect of the assimilation of the S- and Ka-band dual‐wavelength-retrieved water vapor data with radial wind and reflectivity data. The vertical profile of humidity, which provides environmental information before precipitation occurs, was obtained at low levels and thinned into averaged and four-quadrant profiles. Additionally, the following two strategies were examined: 1) assimilation of water vapor data with radar data for the entire 2 h and 2) assimilation of water vapor data in the first hour, and radial velocity and reflectivity data in the second hour. By using the WRF local ensemble transform Kalman filter data assimilation system, three real cases of the Dynamics of the Madden–Julian Oscillation experiment were examined through a series of experiments. The analysis results revealed that assimilating additional water vapor data more markedly improved the analysis at the convective scale than assimilating radial wind and reflectivity data alone. In addition, the strategy of assimilating only retrieved water vapor data in the first hour and radial wind and reflectivity data in the second hour achieved the optimal analysis and subsequent very short-term forecast. The evaluation of quantitative precipitation forecasting demonstrated that assimilating additional retrieved water vapor data distinctly improved the rain forecast compared with assimilating radar data only. When moisture data were assimilated, improved nowcasting could be extended up to 4 h. Furthermore, assimilating moisture profiles into four quadrants achieved more accurate analysis and forecast. Overall, our study demonstrated that the humidify information in nonprecipitation areas is critical for further improving the analysis and forecast of convective weather systems.

Impact of Data Assimilation on Thermal Tides in the Case of Venus Express Wind Observation

AbstractImpacts of horizontal winds assimilation on thermal tides are investigated by using the Venus atmospheric general circulation model for the Earth Simulator local ensemble transform Kalman filter data assimilation system. The assimilated data are horizontal winds derived from Venus ultraviolet images taken by the Venus Monitoring Camera onboard the Venus Express orbiter. The results show that three‐dimensional structures of the thermal tides are significantly improved not only in the horizontal winds but also in the temperature field, even though the observations are available only at the cloud top level on the southern dayside hemisphere approximately once an Earth day. The reproduced temperature fields agree well with recent radio occultation measurements of the Venus Climate Orbiter Akatsuki. The zonal mean fields of the zonal wind and temperature are also improved globally. This study would enable reanalysis of past Venus observations.

Application of ensemble Kalman filter in WRF model to forecast rainfall on monsoon onset period in South Vietnam

Application of ensemble Kalman filter in WRF model to forecast rainfall on monsoon onset period in South Vietnam

Local ensemble transform Kalman filter data assimilation system for the global semi-Lagrangian atmospheric model

Recently, ensemble Kalman filters have come into practical data assimilation for numerical weather prediction models. We give an overview of ensemble Kalman filters and problems that arise with practical implementation of ensemble methods. We present our implementation of the local ensemble transform Kalman filter, one of ensemble square root filters using observation localization. Multiplicative and additive inflations are used to prevent filter divergence and to account for the model error. The implemented assimilation system is tested with the global semi-Lagrangian atmospheric model SL-AV using real observations for 2 months of cyclic assimilation (August and September 2012). The system works stably. Application of the ensemble filter significantly reduces first guess (background) errors and corrects the forecast biases.

A local ensemble transform Kalman filter data assimilation system for the NCEP global model

The accuracy and computational efficiency of a parallel computer implementation of the Local Ensemble Transform Kalman Filter (LETKF) data assimilation scheme on the model component of the 2004 version of the Global Forecast System (GFS) of the National Centers for Environmental Prediction (NCEP) is investigated.Numerical experiments are carried out at model resolution T62L28. All atmospheric observations that were operationally assimilated by NCEP in 2004, except for satellite radiances, are assimilated with the LETKF. The accuracy of the LETKF analyses is evaluated by comparing it to that of the Spectral Statistical Interpolation (SSI), which was the operational global data assimilation scheme of NCEP in 2004. For the selected set of observations, the LETKF analyses are more accurate than the SSI analyses in the Southern Hemisphere extratropics and are comparably accurate in the Northern Hemisphere extratropics and in the Tropics.The computationalwall-clock times achieved on a Beowulf cluster of 3.6 GHz Xeon processors make our implementation of the LETKF on the NCEP GFS a widely applicable analysis-forecast system, especially for research purposes. For instance, the generation of four daily analyses at the resolution of the NCAR-NCEP reanalysis (T62L28) for a full season (90 d), using 40 processors, takes less than 4 d of wall-clock time.

A local ensemble transform Kalman filter data assimilation system for the NCEP global model

A local ensemble transform Kalman filter data assimilation system for the NCEP global model