Low-rank Kalman filtering for efficient state estimation of subsurface advective contaminant transport models
KAUST DepartmentEarth Science and Engineering Program
Applied Mathematics and Computational Science Program
Physical Sciences and Engineering (PSE) Division
Environmental Science and Engineering Program
Earth Sciences and Engineering Program
Earth Fluid Modeling and Prediction Group
Computational Transport Phenomena Lab
Permanent link to this recordhttp://hdl.handle.net/10754/562146
MetadataShow full item record
AbstractAccurate knowledge of the movement of contaminants in porous media is essential to track their trajectory and later extract them from the aquifer. A two-dimensional flow model is implemented and then applied on a linear contaminant transport model in the same porous medium. Because of different sources of uncertainties, this coupled model might not be able to accurately track the contaminant state. Incorporating observations through the process of data assimilation can guide the model toward the true trajectory of the system. The Kalman filter (KF), or its nonlinear invariants, can be used to tackle this problem. To overcome the prohibitive computational cost of the KF, the singular evolutive Kalman filter (SEKF) and the singular fixed Kalman filter (SFKF) are used, which are variants of the KF operating with low-rank covariance matrices. Experimental results suggest that under perfect and imperfect model setups, the low-rank filters can provide estimates as accurate as the full KF but at much lower computational effort. Low-rank filters are demonstrated to significantly reduce the computational effort of the KF to almost 3%. © 2012 American Society of Civil Engineers.
SponsorsThis publication utilized work supported in part by funds from the KAUST GCR Collaborative Fellow program.
Showing items related by title, author, creator and subject.
Reply to "comment on 'Ensemble Kalman filter with the unscented transform'"Luo, Xiaodong; Moroz, Irene M.; Hoteit, Ibrahim (Physica D: Nonlinear Phenomena, Elsevier BV, 2010-09) [Article]This is a reply to the comment of Dr. Sakov on the work "Ensemble Kalman filter with the unscented transform" of Luo and Moroz (2009) . © 2010 Elsevier B.V. All rights reserved.
Ensemble Kalman filtering with residual nudgingLuo, X.; Hoteit, Ibrahim (Tellus A, Co-Action Publishing, 2012-10-03) [Article]Covariance inflation and localisation are two important techniques that are used to improve the performance of the ensemble Kalman filter (EnKF) by (in effect) adjusting the sample covariances of the estimates in the state space. In this work, an additional auxiliary technique, called residual nudging, is proposed to monitor and, if necessary, adjust the residual norms of state estimates in the observation space. In an EnKF with residual nudging, if the residual norm of an analysis is larger than a pre-specified value, then the analysis is replaced by a new one whose residual norm is no larger than a pre-specified value. Otherwise, the analysis is considered as a reasonable estimate and no change is made. A rule for choosing the pre-specified value is suggested. Based on this rule, the corresponding new state estimates are explicitly derived in case of linear observations. Numerical experiments in the 40-dimensional Lorenz 96 model show that introducing residual nudging to an EnKF may improve its accuracy and/or enhance its stability against filter divergence, especially in the small ensemble scenario.
A Two-update Ensemble Kalman Filter for Land Hydrological Data Assimilation with an Uncertain ConstraintKhaki, M.; Ait-El-Fquih, Boujemaa; Hoteit, Ibrahim; Forootan, E.; Awange, J.; Kuhn, M. (Journal of Hydrology, Elsevier BV, 2017-10-25) [Article]Assimilating Gravity Recovery And Climate Experiment (GRACE) data into land hydrological models provides a valuable opportunity to improve the models’ forecasts and increases our knowledge of terrestrial water storages (TWS). The assimilation, however, may harm the consistency between hydrological water fluxes, namely precipitation, evaporation, discharge, and water storage changes. To address this issue, we propose a weak constrained ensemble Kalman filter (WCEnKF) that maintains estimated water budgets in balance with other water fluxes. Therefore, in this study, GRACE terrestrial water storages data are assimilated into the World-Wide Water Resources Assessment (W3RA) hydrological model over the Earth’s land areas covering 2002 – 2012. Multi-mission remotely sensed precipitation measurements from the Tropical Rainfall Measuring Mission (TRMM) and evaporation products from the Moderate Resolution Imaging Spectroradiometer (MODIS), as well as ground-based water discharge measurements are applied to close the water balance equation. The proposed WCEnKF contains two update steps; first, it incorporates observations from GRACE to improve model simulations of water storages, and second, uses the additional observations of precipitation, evaporation, and water discharge to establish the water budget closure. These steps are designed to account for error information associated with the included observation sets during the assimilation process. In order to evaluate the assimilation results, in addition to monitoring the water budget closure errors, in-situ groundwater measurements over the Mississippi River Basin in the US and the Murray-Darling Basin in Australia are used. Our results indicate approximately 24% improvement in the WCEnKF groundwater estimates over both basins compared to the use of (constraint-free) EnKF. WCEnKF also further reduces imbalance errors by approximately 82.53% (on average) and at the same time increases the correlations between the assimilation solutions and the water fluxes.