scholarly journals Multi-model ensemble hydrologic prediction and uncertainties analysis

2014 ◽  
Vol 364 ◽  
pp. 249-254 ◽  
Author(s):  
S. Jiang ◽  
L. Ren ◽  
X. Yang ◽  
M. Ma ◽  
Y. Liu
Keyword(s):  
Parameter Optimization ◽  
Optimization Algorithms ◽  
Model Parameters ◽  
Model Ensemble ◽  
Streamflow Prediction ◽  
Parameter Uncertainties ◽  
Input Uncertainty ◽  
Daily Streamflow ◽  
Shuffled Complex Evolution ◽  
Modelling Uncertainties

Abstract. Modelling uncertainties (i.e. input errors, parameter uncertainties and model structural errors) inevitably exist in hydrological prediction. A lot of recent attention has focused on these, of which input error modelling, parameter optimization and multi-model ensemble strategies are the three most popular methods to demonstrate the impacts of modelling uncertainties. In this paper the Xinanjiang model, the Hybrid rainfall–runoff model and the HYMOD model were applied to the Mishui Basin, south China, for daily streamflow ensemble simulation and uncertainty analysis. The three models were first calibrated by two parameter optimization algorithms, namely, the Shuffled Complex Evolution method (SCE-UA) and the Shuffled Complex Evolution Metropolis method (SCEM-UA); next, the input uncertainty was accounted for by introducing a normally-distributed error multiplier; then, the simulation sets calculated from the three models were combined by Bayesian model averaging (BMA). The results show that both these parameter optimization algorithms generate good streamflow simulations; specifically the SCEM-UA can imply parameter uncertainty and give the posterior distribution of the parameters. Considering the precipitation input uncertainty, the streamflow simulation precision does not improve very much. While the BMA combination not only improves the streamflow prediction precision, it also gives quantitative uncertainty bounds for the simulation sets. The SCEM-UA calculated prediction interval is better than the SCE-UA calculated one. These results suggest that considering the model parameters' uncertainties and doing multi-model ensemble simulations are very practical for streamflow prediction and flood forecasting, from which more precision prediction and more reliable uncertainty bounds can be generated.

scholarly journals Bayesian estimation of parameters in a regional hydrological model

2002 ◽  
Vol 6 (5) ◽  
pp. 883-898 ◽  
Author(s):  
K. Engeland ◽  
L. Gottschalk
Keyword(s):  
Monte Carlo ◽  
Markov Chain ◽  
Markov Chain Monte Carlo ◽  
Hydrological Model ◽  
Parameter Estimates ◽  
Model Parameters ◽  
Parameter Uncertainties ◽  
Daily Streamflow ◽  
Likelihood Model ◽  
Simulation Errors

Abstract. This study evaluates the applicability of the distributed, process-oriented Ecomag model for prediction of daily streamflow in ungauged basins. The Ecomag model is applied as a regional model to nine catchments in the NOPEX area, using Bayesian statistics to estimate the posterior distribution of the model parameters conditioned on the observed streamflow. The distribution is calculated by Markov Chain Monte Carlo (MCMC) analysis. The Bayesian method requires formulation of a likelihood function for the parameters and three alternative formulations are used. The first is a subjectively chosen objective function that describes the goodness of fit between the simulated and observed streamflow, as defined in the GLUE framework. The second and third formulations are more statistically correct likelihood models that describe the simulation errors. The full statistical likelihood model describes the simulation errors as an AR(1) process, whereas the simple model excludes the auto-regressive part. The statistical parameters depend on the catchments and the hydrological processes and the statistical and the hydrological parameters are estimated simultaneously. The results show that the simple likelihood model gives the most robust parameter estimates. The simulation error may be explained to a large extent by the catchment characteristics and climatic conditions, so it is possible to transfer knowledge about them to ungauged catchments. The statistical models for the simulation errors indicate that structural errors in the model are more important than parameter uncertainties. Keywords: regional hydrological model, model uncertainty, Bayesian analysis, Markov Chain Monte Carlo analysis

scholarly journals Quantifying multi-source uncertainties in multi-model predictions using the Bayesian model averaging scheme

2017 ◽  
Vol 49 (3) ◽  
pp. 954-970 ◽  
Author(s):  
Shanhu Jiang ◽  
Liliang Ren ◽  
Chong-Yu Xu ◽  
Shuya Liu ◽  
Fei Yuan ◽  
...  
Keyword(s):  
Uncertainty Analysis ◽  
Bayesian Model ◽  
Bayesian Model Averaging ◽  
Prediction Interval ◽  
Model Averaging ◽  
Streamflow Prediction ◽  
Input Uncertainty ◽  
Shuffled Complex Evolution ◽  
Model Predictions ◽  
Combination Approach

Abstract This study focuses on a quantitative multi-source uncertainty analysis of multi-model predictions. Three widely used hydrological models, i.e., Xinanjiang (XAJ), hybrid rainfall–runoff (HYB), and HYMOD (HYM), were calibrated by two parameter optimization algorithms, namely, shuffled complex evolution (SCE-UA) method and shuffled complex evolution metropolis (SCEM-UA) method on the Mishui basin, south China. The input uncertainty was quantified by utilizing a normally distributed error multiplier. The ensemble simulation sets calculated from the three models were combined using the Bayesian model averaging (BMA) method. Results indicate the following. (1) Both SCE-UA and SCEM-UA resulted in good and comparable streamflow simulations. Specifically, the SCEM-UA implied parameter uncertainty and provided the posterior distribution of the parameters. (2) In terms of the precipitation input uncertainty, precision of streamflow simulations did not improve remarkably. (3) The BMA combination not only improved the precision of streamflow prediction, but also quantified the uncertainty bounds of the simulation. (4) The prediction interval calculated using the SCEM-UA-based BMA combination approach appears superior to that calculated using the SCE-UA-based BMA combination for both high flows and low flows. Results suggest that the comprehensive uncertainty analysis by using the SCEM-UA algorithm and BMA method is superior for streamflow predictions and flood forecasting.

scholarly journals A Review of Tank Model and Its Applicability to Various Korean Catchment Conditions

Water ◽  
2020 ◽  
Vol 12 (12) ◽  
pp. 3588
Author(s):  
Jong Wook Lee ◽  
Sun Dong Chegal ◽  
Seung Oh Lee
Keyword(s):  
Water Resource ◽  
Optimization Algorithms ◽  
Flow Characteristics ◽  
Climatic Conditions ◽  
Model Parameters ◽  
Transfer Model ◽  
Ungauged Catchments ◽  
Modelling Framework ◽  
Shuffled Complex Evolution ◽  
Resource Modelling

This paper reviews a conceptual rainfall-runoff model called Tank which has been widely used over the last 20 years in Korea as a part of a water resource modelling framework for assessing and developing long-term water resource polices. In order to examine the uncertainty of model predictions and the sensitivity of model’s parameters, Monte Carlos and Markov chain-based approaches are applied to five catchments of various Korean geographical and climatic conditions where the catchment sizes are ranged from 83 to 4786 km2. In addition, three optimization algorithms—dynamically dimensioned search (DDS), robust parameter estimation (ROPE), and shuffled complex evolution (SCE)—are selected to test whether the model parameters can be optimized consistently within a narrower range than the uncertainty bounds. From the uncertainty analysis, it is found that there is limited success in refining the priori distributions of the model parameters, indicating there is a high degree of equifinality for some parameters or at least there are large numbers of parameter combinations leading to good solutions within model’s uncertainty bounds. Out of the three optimization algorithms, SCE meets the criteria of the consistency best. It is also found that there are still some parameters that even the SCE method struggles to refine the priori distributions. It means that their contribution to model results is minimal and can take a value within a reasonable range. It suggests that the model may be reconceptualized to be parsimonious and to rationalize some parameters without affecting model’s capacity to replicate historical flow characteristics. Cross-validation indicates that sensitive parameters to catchment characteristics can be transferred when geophysical similarity exists between two catchments. Regionalization can be further improved by using a regression or geophysical similarity-based approach to transfer model parameters to ungauged catchments. It may be beneficial to categorize the model parameters depending on the level of their sensitivities, and a different approach to each category may be applied to regionalize the calibrated parameters.

scholarly journals Parameter Optimization for Real-World ENSO Forecast in an Intermediate Coupled Model

2019 ◽  
Vol 147 (5) ◽  
pp. 1429-1445 ◽  
Author(s):  
Yuchu Zhao ◽  
Zhengyu Liu ◽  
Fei Zheng ◽  
Yishuai Jin
Keyword(s):  
Data Assimilation ◽  
Parameter Optimization ◽  
Real World ◽  
Climate Models ◽  
Coupled Model ◽  
Prediction Skill ◽  
Model Parameters ◽  
Real World Data ◽  
The Real ◽  
Advection Term

Abstract We performed parameter estimation in the Zebiak–Cane model for the real-world scenario using the approach of ensemble Kalman filter (EnKF) data assimilation and the observational data of sea surface temperature and wind stress analyses. With real-world data assimilation in the coupled model, our study shows that model parameters converge toward stable values. Furthermore, the new parameters improve the real-world ENSO prediction skill, with the skill improved most by the parameter of the highest climate sensitivity (gam2), which controls the strength of anomalous upwelling advection term in the SST equation. The improved prediction skill is found to be contributed mainly by the improvement in the model dynamics, and second by the improvement in the initial field. Finally, geographic-dependent parameter optimization further improves the prediction skill across all the regions. Our study suggests that parameter optimization using ensemble data assimilation may provide an effective strategy to improve climate models and their real-world climate predictions in the future.

scholarly journals Regionalizing rainfall–runoff model parameters to predict the daily streamflow of ungauged catchments in the dry tropics

2009 ◽  
Vol 40 (5) ◽  
pp. 433-444 ◽  
Author(s):  
David A. Post
Keyword(s):  
Input Data ◽  
Water Yield ◽  
Model Parameters ◽  
Rainfall Runoff ◽  
Wet Season ◽  
Daily Streamflow ◽  
Total Length ◽  
Dry Tropics ◽  
Rainfall Runoff Model ◽  
Runoff Model

A methodology has been derived which allows an estimate to be made of the daily streamflow at any point within the Burdekin catchment in the dry tropics of Australia. The input data requirements are daily rainfall (to drive the rainfall–runoff model) and mean average wet season rainfall, total length of streams, percent cropping and percent forest in the catchment (to regionalize the parameters of the rainfall–runoff model). The method is based on the use of a simple, lumped parameter rainfall–runoff model, IHACRES (Identification of unit Hydrographs And Component flows from Rainfall, Evaporation and Streamflow data). Of the five parameters in the model, three have been set to constants to reflect regional conditions while the other two have been related to physio-climatic attributes of the catchment under consideration. The parameter defining total catchment water yield (c) has been estimated based on the mean average wet season rainfall, while the streamflow recession time constant (τ) has been estimated based on the total length of streams, percent cropping and percent forest in the catchment. These relationships have been shown to be applicable over a range of scales from 68–130,146 km2. However, three separate relationships were required to define c in the three major physiographic regions of the Burdekin: the upper Burdekin, Bowen and Suttor/lower Burdekin. The invariance of the relationships with scale indicates that the dominant processes may be similar across a range of scales. The fact that different relationships were required for each of the three major regions indicates the geographic limitations of this regionalization approach. For most of the 24 gauged catchments within the Burdekin the regionalized rainfall–runoff models were nearly as good as or better than the rainfall–runoff models calibrated to the observed streamflow. In addition, models often performed better over the simulation period than the calibration period. This indicates that future improvements in regionalization should focus on improving the quality of input data and rainfall–runoff model conceptualization rather than on the regionalization procedure per se.

scholarly journals Calibration of Conceptual Rainfall-Runoff Models Using Global Optimization

2015 ◽  
Vol 2015 ◽  
pp. 1-12 ◽  
Author(s):  
Chao Zhang ◽  
Ru-bin Wang ◽  
Qing-xiang Meng
Keyword(s):  
Global Optimization ◽  
Parameter Optimization ◽  
Hydrological Model ◽  
Optimization Method ◽  
Difficult Problem ◽  
Rainfall Runoff ◽  
True Parameter ◽  
Shuffled Complex Evolution ◽  
Structure Lead ◽  
Parameter Values

Parameter optimization for the conceptual rainfall-runoff (CRR) model has always been the difficult problem in hydrology since watershed hydrological model is high-dimensional and nonlinear with multimodal and nonconvex response surface and its parameters are obviously related and complementary. In the research presented here, the shuffled complex evolution (SCE-UA) global optimization method was used to calibrate the Xinanjiang (XAJ) model. We defined the ideal data and applied the method to observed data. Our results show that, in the case of ideal data, the data length did not affect the parameter optimization for the hydrological model. If the objective function was selected appropriately, the proposed method found the true parameter values. In the case of observed data, we applied the technique to different lengths of data (1, 2, and 3 years) and compared the results with ideal data. We found that errors in the data and model structure lead to significant uncertainties in the parameter optimization.

scholarly journals Modeling the impact of drought on canopy carbon and water fluxes through parameter optimization using an ensemble Kalman filter

2009 ◽  
Vol 6 (4) ◽  
pp. 8279-8309 ◽  
Author(s):  
W. Ju ◽  
S. Wang ◽  
G. Yu ◽  
Y. Zhou ◽  
H. Wang
Keyword(s):  
Kalman Filter ◽  
Parameter Optimization ◽  
Ensemble Kalman Filter ◽  
Terrestrial Ecosystems ◽  
Model Parameters ◽  
Interannual Variations ◽  
Atmospheric Water ◽  
Model Parameterization ◽  
Coniferous Plantation ◽  
The Impact

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.

scholarly journals Linking biogeochemistry to hydro-geometrical variability in tidal estuaries: a generic modeling approach

2015 ◽  
Vol 12 (7) ◽  
pp. 6351-6435
Author(s):  
C. Volta ◽  
G. G. Laruelle ◽  
S. Arndt ◽  
P. Regnier
Keyword(s):  
Sensitivity Analysis ◽  
World Ocean ◽  
Reaction Network ◽  
Co2 Exchange ◽  
Model Parameters ◽  
Parameter Uncertainties ◽  
Tidal Estuaries ◽  
Biogeochemical Indicators ◽  
C And N ◽  
Co2 Sink

Abstract. This study applies the Carbon-Generic Estuary Model (C-GEM) modeling platform to simulate the estuarine biogeochemical dynamics – in particular the air-water CO2 exchange – in three idealized end-member systems covering the main features of tidal alluvial estuaries. C-GEM uses a generic biogeochemical reaction network and a unique set of model parameters extracted from a comprehensive literature survey to perform steady-state simulations representing average conditions for temperate estuaries worldwide. Climate and boundary conditions are extracted from published global databases (e.g. World Ocean Atlas, GLORICH) and catchment model outputs (GlobalNEWS2). The whole-system biogeochemical indicators Net Ecosystem Metabolism (NEM), C and N filtering capacities (FCTC and FCTN, respectively) and CO2 gas exchanges (FCO2) are calculated across the three end-member systems and are related to their main hydrodynamic and transport characteristics. A sensitivity analysis, which propagates the parameter uncertainties, is also carried out, followed by projections of changes in the biogeochemical indicators for the year 2050. Results show that the average C filtering capacities for baseline conditions are 40, 30 and 22% for the marine, mixed and riverine estuary, respectively. This translates into a first-order, global CO2 outgassing flux for tidal estuaries between 0.04 and 0.07 Pg C yr−1. N filtering capacities, calculated in similar fashion, range from 22% for the marine estuary to 18 and 15% for the mixed and the riverine estuary, respectively. Sensitivity analysis performed by varying the rate constants for aerobic degradation, denitrification and nitrification over the range of values reported in the literature significantly widens these ranges for both C and N. Simulations for the year 2050 indicate that all end-member estuaries will remain net heterotrophic and while the riverine and mixed systems will only marginally be affected by river load changes and increase in atmospheric pCO2, the marine estuary is likely to become a significant CO2 sink in its downstream section. In the decades to come, such change of behavior might strengthen the overall CO2 sink of the estuary-coastal ocean continuum.

Nonintrusive Global Sensitivity Analysis for Linear Systems With Process Noise

2019 ◽  
Vol 14 (2) ◽  
Author(s):  
Souransu Nandi ◽  
Tarunraj Singh
Keyword(s):  
Sensitivity Analysis ◽  
Linear Systems ◽  
Global Sensitivity Analysis ◽  
Benchmark Problems ◽  
Model Parameters ◽  
Parameter Uncertainties ◽  
Sobol Indices ◽  
Global Sensitivity ◽  
Time Invariant ◽  
The Impact

The focus of this paper is on the global sensitivity analysis (GSA) of linear systems with time-invariant model parameter uncertainties and driven by stochastic inputs. The Sobol' indices of the evolving mean and variance estimates of states are used to assess the impact of the time-invariant uncertain model parameters and the statistics of the stochastic input on the uncertainty of the output. Numerical results on two benchmark problems help illustrate that it is conceivable that parameters, which are not so significant in contributing to the uncertainty of the mean, can be extremely significant in contributing to the uncertainty of the variances. The paper uses a polynomial chaos (PC) approach to synthesize a surrogate probabilistic model of the stochastic system after using Lagrange interpolation polynomials (LIPs) as PC bases. The Sobol' indices are then directly evaluated from the PC coefficients. Although this concept is not new, a novel interpretation of stochastic collocation-based PC and intrusive PC is presented where they are shown to represent identical probabilistic models when the system under consideration is linear. This result now permits treating linear models as black boxes to develop intrusive PC surrogates.

Sign in / Sign up

Export Citation Format

Share Document