scholarly journals Global Sensitivity Analysis of Parameter Uncertainty in Landscape Evolution Models

Author(s):  
Christopher J. Skinner ◽  
Tom J. Coulthard ◽  
Wolfgang Schwanghart ◽  
Marco J. Van De Wiel ◽  
Greg Hancock

Abstract. Landscape Evolution Models have a long history of use as exploratory models, providing greater understanding of the role large scale processes have on the long-term development of the Earth’s surface. As computational power has advanced so has the development and sophistication of these models. This has seen them applied at increasingly smaller scale and shorter-term simulations at greater detail. However, this has not gone hand-in-hand with more rigorous verifications that are commonplace in the applications of other types of environmental models- for example Sensitivity Analyses. This can be attributed to a paucity of data and methods available in order to calibrate, validate and verify the models, and also to the extra complexity Landscape Evolution Models represent – without these it is not possible to produce a reliable Objective Function against which model performance can be judged. To overcome this deficiency, we present a set of Model Functions – each representing an aspect of model behaviour – and use these to assess the relative sensitivity of a Landscape Evolution Model (CAESAR-Lisflood) to a large set of parameters via a global Sensitivity Analysis using the Morris Method. This novel combination of behavioural Model Functions and the Morris Method provides insight into which parameters are the greatest source of uncertainty in the model, and which have the greatest influence over different model behaviours. The method was repeated over two different catchments, showing that across both catchments and across most model behaviours the choice of Sediment Transport formula was the dominate source of uncertainty in the CAESAR-Lisflood model, although there were some differences between the two catchments. Crucially, different parameters influenced the model behaviours in different ways, with Model Functions related to internal geomorphic changes responding in different ways to those related to sediment yields from the catchment outlet. This method of behavioural sensitivity analysis provides a useful method of assessing the performance of Landscape Evolution Models in the absence of data and methods for an Objective Function approach.

2018 ◽  
Vol 11 (12) ◽  
pp. 4873-4888 ◽  
Author(s):  
Christopher J. Skinner ◽  
Tom J. Coulthard ◽  
Wolfgang Schwanghart ◽  
Marco J. Van De Wiel ◽  
Greg Hancock

Abstract. The evaluation and verification of landscape evolution models (LEMs) has long been limited by a lack of suitable observational data and statistical measures which can fully capture the complexity of landscape changes. This lack of data limits the use of objective function based evaluation prolific in other modelling fields, and restricts the application of sensitivity analyses in the models and the consequent assessment of model uncertainties. To overcome this deficiency, a novel model function approach has been developed, with each model function representing an aspect of model behaviour, which allows for the application of sensitivity analyses. The model function approach is used to assess the relative sensitivity of the CAESAR-Lisflood LEM to a set of model parameters by applying the Morris method sensitivity analysis for two contrasting catchments. The test revealed that the model was most sensitive to the choice of the sediment transport formula for both catchments, and that each parameter influenced model behaviours differently, with model functions relating to internal geomorphic changes responding in a different way to those relating to the sediment yields from the catchment outlet. The model functions proved useful for providing a way of evaluating the sensitivity of LEMs in the absence of data and methods for an objective function approach.


2020 ◽  
Vol 34 (11) ◽  
pp. 1813-1830
Author(s):  
Daniel Erdal ◽  
Sinan Xiao ◽  
Wolfgang Nowak ◽  
Olaf A. Cirpka

Abstract Ensemble-based uncertainty quantification and global sensitivity analysis of environmental models requires generating large ensembles of parameter-sets. This can already be difficult when analyzing moderately complex models based on partial differential equations because many parameter combinations cause an implausible model behavior even though the individual parameters are within plausible ranges. In this work, we apply Gaussian Process Emulators (GPE) as surrogate models in a sampling scheme. In an active-training phase of the surrogate model, we target the behavioral boundary of the parameter space before sampling this behavioral part of the parameter space more evenly by passive sampling. Active learning increases the subsequent sampling efficiency, but its additional costs pay off only for a sufficiently large sample size. We exemplify our idea with a catchment-scale subsurface flow model with uncertain material properties, boundary conditions, and geometric descriptors of the geological structure. We then perform a global-sensitivity analysis of the resulting behavioral dataset using the active-subspace method, which requires approximating the local sensitivities of the target quantity with respect to all parameters at all sampled locations in parameter space. The Gaussian Process Emulator implicitly provides an analytical expression for this gradient, thus improving the accuracy of the active-subspace construction. When applying the GPE-based preselection, 70–90% of the samples were confirmed to be behavioral by running the full model, whereas only 0.5% of the samples were behavioral in standard Monte-Carlo sampling without preselection. The GPE method also provided local sensitivities at minimal additional costs.


Geofluids ◽  
2017 ◽  
Vol 2017 ◽  
pp. 1-16 ◽  
Author(s):  
Yuan Wang ◽  
Jie Ren ◽  
Shaobin Hu ◽  
Di Feng

Salt precipitation is generated near the injection well when dry supercritical carbon dioxide (scCO2) is injected into saline aquifers, and it can seriously impair the CO2 injectivity of the well. We used solid saturation (Ss) to map CO2 injectivity. Ss was used as the response variable for the sensitivity analysis, and the input variables included the CO2 injection rate (QCO2), salinity of the aquifer (XNaCl), empirical parameter m, air entry pressure (P0), maximum capillary pressure (Pmax), and liquid residual saturation (Splr and Sclr). Global sensitivity analysis methods, namely, the Morris method and Sobol method, were used. A significant increase in Ss was observed near the injection well, and the results of the two methods were similar: XNaCl had the greatest effect on Ss; the effect of P0 and Pmax on Ss was negligible. On the other hand, with these two methods, QCO2 had various effects on Ss: QCO2 had a large effect on Ss in the Morris method, but it had little effect on Ss in the Sobol method. We also found that a low QCO2 had a profound effect on Ss but that a high QCO2 had almost no effect on the Ss value.


2020 ◽  
Author(s):  
Gabriele Baroni ◽  
Till Francke

<p>Global sensitivity analysis has been recognized as a fundamental tool to assess the input-output model response and evaluate the role of different sources of uncertainty. Among the different methods, variance- and distribution-based (or also called moment-independent) methods have mostly been applied. The first method relies on variance decomposition while the second method compares the entire distributions. The combination of both methods has also been recognized to provide possibly a better assessment. However, the methods rely on different assumptions and the comparison of indices is not straightforward. For these reasons, the methods are commonly not integrated or even considered as alternative solutions. </p><p>In the present contribution, we show a new strategy to combine the two methods in an effective way to perform a comprehensive global sensitivity analysis based on a generic sampling design. The strategy is tested on three commonly-used analytic functions and one hydrological model. The strategy is compared to the state-of-the-art Jansen/Saltelli approach.</p><p>The results show that the new strategy quantifies main effect and interactions consistently. It also outperforms current best practices by converging with a lower number of model runs. For these reasons, the new strategy can be considered as a new and simple approach to perform global sensitivity analysis that can be easily integrated in any environmental models.</p>


2020 ◽  
Author(s):  
Alena Miftakhova

<p><span>A major tool that supports climate policy decisions, integrated assessment models are highly vulnerable to their initial assumptions and calibrations. Despite the broad literature rich in both single-model and multi-model sensitivity analyses, universal, well-established practices are still missing in this field. This paper endorses structured global sensitivity analysis (GSA) as an indispensable routine in climate–economic modeling. An application of a high-efficiency GSA method based on polynomial chaos expansions to DICE provides two insights. First, only global and comprehensive—as opposed to local or selective—sensitivity analysis delivers a trustworthy picture of the uncertainty propagated through the model. Second, careful treatment of the model’s structure throughout the analysis reconciles the results with established analytical insights—enhancing these insights with more details. The efficient GSA method provides a comprehensive decomposition of the uncertainty in a model’s output while minimizing computational costs, and is hence potentially applicable to models of higher complexity.</span></p>


Sign in / Sign up

Export Citation Format

Share Document