First-order uncertainty analysis using Algorithmic Differentiation of morphodynamic models. (May 2016)

• Record Type: Journal Article
• Authors: Villaret, Catherine; Kopmann, Rebekka; Wyncoll, David; Riehme, Jan; Merkel, Uwe; Naumann, Uwe
• Is Part Of: Computers & geosciences. Volume 90(2016)Part B
• Page Start: 144
• Page End: 151
• Publication Date: 2016-05
• Subjects: Morphodynamics -- Process-based modelling -- Uncertainty analysis -- Sensitivity analysis -- Monte Carlo simulation -- Algorithmic Differentiation; Environmental policy -- Periodicals; 550.5
• DOI: 10.1016/j.cageo.2015.10.012
• Languages: English
• ISSNs: 0098-3004

Abstract

Abstract: We present here an efficient first-order second moment method using Algorithmic Differentiation (FOSM/AD) which can be applied to quantify uncertainty/sensitivities in morphodynamic models. Changes with respect to variable flow and sediment input parameters are estimated with machine accuracy using the technique of Algorithmic Differentiation (AD). This method is particularly attractive for process-based morphodynamic models like the Telemac-2D/Sisyphe model considering the large number of input parameters and CPU time associated to each simulation. The FOSM/AD method is applied to identify the relevant processes in a trench migration experiment (van Rijn, 1987 ). A Tangent Linear Model (TLM) of the Telemac-2D/Sisyphe morphodynamic model (release 6.2) was generated using the AD-enabled NAG Fortran compiler. One single run of the TLM is required per variable input parameter and results are then combined to calculate the total uncertainty. The limits of the FOSM/AD method have been assessed by comparison with Monte Carlo (MC) simulations. Similar results were obtained assuming small standard deviation of the variable input parameters. Both settling velocity and grain size have been identified as the most sensitive input parameters and the uncertainty as measured by the standard deviation of the calculated bed evolution increases with time. Highlights: A first-order second moment method (FOSM) is applied to quantify uncertainty. This method uses Algorithmic Differentiation (AD) and a Tangent Linear Model (TLM). The method is compared with Monte Carlo analysis in a trench migration test case. A TLM of the Telemac-2d/Sisyphe morphodynamic model has been applied. The FOSM/AD method is an efficient alternative to Monte Carlo simulations.