共查询到20条相似文献，搜索用时 562 毫秒
1.
Data assimilation for distributed hydrological catchment modeling via ensemble Kalman filter 总被引：4，自引：0，他引：4
Catchment scale hydrological models are critical decision support tools for water resources management and environment remediation. However, the reliability of hydrological models is inevitably affected by limited measurements and imperfect models. Data assimilation techniques combine complementary information from measurements and models to enhance the model reliability and reduce predictive uncertainties. As a sequential data assimilation technique, the ensemble Kalman filter (EnKF) has been extensively studied in the earth sciences for assimilating in-situ measurements and remote sensing data. Although the EnKF has been demonstrated in land surface data assimilations, there are no systematic studies to investigate its performance in distributed modeling with high dimensional states and parameters. In this paper, we present an assessment on the EnKF with state augmentation for combined state-parameter estimation on the basis of a physical-based hydrological model, Soil and Water Assessment Tool (SWAT). Through synthetic simulation experiments, the capability of the EnKF is demonstrated by assimilating the runoff and other measurements, and its sensitivities are analyzed with respect to the error specification, the initial realization and the ensemble size. It is found that the EnKF provides an efficient approach for obtaining a set of acceptable model parameters and satisfactory runoff, soil water content and evapotranspiration estimations. The EnKF performance could be improved after augmenting with other complementary data, such as soil water content and evapotranspiration from remote sensing retrieval. Sensitivity studies demonstrate the importance of consistent error specification and the potential with small ensemble size in the data assimilation system. 相似文献
2.
Real-time assimilation of streamflow observations into a hydrological routing model: effects of model structures and updating methods 总被引：1，自引：1，他引：0
Maurizio Mazzoleni Haksu Lee Yuqiong Liu Dong-Jun Seo Alessandro Amaranto 《水文科学杂志》2018,63(3):386-407
This paper comparatively assesses the performance of five data assimilation techniques for three-parameter Muskingum routing with a spatially lumped or distributed model structure. The assimilation techniques used include direct insertion (DI), nudging scheme (NS), Kalman filter (KF), ensemble Kalman filter (EnKF) and asynchronous ensemble Kalman filter (AEnKF), which are applied to river reaches in Texas and Louisiana, USA. For both lumped and distributed routing, results from KF, EnKF and AEnKF are sensitive to the error specification. As expected, DI outperformed the other models in the case of lumped modelling, while in distributed routing, KF approaches, particularly AEnKF and EnKF, performed better than DI or nudging, reflecting the benefit of updating distributed states through error covariance modelling in KF approaches. The results of this work would be useful in setting up data assimilation systems that employ increasingly abundant real-time observations using distributed hydrological routing models. 相似文献
3.
Qiang Shu Mariush W. Kemblowski Mac McKee 《Stochastic Environmental Research and Risk Assessment (SERRA)》2005,19(5):361-374
Data assimilation method provides a framework to decrease the uncertainty of hydrological modeling by sequentially incorporating
observations into numerical model. Such a process involves estimating statistical moments of different order based on the
evolution of conditional probability distribution function. Because of the nonlinearity of many hydrological dynamics, explicit
and analytical solutions for moments of state distribution are often impossible. Evensen [J Geophys Res 99(c5): 10143–10162
(1994)] introduced Ensemble Kalman Filtering (EnKF) method to address such problems. We test and evaluate the performance
of EnKF in fusing model predictions and observations for a saturated–unsaturated integral-balance subsurface model. We find
EnKF improve the model predictions, and also we conclude a good estimate of state variance is essential for the success of
EnKF. 相似文献
4.
The performance of an inexpensive, ensemble-based optimal interpolation (EnOI) scheme that uses a stationary ensemble of model
anomalies to approximate forecast error covariances, is compared with that of an ensemble Kalman filter (EnKF). The model
to which the methods are applied is a pair of “perfect”, one-dimensional, linear advection equations for two related variables.
While EnOI is sub-optimal, it can give results that are comparable to those of the EnKF. The computational cost of EnOI is
typically about times less than that of EnKF, where is the ensemble size. We suggest that EnOI may provide a practical and cost-effective alternative to the EnKF for some applications
where computational cost is a limiting factor. We demonstrate that when the ensemble size is smaller than the dimension of
the model’s sub-space, both the EnKF and EnOI may require localisation around each observation to eliminate effects of sampling
error and to increase the effective number of independent ensemble members used to construct an analysis. However, localisation
can degrade an analysis if the length-scales of the localising function are too short. We demonstrate that, as the length-scale
of the localising function is decreased, localisation can significantly compromise the model’s dynamical balances. We also
find that localisation artificially amplifies high frequencies for applications of the EnKF. Based on our experiments, for
applications where localisation is necessary, the length-scales of the localisation should be larger than the decorrelation
length-scales of the variables being updated. 相似文献
5.
Getnet Y. Muluye 《水文科学杂志》2013,58(7):1118-1128
Abstract There is a continuing effort to advance the skill of long-range hydrological forecasts to support water resources decision making. The present study investigates the potential of an extended Kalman filter approach to perform supervised training of a recurrent multilayer perceptron (RMLP) to forecast up to 12-month-ahead lake water levels and streamflows in Canada. The performance of the RMLP was compared with the conventional multilayer perceptron (MLP) using suites of diagnostic measures. The results of the forecasting experiment showed that the RMLP model was able to provide a robust modelling framework capable of describing complex dynamics of the hydrological processes, thereby yielding more accurate and realistic forecasts than the MLP model. The performance of the method in the present study is very promising; however, further investigation is required to ascertain the versatility of the approach in characterizing different water resources and environmental problems. Citation Muluye, G. Y. (2011) Improving long-range hydrological forecasts with extended Kalman filters. Hydrol. Sci. J. 56(7), 1118–1128. 相似文献
6.
In this work the ensemble Kalman filter (EnKF) is applied to investigate the flow and transport processes at the macro-dispersion experiment (MADE) site in Columbus, MS. The EnKF is a sequential data assimilation approach that adjusts the unknown model parameter values based on the observed data with time. The classic advection–dispersion (AD) and the dual-domain mass transfer (DDMT) models are employed to analyze the tritium plume during the second MADE tracer experiment. The hydraulic conductivity (K), longitudinal dispersivity in the AD model, and mass transfer rate coefficient and mobile porosity ratio in the DDMT model, are estimated in this investigation. Because of its sequential feature, the EnKF allows for the temporal scaling of transport parameters during the tritium concentration analysis. Inverse simulation results indicate that for the AD model to reproduce the extensive spatial spreading of the tritium observed in the field, the K in the downgradient area needs to be increased significantly. The estimated K in the AD model becomes an order of magnitude higher than the in situ flowmeter measurements over a large portion of media. On the other hand, the DDMT model gives an estimation of K that is much more comparable with the flowmeter values. In addition, the simulated concentrations by the DDMT model show a better agreement with the observed values. The root mean square (RMS) between the observed and simulated tritium plumes is 0.77 for the AD model and 0.45 for the DDMT model at 328 days. Unlike the AD model, which gives inconsistent K estimates at different times, the DDMT model is able to invert the K values that consistently reproduce the observed tritium concentrations through all times. 相似文献
7.
8.
Hydrological model and observation errors are often non-Gaussian and/or biased, and the statistical properties of the errors are often unknown or not fully known. Thus, determining the true error covariance matrices is a challenge for data assimilation approaches such as the most widely used Kalman filter (KF) and its extensions, which assume Gaussian error nature and need fully known error statistics. This paper introduces H-infinite filter (HF) to hydrological modeling and compares HF with KF under various model and observation error conditions. HF is basically a robust version of KF. When model performance is not well known, or changes unpredictably, HF may be preferred over KF. HF is especially suitable for the cases where the estimation performance in the worst error case needs to be guaranteed. Through the application of HF to a hypothetical hydrologic model, this paper shows that HF is less sensitive to the uncertainty in the initial condition, corrects system bias more effectively, and converges to true state faster after interruptions than KF. In particular, HF performs better in dealing with instant human inputs (irrigation is used as an example), which are characterized by non-stationary, non-Gaussian and not fully known errors. However HF design can be more difficult than KF design due to the sensitivity of HF performance to design parameters (weights for model and observation error terms). Through sensitivity analysis, this paper shows the existence of a certain range of those parameters, in which the “best” value of the parameters is located. The tuning of HF design parameters, which can be based on users’ prior knowledge on the nature of model and observation errors, is critical for the implementation of HF. 相似文献
9.
An approach to handling non-Gaussianity of parameters and state variables in ensemble Kalman filtering 总被引：4，自引：0，他引：4
The ensemble Kalman filter (EnKF) is a commonly used real-time data assimilation algorithm in various disciplines. Here, the EnKF is applied, in a hydrogeological context, to condition log-conductivity realizations on log-conductivity and transient piezometric head data. In this case, the state vector is made up of log-conductivities and piezometric heads over a discretized aquifer domain, the forecast model is a groundwater flow numerical model, and the transient piezometric head data are sequentially assimilated to update the state vector. It is well known that all Kalman filters perform optimally for linear forecast models and a multiGaussian-distributed state vector. Of the different Kalman filters, the EnKF provides a robust solution to address non-linearities; however, it does not handle well non-Gaussian state-vector distributions. In the standard EnKF, as time passes and more state observations are assimilated, the distributions become closer to Gaussian, even if the initial ones are clearly non-Gaussian. A new method is proposed that transforms the original state vector into a new vector that is univariate Gaussian at all times. Back transforming the vector after the filtering ensures that the initial non-Gaussian univariate distributions of the state-vector components are preserved throughout. The proposed method is based in normal-score transforming each variable for all locations and all time steps. This new method, termed the normal-score ensemble Kalman filter (NS-EnKF), is demonstrated in a synthetic bimodal aquifer resembling a fluvial deposit, and it is compared to the standard EnKF. The proposed method performs better than the standard EnKF in all aspects analyzed (log-conductivity characterization and flow and transport predictions). 相似文献
10.
Data assimilation techniques have been proven as an effective tool to improve model forecasts by combining information about observed variables in many areas. This article examines the potential of assimilating surface soil moisture observations into a field‐scale hydrological model, the Root Zone Water Quality Model, to improve soil moisture estimation. The Ensemble Kalman Filter (EnKF), a popular data assimilation technique for nonlinear systems, was applied and compared with a simple direct insertion method. In situ soil moisture data at four different depths (5, 20, 40, and 60 cm) from two agricultural fields (AS1 and AS2) in northeastern Indiana were used for assimilation and validation purposes. Through daily update, the EnKF improved soil moisture estimation compared with the direct insertion method and model results without assimilation, having more distinct improvement at the 5 and 20 cm depths than for deeper layers (40 and 60 cm). Local vertical soil property heterogeneity in AS1 deteriorated soil moisture estimates with the EnKF. Removal of systematic bias in the forecast model was found to be critical for more successful soil moisture data assimilation studies. This study also demonstrates that a more frequent update generally contributes in enhancing the open loop simulation; however, large forecasting error can prevent more frequent update from providing better results. In addition, results indicate that various ensemble sizes make little difference in the assimilation results. An ensemble of 100 members produced results that were comparable with results obtained from larger ensembles. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
11.
Evensen (2003) presents a modification of the Ensemble Kalman Filter (EnKF), in which the observation-error and background-error covariance matrices are both represented by ensembles, in contrast to the usual practice, where only the background error is so represented. It is shown that this modification can cause the ensemble to collapse to a single member, in the common situation where the number of observations is more than twice the number of ensemble members, and to be rank-deficient when the number of observations is greater than or equal to the ensemble size. It is also shown that some further modifications to the scheme, presented by Evensen as offering numerical efficiencies, can prevent this collapse. However, these latter modifications are shown in some simple numerical examples to require tuning to produce acceptable results, which are nevertheless inferior to those of the standard EnKF.Acknowledgements The author acknowledges useful discussions with Peter Steinle, and other participants at the EnKF workshop held in BMRC in November, 2003. 相似文献
12.
Data assimilation in groundwater modelling: ensemble Kalman filter versus ensemble smoothers 
下载免费PDF全文
下载免费PDF全文 Groundwater modelling calls for an effective and robust data integrating method to fill the gap between the model and observation data. The ensemble Kalman filter (EnKF), a real‐time data assimilation method, has been increasingly applied in multiple disciplines such as petroleum engineering and hydrogeology. In this approach, a groundwater model is updated sequentially with measured data such as hydraulic head and concentration. As an alternative to the EnKF, the ensemble smoother (ES) has been proposed for updating groundwater models using all the data together, with much less computational cost. To further improve the performance of the ES, an iterative ES has been proposed for continuously updating the model by assimilating measurements together. In this work, we compare the performance of the EnKF, the ES, and the iterative ES using a synthetic example in groundwater modelling. Hydraulic head data modelled on the basis of the reference conductivity field are used to inversely estimate conductivities at unsampled locations. Results are evaluated in terms of the characterization of conductivity and groundwater flow predictions. It is concluded that (a) the iterative ES works better than the standard ES because of its continuous updating and (b) the iterative ES could achieve results comparable with those of the EnKF, with less computational cost. These findings show that the iterative ES should be paid much more attention for data assimilation in groundwater modelling. 相似文献
13.
The purpose of this paper is to provide a comprehensive presentation and interpretation of the Ensemble Kalman Filter (EnKF) and its numerical implementation. The EnKF has a large user group, and numerous publications have discussed applications and theoretical aspects of it. This paper reviews the important results from these studies and also presents new ideas and alternative interpretations which further explain the success of the EnKF. In addition to providing the theoretical framework needed for using the EnKF, there is also a focus on the algorithmic formulation and optimal numerical implementation. A program listing is given for some of the key subroutines. The paper also touches upon specific issues such as the use of nonlinear measurements, in situ profiles of temperature and salinity, and data which are available with high frequency in time. An ensemble based optimal interpolation (EnOI) scheme is presented as a cost-effective approach which may serve as an alternative to the EnKF in some applications. A fairly extensive discussion is devoted to the use of time correlated model errors and the estimation of model bias.Responsible Editor: Jörg-Olaf Wolff 相似文献
14.
15.
This paper presents a coupling of an ensemble Kalman filter (EnKF) with a discontinuous Galerkin-based, two-dimensional circulation model (DG ADCIRC-2DDI) to improve the state estimation of tidal hydrodynamics including water surface elevations and depth-integrated velocities. The methodology in this paper using EnKF perturbs the modeled hydrodynamics and bottom friction parameterization in the model while assimilating data with inherent error, and demonstrates a capability to apply EnKF within DG ADCIRC-2DDI for data assimilation. Parallel code development presents a unique aspect of the approach taken and is briefly described in the paper, followed by an application to a real estuarine system, the lower St. Johns River in north Florida, for the state estimation of tidal hydrodynamics. To test the value of gauge observations for improving state estimation, a tide modeling case study is performed for the lower St. Johns River successively using one of the four available tide gauging stations in model-data comparison. The results are improved simulations of water surface elevations and depth-integrated velocities using DG ADCIRC-2DDI with EnKF, both locally where data are available and non-locally where data are not available. The methodology, in general, is extensible to other modeling and data applications, for example, the use of remote sensing data, and specifically, can be readily applied as is to study other tidal systems. 相似文献
16.
Martyn P. Clark David E. Rupp Ross A. Woods Xiaogu Zheng Richard P. Ibbitt Andrew G. Slater Jochen Schmidt Michael J. Uddstrom 《Advances in water resources》2008
This paper describes an application of the ensemble Kalman filter (EnKF) in which streamflow observations are used to update states in a distributed hydrological model. We demonstrate that the standard implementation of the EnKF is inappropriate because of non-linear relationships between model states and observations. Transforming streamflow into log space before computing error covariances improves filter performance. We also demonstrate that model simulations improve when we use a variant of the EnKF that does not require perturbed observations. Our attempt to propagate information to neighbouring basins was unsuccessful, largely due to inadequacies in modelling the spatial variability of hydrological processes. New methods are needed to produce ensemble simulations that both reflect total model error and adequately simulate the spatial variability of hydrological states and fluxes. 相似文献
17.
Ensemble Kalman filter, EnKF, as a Monte Carlo sequential data assimilation method has emerged promisingly for subsurface media characterization during past decade. Due to high computational cost of large ensemble size, EnKF is limited to small ensemble set in practice. This results in appearance of spurious correlation in covariance structure leading to incorrect or probable divergence of updated realizations. In this paper, a universal/adaptive thresholding method is presented to remove and/or mitigate spurious correlation problem in the forecast covariance matrix. This method is, then, extended to regularize Kalman gain directly. Four different thresholding functions have been considered to threshold forecast covariance and gain matrices. These include hard, soft, lasso and Smoothly Clipped Absolute Deviation (SCAD) functions. Three benchmarks are used to evaluate the performances of these methods. These benchmarks include a small 1D linear model and two 2D water flooding (in petroleum reservoirs) cases whose levels of heterogeneity/nonlinearity are different. It should be noted that beside the adaptive thresholding, the standard distance dependant localization and bootstrap Kalman gain are also implemented for comparison purposes. We assessed each setup with different ensemble sets to investigate the sensitivity of each method on ensemble size. The results indicate that thresholding of forecast covariance yields more reliable performance than Kalman gain. Among thresholding function, SCAD is more robust for both covariance and gain estimation. Our analyses emphasize that not all assimilation cycles do require thresholding and it should be performed wisely during the early assimilation cycles. The proposed scheme of adaptive thresholding outperforms other methods for subsurface characterization of underlying benchmarks. 相似文献
18.
C. LOUMAGNE M. NORMAND M. RIFFARD A. WEISSE A. QUESNEY S. LE HÉGARAT-MASCLE 《水文科学杂志》2013,58(1):89-102
Abstract The purpose of this paper is to present the methodology set up to derive catchment soil moisture from Earth Observation (EO) data using microwave spaceborne Synthetic Aperture Radar (SAR) images from ERS satellites and to study the improvements brought about by an assimilation of this information into hydrological models. The methodology used to derive EO data is based on the appropriate selection of land cover types for which the radar signal is mainly sensitive to soil moisture variations. Then a hydrological model is chosen, which can take advantage of the new information brought by remote sensing. The assimilation of soil moisture deduced from EO data into hydrological models is based principally on model parameter updating. The main assumption of this method is that the better the model simulates the current hydrological system, the better the following forecast will be. Another methodology used is a sequential one based on Kalman filtering. These methods have been put forward for use in the European AIMWATER project on the Seine catchment upstream of Paris (France) where dams are operated to alleviate floods in the Paris area. 相似文献
19.
The objective of the study is to evaluate the potential of a data assimilation system for real-time flash flood forecasting over small watersheds by updating model states. To this end, the Ensemble Square-Root-Filter (EnSRF) based on the Ensemble Kalman Filter (EnKF) technique was coupled to a widely used conceptual rainfall-runoff model called HyMOD. Two small watersheds susceptible to flash flooding from America and China were selected in this study. The modeling and observational errors were considered in the framework of data assimilation, followed by an ensemble size sensitivity experiment. Once the appropriate model error and ensemble size was determined, a simulation study focused on the performance of a data assimilation system, based on the correlation between streamflow observation and model states, was conducted. The EnSRF method was implemented within HyMOD and results for flash flood forecasting were analyzed, where the calibrated streamflow simulation without state updating was treated as the benchmark or nature run. Results for twenty-four flash-flood events in total from the two watersheds indicated that the data assimilation approach effectively improved the predictions of peak flows and the hydrographs in general. This study demonstrated the benefit and efficiency of implementing data assimilation into a hydrological model to improve flash flood forecasting over small, instrumented basins with potential application to real-time alert systems. 相似文献
20.
本文采用高效集合卡尔曼滤波(EnKF)算法和背景场热层电离层理论模式NCAR-TIEGCM,开发了热层电离层数据同化系统.基于全球空地基GNSS电离层斜TEC观测、CHAMP和TIMED/GUVI热层参量观测构型设计了系列观测系统模拟实验,对热层参量进行估计.实验结果表明,(1)通过集合卡尔曼滤波算法同化电离层TEC观测能够较好地优化热层参量.(2)中性质量密度优化效果在整个同化阶段均有提升,提升百分比能达到40%.(3)积分氧氮比在同化阶段也能得到较好的优化,但在电子密度水平梯度变化剧烈区域效果较差.最后本文对中性质量密度进行了预报评估,结果表明,由于中性成分优化,在地磁平静条件下其预报时间尺度可长达24h. 相似文献