Copula-based methods for global sensitivity analysis with correlated random variables and stochastic processes under incomplete probability information

Uncertainty and sensitivity analysis of cooling and heating loads for building energy planning

ABSTRACTPROSAIL is a combination of the leaf optical properties spectra (PROSPECT) model and the scattering by arbitrarily inclined leaves (SAIL) canopy bidirectional reflectance model. When modelling forest canopy reflectance using the PROSAIL radiative transfer model, the sensitivities of parameters can affect the modelling accuracy. Traditionally, sensitivities have been assessed using local sensitivity analysis (LSA); however, drawbacks to this approach include a lack of consideration for coupled effects between different parameters. In this study, parameter sensitivities in the PROSAIL model were calculated using two global sensitivity analysis (GSA) methods (the Extended Fourier Amplitude Sensitivity Test (EFAST) method and the Morris method), field measurements, and Landsat 5 Thematic Mapper (TM) data for a Moso bamboo forest. The results of GSA were compared with those of LSA in order to identify the key parameters impacting the Moso bamboo forest canopy reflectance, and to provide a reference for model optimization and vegetation canopy inversion improvement. The results showed that: (1) the sensitivities of six major input parameters of the PROSAIL model were generally consistent with the sorting orders of the two GSA methods, but were not in accordance with those from the LSA method, especially in the mid-infrared band; (2) coupled effects among parameters acting on reflectance simulation in visible light bands were greater than those in infrared bands; (3) the simulated canopy reflectance was evaluated using Landsat 5 TM data, and the results simulated based on LSA analysis showed higher error than those based on GSA analysis, because the LSA method ignored the influence of some parameters on canopy reflectance, e.g. leaf mesophyll structure (N), average leaf angle (ALA), leaf water content (Cw), and leaf dry matter content (Cm). However, GSA was able to fully consider the coupled effects among parameters, and thus identified the sensitive parameters impacting on reflectance more accurately.

Physiologically based pharmacokinetic (PBPK) models describe the absorption, distribution, metabolism, and excretion of chemicals. Probabilistic PBPK models can be used to produce distributional estimates of human equivalent doses (HEDs), which are external measures of human exposure predicted to result in a target internal dose, generated with Monte Carlo sampling and reverse dosimetry calculations. Very low HED percentiles represent individuals that are more sensitive to possible adverse effects of chemical exposures and are therefore frequently used in risk evaluation. Details of the parameter distributions used in probabilistic PBPK models impact HED distributions, but obtaining precise distributional estimates for all the parameters would be challenging. Therefore, we sought to determine methods that can identify the extent of model parameter influence on extreme HED percentiles. We first analyzed published PBPK models for dichloromethane and chloroform (for inhalation and oral exposures given 3 internal target levels each) by identifying the overall relative importance of parameters using global sensitivity analysis (GSA) methods. Then, we used a novel yet computationally expensive method to analyze the stability and sensitivity of extreme HED percentiles to input parameter distributions. Applying the traditional GSA methods allowed us to identify subsets of parameters most influential for accurately determining 1st and 99th percentiles, but the specific parameters included in those subsets varied for different models and exposure scenarios. Our results demonstrate that better characterizing PBPK model uncertainty by using precise distributional details for influential parameters informed by GSA methods may improve confidence in estimates of extreme percentiles of HEDs.

Simulation of fate and transport mechanisms in the porous systems from contaminant transport models are affected by uncertainty associated with input model parameters. Sensitivity analysis (SA) provides tools to quantify the sources of uncertainty to the variations in the model output metrics. To quantify the uncertainty associated with physical non-equilibrium contaminant transport model, global SA was conducted for the problem mimicking reactive transport in the saturated soil column conditions. Five global SA methods, namely Morris, RSA (Regionalized Sensitivity Analysis), Sobol, FAST (Fourier Amplitude Sensitivity Test), and PAWN, were tested based on the temporal moments of contaminant concentrations for two output metrics: zeroth temporal moment (ZTM) and mean residence time (MRT). The ranking order of ten input model parameters from global SA methods for two output metrics was compared. Morris SA implied that the ZTM at outlet of the soil column is highly sensitive to sorption distribution coefficient in the mobile and immobile regions and less sensitive to dispersion coefficient and degradation rate constant in the immobile region. The mass-transfer coefficient showed highest non-linear interactions with other flow and transport parameters based on the highest value of standard deviation of elementary effects (EEs) for ZTM output metric. The sorption distribution coefficient in the mobile region and mass-transfer coefficient showed highest sensitivity and non-linear effect toward MRT based on the Morris method. The comparison of global SA methods revealed that the top two sensitive parameters affecting ZTM and MRT were the same from all the considered methods. However, large difference in the ranking order of bottom three sensitive parameters was observed. The sorption distribution coefficient of mobile region and mass-transfer coefficient were observed as the most sensitive parameters affecting ZTM and MRT based on the comparison of all five global SA methods. Overall results suggest that the non-linear contaminant transport model should be examined based on multiple sensitivity output metrics via a multi-model approach. The multi-model global SA approach implemented in this study highlighted its significance in quantifying the interplay of non-linear model parameters.

Global sensitivity analysis and uncertainty quantification of soot formation in an n-dodecane spray flame

Intelligent sensitivity analysis framework based on machine learning for spacecraft thermal design

Development of a behaviour-pattern based global sensitivity analysis procedure for coupled socioeconomic and environmental models

This paper presents an application of the Global Sensitivity Analysis (GSA) method in evaluating performance of impedance measurement algorithms implemented in modern Intelligent Electronic Devices (IEDs) when affected by uncertainty. To be specific, the GSA is applied in testing an algorithm for distance protection of a single transmission line in a system with two lines on the same towers. Factors contributing to uncertainty are discussed in the paper. The GSA method decomposes variance of the estimation error according to the factors and their interactions, and in this way identifies the factors that have the highest impact on the estimation error. Computation of the GSA sensitivity indices relies on the Quasi-Monte Carlo (QMC) technique with Sobol's quasi-random sampling. This is time-consuming for large number of factors therefore, fast pre-screening process based on the Morris method is implemented to remove non-influential parameters and reduce the factor space before running the QMC. The methodology is implemented using fault simulation software DIgSILENT PowerFactory to calculate fault impedance as a function of uncertain parameters (i.e. factors) as well as voltage and current signals that are inputs to an IED under test. In this simulation, factor space is sampled using Sobol's quasi-random sequence generated with the standard GSA software tool called SIMLAB. All simulations of the transmission line faults for various parameter values and calculation of the performance indices (i.e. sensitivities) is automated using the program developed in DIgSILENT Programming Language (DPL). This testing technique is applied to study the fault impedance measurement implemented in SEL-421. (6 pages)

Systems with random variables and random excitations exist widely in various engineering problems. Extending the traditional global reliability sensitivity to this double-stochastic system has important guiding significance for its design optimization. However, because there is a certain coupling between the randomness of variables and the randomness of excitation, this coupling mechanism is difficult to determine in practical projects. Therefore, it is difficult to extend the traditional reliability sensitivity analysis method to this double-stochastic system. In this research, it is assumed that there is no correlation between variables and excitations. Then, combining the first-passage method–based dynamic strength formula and the variance-based sensitivity analysis method, an approximate global reliability sensitivity analysis method for this double-stochastic system is proposed. In order to improve the computational efficiency, a nested loop method based on seven-point estimation is proposed for reliability sensitivity analysis. In order to verify the accuracy and efficiency of the proposed method, a Monte Carlo simulation is given as a reference. Three examples are studied and discussed to illustrate the practicality and feasibility of the proposed method.

Bayesian calibration and uncertainty analysis of an agroecosystem model under different N management practices

Global sensitivity analysis (GSA) of distribution system with respect to stochastic PV variations plays an important role in designing optimal voltage control schemes. This paper proposes a Kriging, i.e., Gaussian process modeling enabled data-driven GSA method. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol index can be conveniently calculated to assess the global sensitivity of voltage to various power injection variations. Comparison results with other model-based GSA methods on the IEEE 37-bus feeder, such as the polynomial chaos expansion and the Monte Carlo approaches demonstrate that the proposed method can achieve accurate GSA outcomes while maintaining high computational efficiency.

SUMMARYThe reliability of model predictions is affected by multiple sources of uncertainty; therefore, most of the efforts for modeling biological systems include a sensitivity analysis step aiming to identify the key contributors to uncertainty. This generates insight about the robustness of the model to variations in environmental conditions, kinetic parameters, initial concentration of the species, or any other source of uncertainty. Local sensitivities measure the robustness of the model to small perturbations on the inputs around their nominal value. There are several numerical methods for the calculation of local sensitivities, but the calculated values should be identical within the numerical accuracy of the method used. In contrast, as will be shown in this contribution, the results of different global sensitivity analysis methods can be very different and highly dependent on the distribution considered for the inputs under evaluation. In this work, derivative‐based global sensitivities are extended to be able to consider an accurate probability density function for the parameters based on the likelihood function. This strategy enforces the areas of the parameter space most likely to reproduce the desired behavior, minimizing the importance of parameter sets with low probability of being optimal to dominate the sensitivity ranking. A model of a biochemical pathway with three enzymatic steps is used here to illustrate the performance of several relevant global sensitivity analysis methods considering different probability density functions for the parameters and revealing important hints about which method and distribution to choose for each type of model and purpose of the analysis. Copyright © 2012 John Wiley & Sons, Ltd.

BackgroundMathematical modeling is being applied to increasingly complex biological systems and datasets; however, the process of analyzing and calibrating against experimental data is often challenging and a rate limiting step in model development. To address this problem, we developed a systematic methodology for calibrating quantitative models of dynamic biological processes and illustrate its utility by validating a model of TRAIL (Tumor necrosis factor Related Apoptosis-Inducing Ligand)-induced cell death.ResultsWe propose a serial framework integrating analysis and calibration modules and we compare various methods for global sensitivity analysis and global parameter estimation. First, adequacy of the network structure is checked by global sensitivity analysis to changes in concentrations of molecular species, validating that the model can reproduce qualitative features of the system behavior derived from experiments or literature surveys. Second, rate parameters are ranked by importance using gradient-based and variance-based sensitivity indices, and we systematically determine the optimal number of parameters to include in model calibration. Third, deterministic, stochastic and hybrid algorithms for global optimization are applied to estimate the values of the most important parameters by fitting to time series data. We compare the performance of these three optimization algorithms.ConclusionsOur proposed framework covers the entire process from validating a proto-model to establishing a realistic model for in silico experiments and thereby provides a generalized workflow for the construction of predictive models of complex network systems.

Four methods for global sensitivity analysis are studied with respect to their applicability in the finite element analysis of composite structures: (1) sensitivities based on the variance conditional expectation involving experimental design, (2) Morris' elementary effects method (including recent improvements to effective trajectory selection), (3) Sobol's main and total sensitivity indices, and (4) the root mean square of approximate derivatives at sampled points. The sensitivities of the first ten natural frequencies to the composite material properties as well as to the ply stack are determined. The methods are compared and an assessment of their applicability is made in the case of a fairly large number (d > 25) of uncertain parameters.

Mixed kernel function support vector regression for global sensitivity analysis