scholarly journals Hierarchical Sensitivity Analysis for Large-scale Process-based Hydrological Modeling with Application in an Amazonian Watershed

2020 ◽  
Author(s):  
Haifan Liu ◽  
Heng Dai ◽  
Jie Niu ◽  
Bill X. Hu ◽  
Han Qiu ◽  
...  
Keyword(s):  
Sensitivity Analysis ◽  
Vadose Zone ◽  
Model Calibration ◽  
Large Scale ◽  
Global Sensitivity Analysis ◽  
Hydrologic Model ◽  
Global Sensitivity ◽  
Uncertainty Sources ◽  
Sensitivity Indices ◽  
Zone Parameters

Abstract. Sensitivity analysis is an effective tool for identifying important uncertainty sources and improving model calibration and predictions, especially for integrated systems with heterogeneous parameter inputs and complex process coevolution. In this work, an advanced hierarchical global sensitivity analysis framework, which integrates the concept of variance-based global sensitivity analysis with a hierarchical uncertainty framework, was implemented to quantitatively analyse several uncertainties associated with a three-dimensional, process-based hydrologic model (PAWS). The uncertainty sources considered include model parameters (three vadose zone parameters, two groundwater parameters, and one overland flow parameter), model structure (different thicknesses to represent unconfined and confined aquifer layers) and climate scenarios. We apply the approach to an ~ 9,000 km2 Amazon catchment modeled at 1 km resolution to provide a demonstration of multiple uncertainty source quantification using a large-scale process-based hydrologic model. The sensitivity indices are assessed based on two important hydrologic outputs: evapotranspiration (ET) and groundwater contribution to streamflow (QG). It was found that, in general, parameters, especially the vadose zone parameters, are the most important uncertainty contributors for all sensitivity indices. In addition, the influence of climate scenarios on ET predictions is also nonignorable. Furthermore, the thickness of the aquifers along the river grid cells is important for both ET and QG. These results can assist in model calibration and provide modelers with a better understanding of the general sources of uncertainty in predictions associated with complex hydrological systems in Amazonia. We demonstrated a pilot example of comprehensive global sensitivity analysis of large-scale, complex hydrological and environmental models in this study. The hierarchical sensitivity analysis methodology used is mathematically rigorous and capable of being implemented in a variety of large-scale hydrological models with various sources of uncertainty.

scholarly journals Hierarchical Sensitivity Analysis for Large Scale Process-based Hydrological Modeling with Application in an Amazonian Watershed

2019 ◽  
Author(s):  
Haifan Liu ◽  
Heng Dai ◽  
Jie Niu ◽  
Bill X. Hu ◽  
Han Qiu ◽  
...  
Keyword(s):  
Sensitivity Analysis ◽  
Large Scale ◽  
Overland Flow ◽  
Global Sensitivity Analysis ◽  
Hydrologic Model ◽  
Model Parameters ◽  
Hydrological Models ◽  
Global Sensitivity ◽  
Uncertainty Sources ◽  
Sensitivity Indices

Abstract. Sensitivity analysis is an effective tool for identifying important uncertainty sources and improving model calibration and predictions, especially for integrated systems with heterogeneous parameter inputs and complex processes coevolution. In this work, an advanced hierarchical global sensitivity analysis framework, which integrates a hierarchical uncertainty framework and a variance-based global sensitivity analysis, was implemented to quantitatively analyze several uncertainties of a three-dimensional, process-based hydrologic model (PAWS). The uncertainty sources considered include model parameters, model structures (with/without overland flow module), and climate forcing. We apply the approach in a ~ 9000 km2 Amazon catchment modeled at 1 km resolution to provide a demonstration of multiple uncertainty source quantification using a large-scale process-based hydrologic model. The sensitivity indices are assessed based on three important hydrologic outputs: evapotranspiration (ET), ground evaporation (EG), and groundwater contribution to streamflow (QG). It is found that, in general, model parameters (especially those within the streamside model grid cells) are the most important uncertainty contributor for all sensitivity indices. In addition, the overland flow module significantly contributes to model predictive uncertainty. These results can assist model calibration and provide modelers a better understanding of the general sources of uncertainty in predictions of complex hydrological systems in Amazonia. We demonstrated a pilot example for comprehensive global sensitivity analysis of large-scale complex hydrological models in this research. The hierarchical sensitivity analysis methodology used is mathematically rigorous and can be applied to a wide range of large-scale hydrological models with various sources of uncertainty.

scholarly journals Hierarchical sensitivity analysis for a large-scale process-based hydrological model applied to an Amazonian watershed

2020 ◽  
Vol 24 (10) ◽  
pp. 4971-4996
Author(s):  
Haifan Liu ◽  
Heng Dai ◽  
Jie Niu ◽  
Bill X. Hu ◽  
Dongwei Gui ◽  
...  
Keyword(s):  
Sensitivity Analysis ◽  
Large Scale ◽  
Global Sensitivity Analysis ◽  
Hydrological Models ◽  
Climate Scenarios ◽  
Analysis Method ◽  
Global Sensitivity ◽  
Uncertainty Sources ◽  
Sensitivity Indices ◽  
Sensitivity Analysis Method

Abstract. Sensitivity analysis methods have recently received much attention for identifying important uncertainty sources (or uncertain inputs) and improving model calibrations and predictions for hydrological models. However, it is still challenging to apply the quantitative and comprehensive global sensitivity analysis method to complex large-scale process-based hydrological models (PBHMs) because of its variant uncertainty sources and high computational cost. Therefore, a global sensitivity analysis method that is capable of simultaneously analyzing multiple uncertainty sources of PBHMs and providing quantitative sensitivity analysis results is still lacking. In an effort to develop a new tool for overcoming these weaknesses, we improved the hierarchical sensitivity analysis method by defining a new set of sensitivity indices for subdivided parameters. A new binning method and Latin hypercube sampling (LHS) were implemented for estimating these new sensitivity indices. For test and demonstration purposes, this improved global sensitivity analysis method was implemented to quantify three different uncertainty sources (parameters, models, and climate scenarios) of a three-dimensional large-scale process-based hydrologic model (Process-based Adaptive Watershed Simulator, PAWS) with an application case in an ∼ 9000 km2 Amazon catchment. The importance of different uncertainty sources was quantified by sensitivity indices for two hydrologic outputs of interest: evapotranspiration (ET) and groundwater contribution to streamflow (QG). The results show that the parameters, especially the vadose zone parameters, are the most important uncertainty contributors for both outputs. In addition, the influence of climate scenarios on ET predictions is also important. Furthermore, the thickness of the aquifers is important for QG predictions, especially in main stream areas. These sensitivity analysis results provide useful information for modelers, and our method is mathematically rigorous and can be applied to other large-scale hydrological models.

Sensitivity analysis and the challenges posed by multiple approaches: a multifaceted mess

2020 ◽  
Author(s):  
Monica Riva ◽  
Aronne Dell'Oca ◽  
Alberto Guadagnini
Keyword(s):  
Sensitivity Analysis ◽  
Model Calibration ◽  
Global Sensitivity Analysis ◽  
System Model ◽  
Model Parameters ◽  
Model Output ◽  
Multiple Model ◽  
Total Uncertainty ◽  
Global Sensitivity ◽  
Sensitivity Indices

<p>Modern models of environmental and industrial systems have reached a relatively high level of complexity. The latter aspect could hamper an unambiguous understanding of the functioning of a model, i.e., how it drives relationships and dependencies among inputs and outputs of interest. Sensitivity Analysis tools can be employed to examine this issue.</p><p>Global sensitivity analysis (GSA) approaches rest on the evaluation of sensitivity across the entire support within which system model parameters are supposed to vary. In this broad context, it is important to note that the definition of a sensitivity metric must be linked to the nature of the question(s) the GSA is meant to address. These include, for example: (i) which are the most important model parameters with respect to given model output(s)?; (ii) could we set some parameter(s) (thus assisting model calibration) at prescribed value(s) without significantly affecting model results?; (iii) at which space/time locations can one expect the highest sensitivity of model output(s) to model parameters and/or knowledge of which parameter(s) could be most beneficial for model calibration?</p><p>The variance-based Sobol’ Indices (e.g., Sobol, 2001) represent one of the most widespread GSA metrics, quantifying the average reduction in the variance of a model output stemming from knowledge of the input. Amongst other techniques, Dell’Oca et al. [2017] proposed a moment-based GSA approach which enables one to quantify the influence of uncertain model parameters on the (statistical) moments of a target model output.</p><p>Here, we embed in these sensitivity indices the effect of uncertainties both in the system model conceptualization and in the ensuing model(s) parameters. The study is grounded on the observation that physical processes and natural systems within which they take place are complex, rendering target state variables amenable to multiple interpretations and mathematical descriptions. As such, predictions and uncertainty analyses based on a single model formulation can result in statistical bias and possible misrepresentation of the total uncertainty, thus justifying the assessment of multiple model system conceptualizations. We then introduce copula-based sensitivity metrics which allow characterizing the global (with respect to the input) value of the sensitivity and the degree of variability (across the whole range of the input values) of the sensitivity for each value that the prescribed model output can possibly undertake, as driven by a governing model. In this sense, such an approach to sensitivity is global with respect to model input(s) and local with respect to model output, thus enabling one to discriminate the relevance of an input across the entire range of values of the modeling goal of interest. The methodology is demonstrated in the context of flow and reactive transport scenarios.</p><p> </p><p><strong>References</strong></p><p>Sobol, I. M., 2001. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Math. Comput. Sim., 55, 271-280.</p><p>Dell’Oca, A., Riva, M., Guadagnini, A., 2017. Moment-based metrics for global sensitivity analysis of hydrological systems. Hydr. Earth Syst. Sci., 21, 6219-6234.</p>

scholarly journals A Global Sensitivity Analysis Framework for Hybrid Simulation with Stochastic Substructures

2021 ◽  
Vol 7 ◽  
Author(s):  
Nikolaos Tsokanas ◽  
Xujia Zhu ◽  
Giuseppe Abbiati ◽  
Stefano Marelli ◽  
Bruno Sudret ◽  
...  
Keyword(s):  
Sensitivity Analysis ◽  
Hybrid Model ◽  
Surrogate Model ◽  
Global Sensitivity Analysis ◽  
Hybrid Simulation ◽  
Hybrid Models ◽  
Analysis Framework ◽  
Global Sensitivity ◽  
Sensitivity Indices ◽  
Realistic Situation

Hybrid simulation is an experimental method used to investigate the dynamic response of a reference prototype structure by decomposing it to physically-tested and numerically-simulated substructures. The latter substructures interact with each other in a real-time feedback loop and their coupling forms the hybrid model. In this study, we extend our previous work on metamodel-based sensitivity analysis of deterministic hybrid models to the practically more relevant case of stochastic hybrid models. The aim is to cover a more realistic situation where the physical substructure response is not deterministic, as nominally identical specimens are, in practice, never actually identical. A generalized lambda surrogate model recently developed by some of the authors is proposed to surrogate the hybrid model response, and Sobol’ sensitivity indices are computed for substructure quantity of interest response quantiles. Normally, several repetitions of every single sample of the inputs parameters would be required to replicate the response of a stochastic hybrid model. In this regard, a great advantage of the proposed framework is that the generalized lambda surrogate model does not require repeated evaluations of the same sample. The effectiveness of the proposed hybrid simulation global sensitivity analysis framework is demonstrated using an experiment.

scholarly journals A Global Sensitivity Analysis Framework for Hybrid Simulation

2021 ◽  
Author(s):  
Giuseppe Abbiati ◽  
Stefano Marelli ◽  
Nikolaos Tsokanas ◽  
Bruno Sudret ◽  
Bozidar Stojadinovic
Keyword(s):  
Sensitivity Analysis ◽  
Hybrid Model ◽  
Polynomial Chaos ◽  
Global Sensitivity Analysis ◽  
Hybrid Simulation ◽  
Polynomial Chaos Expansion ◽  
Chaos Expansion ◽  
Analysis Framework ◽  
Global Sensitivity ◽  
Sensitivity Indices

Hybrid Simulation is a dynamic response simulation paradigm that merges physical experiments and computational models into a hybrid model. In earthquake engineering, it is used to investigate the response of structures to earthquake excitation. In the context of response to extreme loads, the structure, its boundary conditions, damping, and the ground motion excitation itself are all subjected to large parameter variability. However, in current seismic response testing practice, Hybrid Simulation campaigns rely on a few prototype structures with fixed parameters subjected to one or two ground motions of different intensity. While this approach effectively reveals structural weaknesses, it does not reveal the sensitivity of structure's response. This thus far missing information could support the planning of further experiments as well as drive modeling choices in subsequent analysis and evaluation phases of the structural design process.This paper describes a Global Sensitivity Analysis framework for Hybrid Simulation. This framework, based on Sobol' sensitivity indices, is used to quantify the sensitivity of the response of a structure tested using the Hybrid Simulation approach due to the variability of the prototype structure and the excitation parameters. Polynomial Chaos Expansion is used to surrogate the hybrid model response. Thereafter, Sobol' sensitivity indices are obtained as a by-product of polynomial coefficients, entailing a reduced number of Hybrid Simulations compared to a crude Monte Carlo approach. An experimental verification example highlights the excellent performance of Polynomial Chaos Expansion surrogates in terms of stable estimates of Sobol' sensitivity indices in the presence of noise caused by random experimental errors.

Comprehensive global sensitivity analysis of a repository model using different types of transformations and metamodeling techniques

2021 ◽  
Author(s):  
Sabine M. Spiessl ◽  
Dirk-A. Becker ◽  
Sergei Kucherenko
Keyword(s):  
Sensitivity Analysis ◽  
Global Sensitivity Analysis ◽  
Higher Order ◽  
Model Parameters ◽  
High Dimensional Model Representation ◽  
Model Output ◽  
Global Sensitivity ◽  
Sparse Polynomial ◽  
Sensitivity Indices ◽  
Highly Nonlinear

<p>Due to their highly nonlinear, non-monotonic or even discontinuous behavior, sensitivity analysis of final repository models can be a demanding task. Most of the output of repository models is typically distributed over several orders of magnitude and highly skewed. Many values of a probabilistic investigation are very low or even zero. Although this is desirable in view of repository safety it can distort the evidence of sensitivity analysis. For the safety assessment of the system, the highest values of outputs are mainly essential and if those are only a few, their dependence on specific parameters may appear insignificant. By applying a transformation, different model output values are differently weighed, according to their magnitude, in sensitivity analysis. Probabilistic methods of higher-order sensitivity analysis, applied on appropriately transformed model output values, provide a possibility for more robust identification of relevant parameters and their interactions. This type of sensitivity analysis is typically done by decomposing the total unconditional variance of the model output into partial variances corresponding to different terms in the ANOVA decomposition. From this, sensitivity indices of increasing order can be computed. The key indices used most often are the first-order index (SI1) and the total-order index (SIT). SI1 refers to the individual impact of one parameter on the model and SIT represents the total effect of one parameter on the output in interactions with all other parameters. The second-order sensitivity indices (SI2) describe the interactions between two model parameters.</p><p>In this work global sensitivity analysis has been performed with three different kinds of output transformations (log, shifted and Box-Cox transformation) and two metamodeling approaches, namely the Random-Sampling High Dimensional Model Representation (RS-HDMR) [1] and the Bayesian Sparse PCE (BSPCE) [2] approaches. Both approaches are implemented in the SobolGSA software [3, 4] which was used in this work. We analyzed the time-dependent output with two approaches for sensitivity analysis, i.e., the pointwise and generalized approaches. With the pointwise approach, the output at each time step is analyzed independently. The generalized approach considers averaged output contributions at all previous time steps in the analysis of the current step. Obtained results indicate that robustness can be improved by using appropriate transformations and choice of coefficients for the transformation and the metamodel.</p><p>[1] M. Zuniga, S. Kucherenko, N. Shah (2013). Metamodelling with independent and dependent inputs. Computer Physics Communications, 184 (6): 1570-1580.</p><p>[2] Q. Shao, A. Younes, M. Fahs, T.A. Mara (2017). Bayesian sparse polynomial chaos expansion for global sensitivity analysis. Computer Methods in Applied Mechanics and Engineering, 318: 474-496.</p><p>[3] S. M. Spiessl, S. Kucherenko, D.-A. Becker, O. Zaccheus (2018). Higher-order sensitivity analysis of a final repository model with discontinuous behaviour. Reliability Engineering and System Safety, doi: https://doi.org/10.1016/j.ress.2018.12.004.</p><p>[4] SobolGSA software (2021). User manual https://www.imperial.ac.uk/process-systems-engineering/research/free-software/sobolgsa-software/.</p>

scholarly journals Global Sensitivity Analysis of Fuzzy Distribution Parameter on Failure Probability and Its Single-Loop Estimation

2014 ◽  
Vol 2014 ◽  
pp. 1-11 ◽  
Author(s):  
Lei Cheng ◽  
Zhenzhou Lu ◽  
Luyi Li
Keyword(s):  
Sensitivity Analysis ◽  
Failure Probability ◽  
Global Sensitivity Analysis ◽  
Computational Cost ◽  
Single Loop ◽  
Global Sensitivity ◽  
Sensitivity Indices ◽  
Distribution Parameters ◽  
Input Variables ◽  
Sampling Probability

An extending Borgonovo’s global sensitivity analysis is proposed to measure the influence of fuzzy distribution parameters on fuzzy failure probability by averaging the shift between the membership functions (MFs) of unconditional and conditional failure probability. The presented global sensitivity indices can reasonably reflect the influence of fuzzy-valued distribution parameters on the character of the failure probability, whereas solving the MFs of unconditional and conditional failure probability is time-consuming due to the involved multiple-loop sampling and optimization operators. To overcome the large computational cost, a single-loop simulation (SLS) is introduced to estimate the global sensitivity indices. By establishing a sampling probability density, only a set of samples of input variables are essential to evaluate the MFs of unconditional and conditional failure probability in the presented SLS method. Significance of the global sensitivity indices can be verified and demonstrated through several numerical and engineering examples.

Sign in / Sign up

Export Citation Format

Share Document