Global and local sensitivity analysis of the Emission Dispersion Model input parameters

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.

A review of sensitivity analysis methods in building energy analysis

Local and global sensitivity analysis (SA) methods were compared to demonstrate how the choice of SA method leads to differences in the quantification of the relative importance of model parameters in driving canopy reflectance variability. The global SA method used was the Extended Fourier Amplitude Sensitivity Test (EFAST), while a series of basecase simulations were defined for the local SA, and model parameters were perturbed sequentially. Canopy reflectance was simulated using the ProGOSAIL model. Some pronounced differences were found between the two methods in terms of the magnitude of the importance, the rank importance, and the wavelengths at which model parameters are most important in explaining reflectance variability.

There is uncertainty in the results of any mathematical model due to different reasons. It is important to estimate this uncertainty. Sensitivity analysis is commonly used to estimate how the changes in the input parameters affect the solutions of the model. In this paper, we discuss different ways of performing local and global sensitivity analyses and apply them to two models: an epidemic model and a new myocardial infarction model, both based on ordinary differential equations. The first model is a simple model used to explain the ideas, while the second one shows how to apply them to a model with more state variables and parameters. We find that if the parameters are not accurately known, local sensitivity analysis can be misleading and that global sensitivity methods that sample the whole parameter space, varying all the values of the parameters at the same time, are the most reliable. We also show how the sensitivity analysis results can be used to determine the uncertainty in the results of the model. We present numerical simulations.

Dimensionality reduction for surrogate model construction for global sensitivity analysis: Comparison between active subspace and local sensitivity analysis

In this paper, we study the dynamical effects of timely and delayed diagnosis on the spread of COVID-19 in Ghana during its initial phase by using reported data from March 12 to June 19, 2020. The estimated basic reproduction number, ℛ0, for the proposed model is 1.04. One of the main focus of this study is global stability results. Theoretically and numerically, we show that the disease persistence depends on ℛ0. We carry out a local and global sensitivity analysis. The local sensitivity analysis shows that the most positive sensitive parameter is the recruitment rate, followed by the relative transmissibility rate from the infectious with delayed diagnosis to the susceptible individuals. And that the most negative sensitive parameters are: self-quarantined, waiting time of the infectious for delayed diagnosis and the proportion of the infectious with timely diagnosis. The global sensitivity analysis using the partial rank correlation coefficient confirms the directional flow of the local sensitivity analysis. For public health benefit, our analysis suggests that, a reduction in the inflow of new individuals into the country or a reduction in the inter community inflow of individuals will reduce the basic reproduction number and thereby reduce the number of secondary infections (multiple peaks of the infection). Other recommendations for controlling the disease from the proposed model are provided in Section 7.

The aim of sensitivity analysis (SA) is to ascertain how much the uncertainty in the output of a model is influenced by the uncertainty in its input factors. An SA can be performed using different methods, which are classified according to various criteria. One possible classification is that in which global and local approaches are identified. This paper strengthens the role of global SA methods and suggests their use in the context of MCDA. Useful applications of global SA already exist in a variety of fields where numerical models are considered, e.g. in economics, engineering and chemistry. Global SA is quite different in its formulation and application from local SA, which is seen more frequently in the literature. The global sensitivity indices adopted are derived from a variance decomposition scheme: they can be estimated by alternative computational strategies, such as the extended FAST or the method proposed by Sobol’. A truly global SA is capable of apportioning the output uncertainty according to any subgroup of input factors. Hence, the output uncertainty of an MCDA can, for instance, be decomposed into a part due to uncertain model inputs and a part due to poorly defined (or variable) weights attached to the criteria. This information could be useful to the decision maker (DM) since it explains synthetically how much the assessment of an MCDA study is biased by the assessor judgements. Alternative regrouping of the uncertain input elements might shed light on other features of the problem addressed by the DM. Copyright © 1999 John Wiley & Sons, Ltd.

The large number of input factors involved in a sophisticated geotechnical computational model is a challenge in the concept of probabilistic analysis. In the context of model calibration and validation, conducting a sensitivity analysis is substantial as a first step. Sensitivity analysis techniques can determine the key factors which govern the system responses. In this paper, three commonly used sensitivity analysis methods are implemented on a sophisticated geotechnical problem. The computational model of a compressed air energy storage, mined in a rock salt formation, includes many input parameters, each with large amount of uncertainties. Sensitivity measures of different variables involved in the mechanical response of the cavern are computed by different global sensitivity methods, namely, Sobol/Saltelli, Random Balance Design, and Elementary Effect method. Since performing sensitivity analysis requires a large number of model evaluations, the concept of surrogate modelling is utilised to decrease the computational burden. In the following, the accuracy levels of various surrogate techniques are compared. In addition, a comparative study on the applied sensitivity analysis methods shows that the applied sensitivity analysis techniques provide identical parameter importance rankings, although some may also give more information about the system behaviour.

<strong class="journal-contentHeaderColor">Abstract.</strong> The present study focuses on identifying the parameters from the Weather Research and Forecasting (WRF) model that strongly influence the simulation of tropical cyclones over the Bay of Bengal (BoB) region. Three global sensitivity analysis (SA) methods, namely, the Morris One-at-A-Time (MOAT), multivariate adaptive regression splines (MARS), and surrogate-based Sobol', are employed to identify the most sensitive parameters out of 24 tunable parameters corresponding to seven parameterization schemes of the WRF model. Ten tropical cyclones across different categories, such as cyclonic storms, severe cyclonic storms, and very severe cyclonic storms over BoB between 2011 and 2018, are selected in this study. The sensitivity scores of 24 parameters are evaluated for eight meteorological variables. The parameter sensitivity results are consistent across three SA methods for all the variables, and 8 out of the 24 parameters contribute 80 %–90 % to the overall sensitivity scores. It is found that the Sobol' method with Gaussian progress regression as a surrogate model can produce reliable sensitivity results when the available samples exceed 200. The parameters with which the model simulations have the least RMSE values when compared with the observations are considered the optimal parameters. Comparing observations and model simulations with the default and optimal parameters shows that simulations with the optimal set of parameters yield a 16.74 % improvement in the 10 m wind speed, 3.13 % in surface air temperature, 0.73 % in surface air pressure, and 9.18 % in precipitation simulations compared to the default set of parameters.

Global Sensitivity Analysis of Arrhythmic Risk Biomarkers in Cardiac Single Cell Models

ABSTRACTHydrological models demand large numbers of input parameters, which are to be optimally identified for better simulation of various hydrological processes. Identifying the most relevant parameters and their values using efficient sensitivity analysis methods helps to better understand model performance. In this study, the physically-based distributed model SHETRAN is used for hydrological simulation on the Netravathi River Basin in south India and the most important parameters are identified using the Morris screening method. Further, the influence of a particular model parameter on streamflow is quantified using local sensitivity analysis and optimal parameters are obtained for calibration of the SHETRAN model. The results demonstrate the capability of two-stage sensitivity analysis, combining qualitative and quantitative methods in the initial screening-out of insignificant model parameters, identifying parameter interactions and quantifying the contribution of each model parameter to the streamflow. The results of the sensitivity analysis simplified the calibration procedure of SHETRAN for the study area.

The uncertainty of parameters will have a significant impact on slope stability, where sensitivity analysis is a commonly used method in uncertainty research. However, traditional sensitivity analysis method costs much computation time. When calculating the sensitivity index of one parameter, all other parameters are taken as fixed values, and the uncertainty of all parameters cannot be considered simultaneously. Therefore, the variance-based and the moment-independent global sensitivity analysis (GSA) methods are both introduced to determine the influence of geotechnical parameters on slope stability in this study. To solve the importance index of GSA, the least angle regression algorithm, the kernel density estimation, and orthogonal polynomial estimation methods are developed to obtain variance-based importance index and the moment-independent importance index, respectively. The proposed methods allow all variables to change simultaneously within their variation range and have high computational efficiency. The results are in good at with those obtained by the variance-based Monte Carlo simulation method, which is considered as the exact solution forobtaining the importance index. The influence of the correlation between the shear strength parameters (c and φ) on the importance index is also studied, which indicates that the negative correlation will have a great impact on the importance index, which in turn affects the safety assessment of slope. Three engineering cases have been studied for engineering application, and the compared results indicate that the impact of the geotechnical parameters uncertainty on the safety factor (Fs) and failure probability (pf) are different. Therefore, the approaches based on GSA which can integrate the Fs with pf will be a promising approach for slope stability evaluation.

Performing uncertainty quantification (UQ) and sensitivity analysis (SA) is vital when developing a patient-specific physiological model because it can quantify model output uncertainty and estimate the effect of each of the model's input parameters on the mathematical model. By providing this information, UQ and SA act as diagnostic tools to evaluate model fidelity and compare model characteristics with expert knowledge and real world observation. Computational efficiency is an important part of UQ and SA methods and thus optimization is an active area of research. In this work, we investigate a new efficient sampling method for least-squares polynomial approximation, weighted approximate Fekete points (WAFP). We analyze the performance of this method by demonstrating its utility in stochastic analysis of a cardiovascular model that estimates changes in oxyhemoglobin saturation response. Polynomial chaos (PC) expansion using WAFP produced results similar to the more standard Monte Carlo in quantifying uncertainty and identifying the most influential model inputs (including input interactions) when modeling oxyhemoglobin saturation, PC expansion using WAFP was far more efficient. These findings show the usefulness of using WAFP based PC expansion to quantify uncertainty and analyze sensitivity of a oxyhemoglobin dissociation response model. Applying these techniques could help analyze the fidelity of other relevant models in preparation for clinical application.

This paper is concerned with the shape sensitivity analysis of a class of boundary control, constrained, optimal control problems for parabolic systems. The notion of Euler and Lagrange derivatives of a boundary optimal control in the direction of a vector field is introduced. The derivatives are obtained in the form of optimal solutions of auxiliary optimal control problems. The method of sensitivity analysis used in this paper is based on related results on the differential stability of metric projections in Hilbert space onto a convex, closed subset, combined with the material derivative method of shape sensitivity analysis. Parabolic initial-boundary value problems with Dirichlet and Neumann boundary conditions are considered in this paper.

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