Sensitivity analysis of an integrated numerical flow model output to key parameters used in common qualitative vulnerability assessment methods

Author(s):  
Joanna Doummar ◽  
Assaad H. Kassem

<p>Qualitative vulnerability assessment methods applied in karst aquifers rely on key factors in the hydrological compartments usually assigned different weights according to their estimated impact on groundwater vulnerability. Based on an integrated numerical groundwater model on a snow-governed karst catchment area (Assal Spring- Lebanon), the aim of this work is to quantify the importance of the most influential parameters on recharge and spring discharge and outline potential parameters that are not accounted for in standard methods, when in fact they do play a role in the intrinsic vulnerability of a system. The assessment of the model sensitivity and the ranking of parameters are conducted using an automatic calibration tool for local sensitivity analysis in addition to a variance-based local sensitivity assessment of model output time series (recharge and discharge)  for two consecutive years (2016-2017) to various model parameters. The impact of each parameter was normalized to estimate standardized weights for each of the process based key-controlling parameters. Parameters to which model was sensitive were factors related to soil, 2) fast infiltration (bypass function) typical of karst aquifers, 3) climatic parameters (melting temperature and degree day coefficient) and 4) aquifer hydraulic properties that play a major role in groundwater vulnerability inducing a temporal effect and varied recession. Other less important parameters play different roles according to different assigned weights proportional to their ranking. Additionally, the effect of slope/geomorphology (e.g., dolines) was further investigated.  In general, this study shows that the weighting coefficients assigned to key vulnerability factors in the qualitative assessment methods can be reevaluated based on this process-based approach.</p><p> </p><p> </p><p> </p>

2017 ◽  
Vol 17 (1) ◽  
pp. 18-30
Author(s):  
A. Ewusi ◽  
I. Ahenkorah ◽  
J. S. Y. Kuma

Groundwater vulnerability assessment to delineate areas that are susceptible to contamination from mining and anthropogenic activities has become an important element for resource management and landuse planning. In view of the extensive mining in the Tarkwa area, quality of groundwater has become an important issue. This study estimates aquifer vulnerability by applying the SINTACS model which uses seven environmental parameters to evaluate aquifer vulnerability and geographical information system (GIS) in the Tarkwa mining area. Sensitivity analysis has also been carried out to evaluate the relative importance of the model parameters for aquifer vulnerability. The SINTACS model results show that the intrusive rocks within the Tarkwaian and the Birimian rocks are dominated by very high vulnerability classes while the Banket Series is characterised by high vulnerability class. The Huni Sandstones have moderately high vulnerability. In addition, the Kawere Group and the Tarkwa Phyllites displayed medium vulnerability. Analysis from the variogram model shows that all parameters used in the SINTACS model have a strong spatial structure. From statistical analysis, depth to water parameter inflicted the highest impact on the vulnerability of the aquifer followed by effective infiltration, vadose zone media, soil media, aquifer media, topography and hydraulic conductivity in the order of decreasing impact. Sensitivity analysis indicated that the aquifer media, hydraulic characteristics and topography cause large variation in vulnerability index. Depth to water and effective infiltration were found to be more effective in assessing aquifer vulnerability. Keywords: Groundwater, Vulnerability, Tarkwa, SINTACS, GIS


2013 ◽  
Vol 17 (2) ◽  
pp. 461-478 ◽  
Author(s):  
L. Loosvelt ◽  
H. Vernieuwe ◽  
V. R. N. Pauwels ◽  
B. De Baets ◽  
N. E. C. Verhoest

Abstract. Compositional data, such as soil texture, are hard to deal with in the geosciences as standard statistical methods are often inappropriate to analyse this type of data. Especially in sensitivity analysis, the closed character of the data is often ignored. To that end, we developed a method to assess the local sensitivity of a model output with resect to a compositional model input. We adapted the finite difference technique such that the different parts of the input are perturbed simultaneously while the closed character of the data is preserved. This method was applied to a hydrologic model and the sensitivity of the simulated soil moisture content to local changes in soil texture was assessed. Based on a high number of model runs, in which the soil texture was varied across the entire texture triangle, we identified zones of high sensitivity in the texture triangle. In such zones, the model output uncertainty induced by the discrepancy between the scale of measurement and the scale of model application, is advised to be reduced through additional data collection. Furthermore, the sensitivity analysis provided more insight into the hydrologic model behaviour as it revealed how the model sensitivity is related to the shape of the soil moistureretention curve.


2021 ◽  
Author(s):  
Sabine Bauer ◽  
Ivanna Kramer

The knowledge about the impact of structure-specific parameters on the biomechanical behavior of a computer model has an essential meaning for the realistic modeling and system improving. Especially the biomechanical parameters of the intervertebral discs, the ligamentous structures and the facet joints are seen in the literature as significant components of a spine model, which define the quality of the model. Therefore, it is important to understand how the variations of input parameters for these components affect the entire model and its individual structures. Sensitivity analysis can be used to gain the required knowledge about the correlation of the input and output variables in a complex spinal model. The present study analyses the influence of the biomechanical parameters of the intervertebral disc using different sensitivity analysis methods to optimize the spine model parameters. The analysis is performed with a multi-body simulation model of the cervical functional spinal unit C6-C7.


Author(s):  
Souransu Nandi ◽  
Tarunraj Singh

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.


2021 ◽  
Author(s):  
Sabine M. Spiessl ◽  
Dirk-A. Becker ◽  
Sergei Kucherenko

<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>


2011 ◽  
Vol 11 (9) ◽  
pp. 2567-2582 ◽  
Author(s):  
H. Roux ◽  
D. Labat ◽  
P.-A. Garambois ◽  
M.-M. Maubourguet ◽  
J. Chorda ◽  
...  

Abstract. A spatially distributed hydrological model, dedicated to flood simulation, is developed on the basis of physical process representation (infiltration, overland flow, channel routing). Estimation of model parameters requires data concerning topography, soil properties, vegetation and land use. Four parameters are calibrated for the entire catchment using one flood event. Model sensitivity to individual parameters is assessed using Monte-Carlo simulations. Results of this sensitivity analysis with a criterion based on the Nash efficiency coefficient and the error of peak time and runoff are used to calibrate the model. This procedure is tested on the Gardon d'Anduze catchment, located in the Mediterranean zone of southern France. A first validation is conducted using three flood events with different hydrometeorological characteristics. This sensitivity analysis along with validation tests illustrates the predictive capability of the model and points out the possible improvements on the model's structure and parameterization for flash flood forecasting, especially in ungauged basins. Concerning the model structure, results show that water transfer through the subsurface zone also contributes to the hydrograph response to an extreme event, especially during the recession period. Maps of soil saturation emphasize the impact of rainfall and soil properties variability on these dynamics. Adding a subsurface flow component in the simulation also greatly impacts the spatial distribution of soil saturation and shows the importance of the drainage network. Measures of such distributed variables would help discriminating between different possible model structures.


2018 ◽  
Vol 24 (3) ◽  
pp. 293-304 ◽  
Author(s):  
Ismail Chenini ◽  
Adel Zghibi ◽  
Mohamed Haythem Msaddek ◽  
Mahmoud Dlala

Abstract The groundwater vulnerability assessment is normally applied to rural watersheds. However, urbanization modifies the hydrogeological processes. A modified DRASTIC model was adopted to establish a groundwater vulnerability map in an urbanized watershed. The modified DRASTIC model incorporated a land-use map, and net recharge was calculated taking into account the specificity of the urban hydrogeological system. The application of the proposed approach to the Mannouba watershed demonstrates that the groundwater vulnerability indexes range from 80 to 165. The study's results shows that 30 percent of the Mannouba watershed area has a high vulnerability index, 45 percent of the area has a medium index, and 25 percent of the study area has a low vulnerability index. To specify the effect of each DRASTIC factor on the calculated vulnerability index, sensitivity analyses were performed. Land use, topography, and soil media have an important theoretical weight greater than the effective weight. The impact of the vadose zone factor has the most important effective weight and affects the vulnerability index. The sensitivity assessment explored the variation in vulnerability after thematic layer removal. In this analysis, the removal of hydraulic conductivity and impact of vadose zone modified the vulnerability index. Groundwater vulnerability assessment in urbanized watersheds is difficult and has to consider the impact of urbanization in the hydrogeological parameters.


2020 ◽  
Author(s):  
Monica Riva ◽  
Aronne Dell'Oca ◽  
Alberto Guadagnini

<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>


Sign in / Sign up

Export Citation Format

Share Document