scholarly journals Global sensitivity analysis for the Rothermel model based on high-dimensional model representation

2015 ◽  
Vol 45 (11) ◽  
pp. 1474-1479 ◽  
Author(s):  
Yaning Liu ◽  
M. Yousuff Hussaini ◽  
Giray Ökten

Rothermel’s wildland surface fire spread model is widely used in North America. The model outputs depend on a number of input parameters, which can be broadly categorized as fuel model, fuel moisture, terrain, and wind parameters. Due to the inevitable presence of uncertainty in the input parameters, knowing the sensitivity of the model output to a given input parameter can be very useful for understanding and controlling the sources of parametric uncertainty. Instead of obtaining the local sensitivity indices, we perform a global sensitivity analysis that considers the synchronous changes of parameters in their respective ranges. The global sensitivity indices corresponding to different parameter groups are computed by constructing the truncated ANOVA – high dimensional model representation for the model outputs with a polynomial expansion approach. We apply global sensitivity analysis to six standard fuel models, namely short grass, tall grass, chaparral, hardwood litter, timber, and light logging slash. Our sensitivity results show similarities, as well as differences, between fuel models. For example, the sensitivities of the input parameters, i.e., fuel depth, low heat content, and wind, are large in all fuel models and as high as 85% of the total model variance in the fuel model light logging slash. On the other hand, the fuel depth explains around 40% of the total variance in the fuel model light logging slash but only 12% of the total variance in the fuel model short grass. The quantification of the importance of parameters across fuel models helps identify the parameters for which additional resources should be used to lower their uncertainty, leading to effective fire management.

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>


Sign in / Sign up

Export Citation Format

Share Document