scholarly journals Accelerating the EnKF Spinup for Typhoon Assimilation and Prediction

2012 ◽  
Vol 27 (4) ◽  
pp. 878-897 ◽  
Author(s):  
Shu-Chih Yang ◽  
Eugenia Kalnay ◽  
Takemasa Miyoshi

Abstract A mesoscale ensemble Kalman filter (EnKF) for a regional model is often initialized from global analysis products and with initial ensemble perturbations constructed based on the background error covariance used in the three-dimensional variational data assimilation (3DVar) system. Because of the lack of proper mesoscale information, a long spinup period of typically a few days is required for the regional EnKF to reach its asymptotic level of accuracy, and thus, the impact of observations is limited during the EnKF spinup. For the case of typhoon assimilation, such spinup usually corresponds to the stages of generation and development of tropical cyclones, when observations are important but limited over open waters. To improve the analysis quality during the spinup, the “running in place” (RIP) method is implemented within the framework of the local ensemble transform Kalman filter (LETKF) coupled with the Weather Research and Forecasting model (WRF). Results from observing system simulation experiments (OSSEs) for a specific typhoon show that the RIP method is able to accelerate the analysis adjustment of the dynamical structures of the typhoon during the LETKF spinup, and improves both the accuracy of the mean state and the structure of the ensemble-based error covariance. These advantages of the RIP method are found not only in the inner-core structure of the typhoon but also identified in the environmental conditions. As a result, the LETKF-RIP analysis leads to better typhoon prediction, particularly in terms of both track and intensity.

Atmosphere ◽  
2021 ◽  
Vol 12 (8) ◽  
pp. 942
Author(s):  
Benjamin Davis ◽  
Xuguang Wang ◽  
Xu Lu

Six-hourly three-dimensional ensemble variational (3DEnVar) (6H-3DEnVar) data assimilation (DA) assumes constant background error covariance (BEC) during a six-hour DA window and is, therefore, unable to account for temporal evolution of the BEC. This study evaluates the one-hourly 3DEnVar (1H-3DEnVar) and six-hourly 4DEnVar (6H-4DEnVar) DA methods for the analyses and forecasts of hurricanes with rapidly evolving BEC. Both methods account for evolving BEC in a hybrid EnVar DA system. In order to compare these methods, experiments are conducted by assimilating inner core Tail Doppler Radar (TDR) wind for Hurricane Edouard (2014) and by running the Hurricane Weather Research and Forecasting (HWRF) model. In most metrics, 1H-3DEnVar and 6H-4DEnVar analyses and forecasts verify better than 6H-3DEnVar. 6H-4DEnVar produces better thermodynamic analyses than 1H-3DEnVar. Radar reflectivity shows that 1H-3DEnVar produces better structure forecasts. For the first 24–48 h of the intensity forecast, 6H-4DEnVar forecast performs better than 1H-3DEnVar verified against the best track. Degraded 1H-3DEnVar forecasts are found to be associated with background storm center location error as a result of underdispersive ensemble storm center spread. Removing location error in the background improves intensity forecasts of 1H-3DEnVar.


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


2011 ◽  
Vol 2011 ◽  
pp. 1-10 ◽  
Author(s):  
Chien-Ben Chou ◽  
Huei-Ping Huang

This work assesses the effects of assimilating atmospheric infrared sounder (AIRS) observations on typhoon prediction using the three-dimensional variational data assimilation (3DVAR) and forecasting system of the weather research and forecasting (WRF) model. Two major parameters in the data assimilation scheme, the spatial decorrelation scale and the magnitude of the covariance matrix of the background error, are varied in forecast experiments for the track of typhoon Sinlaku over the Western Pacific. The results show that within a wide parameter range, the inclusion of the AIRS observation improves the prediction. Outside this range, notably when the decorrelation scale of the background error is set to a large value, forcing the assimilation of AIRS data leads to degradation of the forecast. This illustrates how the impact of satellite data on the forecast depends on the adjustable parameters for data assimilation. The parameter-sweeping framework is potentially useful for improving operational typhoon prediction.


2013 ◽  
Vol 2013 ◽  
pp. 1-6 ◽  
Author(s):  
Hongze Leng ◽  
Junqiang Song ◽  
Fengshun Lu ◽  
Xiaoqun Cao

This study considers a new hybrid three-dimensional variational (3D-Var) and ensemble Kalman filter (EnKF) data assimilation (DA) method in a non-perfect-model framework, named space-expanded ensemble localization Kalman filter (SELKF). In this method, the localization operation is directly applied to the ensemble anomalies with a Schur Product, rather than to the full error covariance of the state in the EnKF. Meanwhile, the correction space of analysis increment is expanded to a space with larger dimension, and the rank of the forecast error covariance is significantly increased. This scheme can reduce the spurious correlations in the covariance and approximate the full-rank background error covariance well. Furthermore, a deterministic scheme is used to generate the analysis anomalies. The results show that the SELKF outperforms the perturbed EnKF given a relatively small ensemble size, especially when the length scale is relatively long or the observation error covariance is relatively small.


2011 ◽  
Vol 26 (6) ◽  
pp. 868-884 ◽  
Author(s):  
Xuguang Wang

Abstract A hybrid ensemble transform Kalman filter (ETKF)–three-dimensional variational data assimilation (3DVAR) system developed for the Weather Research and Forecasting Model (WRF) was studied for the forecasts of the tracks of two major hurricanes, Ike and Gustav, in 2008 over the Gulf of Mexico. The impacts of the flow-dependent ensemble covariance generated by the ETKF were revealed by comparing the forecasts, analyses, and analysis increments generated by the hybrid data assimilation method with those generated by the 3DVAR that used the static background covariance. The root-mean-square errors of the track forecasts by the hybrid data assimilation (DA) method were smaller than those by the 3DVAR for both Ike and Gustav. Experiments showed that such improvements were due to the use of the flow-dependent covariance provided by the ETKF ensemble in the hybrid DA system. Detailed diagnostics further revealed that the increments produced by the hybrid and the 3DVAR were different for both the analyses of the hurricane itself and its environment. In particular, it was found that the hybrid, using the flow-dependent covariance that gave the hurricane-specific error covariance estimates, was able to systematically adjust the position of the hurricane during the assimilation whereas the 3DVAR was not. The study served as a pilot study to explore and understand the potential of the hybrid method for hurricane data assimilation and forecasts. Caution needs to be taken to extrapolate the results to operational forecast settings.


2009 ◽  
Vol 137 (11) ◽  
pp. 3918-3932 ◽  
Author(s):  
Junjie Liu ◽  
Hong Li ◽  
Eugenia Kalnay ◽  
Eric J. Kostelich ◽  
Istvan Szunyogh

Abstract This study uses the local ensemble transform Kalman filter to assimilate Atmospheric Infrared Sounder (AIRS) specific humidity retrievals with pseudo relative humidity (pseudo-RH) as the observation variable. Three approaches are tested: (i) updating specific humidity with observations other than specific humidity (“passive q”), (ii) updating specific humidity only with humidity observations (“univariate q”), and (iii) assimilating the humidity and the other observations together (“multivariate q”). This is the first time that the performance of the univariate and multivariate assimilation of q is compared within an ensemble Kalman filter framework. The results show that updating the humidity analyses by either AIRS specific humidity retrievals or nonhumidity observations improves both the humidity and wind analyses. The improvement with the multivariate-q experiment is by far the largest for all dynamical variables at both analysis and forecast time, indicating that the interaction between the specific humidity and the other dynamical variables through the background error covariance during data assimilation process yields more balanced analysis fields. In the univariate assimilation of q, the humidity interacts with the other dynamical variables only through the forecast process. The univariate assimilation produces more accurate humidity analyses than those obtained when no humidity observations are assimilated, but it does not improve the accuracy of the zonal wind analyses. The 6-h total column precipitable water forecast also benefits from the improved humidity analyses, with the multivariate q experiment having the largest improvement.


2007 ◽  
Vol 135 (4) ◽  
pp. 1506-1521 ◽  
Author(s):  
Haixia Liu ◽  
Ming Xue ◽  
R. James Purser ◽  
David F. Parrish

Abstract Anisotropic recursive filters are implemented within a three-dimensional variational data assimilation (3DVAR) framework to efficiently model the effect of flow-dependent background error covariance. The background error covariance is based on an estimated error field and on the idea of Riishøjgaard. In the anisotropic case, the background error pattern can be stretched or flattened in directions oblique to the alignment of the grid coordinates and is constructed by applying, at each point, six recursive filters along six directions corresponding, in general, to a special configuration of oblique lines of the grid. The recursive filters are much more efficient than corresponding explicit filters used in an earlier study and are therefore more suitable for real-time numerical weather prediction. A set of analysis experiments are conducted at a mesoscale resolution to examine the effectiveness of the 3DVAR system in analyzing simulated global positioning system (GPS) slant-path water vapor observations from ground-based GPS receivers and observations from collocated surface stations. It is shown that the analyses produced with recursive filters are at least as good as those with corresponding explicit filters. In some cases, the recursive filters actually perform better. The impact of flow-dependent background errors modeled using the anisotropic recursive filters is also examined. The use of anisotropic filters improves the analysis, especially in terms of finescale structures. The analysis system is found to be effective in the presence of typical observational errors. The sensitivity of isotropic and anisotropic recursive-filter analyses to the decorrelation scales is also examined systematically.


2016 ◽  
Vol 144 (12) ◽  
pp. 4849-4865 ◽  
Author(s):  
Keiichi Kondo ◽  
Takemasa Miyoshi

Abstract The ensemble Kalman filter (EnKF) with high-dimensional geophysical systems usually employs up to 100 ensemble members and requires covariance localization to reduce the sampling error in the forecast error covariance between distant locations. The authors’ previous work pioneered implementation of an EnKF with a large ensemble of up to 10 240 members, but this method required application of a relatively broad covariance localization to avoid memory overflow. This study modified the EnKF code to save memory and enabled for the first time the removal of completely covariance localization with an intermediate AGCM. Using the large sample size, this study aims to investigate the analysis and forecast accuracy, as well as the impact of covariance localization when the sampling error is small. A series of 60-day data assimilation cycle experiments with different localization scales are performed under the perfect model scenario to investigate the pure impact of covariance localization. The results show that the analysis and 7-day forecasts are much improved by removing covariance localization and that the long-range covariance between distant locations plays a key role. The eigenvectors of the background error covariance matrix based on the 10 240 ensemble members are explicitly computed and reveal long-range structures.


2014 ◽  
Vol 142 (6) ◽  
pp. 2139-2149 ◽  
Author(s):  
Stephen G. Penny

Abstract Hybrid data assimilation methods combine elements of ensemble Kalman filters (EnKF) and variational methods. While most approaches have focused on augmenting an operational variational system with dynamic error covariance information from an ensemble, this study takes the opposite perspective of augmenting an operational EnKF with information from a simple 3D variational data assimilation (3D-Var) method. A class of hybrid methods is introduced that combines the gain matrices of the ensemble and variational methods, rather than linearly combining the respective background error covariances. A hybrid local ensemble transform Kalman filter (Hybrid-LETKF) is presented in two forms: 1) a traditionally motivated Hybrid/Covariance-LETKF that combines the background error covariance matrices of LETKF and 3D-Var, and 2) a simple-to-implement algorithm called the Hybrid/Mean-LETKF that falls into the new class of hybrid gain methods. Both forms improve analysis errors when using small ensemble sizes and low observation coverage versus either LETKF or 3D-Var used alone. The results imply that for small ensemble sizes, allowing a solution to be found outside of the space spanned by ensemble members provides robustness in both hybrid methods compared to LETKF alone. Finally, the simplicity of the Hybrid/Mean-LETKF design implies that this algorithm can be applied operationally while requiring only minor modifications to an existing operational 3D-Var system.


2016 ◽  
Vol 73 (12) ◽  
pp. 4911-4925 ◽  
Author(s):  
Zhaoxia Pu ◽  
Shixuan Zhang ◽  
Mingjing Tong ◽  
Vijay Tallapragada

Abstract An initial vortex spindown, or strong adjustment to the structure and intensity of a hurricane’s initial vortex, presents a significant problem in hurricane forecasting, as with the NCEP Hurricane Weather Research and Forecasting Model (HWRF), because it can cause significantly degraded intensity forecasts. In this study, the influence of the self-consistent regional ensemble background error covariance on assimilating hurricane inner-core tail Doppler radar (TDR) observations in HWRF is examined with the NCEP gridpoint statistical interpolation (GSI)-based ensemble–three-dimensional variational (3DVAR) hybrid data assimilation system. It is found that the resolution of the background error covariance term, coming from the ensemble forecasts, has notable influence on the assimilation of hurricane inner-core observations and subsequent forecasting results. Specifically, the use of ensemble forecasting at high-resolution native grids results in significant reduction of the vortex spindown problem and thus leads to improved hurricane intensity forecasting. Further diagnoses are conducted to examine the spindown problem with a gradient wind balance. It is found that artificial vortex initialization, performed before data assimilation, can cause strong supergradient winds or imbalance in the vortex inner-core region. Assimilation of hurricane inner-core TDR data can significantly mitigate this imbalance by reducing the supergradient effects. Compared with the use of a global ensemble background error term, application of the self-consistent regional ensemble background covariance to inner-core data assimilation leads to better representation of the mesoscale hurricane inner-core structures. It can also result in more realistic vortex structures in data assimilation even when the observational data are unevenly distributed.


Sign in / Sign up

Export Citation Format

Share Document