Stability criteria for Bayesian calibration of reservoir sedimentation models

Modeling reservoir sedimentation is particularly challenging due to the simultaneous simulation of shallow shores, tributary deltas, and deep waters. The shallow upstream parts of reservoirs, where deltaic avulsion and erosion processes occur, compete with the validity of modeling assumptions used to simulate the deposition of fine sediments in deep waters. We investigate how complex numerical models can be calibrated to accurately predict reservoir sedimentation in the presence of competing model simplifications and identify the importance of calibration parameters for prioritization in measurement campaigns. This study applies Bayesian calibration, a supervised learning technique using surrogate-assisted Bayesian inversion with a Gaussian Process Emulator to calibrate a two-dimensional (2d) hydro-morphodynamic model for simulating sedimentation processes in a reservoir in Albania. Four calibration parameters were fitted to obtain the statistically best possible simulation of bed level changes between 2016 and 2019 through two differently constraining data scenarios. One scenario included measurements from the entire upstream half of the reservoir. Another scenario only included measurements in the geospatially valid range of the numerical model. Model accuracy parameters, Bayesian model evidence, and the variability of the four calibration parameters indicate that Bayesian calibration only converges toward physically meaningful parameter combinations when the calibration nodes are in the valid range of the numerical model. The Bayesian approach also allowed for a comparison of multiple parameters and found that the dry bulk density of the deposited sediments is the most important factor for calibration.