scholarly journals A data assimilation method of the Ensemble Kalman Filter for use in severe dust storm forecasts over China

2007 ◽  
Vol 7 (6) ◽  
pp. 17511-17536
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 severe dust storm episodes occurred in March 2002 over China based on surface observations of dust concentrations to explore its impacts on forecast improvement. A series of sensitivity experiments using our system reveals that the EnKF is an advanced assimilation method to afford better initial conditions with surface observed PM10 in North China and lead to improved forecasts of dust storms, but forecast with large errors can be made by model errors. This result illustrates that it requires identifying and correcting 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.

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.


2016 ◽  
Vol 38 ◽  
pp. 190
Author(s):  
Regis Sperotto de Quadros ◽  
Fabrício Pereira Harter ◽  
Daniela Buske ◽  
Larri Silveira Pereira

Data Assimilation is a procedure to get the initial condition as accurately as possible, through the statistical combination of collected observations and a background field, usually a short-range forecast. In this research a complete assimilation system for the Lorenz equations based on Ensemble Kalman Filter is presented and examined. The Lorenz model is chosen for its simplicity in structure and the dynamic similarities with primitive equations models, such as modern numerical weather forecasting. Based on results, was concluded that, in this implementation, 10 members is the best setting, because there is an overfitting for ensembles with 50 and 100 members. It was also examined if the EnKF is effective to track the control for 20% and 40% of error in the initial conditions. The results show a disagreement between the “truth” and the estimation, especially in the end of integration period, due the chaotic nature of the system.  It was also concluded that EnKF have to be performed sufficiently frequently in order to produce desirable results.


SPE Journal ◽  
2011 ◽  
Vol 16 (02) ◽  
pp. 294-306 ◽  
Author(s):  
Lingzao Zeng ◽  
Haibin Chang ◽  
Dongxiao Zhang

Summary The ensemble Kalman filter (EnKF) has been used widely for data assimilation. Because the EnKF is a Monte Carlo-based method, a large ensemble size is required to reduce the sampling errors. In this study, a probabilistic collocation-based Kalman filter (PCKF) is developed to adjust the reservoir parameters to honor the production data. It combines the advantages of the EnKF for dynamic data assimilation and the polynomial chaos expansion (PCE) for efficient uncertainty quantification. In this approach, all the system parameters and states and the production data are approximated by the PCE. The PCE coefficients are solved with the probabilistic collocation method (PCM). Collocation realizations are constructed by choosing collocation point sets in the random space. The simulation for each collocation realization is solved forward in time independently by means of an existing deterministic solver, as in the EnKF method. In the analysis step, the needed covariance is approximated by the PCE coefficients. In this study, a square-root filter is employed to update the PCE coefficients. After the analysis, new collocation realizations are constructed. With the parameter collocation realizations as the inputs and the state collocation realizations as initial conditions, respectively, the simulations are forwarded to the next analysis step. Synthetic 2D water/oil examples are used to demonstrate the applicability of the PCKF in history matching. The results are compared with those from the EnKF on the basis of the same analysis. It is shown that the estimations provided by the PCKF are comparable to those obtained from the EnKF. The biggest improvement of the PCKF comes from the leading PCE approximation, with which the computational burden of the PCKF can be greatly reduced by means of a smaller number of simulation runs, and the PCKF outperforms the EnKF for a similar computational effort. When the correlation ratio is much smaller, the PCKF still provides estimations with a better accuracy for a small computational effort.


SPE Journal ◽  
2014 ◽  
Vol 20 (01) ◽  
pp. 202-221 ◽  
Author(s):  
Qinzhuo Liao ◽  
Dongxiao Zhang

Summary The ensemble Kalman filter (EnKF) has been widely used for data assimilation. It is challenging, however, when the relation of state and observation is strongly nonlinear. For example, near the flooding front in an immiscible flow, directly updating the saturation by use of the EnKF may lead to nonphysical results. One possible solution, which may be referred to as the restarted EnKF (REnKF), is to update the static state (e.g., permeability and porosity) and rerun the forward model from the initial time to obtain the updated dynamic state (e.g., pressure and saturation). However, it may become time-consuming, especially when the number of assimilation steps is large. In this study, we develop a transformed EnKF (TEnKF), in which the state is represented by displacement as an alternative variable. The displacement is first transformed from the forecasted state, then updated, and finally transformed back to obtain the updated state. Because the relation between displacement and observation is relatively linear, this new method provides a physically meaningful updated state without resolving the forward model. The TEnKF is tested in the history matching of multiphase flow in a 1D homogeneous medium, a 2D heterogeneous reservoir, and a 3D PUNQ-S3 model. The case studies show that the TEnKF produces physical results without the oscillation problem that occurs in the traditional EnKF, whereas the computational effort is reduced compared with the REnKF.


2011 ◽  
Vol 11 (24) ◽  
pp. 12901-12916 ◽  
Author(s):  
X. Tang ◽  
J. Zhu ◽  
Z. F. Wang ◽  
A. Gbaguidi

Abstract. In order to improve the surface ozone forecast over Beijing and surrounding regions, data assimilation method integrated into a high-resolution regional air quality model and a regional air quality monitoring network are employed. Several advanced data assimilation strategies based on ensemble Kalman filter are designed to adjust O3 initial conditions, NOx initial conditions and emissions, VOCs initial conditions and emissions separately or jointly through assimilating ozone observations. As a result, adjusting precursor initial conditions demonstrates potential improvement of the 1-h ozone forecast almost as great as shown by adjusting precursor emissions. Nevertheless, either adjusting precursor initial conditions or emissions show deficiency in improving the short-term ozone forecast at suburban areas. Adjusting ozone initial values brings significant improvement to the 1-h ozone forecast, and its limitations lie in the difficulty in improving the 1-h forecast at some urban site. A simultaneous adjustment of the above five variables is found to be able to reduce these limitations and display an overall better performance in improving both the 1-h and 24-h ozone forecast over these areas. The root mean square errors of 1-h ozone forecast at urban sites and suburban sites decrease by 51% and 58% respectively compared with those in free run. Through these experiments, we found that assimilating local ozone observations is determinant for ozone forecast over the observational area, while assimilating remote ozone observations could reduce the uncertainty in regional transport ozone.


2013 ◽  
Vol 141 (3) ◽  
pp. 900-917 ◽  
Author(s):  
Fuqing Zhang ◽  
Meng Zhang ◽  
Jonathan Poterjoy

Abstract This study examines the performance of a hybrid ensemble-variational data assimilation system (E3DVar) that couples an ensemble Kalman filter (EnKF) with the three-dimensional variational data assimilation (3DVar) system for the Weather Research and Forecasting (WRF) Model. The performance of E3DVar and the component EnKF and 3DVar systems are compared over the eastern United States for June 2003. Conventional sounding and surface observations as well as data from wind profilers, aircraft and ships, and cloud-tracked winds from satellites, are assimilated every 6 h during the experiments, and forecasts are verified using standard sounding observations. Forecasts with 12- to 72-h lead times are found to have noticeably smaller root-mean-square errors when initialized with the E3DVar system, as opposed to the EnKF, especially for the 12-h wind and moisture fields. The E3DVar system demonstrates similar performance as an EnKF, while using less than half the number of ensemble members, and is less sensitive to the use of a multiphysics ensemble to account for model errors. The E3DVar system is also compared with a similar hybrid method that replaces the 3DVar component with the WRF four-dimensional variational data assimilation (4DVar) method (denoted E4DVar). The E4DVar method demonstrated considerable improvements over E3DVar for nearly all model levels and variables at the shorter forecast lead times (12–48 h), but the forecast accuracies of all three ensemble-based methods (EnKF, E3DVar, and E4DVar) converge to similar results at longer lead times (60–72 h). Nevertheless, all methods that used ensemble information produced considerably better forecasts than the two methods that relied solely on static background error covariance (i.e., 3DVar and 4DVar).


2014 ◽  
Vol 142 (8) ◽  
pp. 2915-2934 ◽  
Author(s):  
Hailing Zhang ◽  
Zhaoxia Pu

Abstract A series of numerical experiments are conducted to examine the impact of surface observations on the prediction of landfalls of Hurricane Katrina (2005), one of the deadliest disasters in U.S. history. A specific initial time (0000 UTC 25 August 2005), which led to poor prediction of Hurricane Katrina in several previous studies, is selected to begin data assimilation experiments. Quick Scatterometer (QuikSCAT) ocean surface wind vectors and surface mesonet observations are assimilated with the minimum central sea level pressure and conventional observations from NCEP into an Advanced Research version of the Weather Research and Forecasting Model (WRF) using an ensemble Kalman filter method. Impacts of data assimilation on the analyses and forecasts of Katrina’s track, landfalling time and location, intensity, structure, and rainfall are evaluated. It is found that the assimilation of QuikSCAT and mesonet surface observations can improve prediction of the hurricane track and structure through modifying low-level thermal and dynamical fields such as wind, humidity, and temperature and enhancing low-level convergence and vorticity. However, assimilation of single-level surface observations alone does not ensure reasonable intensity forecasts because of the lack of constraint on the mid- to upper troposphere. When surface observations are assimilated with other conventional data, obvious enhancements are found in the forecasts of track and intensity, realistic convection, and surface wind structures. More importantly, surface data assimilation results in significant improvements in quantitative precipitation forecasts (QPFs) during landfalls.


2007 ◽  
Vol 135 (4) ◽  
pp. 1403-1423 ◽  
Author(s):  
Zhiyong Meng ◽  
Fuqing Zhang

Abstract In Part I of this two-part work, the feasibility of using an ensemble Kalman filter (EnKF) for mesoscale and regional-scale data assimilation through various observing system simulation experiments was demonstrated assuming a perfect forecast model for a winter snowstorm event that occurred on 24–26 January 2000. The current study seeks to explore the performance of the EnKF for the same event in the presence of significant model errors due to physical parameterizations by assimilating synthetic sounding and surface observations with typical temporal and spatial resolutions. The EnKF performance with imperfect models is also examined for a warm-season mesoscale convective vortex (MCV) event that occurred on 10–13 June 2003. The significance of model error in both warm- and cold-season events is demonstrated when the use of different cumulus parameterization schemes within different ensembles results in significantly different forecasts in terms of both ensemble mean and spread. Nevertheless, the EnKF performed reasonably well in most experiments with the imperfect model assumption (though its performance can sometimes be significantly degraded). As in Part I, where the perfect model assumption was utilized, most analysis error reduction comes from larger scales. Results show that using a combination of different physical parameterization schemes in the ensemble forecast can significantly improve filter performance. A multischeme ensemble has the potential to provide better background error covariance estimation and a smaller ensemble bias. There are noticeable differences in the performance of the EnKF for different flow regimes. In the imperfect scenarios considered, the improvement over the reference ensembles (pure ensemble forecasts without data assimilation) after 24 h of assimilation for the winter snowstorm event ranges from 36% to 67%. This is higher than the 26%–45% improvement noted after 36 h of assimilation for the warm-season MCV event. Scale- and flow-dependent error growth dynamics and predictability are possible causes for the differences in improvement. Compared to the power spectrum analyses for the snowstorm, it is found that forecast errors and ensemble spreads in the warm-season MCV event have relatively smaller power at larger scales and an overall smaller growth rate.


Sign in / Sign up

Export Citation Format

Share Document