scholarly journals Research and application of parameter estimation method in hydrological model based on dual ensemble Kalman filter

2021 ◽  
Author(s):  
Mengtian Lu ◽  
Sicheng Lu ◽  
Weihong Liao ◽  
Xiaohui Lei ◽  
Zhaokai Yin ◽  
...  

Abstract Although field measurements and using long hydrological datasets provide a reliable method for parameters' calibration, changes in the underlying basin surface and lack of hydrometeorological data may affect parameter accuracy in streamflow simulation. The ensemble Kalman filter (EnKF) can be used as a real-time parameter correction method to solve this problem. In this study, five representative Xin'anjiang model parameters are selected to study the effects of the initial parameter ensemble distribution and the specific function form of the parameter on EnKF parameter estimation process for both single and multiple parameters. Results indicate: (1) the method of parameter calibration to determine the initial distribution mean can improve the assimilation efficiency; (2) there is mutual interference among the parameters during multiple parameters' estimation which invalidates some conclusions of single-parameter estimation. We applied and evaluated the EnKF method in Jinjiang River Basin, China. Compared to traditional approaches, our method showed a better performance in both basins with long hydrometeorological dataset (an increase of Kling–Gupta efficiency (KGE) from 0.810 to 0.887 and a decrease of bias from −1.08% to −0.74%); and in basins with a lack of hydrometeorological data (an increase of KGE from 0.536 to 0.849 and a decrease of bias from −15.55% to −11.42%).

2011 ◽  
Vol 15 (8) ◽  
pp. 2437-2457 ◽  
Author(s):  
S. Nie ◽  
J. Zhu ◽  
Y. Luo

Abstract. The performance of the ensemble Kalman filter (EnKF) in soil moisture assimilation applications is investigated in the context of simultaneous state-parameter estimation in the presence of uncertainties from model parameters, soil moisture initial condition and atmospheric forcing. A physically based land surface model is used for this purpose. Using a series of identical twin experiments in two kinds of initial parameter distribution (IPD) scenarios, the narrow IPD (NIPD) scenario and the wide IPD (WIPD) scenario, model-generated near surface soil moisture observations are assimilated to estimate soil moisture state and three hydraulic parameters (the saturated hydraulic conductivity, the saturated soil moisture suction and a soil texture empirical parameter) in the model. The estimation of single imperfect parameter is successful with the ensemble mean value of all three estimated parameters converging to their true values respectively in both NIPD and WIPD scenarios. Increasing the number of imperfect parameters leads to a decline in the estimation performance. A wide initial distribution of estimated parameters can produce improved simultaneous multi-parameter estimation performances compared to that of the NIPD scenario. However, when the number of estimated parameters increased to three, not all parameters were estimated successfully for both NIPD and WIPD scenarios. By introducing constraints between estimated hydraulic parameters, the performance of the constrained three-parameter estimation was successful, even if temporally sparse observations were available for assimilation. The constrained estimation method can reduce RMSE much more in soil moisture forecasting compared to the non-constrained estimation method and traditional non-parameter-estimation assimilation method. The benefit of this method in estimating all imperfect parameters simultaneously can be fully demonstrated when the corresponding non-constrained estimation method displays a relatively poor parameter estimation performance. Because all these constraints between parameters were obtained in a statistical sense, this constrained state-parameter estimation scheme is likely suitable for other land surface models even with more imperfect parameters estimated in soil moisture assimilation applications.


2005 ◽  
Vol 12 (3) ◽  
pp. 363-371 ◽  
Author(s):  
J. D. Annan ◽  
D. J. Lunt ◽  
J. C. Hargreaves ◽  
P. J. Valdes

Abstract. We demonstrate the application of an efficient multivariate probabilistic parameter estimation method to a spectral primitive equation atmospheric GCM. The method, which is based on the Ensemble Kalman Filter, is effective at tuning the surface air temperature climatology of the model to both identical twin data and reanalysis data. When 5 parameters were simultaneously tuned to fit the model to reanalysis data, the model errors were reduced by around 35% compared to those given by the default parameter values. However, the precipitation field proved to be insensitive to these parameters and remains rather poor. The model is computationally cheap but chaotic and otherwise realistic, and the success of these experiments suggests that this method should be capable of tuning more sophisticated models, in particular for the purposes of climate hindcasting and prediction. Furthermore, the method is shown to be useful in determining structural deficiencies in the model which can not be improved by tuning, and so can be a useful tool to guide model development. The work presented here is for a limited set of parameters and data, but the scalability of the method is such that it could easily be extended to a more comprehensive parameter set given sufficient observational data to constrain them.


Mathematics ◽  
2021 ◽  
Vol 9 (6) ◽  
pp. 636
Author(s):  
Rabih Ghostine ◽  
Mohamad Gharamti ◽  
Sally Hassrouny ◽  
Ibrahim Hoteit

In this paper, an extended SEIR model with a vaccination compartment is proposed to simulate the novel coronavirus disease (COVID-19) spread in Saudi Arabia. The model considers seven stages of infection: susceptible (S), exposed (E), infectious (I), quarantined (Q), recovered (R), deaths (D), and vaccinated (V). Initially, a mathematical analysis is carried out to illustrate the non-negativity, boundedness, epidemic equilibrium, existence, and uniqueness of the endemic equilibrium, and the basic reproduction number of the proposed model. Such numerical models can be, however, subject to various sources of uncertainties, due to an imperfect description of the biological processes governing the disease spread, which may strongly limit their forecasting skills. A data assimilation method, mainly, the ensemble Kalman filter (EnKF), is then used to constrain the model outputs and its parameters with available data. We conduct joint state-parameters estimation experiments assimilating daily data into the proposed model using the EnKF in order to enhance the model’s forecasting skills. Starting from the estimated set of model parameters, we then conduct short-term predictions in order to assess the predicability range of the model. We apply the proposed assimilation system on real data sets from Saudi Arabia. The numerical results demonstrate the capability of the proposed model in achieving accurate prediction of the epidemic development up to two-week time scales. Finally, we investigate the effect of vaccination on the spread of the pandemic.


2018 ◽  
Vol 140 (12) ◽  
Author(s):  
Kazuya Kusano ◽  
Hironobu Yamakawa ◽  
Kenich Hano

The feasibility of the parameter estimation on the basis of the ensemble Kalman filter (EnKF) for a practical simulation involving model errors was investigated. The three-dimensional flow and thermal simulations for the engine compartment of a test excavator were simulated, and several unknown temperatures used for boundary conditions were estimated with the method. The estimation method was validated in two steps. First, the estimation method was tested with the influence of the model errors removed by virtually creating true values with a simulation. These results showed that the proposed parameter-estimation method can successfully estimate surface temperatures. They also suggested that the appropriate ensemble size can be evaluated from the number of unknown parameters. Second, the estimation method was tested under a practical condition including model errors by using actual measurement data. Model errors were statistically estimated using prior obtained error data concerning other design configurations, and they were added to the observation error in the EnKF. These results showed that taking model errors into account in the EnKF provides more-accurate parameter-estimation results. Moreover, the uncertainty of an estimated parameter can be evaluated with the standard deviation of its distribution.


Author(s):  
Jan Vidar Grindheim ◽  
Inge Revhaug ◽  
Egil Pedersen

A finite difference method (FDM) solving the coupled partial differential equations governing three-dimensional (3D) motions of a towed underwater cable has been implemented in a combined ensemble Kalman filter (EnKF) and ensemble Kalman smoother (EnKS), as a new approach to combined state and parameter estimation for towed underwater cables. A simulation study of the method applied to a seismic streamer has been performed. Cable state variables as well as model parameters are estimated. Parameters estimated are crossline ocean current varying with time as well as cable tangential drag coefficient. The presented results indicate that the method is able to estimate state as well as parameters for seismic streamers.


2021 ◽  
Vol 11 (7) ◽  
pp. 2898
Author(s):  
Humberto C. Godinez ◽  
Esteban Rougier

Simulation of fracture initiation, propagation, and arrest is a problem of interest for many applications in the scientific community. There are a number of numerical methods used for this purpose, and among the most widely accepted is the combined finite-discrete element method (FDEM). To model fracture with FDEM, material behavior is described by specifying a combination of elastic properties, strengths (in the normal and tangential directions), and energy dissipated in failure modes I and II, which are modeled by incorporating a parameterized softening curve defining a post-peak stress-displacement relationship unique to each material. In this work, we implement a data assimilation method to estimate key model parameter values with the objective of improving the calibration processes for FDEM fracture simulations. Specifically, we implement the ensemble Kalman filter assimilation method to the Hybrid Optimization Software Suite (HOSS), a FDEM-based code which was developed for the simulation of fracture and fragmentation behavior. We present a set of assimilation experiments to match the numerical results obtained for a Split Hopkinson Pressure Bar (SHPB) model with experimental observations for granite. We achieved this by calibrating a subset of model parameters. The results show a steady convergence of the assimilated parameter values towards observed time/stress curves from the SHPB observations. In particular, both tensile and shear strengths seem to be converging faster than the other parameters considered.


Water ◽  
2019 ◽  
Vol 11 (7) ◽  
pp. 1520
Author(s):  
Zheng Jiang ◽  
Quanzhong Huang ◽  
Gendong Li ◽  
Guangyong Li

The parameters of water movement and solute transport models are essential for the accurate simulation of soil moisture and salinity, particularly for layered soils in field conditions. Parameter estimation can be achieved using the inverse modeling method. However, this type of method cannot fully consider the uncertainties of measurements, boundary conditions, and parameters, resulting in inaccurate estimations of parameters and predictions of state variables. The ensemble Kalman filter (EnKF) is well-suited to data assimilation and parameter prediction in Situations with large numbers of variables and uncertainties. Thus, in this study, the EnKF was used to estimate the parameters of water movement and solute transport in layered, variably saturated soils. Our results indicate that when used in conjunction with the HYDRUS-1D software (University of California Riverside, California, CA, USA) the EnKF effectively estimates parameters and predicts state variables for layered, variably saturated soils. The assimilation of factors such as the initial perturbation and ensemble size significantly affected in the simulated results. A proposed ensemble size range of 50–100 was used when applying the EnKF to the highly nonlinear hydrological models of the present study. Although the simulation results for moisture did not exhibit substantial improvement with the assimilation, the simulation of the salinity was significantly improved through the assimilation of the salinity and relative solutetransport parameters. Reducing the uncertainties in measured data can improve the goodness-of-fit in the application of the EnKF method. Sparse field condition observation data also benefited from the accurate measurement of state variables in the case of EnKF assimilation. However, the application of the EnKF algorithm for layered, variably saturated soils with hydrological models requires further study, because it is a challenging and highly nonlinear problem.


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.


2001 ◽  
Author(s):  
Jie Xiao ◽  
Bohdan T. Kulakowski

Abstract Vehicle dynamic models include parameters that qualify the dependence of input forces and moments on state and control variables. The accuracy of the model parameter estimates is important for modeling, simulation, and control. In general, the most accurate method for determining values of model parameters is by direct measurement. However, some parameters of vehicle dynamics, such as suspension damping or moments of inertia, are difficult to measure accurately. This study aims at establishing an efficient and accurate parameter estimation method for developing dynamic models for transit buses, such that this method can be easily implemented for simulation and control design purposes. Based on the analysis of robustness, as well as accuracy and efficiency of optimization techniques, a parameter estimation method that integrates Genetic Algorithms and the Maximum Likelihood Estimation is proposed. Choices of output signals and estimation criterion are discussed involving an extensive sensitivity analysis of the predicted output with respect to model parameters. Other experiment-related aspects, such as imperfection of data acquisition, are also considered. Finally, asymptotic Cramer-Rao lower bounds for the covariance of estimated parameters are obtained. Computer simulation results show that the proposed method is superior to gradient-based methods in accuracy, as well as robustness to the initial guesses and measurement uncertainty.


Sign in / Sign up

Export Citation Format

Share Document